tsmall improvements to mcmc_inversion, improve title in histogram plot - 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 054114ed3b91dee3ddb4b77c29a1a4b6c7271f4e
 (DIR) parent 88e977f479b2ecc13689e14219a207e0fdcac247
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Thu, 22 Oct 2015 15:34:30 +0200
       
       small improvements to mcmc_inversion, improve title in histogram plot
       
       Diffstat:
         M matlab/generate_plots.m             |       2 +-
         M matlab/mcmc_inversion.m             |       7 +++++--
       
       2 files changed, 6 insertions(+), 3 deletions(-)
       ---
 (DIR) diff --git a/matlab/generate_plots.m b/matlab/generate_plots.m
       t@@ -230,7 +230,7 @@ title(['Density cross-plots B. Result file =',save_file],'interp','none')
        figure(fh(4)); set(fh(4), 'Visible', show_figures)
        subplot(4,Nwalkers,2)
        %title(['Histograms. Result file =',save_file],'interp','none')
       -title('Distribution of model parameter values','interp','none')
       +title('Distribution of model parameter values')
        
        % position figure windows at certain coordinates on the screen
        if strcmp(show_figures, 'on')
 (DIR) diff --git a/matlab/mcmc_inversion.m b/matlab/mcmc_inversion.m
       t@@ -225,9 +225,12 @@ fs.Sampling = CompleteFsSampling(fs.Sampling);
        fixed_stuff = fs;
        
        fixed_stuff.StartTime = now; %This should allow the program to predict time of finish
       -% ANDERS: consider parfor for parallel computing
       +% ANDERS: consider parfor for parallel computing. However, fixed_stuff and
       +% S structures are incompatible with parfor.
        for iwalk=1:fixed_stuff.Nwalkers
       -    iwalk
       +
       +    disp(['\n#### iwalk = ' iwalk ' ####'])
       +
            fixed_stuff.iwalk = iwalk; %Helps program keep user updated on progress.
            m_starts(:,iwalk) = WalkerStarter(iwalk,fixed_stuff);
            d_starts(:,iwalk) = g(m_starts(:,iwalk),fixed_stuff);