tupdate status text, fix second inversion call - cosmo - front and backend for Markov-Chain Monte Carlo inversion of cosmogenic nuclide concentrations
(HTM) git clone git://src.adamsgaard.dk/cosmo
(DIR) Log
(DIR) Files
(DIR) Refs
(DIR) README
(DIR) LICENSE
---
(DIR) commit 327d29df9476eb52b03f9c9b12269ab2e5316956
(DIR) parent 2d801ff619d5e873f2dcf9a1eb559ddc606d8445
(HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date: Tue, 29 Sep 2015 14:04:18 +0200
update status text, fix second inversion call
Diffstat:
M matlab/m_pakke2014maj11/MetHasLong… | 4 +---
M matlab/mcmc_inversion.m | 3 ++-
2 files changed, 3 insertions(+), 4 deletions(-)
---
(DIR) diff --git a/matlab/m_pakke2014maj11/MetHasLongstep4.m b/matlab/m_pakke2014maj11/MetHasLongstep4.m
t@@ -110,9 +110,7 @@ for it=1:N_run
%%% INSERTED BY ANDERS %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
fid = fopen(statusfile, 'w');
- statusinfo = strcat(...
- 'Elapsed: ', num2str(datestr(ElapsedTime,13)), ...
- ', Remaining: ', num2str(datestr(RemainingTime,13)));
+ statusinfo = strcat(num2str(datestr(RemainingTime,13), ' remaining'));
fprintf(fid, statusinfo);
fclose(fid);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
(DIR) diff --git a/matlab/mcmc_inversion.m b/matlab/mcmc_inversion.m
t@@ -236,7 +236,8 @@ for iwalk=1:fixed_stuff.Nwalkers
seed = fixed_stuff.WalkerSeeds(iwalk);
isBurnIn=0; %<<<<<<<<<<<<<<<<<< ------------------- must be changed when Sampling
[S.ms,S.accepts,S.Qs,S.Qds,S.lump_MetHas]=MetHasLongstep4(...
- mStartSampling,seed,isBurnIn,fixed_stuff);
+ mStartSampling,seed,isBurnIn,fixed_stuff, ...
+ statusfile); %%%% ADDED BY ANDERS
%<<<<< ... End Sampling the posterior
S.fs = fixed_stuff;
S.m_start = m_starts(:,iwalk);