Changeset 5800 in ntrip


Ignore:
Timestamp:
Aug 4, 2014, 6:26:53 PM (8 years ago)
Author:
mervart
Message:
 
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/src/PPP/filter.cpp

    r5799 r5800  
    207207    // ------------------
    208208    //// beg test
    209     DiagonalMatrix PP(iObs+1);
    210     for (int ii = 1; ii < iObs+1; ii++) {
    211       PP(ii,ii) = Sl(ii,ii) * Sl(ii,ii);
    212     }
    213     bncModel::kalman(AA, ll, PP, _QFlt, _xFlt);
     209    DiagonalMatrix PP(iObs+1); PP = 0.0;
     210    for (int ii = 1; ii <= iObs+1; ii++) {
     211      PP(ii,ii) = 1.0 / (Sl(ii,ii) * Sl(ii,ii));
     212    }
     213    ColumnVector dx;
     214    bncModel::kalman(AA, ll, PP, _QFlt, dx);
     215    _xFlt += dx;
    214216    //// end test
    215217
Note: See TracChangeset for help on using the changeset viewer.