Changeset 17846

Timestamp:

04/25/14 11:17:42 (11 years ago)

Author:

Mathieu Morlighem

Message:

NEW: Added Tets

File:

• issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateFaces.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateFaces.cpp

@@ -153 +153 @@
                         elementfaces[cols*4+0] = 4; elementfaces[cols*4+1] = 0;  elementfaces[cols*4+2] = 1; elementfaces[cols*4+3] = 4;  elementfaces[cols*4+4] = 3;
                         break;
+                case TetraEnum:
+                        elementnbv = 4; /*Number of vertices per element*/
+                        elementnbf = 4; /*Number of faces per element*/
+                        facemaxnbv = 3; /*Maximum number of vertices per face*/
+                        cols       = facemaxnbv + 1;
+                        elementfaces         = xNew<int>(elementnbf*cols);
+                        elementfaces_markers = xNew<int>(elementnbf);
+                        /*4 triangles*/
+                        elementfaces_markers[0] = 1;
+                        elementfaces_markers[1] = 1;
+                        elementfaces_markers[2] = 1;
+                        elementfaces[cols*0+0] = 3; elementfaces[cols*0+1] = 0;  elementfaces[cols*0+2] = 1; elementfaces[cols*0+3] = 2;
+                        elementfaces[cols*1+0] = 3; elementfaces[cols*1+1] = 0;  elementfaces[cols*1+2] = 3; elementfaces[cols*1+3] = 1;
+                        elementfaces[cols*2+0] = 3; elementfaces[cols*2+1] = 1;  elementfaces[cols*2+2] = 3; elementfaces[cols*2+3] = 2;
+                        elementfaces[cols*3+0] = 3; elementfaces[cols*3+1] = 0;  elementfaces[cols*3+2] = 2; elementfaces[cols*3+3] = 3;
+                        break;
                 default:
                 _error_("mesh "<< EnumToStringx(iomodel->meshelementtype) <<" not supported");