CosmoMC: seg fault for action=2

Use of Healpix, camb, CLASS, cosmomc, compilers, etc.
Post Reply
Graeme Addison
Posts: 27
Joined: July 17 2014
Affiliation: Johns Hopkins University

CosmoMC: seg fault for action=2

Post by Graeme Addison » August 09 2019

I'm hitting a seg fault when I try running action=2 with the new July 2019 CosmoMC (& new Planck likelihood although it doesn't seem to matter which likelihoods are used). action=0 and action=4 work fine as far as I can tell, so maybe specific to some routine used in the minimizing process? Could I need some extra compiler option etc. that doesn't show up except for action=2?

Pasted last few lines of job output below:

At the return from BOBYQA Number of function values = 156
Least value of F = 3.166520941837106D+02 The corresponding X is:
1.876704D+00 -2.938268D+00 -1.342395D+00 -6.037330D+01 -3.960370D+01
-1.874395D+00 -3.010149D+00
4 Refining minimimum using low temp MCMC
4 Current logLike: 316.652094183711
4 Minimize MCMC with temp 1.000000000000000E-002
3 Stopping as have 140 samples.
slow changes 236 power changes 4
3 MCMC MaxLike = 316.627174630316
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image PC Routine Line Source
cosmomc 000000000069DBBD for__signal_handl Unknown Unknown
libpthread-2.17.s 00002BA777B3D6D0 Unknown Unknown Unknown
cosmomc 00000000006884EF Unknown Unknown Unknown
cosmomc 000000000068938F Unknown Unknown Unknown
cosmomc 000000000068A5AC for_alloc_assign_ Unknown Unknown
cosmomc 00000000005AD1EC Unknown Unknown Unknown
cosmomc 00000000004A1E34 Unknown Unknown Unknown
cosmomc 000000000049E55D Unknown Unknown Unknown
cosmomc 0000000000544CC8 Unknown Unknown Unknown
cosmomc 000000000048D3DF Unknown Unknown Unknown
cosmomc 000000000054A91A Unknown Unknown Unknown
cosmomc 00000000005487D9 Unknown Unknown Unknown
cosmomc 000000000046F8E6 Unknown Unknown Unknown
cosmomc 0000000000553882 Unknown Unknown Unknown
cosmomc 000000000040EEDE Unknown Unknown Unknown
libc-2.17.so 00002BA777D6C445 __libc_start_main Unknown Unknown
cosmomc 000000000040EDE9 Unknown Unknown Unknown

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

Re: CosmoMC: seg fault for action=2

Post by Antony Lewis » August 09 2019

Does cosmomc_debug give something more helpful?

Graeme Addison
Posts: 27
Joined: July 17 2014
Affiliation: Johns Hopkins University

Re: CosmoMC: seg fault for action=2

Post by Graeme Addison » August 09 2019

Output from repeating run with cosmomc_debug:

forrtl: error (73): floating divide by zero
Image PC Routine Line Source
cosmomc_debug 0000000000CAAA8E Unknown Unknown Unknown
libpthread-2.17.s 00002AD02542F6D0 Unknown Unknown Unknown
cosmomc_debug 00000000004CE9C4 powell_constraine 1000 PowellConstrainedMinimize.f90
cosmomc_debug 00000000004C029F powell_constraine 462 PowellConstrainedMinimize.f90
cosmomc_debug 00000000004B69F8 powell_constraine 166 PowellConstrainedMinimize.f90
cosmomc_debug 000000000052C558 minimize_mp_tpowe 253 minimize.f90
cosmomc_debug 0000000000533086 minimize_mp_tpowe 309 minimize.f90
cosmomc_debug 000000000081CDFE MAIN__ 226 driver.F90
cosmomc_debug 00000000004102DE Unknown Unknown Unknown
libc-2.17.so 00002AD025862445 __libc_start_main Unknown Unknown
cosmomc_debug 00000000004101E9 Unknown Unknown Unknown

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

Re: CosmoMC: seg fault for action=2

Post by Antony Lewis » August 10 2019

hmm, seems to be in the minimizer code. That code should not have changed significantly, are you sure the exact example worked with the previous version? If so, please look at the code diff with the version that works for you and let me know if you see anything that changed around the div-by-zero line that could be causing it.

Post Reply