Changeset 84e12af in tol
- Timestamp:
- Feb 12, 2009, 5:06:13 PM (18 years ago)
- Branches:
- 1941-docker-image-to-test, devel, master
- Children:
- 9bfd018
- Parents:
- f4a032f
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/gibbs/_nonLinBlock.tol
rf4a032f r84e12af 146 146 Matrix _.Z := _.W + _.Fi * omega; 147 147 //WriteLn("TRACE"+_MID+"draw.scalar 1.5 _.Z("<<Rows(_.Z)+"x"<<Columns(_.Z)); 148 Set aux = _.segment->ArimaAuxInfo; 149 Set arima = _.segment->Arima; 150 Polyn dif = ARIMAGetDIF(arima); 148 Set segment = _.segment; 149 Set aux = _.segment->ArimaAuxInfo; 150 Set arima = _.segment->Arima; 151 Polyn dif = ARIMAGetDIF(arima); 151 152 // 152 153 VMatrix si = SubRow(SubCol(sigBlk::_.sigmaInv,_.segment->EquIdx),_.segment->EquIdx); … … 196 197 Matrix Zi = EvalFilteredNoise(x); 197 198 VMatrix Zi_ = Mat2VMat(Zi); 198 VMatrix ei_ = arimaBlk::filterJustBlock( _.segment, Zi_);199 VMatrix ei_ = arimaBlk::filterJustBlock(segment, Zi_); 199 200 Matrix ei = VMat2Mat(ei_); 200 201 Real ld = -0.5*MatDat(MtMSqr(ei),1,1); … … 210 211 //WriteLn("TRACE Wi("<<Rows(Wi)+","<<Columns(Wi)+")"); 211 212 VMatrix Wi_ = Mat2VMat(Wi); 212 VMatrix ei_ = arimaBlk::filterJustBlock( _.segment, Wi_);213 VMatrix ei_ = arimaBlk::filterJustBlock(segment, Wi_); 213 214 //WriteLn("TRACE ei_("<<VRows(ei_)+","<<VColumns(ei_)+")"); 214 215 Matrix ei = VMat2Mat(Sub(ei_,d+1,1,n,1));
Note: See TracChangeset
for help on using the changeset viewer.
