Re: [Libmesh-users] mpi split: possible bug

2016-03-11 Thread Roy Stogner
On Fri, 11 Mar 2016, Manav Bhatia wrote: > While we are at it, I have one more observation: > > Consider the global comm with ranks {0, 1, 2, 3}. > > Then, after the following intended split: {0}, which requires rank > {0} to request split with a color and ranks {1, 2, 3} with > MPI_UNDEFINED, th

Re: [Libmesh-users] mpi split: possible bug

2016-03-11 Thread Manav Bhatia
Looks good to me, Roy. While we are at it, I have one more observation: Consider the global comm with ranks {0, 1, 2, 3}. Then, after the following intended split: {0}, which requires rank {0} to request split with a color and ranks {1, 2, 3} with MPI_UNDEFINED, the mpi_communicator on the

Re: [Libmesh-users] mpi split: possible bug

2016-03-11 Thread Roy Stogner
On Fri, 11 Mar 2016, Roy Stogner wrote: > On Fri, 11 Mar 2016, Manav Bhatia wrote: > >> From the current status of the code in parallel.h, it looks like >> mpi_comm_free is called only when _I_duped_it is true, which is >> true only when Communicator::duplicate() is called. >> >> I think the

Re: [Libmesh-users] mpi split: possible bug

2016-03-11 Thread Roy Stogner
On Fri, 11 Mar 2016, Manav Bhatia wrote: > From the current status of the code in parallel.h, it looks like > mpi_comm_free is called only when _I_duped_it is true, which is > true only when Communicator::duplicate() is called. > > I think the communicator created using split() also needs to

[Libmesh-users] mpi split: possible bug

2016-03-10 Thread Manav Bhatia
Hi, From the current status of the code in parallel.h, it looks like mpi_comm_free is called only when _I_duped_it is true, which is true only when Communicator::duplicate() is called. I think the communicator created using split() also needs to be destroyed using comm_free. But since _I