Page 1 of 1

CosmoMC: DoAbort and MPI

Posted: March 13 2006
by Loison Hoi
Hi,

I have a problem about the DoAbort statement at the end of subroutine MCMCsample.

If we use MPI, and the chains do not converge (e.g. sample = 20 to test whether MPI runs well), then cosmomc will hit the DoAbort statement at the end of MCMCsample, and returns an error like:

[0] MPI Abort by user Aborting program !
[0] Aborting program!
p0_10455: p4_error: : 4

The DoAbort statement in MCMCsample is right before the DoStop statement in the driver.F90, which just writes down the running time and the number of slow proposals. Therefore, the chains are finished normally and I think the DoAbort needs not to be there. A naive reason is that the DoAbort statement is inside the if-statement, so setting Feedback=0 can avoid the error. Anyway, I haven't tried a large number of sample to see if the error message still comes up.

I am using July 2005 CosmoMC.

Thanks.

Loison Hoi
12 Mar 2006

Re: CosmoMC: DoAbort and MPI

Posted: March 13 2006
by Antony Lewis
I agree- will delete it. (however presumably pretty harmless since the chains have finished anyway)