Changeset 3393 in ntrip for trunk


Ignore:
Timestamp:
Aug 31, 2011, 11:45:02 AM (13 years ago)
Author:
mervart
Message:
 
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/bncmodel.cpp

    r3392 r3393  
    12291229    }
    12301230
    1231     bool outlierDetected = false;
    1232 
    12331231    // Loop over all Combinations of used Satellites
    12341232    // ---------------------------------------------
    12351233    do {
    1236 
    1237       if (outlierDetected) {
    1238         _log += "TRY WITH PRNs: ";
    1239         for (unsigned ip = 0; ip < usedPrns.size(); ip++) {
    1240           _log += usedPrns[ip] + ' ';
    1241         }
    1242         _log += '\n';
    1243       }
    12441234
    12451235      // Remove Neglected Satellites from epoData
     
    12741264        }
    12751265       
     1266        // Prepare first-design Matrix, vector observed-computed
     1267        // -----------------------------------------------------
    12761268        Matrix          AA(nObs, nPar);  // first design matrix
    12771269        ColumnVector    ll(nObs);        // tems observed-computed
     
    12791271       
    12801272        unsigned iObs = 0;
    1281        
    1282         // GPS
    1283         // ---
    12841273        QMapIterator<QString, t_satData*> it(epoData->satData);
    12851274        while (it.hasNext()) {
     
    13171306        // Check the residuals
    13181307        // -------------------
    1319         outlierDetected = outlierDetection(iPhase, vv, epoData->satData);
    1320 
    1321         if (outlierDetected) {
     1308        if ( outlierDetection(iPhase, vv, epoData->satData) ) {
    13221309          restoreState(epoData);
    13231310          break;
    13241311        }
    13251312
     1313        // Set estimated values
     1314        // --------------------
     1315        else if (_usePhase || iPhase == 1) {
     1316          QVectorIterator<bncParam*> itPar(_params);
     1317          while (itPar.hasNext()) {
     1318            bncParam* par = itPar.next();
     1319            par->xx += dx(par->index);
     1320          }
     1321          return success;
     1322        }
     1323
    13261324      } // for (int iPhase = 0; iPhase <= (_usePhase ? 1 : 0); iPhase++)
    1327 
    1328       // Update Parameters
    1329       // -----------------
    1330       if (!outlierDetected) {
    1331         QVectorIterator<bncParam*> itPar(_params);
    1332         while (itPar.hasNext()) {
    1333           bncParam* par = itPar.next();
    1334           par->xx += dx(par->index);
    1335         }
    1336         return success;
    1337       }
    13381325
    13391326    } while ( next_combination(allPrns.begin(), allPrns.end(),
Note: See TracChangeset for help on using the changeset viewer.