Changeset 0bdc6a7 in tol

Timestamp:

Dec 2, 2008, 11:05:28 AM (18 years ago)

Author:

Víctor de Buen <vdebuen@…>

Branches:

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

Children:

8c8cb65

Parents:

3e1e359

Message:

Fixed bug in resuming

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

File:

• tol/stdlib/math/stat/models/bayesian/bysMcmc/_build.tol (modified) (2 diffs)

tol/stdlib/math/stat/models/bayesian/bysMcmc/_build.tol

@@ -119 +119 @@
           Real m_ = Case(
             resume_>0, Min(m,resume_),
-            resume_<0,  Max(0,m+resume_) );
+            resume_<0, Max(0,m+resume_) );
           If(m_==0,
           {
@@ -128 +128 @@
           {
             Real WriteLn(_MID+"resuming "<<m_+" old simulations");
-            Real If(m_==m,0,
+            Matrix resumed = If(m_==m,
+            {
+              MatReadFile(_.path)
+            },
             {
               Real WriteLn(_MID+"Removing last "<<Real(m-m_)+" old simulations");
-              Matrix resumed = Sub(old.mcmc,1,1,m_,_.n);
-              Matrix MatWriteFile(_.path, resumed);
-              1
+              Matrix aux = MatReadRows(_.path,1,m_,1);
+              Matrix MatWriteFile(_.path, aux);
+              aux
             });
-            Matrix lastRow = Sub(old.mcmc,m_,1,1,_.n);
+            Matrix lastRow = Sub(resumed,m_,1,1,_.n);
           //WriteLn("TRACE "+_MID+"lastRow="<<lastRow);
             Real _.sampler::setStore(lastRow);