Changeset 3085 for inundation/pymetis/metis-4.0/Lib/meshpart.c

Timestamp:

Jun 5, 2006, 5:03:48 PM (19 years ago)

Author:

ole

Message:

Moved pyvolution.mesh.py to pyvolution.test_mesh.py in order to avoid confusion with pmesh.

File:

• inundation/pymetis/metis-4.0/Lib/meshpart.c (modified) (4 diffs)

inundation/pymetis/metis-4.0/Lib/meshpart.c

@@ -31 +31 @@
   esize = esizes[*etype];
 
+  printf("Inside METIS_PartMeshNodal\n");
+  
   if (*numflag == 1)
     ChangeMesh2CNumbering((*ne)*esize, elmnts);
@@ -37 +39 @@
   adjncy = idxmalloc(20*(*nn), "METIS_MESHPARTNODAL: adjncy");
 
+  printf("Calling METIS_MeshToNodal\n");
   METIS_MeshToNodal(ne, nn, elmnts, etype, &pnumflag, xadj, adjncy);
-
+  printf("After METIS_MeshToNodal\n");
+    
   adjncy = realloc(adjncy, xadj[*nn]*sizeof(idxtype));
 
   options[0] = 0;
+  printf("Calling METIS_PartGraphKway\n");  
   METIS_PartGraphKway(nn, xadj, adjncy, NULL, NULL, &wgtflag, &pnumflag, nparts, options, edgecut, npart);
 
   /* OK, now compute an element partition based on the nodal partition npart */
+  printf("Starting loop 1\n");  
   idxset(*ne, -1, epart);
   pwgts = idxsmalloc(*nparts, 0, "METIS_MESHPARTNODAL: pwgts");
@@ -59 +65 @@
   }
 
+  printf("Starting loop 2\n");    
   maxpwgt = 1.03*(*ne)/(*nparts);
   for (i=0; i<*ne; i++) {
@@ -96 +103 @@
   }
 
+  printf("ChangeMesh2FNumbering2\n");    
   if (*numflag == 1)
     ChangeMesh2FNumbering2((*ne)*esize, elmnts, *ne, *nn, epart, npart);
 
+  printf("Free\n");        
   GKfree(&xadj, &adjncy, &pwgts, LTERM);
+  printf("Done\n");      
 
 }