[cig-commits] commit: Fix a bug in StressBC in 3D with tangential stresses

[Mercurial]

changeset:   930:3463bda1bc2a
tag:         tip
user:        Walter Landry <wlandry at caltech.edu>
date:        Sun Apr 29 05:36:42 2012 -0700
files:       Utils/src/StressBC.cxx
description:
Fix a bug in StressBC in 3D with tangential stresses


diff -r d54e97f15481 -r 3463bda1bc2a Utils/src/StressBC.cxx
--- a/Utils/src/StressBC.cxx	Wed Apr 18 00:31:14 2012 -0700
+++ b/Utils/src/StressBC.cxx	Sun Apr 29 05:36:42 2012 -0700
@@ -450,13 +450,13 @@ void _StressBC_AssembleElement( void* fo
                       if(dim==2)
                         {
                           geometric_factor=
-                            jac[surface][(d+1)%dim]*(d==face/2 ? -1 : 1);
+                            jac[surface][(d+1)%dim]*(d==local_norm ? 1 : -1);
                         }
                       else
-                        geometric_factor=jac[surface][(d+1)%dim]
-                          * jac[surface2][(d+2)%dim]
-                          - jac[surface][(d+2)%dim]
-                          * jac[surface2][(d+1)%dim];
+                        geometric_factor=
+                          (jac[surface][(d+1)%dim] * jac[surface2][(d+2)%dim]
+                           - jac[surface][(d+2)%dim] * jac[surface2][(d+1)%dim])
+                          *(d==local_norm ? 1 : -1);
                     }
                   else
                     {