CosmoMC: DoAbort and MPI

Use of Healpix, camb, CLASS, cosmomc, compilers, etc.
Post Reply
Loison Hoi
Posts: 23
Joined: February 21 2006
Affiliation: McGill University

CosmoMC: DoAbort and MPI

Post by Loison Hoi » March 13 2006

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

Antony Lewis
Posts: 1501
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

Re: CosmoMC: DoAbort and MPI

Post by Antony Lewis » March 13 2006

I agree- will delete it. (however presumably pretty harmless since the chains have finished anyway)

Post Reply