- Timestamp:
- Aug 30, 2011, 5:14:09 PM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNC/bncmodel.cpp
r3383 r3384 1238 1238 Tracer tracer("bncModel::update_p"); 1239 1239 1240 // Bancroft Solution 1241 // ----------------- 1242 if (cmpBancroft(epoData) != success) { 1243 return failure; 1244 } 1245 1240 1246 rememberState(epoData); 1241 1247 … … 1249 1255 allPrns.push_back(satData->prn); 1250 1256 } 1251 1252 1257 std::vector<QString> usedPrns; 1253 1258 … … 1275 1280 } 1276 1281 1282 // Remove Neglected Satellites form epoData 1283 // ---------------------------------------- 1277 1284 for (unsigned ip = 0; ip < allPrns.size(); ip++) { 1278 1285 QString prn = allPrns[ip]; … … 1281 1288 } 1282 1289 } 1290 if (epoData->sizeSys('G') < MINOBS) { 1291 continue; 1292 } 1283 1293 1284 1294 // First update using code observations, then phase observations … … 1286 1296 for (int iPhase = 0; iPhase <= (_usePhase ? 1 : 0); iPhase++) { 1287 1297 1288 // Bancroft Solution1289 // -----------------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 1324 1298 // Status Prediction 1325 1299 // -----------------
Note:
See TracChangeset
for help on using the changeset viewer.