/[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 1315 by ksteube, Tue Sep 25 02:41:13 2007 UTC revision 1326 by ksteube, Mon Oct 1 08:10:41 2007 UTC
# Line 57  void Finley_Mesh_distributeByRankOfDOF(F Line 57  void Finley_Mesh_distributeByRankOfDOF(F
57            #pragma omp parallel for private(n) schedule(static)            #pragma omp parallel for private(n) schedule(static)
58            for (n=0;n<len;n++) tmp_node_localDOF_mask[n]=-1;            for (n=0;n<len;n++) tmp_node_localDOF_mask[n]=-1;
59            #pragma omp parallel for private(n) schedule(static)            #pragma omp parallel for private(n) schedule(static)
60            for (n=0;n<self->Nodes->numNodes;n++) tmp_node_localDOF_mask[self->Nodes->globalDegreesOfFreedom[n]-min_id]=n;            for (n=0;n<self->Nodes->numNodes;n++) {
61    #ifdef BOUNDS_CHECK
62                 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); }
63    #endif
64             tmp_node_localDOF_mask[self->Nodes->globalDegreesOfFreedom[n]-min_id]=n;
65          }
66        
67            numDOFs=0;            numDOFs=0;
68            for  (n=0;n<len;n++) {            for  (n=0;n<len;n++) {
69                k=tmp_node_localDOF_mask[n];                k=tmp_node_localDOF_mask[n];
70                if (k>=0) {                if (k>=0) {
71    #ifdef BOUNDS_CHECK
72                     if ((k) >= len || (k) < 0) { printf("BOUNDS_CHECK %s %d k=%d n=%d max_id=%d min_id=%d len=%d\n", __FILE__, __LINE__, k, n, max_id, min_id, len); exit(1); }
73    #endif
74                   tmp_node_localDOF_map[k]=numDOFs;                   tmp_node_localDOF_map[k]=numDOFs;
75                   numDOFs++;                   numDOFs++;
76                }                }

Legend:
Removed from v.1315  
changed lines
  Added in v.1326

  ViewVC Help
Powered by ViewVC 1.1.26