On Wed, Jun 13, 2018 at 10:14 AM, Michael Povolotskyi <mpovo...@purdue.edu>
wrote:
> Thank you John,
>
> I do have some MPI code before LibMeshInit in my real application.
>
> Also, If I do not call MPI_Finalize at the very end I'm getting warning
> messages from MPI.
>
> Thank you for you explanation with .pvtu, I have changed my code to
> libMesh::VTKIO out(mesh); out.write("mesh1.pvtu"); and this worked.
>
> Do you have any suggestions for the Problem #2 that I reported?
>
> It seems to me that the problem is that the communicator of the Mesh
> object is smaller than the communicator of the VTK. This causes a problem
> for me.
>
Hmm, I don't think we have considered this use case in the past, the
vtkMPIController is initialized in libmesh.C, and it must be assuming
MPI_COMM_WORLD...
_vtk_mpi_controller = vtkMPIController::New();
_vtk_mpi_controller->Initialize(&argc, const_cast<char ***>(&argv),
/*initialized_externally=*/1);
_vtk_mpi_controller->SetGlobalController(_vtk_mpi_controller);
In order for the parallel VTK writer to work, the Mesh probably therefore
also needs to be using MPI_COMM_WORLD.
--
John
------------------------------------------------------------------------------
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
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users
