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 found calling prepare_for_use on its own caused a crash, due to a 
Parallel method. So no, I don't think I am.

Cody - yes, each mesh is entirely separate. Unfortunately I don't think 
multiple jobs will do it, because (and I left this out) the serial 
meshes are nested in a "global" mesh, so they need to talk to each other.

On 09/02/15 18:23, Roy Stogner wrote:
>
> 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 size
> 1) Communicator?
> ---
> Roy

-- 
Peter Wijeratne
Research Associate
Centre for Medical Image Computing,
c/o Department of Medical Physics and Bioengineering,
Engineering Front Building,
Malet Place,
University College London,
London, UK
WC1E 6BT

Tel: 020 7679 0177


------------------------------------------------------------------------------
Dive into the World of Parallel Programming. The Go Parallel Website,
sponsored by Intel and developed in partnership with Slashdot Media, is your
hub for all things parallel software development, from weekly thought
leadership blogs to news, videos, case studies, tutorials and more. Take a
look and join the conversation now. http://goparallel.sourceforge.net/
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to