Re: [Libmesh-users] Multiple meshes in parallel

2015-02-17 Thread Peter Wijeratne
Thanks very much Roy, that's done the trick. Peter On 16/02/15 17:30, Roy Stogner wrote: > > On Mon, 16 Feb 2015, Peter Wijeratne wrote: > >> Still having the same problem. I now initialise using the >> Communicator from LibMeshInit::comm(); is this correct, or should I >> be creating my own Commu

Re: [Libmesh-users] Multiple meshes in parallel

2015-02-16 Thread Roy Stogner
On Mon, 16 Feb 2015, Peter Wijeratne wrote: > Still having the same problem. I now initialise using the > Communicator from LibMeshInit::comm(); is this correct, or should I > be creating my own Communicators? No, you'll need to create your own subcommunicators. The communicator from LibMeshIni

Re: [Libmesh-users] Multiple meshes in parallel

2015-02-16 Thread Peter Wijeratne
Still having the same problem. I now initialise using the Communicator from LibMeshInit::comm(); is this correct, or should I be creating my own Communicators? Thanks, Peter On 09/02/15 20:00, Derek Gaston wrote: > If you're trying to solve a different problem on each processor you'll > want t

Re: [Libmesh-users] Multiple meshes in parallel

2015-02-09 Thread Derek Gaston
If you're trying to solve a different problem on each processor you'll want to create an MPI sub-communicator for each processor and then initialize the Mesh with that (Note the Mesh constructor that takes a Parallel::Communicator: http://libmesh.github.io/doxygen/classlibMesh_1_1Mesh.html ). I be

Re: [Libmesh-users] Multiple meshes in parallel

2015-02-09 Thread Peter Wijeratne
Thanks for the quick replies! Roy - I actually do something a bit messy: SerialMesh* rmsh = new SerialMesh(1); GmshIO(*rmsh).read(filename.str().c_str()); rmsh->skip_partitioning(true); rmsh->prepare_for_use(); rmsh->skip_partitioning(false); rmsh->partition(1); I f

Re: [Libmesh-users] Multiple meshes in parallel

2015-02-09 Thread Roy Stogner
On Mon, 9 Feb 2015, Peter Wijeratne wrote: > Bit of a special (as opposed to general) problem here. I'm trying to > create a serial mesh on each of N processors, each to be kept on its > respective processor and solved separately. Are you initializing each Mesh with an appropriate (in this case

Re: [Libmesh-users] Multiple meshes in parallel

2015-02-09 Thread Cody Permann
There's an assumption that all the processors are working on the same problem so each processor gets a part of the mesh. Are each of the meshes truly independent so they don't share any DOFs? If they really are independent have you considered just running thousands of jobs instead? Cody On Mon

[Libmesh-users] Multiple meshes in parallel

2015-02-09 Thread Peter Wijeratne
Hello libMesh community, Bit of a special (as opposed to general) problem here. I'm trying to create a serial mesh on each of N processors, each to be kept on its respective processor and solved separately. This causes problems when initialising a given mesh's EquationSystems; I've narrowed the