Changeset 3384 in ntrip for trunk/BNC


Ignore:
Timestamp:
Aug 30, 2011, 5:14:09 PM (13 years ago)
Author:
mervart
Message:
 
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/bncmodel.cpp

    r3383 r3384  
    12381238  Tracer tracer("bncModel::update_p");
    12391239
     1240  // Bancroft Solution
     1241  // -----------------
     1242  if (cmpBancroft(epoData) != success) {
     1243    return failure;
     1244  }
     1245
    12401246  rememberState(epoData);
    12411247
     
    12491255    allPrns.push_back(satData->prn);
    12501256  }
    1251 
    12521257  std::vector<QString> usedPrns;
    12531258
     
    12751280      }
    12761281
     1282      // Remove Neglected Satellites form epoData
     1283      // ----------------------------------------
    12771284      for (unsigned ip = 0; ip < allPrns.size(); ip++) {
    12781285        QString prn = allPrns[ip];
     
    12811288        }
    12821289      }
     1290      if (epoData->sizeSys('G') < MINOBS) {
     1291        continue;
     1292      }
    12831293
    12841294      // First update using code observations, then phase observations
     
    12861296      for (int iPhase = 0; iPhase <= (_usePhase ? 1 : 0); iPhase++) {
    12871297     
    1288         // Bancroft Solution
    1289         // -----------------
    1290         if (iPhase == 0) {     
    1291           if (cmpBancroft(epoData) != success) {
    1292             restoreState(epoData);
    1293             return failure;
    1294           }
    1295           else {
    1296             if (nNeglected == 0) {
    1297               _epoData_sav->deepCopy(epoData);
    1298             }
    1299           }
    1300         }
    1301         else {
    1302           if (epoData->sizeSys('G') < MINOBS) {
    1303             restoreState(epoData);
    1304             _log += "bncModel::update_p: not enough data\n";
    1305             return failure;
    1306           }
    1307           unsigned numSatNoSlip = 0;
    1308           QVectorIterator<bncParam*> itPar(_params);
    1309           while (itPar.hasNext()) {
    1310             bncParam* par = itPar.next();
    1311             if (par->type == bncParam::AMB_L3 && par->prn[0] == 'G') {
    1312               if (par->numEpo >= 1) {
    1313                 ++numSatNoSlip;
    1314               }
    1315             }
    1316           }
    1317           if (numSatNoSlip > 0 && numSatNoSlip < MINOBS) {
    1318             restoreState(epoData);
    1319             _log += "bncModel::update_p: not enough GPS satellites without cycle-slips\n";
    1320             return failure;
    1321           }
    1322         }
    1323        
    13241298        // Status Prediction
    13251299        // -----------------
Note: See TracChangeset for help on using the changeset viewer.