Changeset 3140 in ntrip


Ignore:
Timestamp:
Mar 25, 2011, 3:41:05 PM (13 years ago)
Author:
mervart
Message:
 
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/combination/bnccomb.cpp

    r3138 r3140  
    604604    const double Pl = 1.0 / (0.05 * 0.05);
    605605
    606     const int nCon = (_firstReg == false) ? 2 + MAXPRN_GPS : 2;
     606    const int nCon = (_firstReg == false) ? 1 + MAXPRN_GPS : 1;
    607607    Matrix         AA(nObs+nCon, nPar);  AA = 0.0;
    608608    ColumnVector   ll(nObs+nCon);        ll = 0.0;
     
    639639    const double Ph = 1.e6;
    640640    int iCond = 1;
    641     PP(nObs+iCond) = Ph;
    642     for (int iPar = 1; iPar <= _params.size(); iPar++) {
    643       cmbParam* pp = _params[iPar-1];
    644       if      (pp->type == cmbParam::AC_offset) {
    645         AA(nObs+iCond, iPar) = 1.0;
    646       }
    647     }
    648 
    649     ++iCond;
    650641    PP(nObs+iCond) = Ph;
    651642    for (int iPar = 1; iPar <= _params.size(); iPar++) {
     
    681672      ColumnVector vv = ll - AA * dx;
    682673
     674      for (int iv = 1; iv < vv.Nrows(); iv++) {
     675        out << vv(iv) << " ";
     676      }
     677      out << endl;
     678
    683679      int    maxResIndex;
    684680      double maxRes = vv.maximum_absolute_value1(maxResIndex);   
Note: See TracChangeset for help on using the changeset viewer.