- Timestamp:
- Mar 25, 2011, 1:40:44 PM (14 years ago)
- Location:
- trunk/BNC/combination
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNC/combination/bnccomb.cpp
r3137 r3138 202 202 } 203 203 } 204 205 // Not yet regularized 206 // ------------------- 207 _firstReg = false; 204 208 } 205 209 … … 600 604 const double Pl = 1.0 / (0.05 * 0.05); 601 605 602 const int nCon = 2 + MAXPRN_GPS;606 const int nCon = (_firstReg == false) ? 2 + MAXPRN_GPS : 2; 603 607 Matrix AA(nObs+nCon, nPar); AA = 0.0; 604 608 ColumnVector ll(nObs+nCon); ll = 0.0; … … 653 657 } 654 658 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 } 663 670 } 664 671 } -
trunk/BNC/combination/bnccomb.h
r3052 r3138 83 83 QVector<cmbParam*> _params; 84 84 SymmetricMatrix _QQ; 85 bool _firstReg; 85 86 QByteArray _log; 86 87 QString _masterAC;
Note:
See TracChangeset
for help on using the changeset viewer.