28 #ifndef _TOMOGRAPHER_MHRWVALUEERRORBINSCONVERGEDCONTROLLER_H 29 #define _TOMOGRAPHER_MHRWVALUEERRORBINSCONVERGEDCONTROLLER_H 66 template<
typename ValueHistogramWithBinningMHRWStatsCollectorType_,
67 typename IterCountIntType_,
68 typename BaseLoggerType_>
76 typedef ValueHistogramWithBinningMHRWStatsCollectorType_
77 ValueHistogramWithBinningMHRWStatsCollectorType;
79 typedef IterCountIntType_ IterCountIntType;
81 typedef BaseLoggerType_ BaseLoggerType;
84 const ValueHistogramWithBinningMHRWStatsCollectorType & value_stats_collector;
86 const IterCountIntType check_frequency_sweeps;
88 IterCountIntType last_forbidden_iter_number;
90 const Eigen::Index max_allowed_unknown;
91 const Eigen::Index max_allowed_unknown_notisolated;
92 const Eigen::Index max_allowed_not_converged;
94 const double max_add_run_iters;
105 const ValueHistogramWithBinningMHRWStatsCollectorType & value_stats_collector_,
106 BaseLoggerType & baselogger_,
107 IterCountIntType check_frequency_sweeps_ = 1024,
108 Eigen::Index max_allowed_unknown_ = 0,
109 Eigen::Index max_allowed_unknown_notisolated_ = 0,
110 Eigen::Index max_allowed_not_converged_ = 0,
111 double max_add_run_iters_ = 1.5
113 : value_stats_collector(value_stats_collector_),
114 check_frequency_sweeps(check_frequency_sweeps_),
115 last_forbidden_iter_number(0),
116 max_allowed_unknown(max_allowed_unknown_),
117 max_allowed_unknown_notisolated(max_allowed_unknown_notisolated_),
118 max_allowed_not_converged(max_allowed_not_converged_),
119 max_add_run_iters(max_add_run_iters_),
120 llogger(
"Tomographer::MHRWValueErrorBinsConvergedAdjuster", baselogger_)
124 template<
typename MHRWParamsType,
typename MHWalker,
typename MHRandomWalkType>
125 inline void init(MHRWParamsType & ,
const MHWalker & ,
126 const MHRandomWalkType & )
const 130 template<
typename MHRWParamsType,
typename MHWalker,
typename CountIntType,
typename MHRandomWalkType>
131 bool allowDoneThermalization(
const MHRWParamsType & ,
const MHWalker & ,
132 CountIntType ,
const MHRandomWalkType & )
const 137 template<
typename MHRWParamsType,
typename MHWalker,
typename MHRandomWalkType>
138 bool allowDoneRuns(
const MHRWParamsType & params,
const MHWalker & ,
139 IterCountIntType iter_k,
const MHRandomWalkType & )
142 if (check_frequency_sweeps == 0) {
149 if (last_forbidden_iter_number > 0 &&
150 (iter_k-last_forbidden_iter_number) < params.n_sweep*check_frequency_sweeps) {
156 stream <<
"params=" << params <<
", iter_k=" << iter_k
157 <<
", max_add_run_iters=" << max_add_run_iters;
162 if (max_add_run_iters > 0 && iter_k > max_add_run_iters * params.n_run * params.n_sweep) {
164 stream <<
"Ending random walk after reaching maximum sweep number " 165 << iter_k/params.n_sweep <<
" (" 166 << 100.0*iter_k/(params.n_sweep*params.n_run) <<
"% of set run length)";
174 const auto& binning_analysis = value_stats_collector.getBinningAnalysis();
175 const auto binmeans = value_stats_collector.binMeans();
177 const auto error_levels = binning_analysis.calcErrorLevels(binmeans);
178 const auto conv_status = binning_analysis.determineErrorConvergence(error_levels);
182 logger.longdebug([&](
std::ostream & stream) { stream <<
"Convergence summary = " << conv_summary; }) ;
184 if (conv_summary.n_not_converged > max_allowed_not_converged ||
185 conv_summary.n_unknown > max_allowed_unknown ||
186 (conv_summary.n_unknown-conv_summary.n_unknown_isolated) > max_allowed_unknown_notisolated) {
195 template<
typename MHRWParamsType,
typename MHWalker,
typename MHRandomWalkType>
196 inline void thermalizingDone(
const MHRWParamsType & ,
const MHWalker & ,
197 const MHRandomWalkType & )
const 201 template<
typename MHRWParamsType,
typename MHWalker,
typename MHRandomWalkType>
202 inline void done(MHRWParamsType & ,
const MHWalker & ,
203 const MHRandomWalkType & )
const 214 template<
typename IterCountIntType_ = int,
216 typename ValueHistogramWithBinningMHRWStatsCollectorType_ = void,
217 typename BaseLoggerType_ =
void>
222 const ValueHistogramWithBinningMHRWStatsCollectorType_ & value_stats_collector_,
223 BaseLoggerType_ & baselogger_,
224 IterCountIntType_ check_frequency_sweeps_ = 1024,
225 Eigen::Index max_allowed_unknown_ = 0,
226 Eigen::Index max_allowed_unknown_notisolated_ = 0,
227 Eigen::Index max_allowed_not_converged_ = 0,
228 double max_add_run_iters = 1.5
234 value_stats_collector_,
236 check_frequency_sweeps_,
237 max_allowed_unknown_,
238 max_allowed_unknown_notisolated_,
239 max_allowed_not_converged_,
Utilities for formatting strings.
static BinningErrorBarConvergenceSummary fromConvergedStatus(const Eigen::Ref< const Eigen::ArrayXi > &converged_status)
Construct a summary object from a list of converged status obtained from BinningAnalysis::determineEr...
A MHRWController Interface which stops the random walk after enough samples have been taken to make a...
Base namespace for the Tomographer project.
#define TOMO_ORIGIN
Use this as argument for a Tomographer::Logger::LocalLogger constructor .
LocalLogger< LocalLogger< BaseLoggerType > > subLogger(std::string new_prefix)
Create a sub-logger.
Never adjust the parameters of the random walk.
MHRWValueErrorBinsConvergedController(const ValueHistogramWithBinningMHRWStatsCollectorType &value_stats_collector_, BaseLoggerType &baselogger_, IterCountIntType check_frequency_sweeps_=1024, Eigen::Index max_allowed_unknown_=0, Eigen::Index max_allowed_unknown_notisolated_=0, Eigen::Index max_allowed_not_converged_=0, double max_add_run_iters_=1.5)
Constructor.
Some C++ utilities, with a tad of C++11 tricks.
Managing the need for specific overrides to operator new() for some types (especially Eigen types) ...
MHRWValueErrorBinsConvergedController< ValueHistogramWithBinningMHRWStatsCollectorType_, IterCountIntType_, BaseLoggerType_ > mkMHRWValueErrorBinsConvergedController(const ValueHistogramWithBinningMHRWStatsCollectorType_ &value_stats_collector_, BaseLoggerType_ &baselogger_, IterCountIntType_ check_frequency_sweeps_=1024, Eigen::Index max_allowed_unknown_=0, Eigen::Index max_allowed_unknown_notisolated_=0, Eigen::Index max_allowed_not_converged_=0, double max_add_run_iters=1.5)
Convenience function to create a MHRWValueErrorBinsConvergedController (using template argument deduc...
Definitions for MHRWStatsCollector Interface's.
Routines for performing a Metropolis-Hastings random walk.
Utilities for logging messages.