Changeset 5e3cfca in tol

Timestamp:

Feb 9, 2010, 1:47:47 PM (17 years ago)

Author:

Víctor de Buen <vdebuen@…>

Branches:

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

Children:

1134230

Parents:

681fb68

Message:

Just traces

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

File:

• tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/_oneNode.tol (modified) (12 diffs)

tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/_oneNode.tol

@@ -83 +83 @@
   {
     Text _MID = "[BysMcmc::Bsr::OneNode::GetCycler]";
-  //WriteLn(_MID+" TRACE 1");
+  //WriteLn(" TRACE "+_MID+" 1");
     //Output matrix
     VMatrix Y = data::_.Y;
-  //WriteLn(_MID+" TRACE 2");
+  //WriteLn("TRACE "+_MID+" 2");
     //Input matrix
     VMatrix X = data::_.X;
 
-  //WriteLn(_MID+" TRACE 3");
+  //WriteLn("TRACE "+_MID+" 3");
     //Observational data length
     Real m = VRows(Y);
@@ -106 +106 @@
     Matrix  domain.A.triplet = Rand(0,3,0,0);
     
-  //WriteLn(_MID+" TRACE 4");
+  //WriteLn("TRACE "+_MID+" 4");
     //Builds linear block and collects prior and domain information
     Set linearBlock = For(1,n, Set(Real j)
@@ -112 +112 @@
       Set pi = (data::_.linearParamInfo)[j];
       Set aux = [[@BSR.LinearBlock(pi->Name, pi->InitValue, True, False)]];
-    //WriteLn(_MID+" TRACE 4.["<<j+"].1");
+    //WriteLn("TRACE "+_MID+" 4.["<<j+"].1");
       Real prior.Average = If(ObjectExist("Real","pi->Prior.Average"),
                               pi->Prior.Average, ?);
@@ -129 +129 @@
         0
       });
-    //WriteLn(_MID+" TRACE 4.["<<j+"].2");
+    //WriteLn("TRACE "+_MID+" 4.["<<j+"].2");
       Real prior.LowerBound = If(ObjectExist("Real","pi->Prior.LowerBound"),
                               pi->Prior.LowerBound, -1/0);
@@ -140 +140 @@
         0
       });
-    //WriteLn(_MID+" TRACE 4.["<<j+"].3");
+    //WriteLn("TRACE "+_MID+" 4.["<<j+"].3");
       Real prior.UpperBound = If(ObjectExist("Real","pi->Prior.UpperBound"),
                               pi->Prior.UpperBound, +1/0);
@@ -157 +157 @@
     //Checks unique name and prepares for hashed searching by name
     Real SetIndexByName(linearBlock);
-  //WriteLn(_MID+" TRACE 5");
+  //WriteLn("TRACE "+_MID+" 5");
     //Domain constraints
     VMatrix domain.A = Triplet(domain.A.triplet,domain.num,n);
@@ -184 +184 @@
     Real m.full = m+prior.num;
 
-  //WriteLn(_MID+" TRACE 6");
+  //WriteLn("TRACE "+_MID+" 6");
     //Observational information
     Real sigmaIsUnknown = IsUnknown(data::_.sigma);
@@ -207 +207 @@
     ]]};
 
-  //WriteLn(_MID+" TRACE 7");
+  //WriteLn("TRACE "+_MID+" 7.1");
     //Prior information
-    VMatrix prior.X = Triplet(prior.X.triplet,prior.num,n);
+    VMatrix prior.X = Convert(Triplet(prior.X.triplet,prior.num,n),"Cholmod.R.Sparse");
+  //WriteLn("TRACE "+_MID+" 7.2");
     Real sigmaIndex = If(sigmaIsUnknown,1,0);
     Set noisePrior = For(1,Card(prior.index),Set(Real j)
     {
+    //If(!j%100, WriteLn("TRACE "+_MID+"  8."<<j));
       Set pi = (data::_.linearParamInfo)[prior.index[j] ];
       Real estimSigma = Case(
@@ -246 +248 @@
         Set NonLinFilters = Copy(Empty)
       )
-   });
-
-  //WriteLn(_MID+" TRACE 8");
+    });
+
+  //WriteLn("TRACE "+_MID+" 9.1");
     //All allowed segments in this simplified API are obs and prior (optionally)
     Set noiseDistrib_ = noiseObs << noisePrior;
@@ -259 +261 @@
     //Checks unique name and prepares for hashed searching by name
     Real SetIndexByName(noiseDistrib);
-  //WriteLn(_MID+" TRACE 9");
+  //WriteLn("TRACE "+_MID+" 9.2");
     //Building generic model definition
     Set modelDef = @BSR.ModelDef (
@@ -275 +277 @@
       VMatrix A = domain.A << order.A
     );
+  //WriteLn("TRACE "+_MID+" 9.3");
     BysMcmc::@Cycler cycler = BysMcmc::Bsr::Gibbs::GetCyclerFromParsed
-   (
+    (
       modelDef,
       BysMcmc::@McmcPartialHandler::Null(0),
@@ -282 +285 @@
       BysMcmc::@Notifier::Null(0)
     );
+  //WriteLn("TRACE "+_MID+" 9.4");
     Set Append(cycler::auxInfo,[[data]]);
-  //WriteLn(_MID+" TRACE 12");
+  //WriteLn("TRACE "+_MID+" 9.5");
     cycler
   };