Changeset 943545c in tol
- Timestamp:
- Mar 11, 2009, 5:58:55 PM (18 years ago)
- Branches:
- 1941-docker-image-to-test, devel, master
- Children:
- 91f3555
- Parents:
- 9e0ac9c
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
tol/stdlib/math/stat/models/bayesian/bysMcmc/_build.tol
r9e0ac9c r943545c 33 33 NameBlock _.sampler = modelSampler; 34 34 NameBlock _.notifier = notifier; 35 NameBlock _.config = config; 35 36 //Stores the sampler name 36 37 Text _.name = _.sampler::get.name(0)+"."+_.sampler::get.session(0); … … 41 42 //Timer member that can be externally consulted 42 43 Real _.time = 0; 43 Real _mcmc.maxTime = config::mcmc.maxTime;44 Real _mcmc.maxTime = _.config::mcmc.maxTime; 44 45 //Setting NameBlock description 45 46 Text _.autodoc.description="MCMC cycling builder for model "+_.name; … … 49 50 //-------------------------------------------------------------------------- 50 51 //The number of simulations cached in RAM before to disc writing 51 Real cacheLength_ = config::mcmc.cacheLength;52 Real cacheLength_ = _.config::mcmc.cacheLength; 52 53 //If resume is true and there was old simulations in _.path, new ones will 53 54 //be added at the end. In other case old simulations will be lost. 54 Real resume_ = config::do.resume;55 Real resume_ = _.config::do.resume; 55 56 //Bayesian inference configuration 56 57 Set bysInf_ = {BysInf.Report.Config 57 58 ( 58 Real raftery.diag.q = config::report.raftery.diag.q,59 Real raftery.diag.r = config::report.raftery.diag.r,60 Real raftery.diag.s = config::report.raftery.diag.s,61 Real raftery.diag.eps = config::report.raftery.diag.eps,62 Real acf.lag = config::report.acf.lag,63 Real histogram.parts = config::report.histogram.parts,64 Real kerDens.points = config::report.kerDens.points59 Real raftery.diag.q = _.config::report.raftery.diag.q, 60 Real raftery.diag.r = _.config::report.raftery.diag.r, 61 Real raftery.diag.s = _.config::report.raftery.diag.s, 62 Real raftery.diag.eps = _.config::report.raftery.diag.eps, 63 Real acf.lag = _.config::report.acf.lag, 64 Real histogram.parts = _.config::report.histogram.parts, 65 Real kerDens.points = _.config::report.kerDens.points 65 66 )}; 66 67 … … 177 178 //WriteLn("TRACE "+_MID+"numSim="<<numSim); 178 179 Real okMls = Case(Or(numSim>1, 179 config::mls.method==BysMcmc::Options::Mls.Method::None),180 _.config::mls.method==BysMcmc::Options::Mls.Method::None), 180 181 { 181 182 1 182 183 }, 183 config::mls.method==BysMcmc::Options::Mls.Method::FunnelGibbs,184 _.config::mls.method==BysMcmc::Options::Mls.Method::FunnelGibbs, 184 185 { 185 186 Real numIter = 1; 186 While(numIter<= config::mls.maxIter,187 While(numIter<=_.config::mls.maxIter, 187 188 { 188 189 VMatrix param = _.sampler::maxLikelihood(numIter); … … 194 195 { 195 196 WriteLn(_MID+"Unknown Maximum Likelihood Step method "+ 196 config::mls.method, "E");197 _.config::mls.method, "E"); 197 198 0 198 199 });
Note: See TracChangeset
for help on using the changeset viewer.
