/[escript]/trunk/finley/src/Mesh_read.c
ViewVC logotype

Diff of /trunk/finley/src/Mesh_read.c

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

revision 1738 by ksteube, Wed Aug 27 05:47:03 2008 UTC revision 1739 by gross, Fri Aug 29 06:19:53 2008 UTC
# Line 106  Finley_Mesh* Finley_Mesh_read(char* fnam Line 106  Finley_Mesh* Finley_Mesh_read(char* fnam
106                      for (i0 = 0; i0 < numEle; i0++) {                      for (i0 = 0; i0 < numEle; i0++) {
107                        fscanf(fileHandle_p, "%d %d", &mesh_p->Elements->Id[i0], &mesh_p->Elements->Tag[i0]);                        fscanf(fileHandle_p, "%d %d", &mesh_p->Elements->Id[i0], &mesh_p->Elements->Tag[i0]);
108                        mesh_p->Elements->Color[i0]=i0;                        mesh_p->Elements->Color[i0]=i0;
109                          mesh_p->Elements->Owner[i0]=0;
110                        for (i1 = 0; i1 < mesh_p->Elements->ReferenceElement->Type->numNodes; i1++) {                        for (i1 = 0; i1 < mesh_p->Elements->ReferenceElement->Type->numNodes; i1++) {
111                             fscanf(fileHandle_p, " %d",                             fscanf(fileHandle_p, " %d",
112                                &mesh_p->Elements->Nodes[INDEX2(i1, i0, mesh_p->Elements->ReferenceElement->Type->numNodes)]);                                &mesh_p->Elements->Nodes[INDEX2(i1, i0, mesh_p->Elements->ReferenceElement->Type->numNodes)]);
# Line 133  Finley_Mesh* Finley_Mesh_read(char* fnam Line 134  Finley_Mesh* Finley_Mesh_read(char* fnam
134                        for (i0 = 0; i0 < numEle; i0++) {                        for (i0 = 0; i0 < numEle; i0++) {
135                          fscanf(fileHandle_p, "%d %d", &mesh_p->FaceElements->Id[i0], &mesh_p->FaceElements->Tag[i0]);                          fscanf(fileHandle_p, "%d %d", &mesh_p->FaceElements->Id[i0], &mesh_p->FaceElements->Tag[i0]);
136                          mesh_p->FaceElements->Color[i0]=i0;                          mesh_p->FaceElements->Color[i0]=i0;
137                            mesh_p->FaceElements->Owner[i0]=0;
138                          for (i1 = 0; i1 < mesh_p->FaceElements->ReferenceElement->Type->numNodes; i1++) {                          for (i1 = 0; i1 < mesh_p->FaceElements->ReferenceElement->Type->numNodes; i1++) {
139                               fscanf(fileHandle_p, " %d",                               fscanf(fileHandle_p, " %d",
140                                  &mesh_p->FaceElements->Nodes[INDEX2(i1, i0, mesh_p->FaceElements->ReferenceElement->Type->numNodes)]);                                  &mesh_p->FaceElements->Nodes[INDEX2(i1, i0, mesh_p->FaceElements->ReferenceElement->Type->numNodes)]);
# Line 160  Finley_Mesh* Finley_Mesh_read(char* fnam Line 162  Finley_Mesh* Finley_Mesh_read(char* fnam
162                        for (i0 = 0; i0 < numEle; i0++) {                        for (i0 = 0; i0 < numEle; i0++) {
163                          fscanf(fileHandle_p, "%d %d", &mesh_p->ContactElements->Id[i0], &mesh_p->ContactElements->Tag[i0]);                          fscanf(fileHandle_p, "%d %d", &mesh_p->ContactElements->Id[i0], &mesh_p->ContactElements->Tag[i0]);
164                          mesh_p->ContactElements->Color[i0]=i0;                          mesh_p->ContactElements->Color[i0]=i0;
165                            mesh_p->ContactElements->Owner[i0]=0;
166                          for (i1 = 0; i1 < mesh_p->ContactElements->ReferenceElement->Type->numNodes; i1++) {                          for (i1 = 0; i1 < mesh_p->ContactElements->ReferenceElement->Type->numNodes; i1++) {
167                              fscanf(fileHandle_p, " %d",                              fscanf(fileHandle_p, " %d",
168                                 &mesh_p->ContactElements->Nodes[INDEX2(i1, i0, mesh_p->ContactElements->ReferenceElement->Type->numNodes)]);                                 &mesh_p->ContactElements->Nodes[INDEX2(i1, i0, mesh_p->ContactElements->ReferenceElement->Type->numNodes)]);
# Line 187  Finley_Mesh* Finley_Mesh_read(char* fnam Line 190  Finley_Mesh* Finley_Mesh_read(char* fnam
190                     for (i0 = 0; i0 < numEle; i0++) {                     for (i0 = 0; i0 < numEle; i0++) {
191                       fscanf(fileHandle_p, "%d %d", &mesh_p->Points->Id[i0], &mesh_p->Points->Tag[i0]);                       fscanf(fileHandle_p, "%d %d", &mesh_p->Points->Id[i0], &mesh_p->Points->Tag[i0]);
192                       mesh_p->Points->Color[i0]=i0;                       mesh_p->Points->Color[i0]=i0;
193                         mesh_p->Points->Owner[i0]=0;
194                       for (i1 = 0; i1 < mesh_p->Points->ReferenceElement->Type->numNodes; i1++) {                       for (i1 = 0; i1 < mesh_p->Points->ReferenceElement->Type->numNodes; i1++) {
195                           fscanf(fileHandle_p, " %d",                           fscanf(fileHandle_p, " %d",
196                              &mesh_p->Points->Nodes[INDEX2(i1, i0, mesh_p->Points->ReferenceElement->Type->numNodes)]);                              &mesh_p->Points->Nodes[INDEX2(i1, i0, mesh_p->Points->ReferenceElement->Type->numNodes)]);
# Line 474  Finley_Mesh* Finley_Mesh_read_MPI(char* Line 478  Finley_Mesh* Finley_Mesh_read_MPI(char*
478        for (i0=0; i0<chunkEle; i0++) {        for (i0=0; i0<chunkEle; i0++) {
479          mesh_p->Elements->Id[i0]    = tempInts[i0*(2+numNodes)+0];          mesh_p->Elements->Id[i0]    = tempInts[i0*(2+numNodes)+0];
480          mesh_p->Elements->Tag[i0]   = tempInts[i0*(2+numNodes)+1];          mesh_p->Elements->Tag[i0]   = tempInts[i0*(2+numNodes)+1];
481                mesh_p->Elements->Owner[i0]  =mpi_info->rank;
482              mesh_p->Elements->Color[i0] = i0;              mesh_p->Elements->Color[i0] = i0;
483          for (i1 = 0; i1 < numNodes; i1++) {          for (i1 = 0; i1 < numNodes; i1++) {
484            mesh_p->Elements->Nodes[INDEX2(i1, i0, numNodes)] = tempInts[i0*(2+numNodes)+2+i1];            mesh_p->Elements->Nodes[INDEX2(i1, i0, numNodes)] = tempInts[i0*(2+numNodes)+2+i1];
# Line 573  Finley_Mesh* Finley_Mesh_read_MPI(char* Line 578  Finley_Mesh* Finley_Mesh_read_MPI(char*
578        for (i0=0; i0<chunkEle; i0++) {        for (i0=0; i0<chunkEle; i0++) {
579          mesh_p->FaceElements->Id[i0]    = tempInts[i0*(2+numNodes)+0];          mesh_p->FaceElements->Id[i0]    = tempInts[i0*(2+numNodes)+0];
580          mesh_p->FaceElements->Tag[i0]   = tempInts[i0*(2+numNodes)+1];          mesh_p->FaceElements->Tag[i0]   = tempInts[i0*(2+numNodes)+1];
581                mesh_p->FaceElements->Owner[i0]  =mpi_info->rank;
582              mesh_p->FaceElements->Color[i0] = i0;              mesh_p->FaceElements->Color[i0] = i0;
583          for (i1 = 0; i1 < numNodes; i1++) {          for (i1 = 0; i1 < numNodes; i1++) {
584            mesh_p->FaceElements->Nodes[INDEX2(i1, i0, numNodes)] = tempInts[i0*(2+numNodes)+2+i1];            mesh_p->FaceElements->Nodes[INDEX2(i1, i0, numNodes)] = tempInts[i0*(2+numNodes)+2+i1];
# Line 672  Finley_Mesh* Finley_Mesh_read_MPI(char* Line 678  Finley_Mesh* Finley_Mesh_read_MPI(char*
678        for (i0=0; i0<chunkEle; i0++) {        for (i0=0; i0<chunkEle; i0++) {
679          mesh_p->ContactElements->Id[i0] = tempInts[i0*(2+numNodes)+0];          mesh_p->ContactElements->Id[i0] = tempInts[i0*(2+numNodes)+0];
680          mesh_p->ContactElements->Tag[i0]    = tempInts[i0*(2+numNodes)+1];          mesh_p->ContactElements->Tag[i0]    = tempInts[i0*(2+numNodes)+1];
681                mesh_p->ContactElements->Owner[i0]  =mpi_info->rank;
682              mesh_p->ContactElements->Color[i0] = i0;              mesh_p->ContactElements->Color[i0] = i0;
683          for (i1 = 0; i1 < numNodes; i1++) {          for (i1 = 0; i1 < numNodes; i1++) {
684            mesh_p->ContactElements->Nodes[INDEX2(i1, i0, numNodes)] = tempInts[i0*(2+numNodes)+2+i1];            mesh_p->ContactElements->Nodes[INDEX2(i1, i0, numNodes)] = tempInts[i0*(2+numNodes)+2+i1];
# Line 771  Finley_Mesh* Finley_Mesh_read_MPI(char* Line 778  Finley_Mesh* Finley_Mesh_read_MPI(char*
778        for (i0=0; i0<chunkEle; i0++) {        for (i0=0; i0<chunkEle; i0++) {
779          mesh_p->Points->Id[i0]  = tempInts[i0*(2+numNodes)+0];          mesh_p->Points->Id[i0]  = tempInts[i0*(2+numNodes)+0];
780          mesh_p->Points->Tag[i0] = tempInts[i0*(2+numNodes)+1];          mesh_p->Points->Tag[i0] = tempInts[i0*(2+numNodes)+1];
781                mesh_p->Points->Owner[i0]  =mpi_info->rank;
782              mesh_p->Points->Color[i0] = i0;              mesh_p->Points->Color[i0] = i0;
783          for (i1 = 0; i1 < numNodes; i1++) {          for (i1 = 0; i1 < numNodes; i1++) {
784            mesh_p->Points->Nodes[INDEX2(i1, i0, numNodes)] = tempInts[i0*(2+numNodes)+2+i1];            mesh_p->Points->Nodes[INDEX2(i1, i0, numNodes)] = tempInts[i0*(2+numNodes)+2+i1];

Legend:
Removed from v.1738  
changed lines
  Added in v.1739

  ViewVC Help
Powered by ViewVC 1.1.26