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

Timestamp:

06/23/14 09:47:02 (11 years ago)

Author:

Mathieu Morlighem

Message:

NEW: added CrouzeixRaviart elements in FS

File:

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

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

@@ -1975 +1975 @@
                         tria_node_ids[8]=iomodel->nodecounter+iomodel->numberofvertices+iomodel->numberofedges+iomodel->elements[3*index+2];
                         break;
+                case CrouzeixRaviartEnum:
+                        numnodes        = 10;
+                        tria_node_ids   = xNew<int>(numnodes);
+                        tria_node_ids[0]=iomodel->nodecounter+iomodel->elements[3*index+0];
+                        tria_node_ids[1]=iomodel->nodecounter+iomodel->elements[3*index+1];
+                        tria_node_ids[2]=iomodel->nodecounter+iomodel->elements[3*index+2];
+                        tria_node_ids[3]=iomodel->nodecounter+iomodel->numberofvertices+iomodel->elementtoedgeconnectivity[3*index+0]+1;
+                        tria_node_ids[4]=iomodel->nodecounter+iomodel->numberofvertices+iomodel->elementtoedgeconnectivity[3*index+1]+1;
+                        tria_node_ids[5]=iomodel->nodecounter+iomodel->numberofvertices+iomodel->elementtoedgeconnectivity[3*index+2]+1;
+                        tria_node_ids[6]=iomodel->nodecounter+iomodel->numberofvertices+iomodel->numberofedges+index+1;
+
+                        tria_node_ids[7]=iomodel->nodecounter+iomodel->numberofvertices+iomodel->numberofedges+iomodel->numberofelements+3*index+1;
+                        tria_node_ids[8]=iomodel->nodecounter+iomodel->numberofvertices+iomodel->numberofedges+iomodel->numberofelements+3*index+2;
+                        tria_node_ids[9]=iomodel->nodecounter+iomodel->numberofvertices+iomodel->numberofedges+iomodel->numberofelements+3*index+3;
+                        break;
                 default:
                         _error_("Finite element "<<EnumToStringx(finiteelement_type)<<" not supported yet");