tadd debug information - 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 77315fb555bef881b3b95155e01dd4a349a13d07
 (DIR) parent 9c54020ea1b3f6bf881e4987ffd8e2a44363884f
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Thu, 22 Oct 2015 15:54:43 +0200
       
       add debug information
       
       Diffstat:
         M matlab/generate_plots.m             |       7 ++++---
         M matlab/mcmc_inversion.m             |       6 ++++--
       
       2 files changed, 8 insertions(+), 5 deletions(-)
       ---
 (DIR) diff --git a/matlab/generate_plots.m b/matlab/generate_plots.m
       t@@ -178,11 +178,12 @@ for i1 = 1:M
        end
        
        
       -% Plot epsilon_rate_int in one histogram per walker
       +% Plot one histogram per model parameter per walker
        fh = [fh;figure('visible', show_figures)];
       -for i1 = 1:M
       -  for iwalk=1:min(4,Nwalkers)
       +for i1 = 1:M % for each model parameter
       +  for iwalk=1:min(4,Nwalkers) % for each walker
            isub = (i1-1)*4 + iwalk;
       +    i1, M, iwalk, Nwalkers
            subplot(4,Nwalkers,isub)
            Nhistc=histc(Ss{iwalk}.ms(i1,:),xbins{i1});
            bar(xbins{i1},Nhistc,'histc')
 (DIR) diff --git a/matlab/mcmc_inversion.m b/matlab/mcmc_inversion.m
       t@@ -228,8 +228,10 @@ fixed_stuff.StartTime = now; %This should allow the program to predict time of f
        % ANDERS: consider parfor for parallel computing. However, fixed_stuff and
        % S structures are incompatible with parfor.
        for iwalk=1:fixed_stuff.Nwalkers
       -
       -    disp(['\n#### iwalk = ' iwalk ' ####'])
       +    
       +    disp(' ')
       +    disp(['#### iwalk = ' num2str(iwalk) '/' ...
       +        num2str(fixed_stuff.Nwalkers) ' ####'])
        
            fixed_stuff.iwalk = iwalk; %Helps program keep user updated on progress.
            m_starts(:,iwalk) = WalkerStarter(iwalk,fixed_stuff);