[cig-commits] commit: Clean up

Mercurial hg at geodynamics.org
Tue Oct 18 20:13:56 PDT 2011


changeset:   892:8d657451f1ce
user:        Walter Landry <wlandry at caltech.edu>
date:        Tue Oct 18 20:11:41 2011 -0700
files:       Utils/src/StressBC.cxx
description:
Clean up


diff -r 06656e8d7824 -r 8d657451f1ce Utils/src/StressBC.cxx
--- a/Utils/src/StressBC.cxx	Tue Oct 18 16:41:42 2011 -0700
+++ b/Utils/src/StressBC.cxx	Tue Oct 18 20:11:41 2011 -0700
@@ -394,9 +394,8 @@ void _StressBC_AssembleElement( void* fo
                 }
               double Ni[num_nodes];
               ElementType_EvaluateShapeFunctionsAt(mesh->feElType,local_coord,Ni);
-              ElementType_Jacobian_AxisIndependent(mesh->feElType,mesh,lElement_I,
-                                                   local_coord,dim,jac,NULL,
-                                                   0,1,2);
+              ElementType_Jacobian(mesh->feElType,mesh,lElement_I,
+                                   local_coord,dim,jac,NULL);
               FeMesh_CoordLocalToGlobal(mesh,lElement_I,local_coord,global_coord);
 
               /* Get the applied stress */
@@ -474,8 +473,7 @@ void _StressBC_AssembleElement( void* fo
                       RegularMeshUtils_Node_1DTo3D
                         (mesh,Mesh_DomainToGlobal(mesh,MT_VERTEX,
                                                   elementNodes[node]),ijk);
-                      double total=stress*weights[i]*geometric_factor
-                        *Ni[node];
+                      double total=stress*weights[i]*geometric_factor*Ni[node];
                       if(dim==3)
                         {
                           total*=weights[j];



More information about the CIG-COMMITS mailing list