Changeset 24393

Timestamp:

11/25/19 11:49:45 (5 years ago)

Author:

Mathieu Morlighem

Message:

CHG: more info for METIS

File:

• issm/trunk-jpl/src/c/toolkits/metis/patches/METIS_PartMeshNodalPatch.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/toolkits/metis/patches/METIS_PartMeshNodalPatch.cpp

@@ -13 +13 @@
         METIS_PartMeshNodal(pnumberofelements,pnumberofnodes, index, petype, pnumflag, pnum_procs, pedgecut, epart, npart); 
         #elif _METIS_VERSION_ == 5
+
         /*This interface is heavily changed. More options, different ways of meshing, etc ...: */ 
-        int i;
 
+        /*Create options*/
         idx_t options[METIS_NOPTIONS];
-        idx_t objval;
-        idx_t* eptr=NULL;
-        idx_t  k=0;
-        real_t* tpwgts=NULL;
-
-        /*setup options: */
         METIS_SetDefaultOptions(options);
 
-        options[METIS_OPTION_PTYPE]   = 1;
-        options[METIS_OPTION_OBJTYPE] = 0;
-        options[METIS_OPTION_CTYPE]   = 1;
-        options[METIS_OPTION_IPTYPE]  = 4;
-        options[METIS_OPTION_RTYPE]   = 1;
-        options[METIS_OPTION_DBGLVL]  = 0;
-        options[METIS_OPTION_UFACTOR] = 30;
-        options[METIS_OPTION_MINCONN] = 0;
-        options[METIS_OPTION_CONTIG]  = 0;
-        options[METIS_OPTION_SEED]    = -1;
-        options[METIS_OPTION_NITER]   = 10;
-        options[METIS_OPTION_NCUTS]   = 1;
+        options[METIS_OPTION_PTYPE]   = METIS_PTYPE_KWAY;     /* partitioning method  */
+        options[METIS_OPTION_OBJTYPE] = METIS_OBJTYPE_VOL;    /* type of objective */
+        options[METIS_OPTION_CTYPE]   = METIS_CTYPE_SHEM;     /* matching scheme to be used during coarsening.*/
+        options[METIS_OPTION_IPTYPE]  = METIS_IPTYPE_METISRB; /* algorithm used during initial partitioning*/
+        options[METIS_OPTION_RTYPE]   = METIS_RTYPE_GREEDY;   /* algorithm used for refinement*/
+        options[METIS_OPTION_DBGLVL]  = 0;                    /* amount of progress/debugging information will be printed */
+        options[METIS_OPTION_UFACTOR] = 30;                   /* maximum allowed load imbalance among the partitions*/
+        options[METIS_OPTION_MINCONN] = 0;                    /* explicitly minimize the maximum connectivity ?*/
+        options[METIS_OPTION_CONTIG]  = 0;                    /* force contiguous partitions?*/
+        options[METIS_OPTION_SEED]    = -1;                   /* seed for the random number generator*/
+        options[METIS_OPTION_NITER]   = 10;                   /* number of iterations for the refinement algorithms*/
+        options[METIS_OPTION_NCUTS]   = 1;                    /* number of different partitionings that it will compute*/
 
         /*create eptr: */
-        eptr=xNew<idx_t>((*pnumberofelements+1));
+        idx_t  k=0;
+        idx_t* eptr=xNew<idx_t>((*pnumberofelements+1));
         eptr[0]=0;
-        for(i=0;i<*pnumberofelements;i++){
+        for(int i=0;i<*pnumberofelements;i++){
                 k+=3;
                 eptr[i+1]=k;
@@ -47 +43 @@
 
         /*create tpwgts: */
-        tpwgts=xNew<real_t>(*pnum_procs);
-        for(i=0;i<*pnum_procs;i++){
+        real_t* tpwgts=xNew<real_t>(*pnum_procs);
+        for(int i=0;i<*pnum_procs;i++){
                 tpwgts[i]=1.0/(*pnum_procs);
         }
 
+        /*Call METIS*/
+        idx_t objval;
         METIS_PartMeshNodal(pnumberofelements,pnumberofnodes, eptr, index,
                         NULL, NULL, pnum_procs, tpwgts, options, &objval,epart, npart);