Changeset 15611 for issm/trunk-jpl/src/c/modules/IoModelToConstraintsx/IoModelToConstraintsx.cpp

Timestamp:

07/25/13 09:32:50 (12 years ago)

Author:

Mathieu Morlighem

Message:

NEW: make the distinction between faces and edges for 3d Quadratic elements (quadratic requires edges, DG requires faces, the format is not the same)

File:

• issm/trunk-jpl/src/c/modules/IoModelToConstraintsx/IoModelToConstraintsx.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/modules/IoModelToConstraintsx/IoModelToConstraintsx.cpp ¶

@@ -57 +57 @@
                         for(i=0;i<iomodel->numberofedges;i++){
                                 if(my_edges[i]){
-                                        v1 = iomodel->edges[4*i+0]-1;
-                                        v2 = iomodel->edges[4*i+1]-1;
+                                        v1 = iomodel->edges[2*i+0]-1;
+                                        v2 = iomodel->edges[2*i+1]-1;
                                         if(!xIsNan<IssmDouble>(spcdata[v1]) && !xIsNan<IssmDouble>(spcdata[v2])){
                                                 constraints->AddObject(new SpcStatic(iomodel->constraintcounter+count+1,iomodel->nodecounter+iomodel->numberofvertices+i+1,
@@ -99 +99 @@
                         for(i=0;i<iomodel->numberofedges;i++){
                                 if(my_edges[i]){
-                                        v1 = iomodel->edges[4*i+0]-1;
-                                        v2 = iomodel->edges[4*i+1]-1;
+                                        v1 = iomodel->edges[2*i+0]-1;
+                                        v2 = iomodel->edges[2*i+1]-1;
                                         values=xNew<IssmDouble>(N);
                                         spcpresent=false;