Changeset 3138 in ntrip for trunk/BNC


Ignore:
Timestamp:
Mar 25, 2011, 1:40:44 PM (14 years ago)
Author:
mervart
Message:
 
Location:
trunk/BNC/combination
Files:
2 edited

Legend:

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

    r3137 r3138  
    202202    }
    203203  }
     204
     205  // Not yet regularized
     206  // -------------------
     207  _firstReg = false;
    204208}
    205209
     
    600604    const double Pl = 1.0 / (0.05 * 0.05);
    601605
    602     const int nCon = 2 + MAXPRN_GPS;
     606    const int nCon = (_firstReg == false) ? 2 + MAXPRN_GPS : 2;
    603607    Matrix         AA(nObs+nCon, nPar);  AA = 0.0;
    604608    ColumnVector   ll(nObs+nCon);        ll = 0.0;
     
    653657    }
    654658
    655     for (int iGps = 1; iGps <= MAXPRN_GPS; iGps++) {
    656       ++iCond;
    657       QString prn = QString("G%1").arg(iGps, 2, 10, QChar('0'));
    658       PP(nObs+1+iGps) = Ph;
    659       for (int iPar = 1; iPar <= _params.size(); iPar++) {
    660         cmbParam* pp = _params[iPar-1];
    661         if (pp->type == cmbParam::Sat_offset && pp->prn == prn) {
    662           AA(nObs+iCond, iPar) = 1.0;
     659    if (!_firstReg) {
     660      _firstReg = true;
     661      for (int iGps = 1; iGps <= MAXPRN_GPS; iGps++) {
     662        ++iCond;
     663        QString prn = QString("G%1").arg(iGps, 2, 10, QChar('0'));
     664        PP(nObs+1+iGps) = Ph;
     665        for (int iPar = 1; iPar <= _params.size(); iPar++) {
     666          cmbParam* pp = _params[iPar-1];
     667          if (pp->type == cmbParam::Sat_offset && pp->prn == prn) {
     668            AA(nObs+iCond, iPar) = 1.0;
     669          }
    663670        }
    664671      }
  • trunk/BNC/combination/bnccomb.h

    r3052 r3138  
    8383  QVector<cmbParam*>    _params;
    8484  SymmetricMatrix       _QQ;
     85  bool                  _firstReg;
    8586  QByteArray            _log;
    8687  QString               _masterAC;
Note: See TracChangeset for help on using the changeset viewer.