Dear Libmesh developers and users,

I have a simple code which runs perfectly on a single MPI process, but fails on more than one MPI process.

The code is as follows:

using namespace libMesh;

int main (int argc, char ** argv)
{
  {
    LibMeshInit init (argc, argv);

    Mesh mesh(init.comm());
    MeshTools::Generation::build_cube (mesh, 5, 5, 5);

    mesh.print_info();

    TecplotIO mesh_output(mesh);

    mesh_output.write("mesh.dat");
  }

  return 0;
}

The code is frozen at some MPI related calls inside mesh.prepare_for_use(). The libmesh was configured as:

./configure PETSC_DIR=/home/mpovolot/photonics/petsc/petsc-3.7.6 MPIHOME=/apps/conte/openmpi/1.10.1/gcc-5.2.0 PETSC_ARCH=arch-linux2-cxx-opt F77=mpif77 CC=mpicc GCC=mpicc CXX=mpiCC --disable-vtk --with-vtk-include=/depot/itap/tsaiwei/apps/vtk/6.3.0_test/include/vtk-6.3/ --with-vtk-lib=/depot/itap/tsaiwei/apps/vtk/6.3.0_test/lib --enable-tetgen --enable-triangle --enable-tecplot --disable-nemesis --disable-strict-lgpl --enable-parmesh --enable-amr --enable-shared=yes --disable-glibcxx-debugging --prefix=/home/mpovolot/photonics/libmesh/libmesh-1.2.0/built

Could you, please, tell me what is wrong with my code?

Thank you,

Michael.



------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to