Changeset 2124 in ntrip
- Timestamp:
- Dec 28, 2009, 9:06:59 AM (15 years ago)
- Location:
- trunk/BNC
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNC/bncmodel.cpp
r2118 r2124 172 172 173 173 if (epoData->size() < MINOBS) { 174 _log += "\nNot enough data"; 174 175 return failure; 175 176 } … … 399 400 t_irc bncModel::update(t_epoData* epoData) { 400 401 402 _log = "Precise Point Positioning"; 403 401 404 SymmetricMatrix QQsav; 402 405 ColumnVector dx; … … 410 413 // ----------------- 411 414 if (cmpBancroft(epoData) != success) { 415 _log += "\nBancroft failed"; 416 emit newMessage(_log, false); 412 417 return failure; 413 418 } 414 419 415 420 if (epoData->size() < MINOBS) { 421 _log += "\nNot enough data"; 422 emit newMessage(_log, false); 416 423 return failure; 417 424 } … … 482 489 // Set Solution Vector 483 490 // ------------------- 491 ostringstream str1; 492 str1.setf(ios::fixed); 484 493 QVectorIterator<bncParam*> itPar(_params); 485 494 while (itPar.hasNext()) { 486 495 bncParam* par = itPar.next(); 487 496 par->xx += dx(par->index); 488 } 497 if (par->type == bncParam::RECCLK) { 498 str1 << "\n clk = " << setw(6) << setprecision(3) << par->xx 499 << " +- " << setw(6) << setprecision(3) 500 << sqrt(_QQ(par->index,par->index)); 501 } 502 else if (par->type == bncParam::AMB_L3) { 503 str1 << "\n amb " << par->prn.toAscii().data() << " = " 504 << setw(6) << setprecision(3) << par->xx 505 << " +- " << setw(6) << setprecision(3) 506 << sqrt(_QQ(par->index,par->index)); 507 } 508 } 509 _log += str1.str().c_str(); 489 510 490 511 // Message (both log file and screen) 491 512 // ---------------------------------- 492 ostringstream str ;493 str .setf(ios::fixed);494 str << " PPP " << _staID.data() << " "495 << epoData->tt.timestr(1) << " " << epoData->size() << " "496 << setw(14) << setprecision(3) << x() << " +- "497 << setw(6) << setprecision(3) << sqrt(_QQ(1,1)) << " "498 << setw(14) << setprecision(3) << y() << " +- "499 << setw(6) << setprecision(3) << sqrt(_QQ(2,2)) << " "500 << setw(14) << setprecision(3) << z() << " +- "501 << setw(6) << setprecision(3) << sqrt(_QQ(3,3));513 ostringstream str2; 514 str2.setf(ios::fixed); 515 str2 << " PPP " << _staID.data() << " " 516 << epoData->tt.timestr(1) << " " << epoData->size() << " " 517 << setw(14) << setprecision(3) << x() << " +- " 518 << setw(6) << setprecision(3) << sqrt(_QQ(1,1)) << " " 519 << setw(14) << setprecision(3) << y() << " +- " 520 << setw(6) << setprecision(3) << sqrt(_QQ(2,2)) << " " 521 << setw(14) << setprecision(3) << z() << " +- " 522 << setw(6) << setprecision(3) << sqrt(_QQ(3,3)); 502 523 if (_estTropo) { 503 str << " " << setw(6) << setprecision(3) << trp() << " +- " 504 << setw(6) << setprecision(3) << sqrt(_QQ(5,5)); 505 } 506 507 emit newMessage(QString(str.str().c_str()).toAscii(), true); 524 str2 << " " << setw(6) << setprecision(3) << trp() << " +- " 525 << setw(6) << setprecision(3) << sqrt(_QQ(5,5)); 526 } 527 528 emit newMessage(_log, false); 529 emit newMessage(QByteArray(str2.str().c_str()), true); 508 530 509 531 return success; … … 548 570 _QQ = QQsav; 549 571 550 QByteArray msg = "Outlier Code " + prn.toAscii() + " " 551 + QByteArray::number(vvMaxCode, 'f', 3); 552 emit newMessage(msg, true); 572 _log += "\nOutlier Code " + prn.toAscii() + " " 573 + QByteArray::number(vvMaxCode, 'f', 3); 553 574 554 575 return 1; … … 561 582 _QQ = QQsav; 562 583 563 QByteArray msg = "Outlier Phase " + prn.toAscii() + " " 564 + QByteArray::number(vvMaxPhase, 'f', 3); 565 emit newMessage(msg, true); 584 _log += "\nOutlier Phase " + prn.toAscii() + " " 585 + QByteArray::number(vvMaxPhase, 'f', 3); 566 586 567 587 return 1; -
trunk/BNC/bncmodel.h
r2113 r2124 82 82 bool _usePhase; 83 83 bool _estTropo; 84 QByteArray _log; 84 85 }; 85 86
Note:
See TracChangeset
for help on using the changeset viewer.