Changeset 584d9b3 in tol


Ignore:
Timestamp:
Apr 2, 2009, 9:53:05 AM (18 years ago)
Author:
Víctor de Buen <vdebuen@…>
Branches:
1941-docker-image-to-test, devel, master
Children:
5937c50
Parents:
6f124ed
Message:

Just traces

git-svn-id: ​https://www.tol-project.org/svn/tolp/trunk@854 be4846ae-f051-483a-831f-b64e7d08527f

Location:
tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/gibbs
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/gibbs/_basicMaster.tol

    r6f124ed r584d9b3  
    232232  //WriteLn("TRACE"+_MID+" draw.sig.lin 1.1");
    233233    Real DoExternalFilter(0);
     234  //WriteLn("TRACE"+_MID+" draw.sig.lin 1.2");
    234235    VMatrix om = If(Not(_.outMis.blk::_.blockDef::_.numCol),Rand(0,1,0,0),
    235236    {
    … …  
    244245                         _.beta, _.YplusV-_.modelDef->X*_.beta)
    245246    });
     247  //WriteLn("TRACE"+_MID+" draw.sig.lin 1.4");
    246248    VMatrix _.XplusU := _.inpMis.blk::filter(_.modelDef->X);
    247249    VMatrix _.filter := _.XplusU*_.beta;
    … …  
    259261    VMatrix si = _.sig.blk::_.sigmaInv;
    260262  //WriteLn("TRACE"+_MID+" draw.sig.lin 7");
    261 /*
    262     VMatrix Y01 = _.arm.blk::filter(si*_.YplusV);
    263   //WriteLn("TRACE"+_MID+" draw.sig.lin 6");
    264     VMatrix X01 = _.arm.blk::filter(si*_.XplusU);
    265   //WriteLn("TRACE"+_MID+" draw.sig.lin 8");
    266     VMatrix _.beta := If(Not(_.mln.blk::_.N),Rand(0,1,0,0),
    267       _.mln.blk::draw(numSim, Y01, X01));
    268 */
     263  //WriteLn("TRACE"+_MID+" draw.sig.lin 9");
    269264    VMatrix _.beta := If(Not(_.mln.blk::_.N),Rand(0,1,0,0),
    270265      _.mln.blk::draw(numSim, si*_.YplusV, si*_.XplusU, _.arm.blk));
    271   //WriteLn("TRACE"+_MID+" draw.sig.lin 9");
     266  //WriteLn("TRACE"+_MID+" draw.sig.lin 10");
    272267    _.beta << sig << arm << im << om
    273268  };

  • tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/gibbs/_mainLinearBlock.tol

    r6f124ed r584d9b3  
    7878  //WriteLn("TRACE "+_MID+" draw.gibbs 1");
    7979    Real _.MH.consecutive := 0;
     80  //WriteLn("TRACE "+_MID+" draw.gibbs 2");
    8081    VMatrix Y01 = arm.blk::filter(Y);
     82  //WriteLn("TRACE "+_MID+" draw.gibbs 3");
    8183    VMatrix X01 = arm.blk::filter(X);
     84  //WriteLn("TRACE "+_MID+" draw.gibbs 4");
    8285    Real _.lin.blk::solve(numSim, Y01, X01);
     86  //WriteLn("TRACE "+_MID+" draw.gibbs 5");
    8387    VMatrix b =_.lin.blk::draw(numSim);
     88  //WriteLn("TRACE "+_MID+" draw.gibbs 6");
    8489    VMatrix ep = Y01-X01*b;
     90  //WriteLn("TRACE "+_MID+" draw.gibbs 7");
    8591    Real _.betaLogDens := -0.5*VMatDat(MtMSqr(ep),1,1);
    8692    b
    … …  
    162168  ////////////////////////////////////////////////////////////////////////////
    163169  {
     170  //WriteLn("TRACE "+_MID+" draw 1");
    164171    If(Not(_.blockDef::enabled), _.blockDef::_.store,
    165172    {
     173    //WriteLn("TRACE "+_MID+" draw 2");
    166174      Real _.calculatedSimulation := _.calculatedSimulation+1;
     175    //WriteLn("TRACE "+_MID+" draw 3");
    167176      If(And(_.calculatedSimulation >_.MH.useAfterIter,
    168177             _.MH.consecutive <= _.MH.maxConsecutive),

  • tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/gibbs/_outputMissingBlock.tol

    r6f124ed r584d9b3  
    204204  {
    205205    Real _.MH.consecutive := 0;
     206  //WriteLn("TRACE "+_MID+" draw.gibbs 1");
    206207    Real _.lin.blk::solve(numSim, ls::Y01, ls::X01);
     208  //WriteLn("TRACE "+_MID+" draw.gibbs 2");
    207209    VMatrix _.blockDef::_.store := _.lin.blk::draw(numSim);
     210  //WriteLn("TRACE "+_MID+" draw.gibbs 3");
    208211    VMatrix e = ls::Y01-ls::X01*_.blockDef::_.store;
     212  //WriteLn("TRACE "+_MID+" draw.gibbs 4");
    209213    Real _.betaLogDens := MatSum(Log(f01(VMat2Mat(e))));
     214  //WriteLn("TRACE "+_MID+" draw.gibbs 5");
    210215    _.blockDef::_.store
    211216  };
    … …  
    290295  ////////////////////////////////////////////////////////////////////////////
    291296  {
     297  //WriteLn("TRACE "+_MID+" draw 1");
    292298    If(Not(_.blockDef::enabled), _.blockDef::_.store,
    293299    {
     300    //WriteLn("TRACE "+_MID+" draw 2");
    294301      NameBlock ls = buildLinearSystem(arima,si,Z);
     302    //WriteLn("TRACE "+_MID+" draw 3");
    295303      Real _.calculatedSimulation := _.calculatedSimulation+1;
     304    //WriteLn("TRACE "+_MID+" draw 4");
    296305      If(And(_.calculatedSimulation >_.MH.useAfterIter,
    297306             _.MH.consecutive <= _.MH.maxConsecutive),
Note: See TracChangeset for help on using the changeset viewer.