Changeset 824d400 in tol


Ignore:
Timestamp:
Aug 25, 2009, 11:35:32 AM (17 years ago)
Author:
Víctor de Buen <vdebuen@…>
Branches:
1941-docker-image-to-test, devel, master
Children:
d8928a5
Parents:
6be9c2a
Message:

Method build must set all to zero or empty when _.N==0

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • tol/stdlib/math/stat/models/bayesian/bysMcmc/_partialSampler.tol

    r6be9c2a r824d400  
    165165    Set fixedParamIdx = Copy(Empty);
    166166    Set reloadParamIdx = Copy(Empty);
    167     If(_.N,
     167    If(!_.N,
     168    {
     169      Matrix _.fixedParamIdx  := Rand(0,1,0,0);
     170      Matrix _.reloadParamIdx := Rand(0,1,0,0);
     171      Matrix _.freeParamIdx   := Rand(0,1,0,0);
     172      Set    _.fixedParam     := Copy(Empty);
     173      Set    _.reloadParam    := Copy(Empty);
     174      Set    _.freeParam      := Copy(Empty);
     175      Real   _.fixedParamNum  := 0;
     176      Real   _.reloadParamNum := 0;
     177      Real   _.freeParamNum   := 0;
     178      1
     179    },
    168180    {
    169181      Real If(Card(_.mcmcPartialHandler::_.fixedParam),
    … …  
    189201        1
    190202      });
     203      Set freeParamIdx = Range(1,_.N,1)-(fixedParamIdx << reloadParamIdx);
     204      Matrix _.freeParamIdx := SetCol(freeParamIdx);
     205      Set _.freeParam := ExtractByIndex(blockParam, freeParamIdx);
     206      Real _.fixedParamNum  := Rows(_.fixedParamIdx);
     207      Real _.reloadParamNum := Rows(_.reloadParamIdx);
     208      Real _.freeParamNum   := Rows(_.freeParamIdx);
     209      VMatrix If(_.fixedParamNum | _.reloadParamNum,
     210      {
     211        getMergedFixedAndReloadedSample(0)
     212      });
    191213      1
    192     });
    193     Set freeParamIdx = Range(1,_.N,1)-(fixedParamIdx << reloadParamIdx);
    194     Matrix _.freeParamIdx := SetCol(freeParamIdx);
    195     Set _.freeParam := ExtractByIndex(blockParam, freeParamIdx);
    196     Real _.fixedParamNum  := Rows(_.fixedParamIdx);
    197     Real _.reloadParamNum := Rows(_.reloadParamIdx);
    198     Real _.freeParamNum   := Rows(_.freeParamIdx);
    199     VMatrix If(_.fixedParamNum | _.reloadParamNum,
    200     {
    201       getMergedFixedAndReloadedSample(0)
    202     });
    203     1
     214    })
    204215  };
    205216
Note: See TracChangeset for help on using the changeset viewer.