Changeset 4c93933 in tol

Timestamp:

Oct 27, 2008, 8:37:27 AM (18 years ago)

Author:

Víctor de Buen <vdebuen@…>

Branches:

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

Children:

dacdfcd8

Parents:

64f4b50

Message:

Storing failed factorization matrix

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

File:

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

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

@@ -112 +112 @@
   ////////////////////////////////////////////////////////////////////////////
   {
+    WriteLn("TRACE"+_MID+" factorize X("<<VRows(X)+"x"+VColumns(X)+")");
     If(And(design.is.constant,VRows(_.L)), 0,
     {
       Real sparse = VNonNullCells(X)/(VRows(X)*VColumns(X));
-
+      Real numErr = Copy(NError);
       VMatrix _.L := If(sparse<.20, 
       {
@@ -125 +126 @@
         VMatrix XtX  = Convert(XtX_,"Cholmod.R.Sparse");
         CholeskiFactor(XtX,"X")
+      });
+      Real If(Copy(NError)>numErr,
+      {
+        Text path = TolAppDataPath+"tmp/"+ToName(_MID)+".oza";
+        WriteLn("Failed factorization matrix will be stored into "+path);
+        Ois.Store([[X]],path)
       });
       VMatrix If(_.R, _.Ls := Convert(_.L, "Cholmod.R.Sparse"));