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

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

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

revision 1388 by trankine, Fri Jan 11 07:45:58 2008 UTC revision 1739 by gross, Fri Aug 29 06:19:53 2008 UTC
# Line 35  void Finley_Mesh_distributeByRankOfDOF(F Line 35  void Finley_Mesh_distributeByRankOfDOF(F
35       mpiRankOfDOF=TMPMEMALLOC(self->Nodes->numNodes,Paso_MPI_rank);       mpiRankOfDOF=TMPMEMALLOC(self->Nodes->numNodes,Paso_MPI_rank);
36       if (!Finley_checkPtr(mpiRankOfDOF)) {       if (!Finley_checkPtr(mpiRankOfDOF)) {
37    
38    
39          Finley_NodeFile_assignMPIRankToDOFs(self->Nodes,mpiRankOfDOF,dof_distribution);          Finley_NodeFile_assignMPIRankToDOFs(self->Nodes,mpiRankOfDOF,dof_distribution);
40          /* first the elements are redistributed according to mpiRankOfDOF */          /* first the elements are redistributed according to mpiRankOfDOF */
41          /* at the input the Node tables refering to a the local labeling of the nodes */          /* at the input the Node tables refering to a the local labeling of the nodes */
# Line 60  void Finley_Mesh_distributeByRankOfDOF(F Line 61  void Finley_Mesh_distributeByRankOfDOF(F
61            for (n=0;n<self->Nodes->numNodes;n++) tmp_node_localDOF_map[n]=-1;            for (n=0;n<self->Nodes->numNodes;n++) tmp_node_localDOF_map[n]=-1;
62            #pragma omp parallel for private(n) schedule(static)            #pragma omp parallel for private(n) schedule(static)
63            for (n=0;n<self->Nodes->numNodes;n++) {            for (n=0;n<self->Nodes->numNodes;n++) {
64  #ifdef BOUNDS_CHECK               #ifdef BOUNDS_CHECK
65               if ((self->Nodes->globalDegreesOfFreedom[n]-min_id) >= len || (self->Nodes->globalDegreesOfFreedom[n]-min_id) < 0) { printf("BOUNDS_CHECK %s %d\n", __FILE__, __LINE__); exit(1); }               if ((self->Nodes->globalDegreesOfFreedom[n]-min_id) >= len || (self->Nodes->globalDegreesOfFreedom[n]-min_id) < 0) { printf("BOUNDS_CHECK %s %d\n", __FILE__, __LINE__); exit(1); }
66  #endif               #endif
67           tmp_node_localDOF_mask[self->Nodes->globalDegreesOfFreedom[n]-min_id]=n;           tmp_node_localDOF_mask[self->Nodes->globalDegreesOfFreedom[n]-min_id]=n;
68        }        }
69        

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

  ViewVC Help
Powered by ViewVC 1.1.26