Changeset 0d654d6 in tol

Timestamp:

Dec 17, 2009, 3:04:08 PM (17 years ago)

Author:

Víctor de Buen <vdebuen@…>

Branches:

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

Children:

ff5d8e2

Parents:

57052ff

Message:

Just traces

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

File:

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

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

@@ -68 +68 @@
   VMatrix _.D  = SetDiag(MatSet(Tra(_.S))[1]);
   Matrix _one = Rand(1,1,1,1);
+  Matrix _zero = _.Y*0;
   Matrix _emptyMatrix = Rand(0,0,0,0);
 
@@ -104 +105 @@
   ///////////////////////////////////////////////////////////////////////////
   {
+  //WriteLn("TRACE"+_MID+"eval 1");
+  //WriteLn("TRACE"+_MID+"eval _.externalInfo[3]=\n"<<Matrix VMat2Mat(_.externalInfo[3],1));
+  //WriteLn("TRACE"+_MID+"eval _.externalInfo[4]=\n"<<Matrix SetRow(_.externalInfo[4]));
     Matrix Xb = VMat2Mat(SubRow(_.externalInfo[3],_.externalInfo[4]));
-    _.Y-RandTruncatedMultNormal(Xb, _one, _.A, _.B)
+  //WriteLn("TRACE"+_MID+"eval Xb=\n"<<Matrix Tra(Xb));
+    Matrix e = RandTruncatedMultNormal(Xb, _one, _.A, _.B);
+  //WriteLn("TRACE"+_MID+"eval e=\n"<<Matrix Tra(e));
+    _.Y-e
   }