Changeset 6f124ed in tol

Timestamp:

Apr 2, 2009, 9:52:37 AM (18 years ago)

Author:

Víctor de Buen <vdebuen@…>

Branches:

1941-docker-image-to-test, devel, master

Children:

584d9b3

Parents:

b4134311

Message:

Cannot assign a system Code to a member as drawer that will be modified

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

File:

• tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/gibbs/_stdLinearBlock.tol (modified) (6 diffs)

tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/gibbs/_stdLinearBlock.tol

@@ -70 +70 @@
   Real _.betaLogDens = -1/0;
   //Stores the associated drawing method
-  Code drawer = UnknownCode;
+  Real noDraw(Real unused) { ? };
+  Code drawer = noDraw;
 
   ////////////////////////////////////////////////////////////////////////////
@@ -242 +243 @@
   ////////////////////////////////////////////////////////////////////////////
   {
+  //WriteLn("TRACE "+_MID+" draw.unconstrained 1");
     VMatrix z = Gaussian(_.N,1,0,1);
     Real _.betaLogDens := -0.5*VMatDat(MtMSqr(z),1,1);
@@ -254 +256 @@
   ////////////////////////////////////////////////////////////////////////////
   {
-  //WriteLn("TRACE "+_MID+" draw.constrained 1.1 Max(A*b0-a)="<<VMatMax(_.A*_.blockDef::_.store - a));
+  //WriteLn("TRACE "+_MID+" draw.constrained 1");
     VMatrix d   = _.a - _.A*nu;
+  //WriteLn("TRACE "+_MID+" draw.constrained 2");
     VMatrix If(changed.X, _.D  := Tra(CholeskiSolve(_.L,_.At,"PtL")));
+  //WriteLn("TRACE "+_MID+" draw.constrained 3");
     VMatrix z0  = Tra(Tra(_.blockDef::_.store-nu)*_.Ls);
-  //WriteLn("TRACE "+_MID+" draw.constrained 1.2 Max(D*z0-d)="<<VMatMax(_.D*z0 - d));
+  //WriteLn("TRACE "+_MID+" draw.constrained 4");
     VMatrix z = TruncStdGaussian(_.D,d,z0,1,
                  _.blockDef::_.config::bsr.truncMNormal.gibbsNumIter);
+  //WriteLn("TRACE "+_MID+" draw.constrained 5");
     Real distBorder = VMatMax(_.D*z-d);
+  //WriteLn("TRACE "+_MID+" draw.constrained 6");
     Real If(distBorder<=0,
     {
+    //WriteLn("TRACE "+_MID+" draw.constrained 7");
       VMatrix _.beta0 := _.blockDef::_.store;
       1
@@ -272 +279 @@
       If(!VRows(_.beta0),0,
       {
+      //WriteLn("TRACE "+_MID+" draw.constrained 8");
         VMatrix z0 := Tra(Tra(_.beta0-nu)*_.Ls);
+      //WriteLn("TRACE "+_MID+" draw.constrained 9");
         VMatrix z  := TruncStdGaussian(_.D,d,z0,1,
                      _.blockDef::_.config::bsr.truncMNormal.gibbsNumIter);
@@ -279 +288 @@
       1
     });
+  //WriteLn("TRACE "+_MID+" draw.constrained 10");
     Real _.betaLogDens := -0.5*VMatDat(MtMSqr(z),1,1);
+  //WriteLn("TRACE "+_MID+" draw.constrained 11");
     VMatrix  (_.blockDef::_.store) := nu + CholeskiSolve(_.L,z,"LtP");
+  //WriteLn("TRACE "+_MID+" draw.constrained 12");
     Copy(_.blockDef::_.store)
   };
@@ -290 +302 @@
   ////////////////////////////////////////////////////////////////////////////
   {
+  //WriteLn("TRACE "+_MID+" draw 1");
     If(_.blockDef::enabled,
     {
+    //WriteLn("TRACE "+_MID+" draw 2");
       VMatrix aux = drawer(numSim, _.solution::nu);
-    //WriteLn("TRACE "+_MID+" draw");
+    //WriteLn("TRACE "+_MID+" draw 3");
       Real _.blockDef::checkNonFiniteValues(numSim);
       aux