[cig-commits] r5375 - in long/3D/Gale/trunk: . src/PICellerator/Utils/src

[walter at geodynamics.org]

Author: walter
Date: 2006-11-29 10:08:10 -0800 (Wed, 29 Nov 2006)
New Revision: 5375

Modified:
   long/3D/Gale/trunk/
   long/3D/Gale/trunk/src/PICellerator/Utils/src/StressBC.c
Log:
 r1120 at earth:  boo | 2006-11-29 01:21:51 -0800
 Make boundary detection more robust in StressBC



Property changes on: long/3D/Gale/trunk
___________________________________________________________________
Name: svk:merge
   - 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:1119
   + 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:1120

Modified: long/3D/Gale/trunk/src/PICellerator/Utils/src/StressBC.c
===================================================================
--- long/3D/Gale/trunk/src/PICellerator/Utils/src/StressBC.c	2006-11-29 18:08:07 UTC (rev 5374)
+++ long/3D/Gale/trunk/src/PICellerator/Utils/src/StressBC.c	2006-11-29 18:08:10 UTC (rev 5375)
@@ -344,14 +344,18 @@
 
 double _StressBC_CalcForce( void* forceTerm, Mesh* mesh, Element_DomainIndex lElement_I, Index node_I ) {
 	StressBC*               self               = (StressBC*) forceTerm;
-        Coord *coord;
-
-        coord=&(Mesh_CoordAt(mesh,mesh->elementNodeTbl[lElement_I][node_I]));
-
-        if((*coord)[1]<0.001)
-          {
+	MeshLayout*		meshLayout = mesh->layout;
+	HexaMD*			decomp = (HexaMD*)meshLayout->decomp;
+        
+        IJK			ijk;
+		
+        RegularMeshUtils_Node_1DTo3D( decomp, mesh->elementNodeTbl[lElement_I][node_I],
+                                      &ijk[0], &ijk[1], &ijk[2] );      
+        if( ijk[1] == 0 ) {
+          Coord *coord;
 /*             return self->force; */
-            return exp(-(0.5-(*coord)[0])*(0.5-(*coord)[0])/0.01)*self->force;
-          }
+          coord=&(Mesh_CoordAt(mesh,mesh->elementNodeTbl[lElement_I][node_I]));
+          return exp(-(0.5-(*coord)[0])*(0.5-(*coord)[0])/0.01)*self->force;
+        }
 	return 0.0;
 }