[cig-commits] r5345 - mc/3D/CitcomS/trunk/lib

[tan2 at geodynamics.org]

Author: tan2
Date: 2006-11-22 13:25:47 -0800 (Wed, 22 Nov 2006)
New Revision: 5345

Modified:
   mc/3D/CitcomS/trunk/lib/Full_parallel_related.c
Log:
Ignoring overlapping nodes when counting # of global nodes

Modified: mc/3D/CitcomS/trunk/lib/Full_parallel_related.c
===================================================================
--- mc/3D/CitcomS/trunk/lib/Full_parallel_related.c	2006-11-22 20:25:40 UTC (rev 5344)
+++ mc/3D/CitcomS/trunk/lib/Full_parallel_related.c	2006-11-22 21:25:47 UTC (rev 5345)
@@ -464,16 +464,17 @@
       }       /* end for m */
     }   /* end for level */
 
-/*  ii=0; */
-/*  for (m=1;m<=E->sphere.caps_per_proc;m++) */
-/*     for (node=1;node<=E->lmesh.nno;node++) */
-/*       if(E->node[m][node] & SKIPS) */
-/*         ii+=1; */
+ /* count # of global nodes, ignoring overlapping nodes */
+ ii=0;
+ for (m=1;m<=E->sphere.caps_per_proc;m++)
+    for (node=1;node<=E->lmesh.nno;node++)
+      if(E->node[m][node] & SKIPS)
+        ++ii;
 
-/*  MPI_Allreduce(&ii, &node  ,1,MPI_INT,MPI_SUM,E->parallel.world); */
+ MPI_Allreduce(&ii, &node , 1, MPI_INT,MPI_SUM, E->parallel.world);
 
-/*  E->mesh.nno = E->lmesh.nno*E->parallel.nproc - node - 2*E->mesh.noz; */
-/*  E->mesh.neq = E->mesh.nno*3; */
+ E->mesh.nno = E->lmesh.nno*E->parallel.nproc - node - 2*E->mesh.noz;
+ E->mesh.neq = E->mesh.nno*3;
 
 if (E->control.verbose) {
  fprintf(E->fp_out,"output_shared_nodes %d \n",E->parallel.me);