Changeset 16328

Timestamp:

10/08/13 09:03:34 (11 years ago)

Author:

Mathieu Morlighem

Message:

CHG: added support for mesh 2d vertical partitioning

File:

• issm/trunk-jpl/src/c/modules/MeshPartitionx/MeshPartitionx.h (modified) (1 diff)

issm/trunk-jpl/src/c/modules/MeshPartitionx/MeshPartitionx.h

@@ -33 +33 @@
         switch(meshtype){
                 case Mesh2DhorizontalEnum:
+                        epart=xNew<int>(numberofelements);
+                        npart=xNew<int>(numberofnodes);
+                        index=xNew<int>(elements_width*numberofelements);
+                        for (i=0;i<numberofelements;i++){
+                                for (j=0;j<elements_width;j++){
+                                        *(index+elements_width*i+j)=(*(elements+elements_width*i+j))-1; //-1 for C indexing in Metis
+                                }
+                        }
+
+                        /*Partition using Metis:*/
+                        if (num_procs>1){
+#ifdef _HAVE_METIS_
+                                METIS_PartMeshNodalPatch(&numberofelements,&numberofnodes, index, &etype, &numflag, &num_procs, &edgecut, epart, npart);
+#else
+                                _error_("metis has not beed installed. Cannot run with more than 1 cpu");
+#endif
+                        }
+                        else if (num_procs==1){
+                                /*METIS does not know how to deal with one cpu only!*/
+                                for (i=0;i<numberofelements;i++) epart[i]=0;
+                                for (i=0;i<numberofnodes;i++)    npart[i]=0;
+                        }
+                        else _error_("At least one processor is required");
+                        break;
+                case Mesh2DverticalEnum:
                         epart=xNew<int>(numberofelements);
                         npart=xNew<int>(numberofnodes);