Changeset 3140 in ntrip
- Timestamp:
- Mar 25, 2011, 3:41:05 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNC/combination/bnccomb.cpp
r3138 r3140 604 604 const double Pl = 1.0 / (0.05 * 0.05); 605 605 606 const int nCon = (_firstReg == false) ? 2 + MAXPRN_GPS : 2;606 const int nCon = (_firstReg == false) ? 1 + MAXPRN_GPS : 1; 607 607 Matrix AA(nObs+nCon, nPar); AA = 0.0; 608 608 ColumnVector ll(nObs+nCon); ll = 0.0; … … 639 639 const double Ph = 1.e6; 640 640 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;650 641 PP(nObs+iCond) = Ph; 651 642 for (int iPar = 1; iPar <= _params.size(); iPar++) { … … 681 672 ColumnVector vv = ll - AA * dx; 682 673 674 for (int iv = 1; iv < vv.Nrows(); iv++) { 675 out << vv(iv) << " "; 676 } 677 out << endl; 678 683 679 int maxResIndex; 684 680 double maxRes = vv.maximum_absolute_value1(maxResIndex);
Note:
See TracChangeset
for help on using the changeset viewer.