Hi,
I have a few questions about reading data in parallel using libmesh. The
program is running on n processes, and I wish to have n partitions on this mesh
which means each data partition is assigned to one process.
1. How to make sure that each process just contains one part of the entire mesh
and the corresponding nodal values?
2. How can each process maintain the data values of ghost nodes since its local
elements may contains the nodes that are not in its local nodes?
Here is the initialization code. It seems that each process can access the
entire mesh and the entire nodal solution. Does it mean that each process hold
the entire data?
libMesh::UnstructuredMesh *_mesh;
libMesh::ExodusII_IO *_exio;
libMesh::EquationSystems *_eqsys;
libMesh::System *_asys;
libMesh::PointLocatorTree *_locator;
_mesh = new Mesh(comm());
_exio = new ExodusII_IO(*_mesh);
_exio->read(filename);
_mesh->allow_renumbering(false);
_mesh->prepare_for_use();
unsigned int _Vx_var, _Vy_var, _Vz_var;
_eqsys = new EquationSystems(*_mesh);
_asys = &(_eqsys->add_system<System>("Velocity"));
_Vx_var = _asys->add_variable("Vx", FIRST, LAGRANGE);
_Vy_var = _asys->add_variable("Vy", FIRST, LAGRANGE);
_Vz_var = _asys->add_variable("Vz", FIRST, LAGRANGE);
_eqsys->init();
_locator = new PointLocatorTree(*_mesh);
_exio->copy_nodal_solution(*_asys, "Vx", varnames[0]); // velocityX
_exio->copy_nodal_solution(*_asys, "Vy", varnames[1]); // velocityY
_exio->copy_nodal_solution(*_asys, "Vz", varnames[2]); // velocityZ
libMesh::UniquePtr<libMesh::NumericVector<libMesh::Number> > _as =
_asys->current_local_solution->clone();
Best regards,
------------------------------------------------------------------------------
Attend Shape: An AT&T Tech Expo July 15-16. Meet us at AT&T Park in San
Francisco, CA to explore cutting-edge tech and listen to tech luminaries
present their vision of the future. This family event has something for
everyone, including kids. Get more information and register today.
http://sdm.link/attshape
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
