CosmoCoffee

Xin Wang · Post by **Xin Wang** » May 07 2013

Hi, all,

Could anyone please comment on how the .chk files come into being? Recently I've hacked the cosmomc to add one module, which is capable of contributing a very large -ln(likelihood) term (say up to 10^8) to the total 0.5*chi2. As a consequence, checkpoint never works and meanwhile chains are far from converging, no matter how much time it takes or how many points they accumulate.

Thanks a million!

Post by **Antony Lewis** » May 14 2013

Checkpointing was accidentally broken in the April version. To fix, in MCMC.f90 move

Code: Select all

        if &#40;checkpoint_burn/=0&#41; checkpoint_burn = checkpoint_burn-1

down until outside the IF (thanks Silvia for reporting).

Xin Wang · Post by **Xin Wang** » May 16 2013

Thanks, Antony! You mean to revise the body of subroutine MoveDone into this?

Code: Select all

    if &#40;accpt&#41; then
        num_accept = num_accept + 1
        thin=1
        if &#40;present&#40;thin_fac&#41;&#41; thin=thin_fac
        want= num_accept> burn_in .and. checkpoint_burn ==0  .and. CurLike /= LogZero
        if &#40;want&#41; then
            acc = acc + mult
            indep_acc= indep_acc + mult
        end if
        if &#40;acc >= thin&#41; then
            if &#40;want&#41; then
                output_lines = output_lines +1
                call WriteParams&#40;CurParams,real&#40;acc/thin,mcp&#41;,CurLike&#41;
            end if
            acc = mod&#40;acc, thin&#41;
        end if
if &#40;checkpoint_burn/=0&#41; checkpoint_burn = checkpoint_burn-1		! checkpoint_burn fixed
        if &#40;indep_sample /= 0 .and. indep_acc >= indep_sample*thin .and. want&#41; then
            call WriteIndepSample&#40;CurParams, CurLike,real&#40;indep_acc/&#40;indep_sample*thin&#41;,mcp&#41;&#41;
            indep_acc = mod&#40;indep_acc, indep_sample*thin&#41;
        end if
        if &#40;CurLike < MaxLike&#41; MaxLike = CurLike
        mult=1
    else
        mult = mult + 1
    end if

Antony Lewis wrote:Checkpointing was accidentally broken in the April version. To fix, in MCMC.f90 move
Code: Select all
        if &#40;checkpoint_burn/=0&#41; checkpoint_burn = checkpoint_burn-1
down until outside the IF (thanks Silvia for reporting).

CosmoCoffee

On checkpoint functionality and chains' convergence

On checkpoint functionality and chains'' convergence

Re: On checkpoint functionality and chains'''''''' convergen

Re: On checkpoint functionality and chains'''''''''''''''' c