Changeset 17855

Timestamp:

04/28/14 11:48:49 (11 years ago)

Author:

seroussi

Message:

BUG: fixed NormalBase and NormalTop for flowband models

File:

• issm/trunk-jpl/src/c/classes/Elements/Tria.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/classes/Elements/Tria.cpp

@@ -1665 +1665 @@
 
         /*Build unit outward pointing vector*/
+        int index1,index2;
         IssmDouble vector[2];
         IssmDouble norm;
 
-        vector[0]=xyz_list[1*3+0] - xyz_list[0*3+0];
-        vector[1]=xyz_list[1*3+1] - xyz_list[0*3+1];
+        this->EdgeOnBaseIndices(&index1,&index2);
+        vector[0]=xyz_list[index2*3+0] - xyz_list[index1*3+0];
+        vector[1]=xyz_list[index2*3+1] - xyz_list[index1*3+1];
 
         norm=sqrt(vector[0]*vector[0] + vector[1]*vector[1]);
@@ -1681 +1683 @@
 
         /*Build unit outward pointing vector*/
+        int index1,index2;
         IssmDouble vector[2];
         IssmDouble norm;
 
-        vector[0]=xyz_list[1*3+0] - xyz_list[0*3+0];
-        vector[1]=xyz_list[1*3+1] - xyz_list[0*3+1];
+        this->EdgeOnSurfaceIndices(&index1,&index2);
+        vector[0]=xyz_list[index2*3+0] - xyz_list[index1*3+0];
+        vector[1]=xyz_list[index2*3+1] - xyz_list[index1*3+1];
 
         norm=sqrt(vector[0]*vector[0] + vector[1]*vector[1]);