Changeset 7297 in ntrip for trunk/BNC/src/combination


Ignore:
Timestamp:
Sep 22, 2015, 4:42:55 PM (10 years ago)
Author:
stuerze
Message:

some re-naming

File:
1 edited

Legend:

Unmodified
Added
Removed

  • TabularUnified trunk/BNC/src/combination/bnccomb.cpp

    r7132 r7297  
    1111 * Created:    22-Jan-2011
    1212 *
    13  * Changes:   
     13 * Changes:
    1414 *
    1515 * -----------------------------------------------------------------------*/
    … …  
    7878////////////////////////////////////////////////////////////////////////////
    7979double bncComb::cmbParam::partial(const QString& AC_, const QString& prn_) {
    80  
     80
    8181  if      (type == offACgps) {
    8282    if (AC == AC_ && prn_[0] == 'G') {
    … …  
    103103}
    104104
    105 // 
     105//
    106106////////////////////////////////////////////////////////////////////////////
    107107QString bncComb::cmbParam::toString() const {
    108108
    109109  QString outStr;
    110  
     110
    111111  if      (type == offACgps) {
    112112    outStr = "AC offset GPS " + AC;
    … …  
    131131  bncSettings settings;
    132132
    133   QStringList combineStreams = settings.value("combineStreams").toStringList();
     133  QStringList cmbStreams = settings.value("cmbStreams").toStringList();
    134134
    135135  _cmbSampl = settings.value("cmbSampl").toInt();
    … …  
    140140  _masterMissingEpochs = 0;
    141141
    142   if (combineStreams.size() >= 1 && !combineStreams[0].isEmpty()) {
    143     QListIterator<QString> it(combineStreams);
     142  if (cmbStreams.size() >= 1 && !cmbStreams[0].isEmpty()) {
     143    QListIterator<QString> it(cmbStreams);
    144144    while (it.hasNext()) {
    145145      QStringList hlp = it.next().split(" ");
    … …  
    157157  _rtnetDecoder = 0;
    158158
    159   connect(this,     SIGNAL(newMessage(QByteArray,bool)), 
     159  connect(this,     SIGNAL(newMessage(QByteArray,bool)),
    160160          BNC_CORE, SLOT(slotMessage(const QByteArray,bool)));
    161161
    … …  
    197197      for (unsigned iGps = 1; iGps <= t_prn::MAXPRN_GPS; iGps++) {
    198198        QString prn = QString("G%1_0").arg(iGps, 2, 10, QChar('0'));
    199         _params.push_back(new cmbParam(cmbParam::offACSat, ++nextPar, 
     199        _params.push_back(new cmbParam(cmbParam::offACSat, ++nextPar,
    200200                                       AC->name, prn));
    201201      }
    … …  
    204204        for (unsigned iGlo = 1; iGlo <= t_prn::MAXPRN_GLONASS; iGlo++) {
    205205          QString prn = QString("R%1_0").arg(iGlo, 2, 10, QChar('0'));
    206           _params.push_back(new cmbParam(cmbParam::offACSat, ++nextPar, 
     206          _params.push_back(new cmbParam(cmbParam::offACSat, ++nextPar,
    207207                                         AC->name, prn));
    208208        }
    … …  
    219219      }
    220220    }
    221    
     221
    222222    // Initialize Variance-Covariance Matrix
    223223    // -------------------------------------
    … …  
    353353      continue;
    354354    }
    355  
     355
    356356    // Create new correction
    357357    // ---------------------
    … …  
    424424}
    425425
    426 // Change the correction so that it refers to last received ephemeris 
     426// Change the correction so that it refers to last received ephemeris
    427427////////////////////////////////////////////////////////////////////////////
    428428void bncComb::switchToLastEph(t_eph* lastEph, cmbCorr* corr) {
    … …  
    471471  QTextStream out(&_log, QIODevice::WriteOnly);
    472472
    473   out << endl <<           "Combination:" << endl 
     473  out << endl <<           "Combination:" << endl
    474474      << "------------------------------" << endl;
    475475
    … …  
    515515        if (AC->numObs > 0) {
    516516          out << "Switching Master AC "
    517               << _masterOrbitAC.toAscii().data() << " --> " 
    518               << AC->name.toAscii().data()   << " " 
    519               << _resTime.datestr().c_str()    << " " 
     517              << _masterOrbitAC.toAscii().data() << " --> "
     518              << AC->name.toAscii().data()   << " "
     519              << _resTime.datestr().c_str()    << " "
    520520              << _resTime.timestr().c_str()    << endl;
    521521          _masterOrbitAC = AC->name;
    … …  
    550550        }
    551551      }
    552       out << _resTime.datestr().c_str() << " " 
     552      out << _resTime.datestr().c_str() << " "
    553553          << _resTime.timestr().c_str() << " ";
    554554      out.setRealNumberNotation(QTextStream::FixedNotation);
    … …  
    618618
    619619    int     maxResIndex;
    620     double  maxRes = vv.maximum_absolute_value1(maxResIndex);   
     620    double  maxRes = vv.maximum_absolute_value1(maxResIndex);
    621621    out.setRealNumberNotation(QTextStream::FixedNotation);
    622     out.setRealNumberPrecision(3); 
     622    out.setRealNumberPrecision(3);
    623623    out << _resTime.datestr().c_str() << " " << _resTime.timestr().c_str()
    624624        << " Maximum Residuum " << maxRes << ' '
    … …  
    627627      for (int iPar = 1; iPar <= _params.size(); iPar++) {
    628628        cmbParam* pp = _params[iPar-1];
    629         if (pp->type == cmbParam::offACSat            && 
     629        if (pp->type == cmbParam::offACSat            &&
    630630            pp->AC   == corrs()[maxResIndex-1]->_acName &&
    631631            pp->prn  == corrs()[maxResIndex-1]->_prn.mid(0,3)) {
    … …  
    675675      eph->getCrd(_resTime, xc, vv, false);
    676676
    677       out << _resTime.datestr().c_str() << " " 
     677      out << _resTime.datestr().c_str() << " "
    678678          << _resTime.timestr().c_str() << " ";
    679679      out.setFieldWidth(3);
    … …  
    704704  _resTime.civil_time(hour, minute, sec);
    705705
    706   outLines.sprintf("*  %4d %2d %2d %d %d %12.8f\n", 
     706  outLines.sprintf("*  %4d %2d %2d %d %d %12.8f\n",
    707707                   year, month, day, hour, minute, sec);
    708708
    … …  
    797797////////////////////////////////////////////////////////////////////////////
    798798t_irc bncComb::createAmat(Matrix& AA, ColumnVector& ll, DiagonalMatrix& PP,
    799                           const ColumnVector& x0, 
     799                          const ColumnVector& x0,
    800800                          QMap<QString, cmbCorr*>& resCorr) {
    801801
    802802  unsigned nPar = _params.size();
    803   unsigned nObs = corrs().size(); 
     803  unsigned nObs = corrs().size();
    804804
    805805  if (nObs == 0) {
    … …  
    858858        cmbParam* pp = _params[iPar-1];
    859859        if ( AA.Column(iPar).maximum_absolute_value() > 0.0 &&
    860              pp->type == cmbParam::offACSat                 && 
     860             pp->type == cmbParam::offACSat                 &&
    861861             pp->prn == prn) {
    862862          AA(nObs+iCond, iPar) = 1.0;
    … …  
    872872//          cmbParam* pp = _params[iPar-1];
    873873//          if ( AA.Column(iPar).maximum_absolute_value() > 0.0 &&
    874 //               pp->type == cmbParam::offACSat                 && 
     874//               pp->type == cmbParam::offACSat                 &&
    875875//               pp->prn == prn) {
    876876//            AA(nObs+iCond, iPar) = 1.0;
    … …  
    899899  // ----------------------
    900900  while (true) {
    901    
     901
    902902    // Remove Satellites that are not in Master
    903903    // ----------------------------------------
    … …  
    922922      }
    923923    }
    924    
     924
    925925    // Count Number of Observations per Satellite and per AC
    926926    // -----------------------------------------------------
    … …  
    945945      }
    946946    }
    947    
     947
    948948    // Clean-Up the Paramters
    949949    // ----------------------
    … …  
    952952    }
    953953    _params.clear();
    954    
     954
    955955    // Set new Parameters
    956956    // ------------------
    957957    int nextPar = 0;
    958    
     958
    959959    QMapIterator<QString, int> itAC(numObsAC);
    960960    while (itAC.hasNext()) {
    … …  
    968968        }
    969969      }
    970     } 
    971    
     970    }
     971
    972972    QMapIterator<QString, int> itPrn(numObsPrn);
    973973    while (itPrn.hasNext()) {
    … …  
    978978        _params.push_back(new cmbParam(cmbParam::clkSat, ++nextPar, "", prn));
    979979      }
    980     } 
    981    
     980    }
     981
    982982    int nPar = _params.size();
    983     ColumnVector x0(nPar); 
     983    ColumnVector x0(nPar);
    984984    x0 = 0.0;
    985    
     985
    986986    // Create First-Design Matrix
    987987    // --------------------------
    … …  
    992992      return failure;
    993993    }
    994    
     994
    995995    ColumnVector vv;
    996996    try {
    … …  
    10081008
    10091009    int     maxResIndex;
    1010     double  maxRes = vv.maximum_absolute_value1(maxResIndex);   
     1010    double  maxRes = vv.maximum_absolute_value1(maxResIndex);
    10111011    out.setRealNumberNotation(QTextStream::FixedNotation);
    1012     out.setRealNumberPrecision(3); 
     1012    out.setRealNumberPrecision(3);
    10131013    out << _resTime.datestr().c_str() << " " << _resTime.timestr().c_str()
    10141014        << " Maximum Residuum " << maxRes << ' '
    … …  
    10231023      out << "  OK" << endl;
    10241024      out.setRealNumberNotation(QTextStream::FixedNotation);
    1025       out.setRealNumberPrecision(3); 
     1025      out.setRealNumberPrecision(3);
    10261026      for (int ii = 0; ii < vv.Nrows(); ii++) {
    10271027        const cmbCorr* corr = corrs()[ii];
    1028         out << _resTime.datestr().c_str() << ' ' 
     1028        out << _resTime.datestr().c_str() << ' '
    10291029            << _resTime.timestr().c_str() << " "
    10301030            << corr->_acName << ' ' << corr->_prn.mid(0,3);
    … …  
    10921092        meanRao[prn].ReSize(4);
    10931093        meanRao[prn].Rows(1,3) = corr->_orbCorr._xr;
    1094         meanRao[prn](4)        = 1; 
     1094        meanRao[prn](4)        = 1;
    10951095      }
    10961096      else {
    10971097        meanRao[prn].Rows(1,3) += corr->_orbCorr._xr;
    1098         meanRao[prn](4)        += 1; 
     1098        meanRao[prn](4)        += 1;
    10991099      }
    11001100      if (numCorr.find(prn) == numCorr.end()) {
    … …  
    11051105      }
    11061106    }
    1107    
     1107
    11081108    // Compute Differences wrt Mean, find Maximum
    11091109    // ------------------------------------------
    … …  
    11271127          maxDiff[prn] = corr;
    11281128        }
    1129       } 
    1130     }
    1131    
     1129      }
     1130    }
     1131
    11321132    if (_ACs.size() == 1) {
    11331133      break;
    … …  
    11501150          out << _resTime.datestr().c_str()    << " "
    11511151              << _resTime.timestr().c_str()    << " "
    1152               << "Orbit Outlier: " 
    1153               << corr->_acName.toAscii().data() << " " 
     1152              << "Orbit Outlier: "
     1153              << corr->_acName.toAscii().data() << " "
    11541154              << prn.mid(0,3).toAscii().data()           << " "
    1155               << corr->_iod                     << " " 
     1155              << corr->_iod                     << " "
    11561156              << norm                           << endl;
    11571157          delete corr;
    … …  
    11611161      }
    11621162    }
    1163    
     1163
    11641164    if (!removed) {
    11651165      break;
    … …  
    11701170}
    11711171
    1172 // 
     1172//
    11731173////////////////////////////////////////////////////////////////////////////
    11741174void bncComb::slotProviderIDChanged(QString mountPoint) {
Note: See TracChangeset for help on using the changeset viewer.