/[escript]/trunk-mpi-branch/finley/src/Mesh_readGmsh.c
ViewVC logotype

Diff of /trunk-mpi-branch/finley/src/Mesh_readGmsh.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1306 by ksteube, Tue Sep 18 05:51:09 2007 UTC revision 1307 by gross, Tue Sep 18 08:53:11 2007 UTC
# Line 235  Finley_Mesh* Finley_Mesh_readGmsh(char* Line 235  Finley_Mesh* Finley_Mesh_readGmsh(char*
235         }         }
236                if (! Finley_noError()) break;                if (! Finley_noError()) break;
237                for(j = 0; j < numNodesPerElement; j++) fscanf(fileHandle_p, "%d", &vertices[INDEX2(j,e,MAX_numNodes_gmsh)]);                for(j = 0; j < numNodesPerElement; j++) fscanf(fileHandle_p, "%d", &vertices[INDEX2(j,e,MAX_numNodes_gmsh)]);
238                  /* for tet10 the last two nodes need to be swapped */
239                  if (element_type[e]==Tet10) {
240                       itmp=vertices[INDEX2(9,e,MAX_numNodes_gmsh)];
241                       vertices[INDEX2(9,e,MAX_numNodes_gmsh)]=vertices[INDEX2(8,e,MAX_numNodes_gmsh)];
242                       vertices[INDEX2(8,e,MAX_numNodes_gmsh)]=itmp;
243                  }
244              }              }
245              /* all elements have been read, now we have to identify the elements for finley */              /* all elements have been read, now we have to identify the elements for finley */
246                    
# Line 283  Finley_Mesh* Finley_Mesh_readGmsh(char* Line 289  Finley_Mesh* Finley_Mesh_readGmsh(char*
289                              mesh_p->Elements->Id[numElements]=id[e];                              mesh_p->Elements->Id[numElements]=id[e];
290                              mesh_p->Elements->Tag[numElements]=tag[e];                              mesh_p->Elements->Tag[numElements]=tag[e];
291                              mesh_p->Elements->Color[numElements]=numElements;                              mesh_p->Elements->Color[numElements]=numElements;
292                              for (j = 0; j<  mesh_p->Elements->ReferenceElement->Type->numNodes; ++j)                              for (j = 0; j<  mesh_p->Elements->ReferenceElement->Type->numNodes; ++j)  {
293                                       mesh_p->Elements->Nodes[INDEX2(j, numElements, mesh_p->Elements->ReferenceElement->Type->numNodes)]=                                    mesh_p->Elements->Nodes[INDEX2(j, numElements, mesh_p->Elements->ReferenceElement->Type->numNodes)]=vertices[INDEX2(j,e,MAX_numNodes_gmsh)];
294                                                                                                       vertices[INDEX2(j,e,MAX_numNodes_gmsh)];                              }
295                              numElements++;                              numElements++;
296                           } else if (element_type[e] == final_face_element_type) {                           } else if (element_type[e] == final_face_element_type) {
297                              mesh_p->FaceElements->Id[numFaceElements]=id[e];                              mesh_p->FaceElements->Id[numFaceElements]=id[e];
298                              mesh_p->FaceElements->Tag[numFaceElements]=tag[e];                              mesh_p->FaceElements->Tag[numFaceElements]=tag[e];
299                              mesh_p->FaceElements->Color[numFaceElements]=numFaceElements;                              mesh_p->FaceElements->Color[numFaceElements]=numFaceElements;
300                              for (j = 0; j<  mesh_p->FaceElements->ReferenceElement->Type->numNodes; ++j)                              for (j = 0; j<  mesh_p->FaceElements->ReferenceElement->Type->numNodes; ++j) {
301                                       mesh_p->FaceElements->Nodes[INDEX2(j, numFaceElements, mesh_p->FaceElements->ReferenceElement->Type->numNodes)]=                                       mesh_p->FaceElements->Nodes[INDEX2(j, numFaceElements, mesh_p->FaceElements->ReferenceElement->Type->numNodes)]=vertices[INDEX2(j,e,MAX_numNodes_gmsh)];
302                                                                                                       vertices[INDEX2(j,e,MAX_numNodes_gmsh)];                              }
303                              numFaceElements++;                              numFaceElements++;
304                           }                           }
305                        }                        }

Legend:
Removed from v.1306  
changed lines
  Added in v.1307

  ViewVC Help
Powered by ViewVC 1.1.26