The errors do not occur at all with the old Metis 4.0.3. So there must be some incompatibility between Getfem and Metis 5
- Torquil On 18 November 2013 15:33, Torquil Macdonald Sørensen <[email protected]>wrote: > Ok, I shall take a look at the available MPI-related functions. Thanks for > the pointer. > > When I try the method mesh.get_mpi_region() in one of my simple test > programs, or when I run tests/elastostatic.cc, with mpiexec -n N, for N>=2, > I get error messages: > > "Input Error: Incorrect ncuts." > > E.g. I get it from this test program (from the line that calls > get_mpi_region): > > #include "getfem/getfem_mesh.h" > #include "getfem/getfem_regular_meshes.h" > #include "iostream" > int main() > { > MPI_Init(0,0); > getfem::mesh mesh; > getfem::regular_mesh(mesh, "GT='GT_PK(2,1)'; NSUBDIV=[4,4]; ORG=[0,0]; > SIZES=[1,1]; NOISED=0"); > std::cout << mesh.get_mpi_region() << std::endl; > MPI_Finalize(); > return 0; > } > > Running with mpiexec -n 1 is fine. The program then just outputs > "ALL_CONVEXES". > > My Metis is version 5.1.0. Something similar was discussed here: > > http://glaros.dtc.umn.edu/gkhome/node/853 > > Maybe Getfem needs the older Metis 4? > > - Torquil > > > > On 18 November 2013 13:41, Yves Renard <[email protected]> wrote: > >> >> Dear Torquil, >> >> With --enable-paralevel=2, metis is called to give a partition of each >> mesh. >> The model bricks automatically make the assembly only on the part of the >> mesh corresponding to the process number. >> This is done by a call to the following method of getfem_mesh: >> mesh().intersect_with_mpi_region(rg) >> >> But this does not affect the numbering of nodes. In particular, each >> process have a complete mesh. >> >> Yves. >> >> >> >> >> >> Le 18/11/2013 12:20, Torquil Macdonald Sørensen a écrit : >> >> Hi! >> >> Should I see a difference in the mesh node enumeration, after having >> compiled Getfem with --enable-paralevel=2 and Metis support, and my >> program with -DGETFEM_PARA_LEVEL=2 -DGMM_USES_MPI=1? >> >> According to the documentation, GetFEM should then use Metis to >> "implicitly partition" the mesh. >> >> Does that mean that the Metis based partitioning goes on behind the >> scenes and no changes would be seen in mesh.points() when switching >> between "mpiexec -n 1 ./prog" and "mpiexec -n 4 ./prog" >> >> Best regards >> Torquil Sørensen >> >> >> >> _______________________________________________ >> Getfem-users mailing >> [email protected]https://mail.gna.org/listinfo/getfem-users >> >> >> >> -- >> >> Yves Renard ([email protected]) tel : (33) 04.72.43.87.08 >> Pole de Mathematiques, INSA-Lyon fax : (33) 04.72.43.85.29 >> 20, rue Albert Einstein >> 69621 Villeurbanne Cedex, FRANCE >> http://math.univ-lyon1.fr/~renard >> >> --------- >> >> >> _______________________________________________ >> Getfem-users mailing list >> [email protected] >> https://mail.gna.org/listinfo/getfem-users >> >> >
_______________________________________________ Getfem-users mailing list [email protected] https://mail.gna.org/listinfo/getfem-users
