Changeset a9a5553 in tol
- Timestamp:
- Sep 22, 2009, 1:24:19 PM (17 years ago)
- Branches:
- 1941-docker-image-to-test, devel, master
- Children:
- 5c7ee7a
- Parents:
- cb115ec
- Location:
- tol_tests/tol/BSR
- Files:
-
- 5 edited
-
test_03/test.tol (modified) (1 diff)
-
test_12/test.tol (modified) (2 diffs)
-
test_13/test.tol (modified) (1 diff)
-
test_15/test.tol (modified) (4 diffs)
-
test_16/test.tol (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
tol_tests/tol/BSR/test_03/test.tol
rcb115ec ra9a5553 17 17 Real numErr0 = Copy(NError); 18 18 19 Real PutRandomSeed( 208774589);19 Real PutRandomSeed(0); 20 20 21 21 //Random seed setting -
tol_tests/tol/BSR/test_12/test.tol
rcb115ec ra9a5553 177 177 Code bsr.scalarSampler.nonLinBlock = 178 178 BysMcmc::Options::Scalar.Sampler::Slice; 179 //Metodo de aproximacion del arima 180 Text bsr.arimaFilter = BysMcmc::Options::Arima.Filter::FastCholSea; 179 181 //Report configuration 180 182 Real report.raftery.diag.q = 0.25; … … 186 188 Real report.kerDens.points = 0; 187 189 //Generic flags 188 Real do.resume = 0; 189 Real do.report = True; 190 Real do.eval = True; 191 Real do.linear.effects = True; 192 //Metodo de aproximacion del arima 193 Text bsr.arimaFilter = BysMcmc::Options::Arima.Filter::FastCholSea//; 194 //Real bsr.arimaSkipIter = 100 190 Real do.resume = 0; 191 Real do.report = True; 192 Real do.eval = True; 193 Real do.linear.effects = True; 194 Real do.mcmc.profiler = True 195 195 ]]; 196 197 196 198 197 BysMcmc::Cycler cycler = BysMcmc::Bsr::DynHlm::DBApi::GetCycler -
tol_tests/tol/BSR/test_13/test.tol
rcb115ec ra9a5553 156 156 157 157 Set bsr.parsed = BSR.ModelDef( 158 "joint", 158 159 parsed->DocInfo, 159 160 parsed->LinearBlock, -
tol_tests/tol/BSR/test_15/test.tol
rcb115ec ra9a5553 33 33 Text bsr.arimaFilter = BysMcmc::Options::Arima.Filter::FastChol; 34 34 Real bsr.arimaSkipIter = 1; 35 35 /* 36 36 Real bsr.linBlk.MH.useAfterIter = 10; 37 37 Real bsr.linBlk.MH.maxRejected = 3; 38 38 Real bsr.linBlk.MH.maxConsecutive = 100; 39 39 */ 40 40 Real bsr.linBlk.X.drop = 1.E-15; 41 41 Real bsr.linBlk.usePreconditioning = 0; … … 64 64 Set LinearBlk = Select(recover::_.colNames, Real(Text param) 65 65 { 66 If(!TextBeginWith(param,"LinearBlk::"),0,Rand(0,1)< .98)66 If(!TextBeginWith(param,"LinearBlk::"),0,Rand(0,1)<1) 67 67 }); 68 68 … … 77 77 }); 78 78 79 Set ArimaBlk = Select(recover::_.colNames, Real(Text param) 80 { 81 If(!TextBeginWith(param,"ArimaBlk::"),0,Rand(0,1)<1) 82 }); 83 79 84 Set fixedParam = Copy(Empty); 80 85 VMatrix fixedValues = Rand(0,1,0,0); 81 Set reloadParam = LinearBlk << MissingBlk << SigmaBlk ;86 Set reloadParam = LinearBlk << MissingBlk << SigmaBlk << ArimaBlk; 82 87 83 88 BysMcmc::McmcPartialHandler partialHandler = … … 111 116 112 117 Real cycler::_.sampler::_.arm.blk::refreshFrequency := 0; 113 114 118 //Real cycler::initialize(0); 115 116 119 BysMcmc::Estimation estim = BysMcmc::Estimation::New(cycler); 117 118 Real TolOprProfiler.Enabled := True;119 120 Real estim::Run(0); 120 Real TolOprProfiler.Dump("ses.01");121 122 121 Real FClose(recover::_.mcmc.handler); 123 122 -
tol_tests/tol/BSR/test_16/test.tol
rcb115ec ra9a5553 29 29 Real do.mcmc.profiler = True; 30 30 //Partial sampling 31 Real partialSampling.Sequential = False31 Real partialSampling.Sequential = True 32 32 ]]; 33 34 35 33 36 34 BysMcmc::McmcRecover recover = BysMcmc::McmcRecover::Load( … … 41 39 { 42 40 If(!TextBeginWith(param,"Node")| 43 TextEndAt(param,"_Sigma2"),0,Rand(0,1)< .98)41 TextEndAt(param,"_Sigma2"),0,Rand(0,1)<0.98) 44 42 }); 45 43 46 44 Set MissingBlk = Select(recover::_.colNames, Real(Text param) 47 45 { 48 If(!TextBeginWith(param,"MissInp."),0,Rand(0,1)< .98)46 If(!TextBeginWith(param,"MissInp."),0,Rand(0,1)<0.98) 49 47 }); 50 48 51 49 Set SigmaBlk = Select(recover::_.colNames, Real(Text param) 52 50 { 53 If(!TextEndAt(param,"_Sigma2"),0,Rand(0,1)<0. 5)51 If(!TextEndAt(param,"_Sigma2"),0,Rand(0,1)<0.98) 54 52 }); 53 54 Set ArimaBlk = Select(recover::_.colNames, Real(Text param) 55 { 56 If(!TextBeginWith(param,"ArmaBlk"),0,Rand(0,1)<.98) 57 }); 58 55 59 56 60 Set fixedParam = Copy(Empty); 57 61 VMatrix fixedValues = Rand(0,1,0,0); 58 Set reloadParam = LinearBlk << MissingBlk << SigmaBlk; 62 Set reloadParam = LinearBlk << MissingBlk << SigmaBlk;// << ArimaBlk; 59 63 60 64 BysMcmc::McmcPartialHandler partialHandler =
Note: See TracChangeset
for help on using the changeset viewer.
