Changeset 3439 in ntrip


Ignore:
Timestamp:
Sep 22, 2011, 10:36:25 AM (13 years ago)
Author:
mervart
Message:
 
Location:
trunk/BNC/combination
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/combination/bnccomb.cpp

    r3438 r3439  
    3333
    3434const double sig0_offAC    = 1000.0;
    35 const double sigP_offAC    = 1000.0;
    3635const double sig0_offACSat =  100.0;
    3736const double sigP_offACSat =    0.0;
     
    5554
    5655  if      (type == offAC) {
    57     sig_0 = sig0_offAC;
    58     sig_P = sigP_offAC;
     56    epoSpec = true;
     57    sig0    = sig0_offAC;
     58    sigP    = sig0;
    5959  }
    6060  else if (type == offACSat) {
    61     sig_0 = sig0_offACSat;
    62     sig_P = sigP_offACSat;
     61    epoSpec = false;
     62    sig0    = sig0_offACSat;
     63    sigP    = sigP_offACSat;
    6364  }
    6465  else if (type == clkSat) {
    65     sig_0 = sig0_clkSat;
    66     sig_P = sigP_clkSat;
     66    epoSpec = false;
     67    sig0    = sig0_clkSat;
     68    sigP    = sigP_clkSat;
    6769  }
    6870}
     
    165167  for (int iPar = 1; iPar <= _params.size(); iPar++) {
    166168    cmbParam* pp = _params[iPar-1];
    167     _QQ(iPar,iPar) = pp->sig_0 * pp->sig_0;
     169    _QQ(iPar,iPar) = pp->sig0 * pp->sig0;
    168170  }
    169171
     
    388390  for (int iPar = 1; iPar <= _params.size(); iPar++) {
    389391    cmbParam* pp = _params[iPar-1];
    390     _QQ(iPar,iPar) += pp->sig_P * pp->sig_P;
     392    if (pp->epoSpec) {
     393      pp->xx = 0.0;
     394      _QQ.Row(iPar)    = 0.0;
     395      _QQ.Column(iPar) = 0.0;
     396      _QQ(iPar,iPar) = pp->sig0 * pp->sig0;
     397    }
     398    else {
     399      _QQ(iPar,iPar) += pp->sigP * pp->sigP;
     400    }
    391401    x0(iPar) = pp->xx;
    392402  }
     
    477487          QQ_sav.Row(iPar)    = 0.0;
    478488          QQ_sav.Column(iPar) = 0.0;
    479           QQ_sav(iPar,iPar)   = pp->sig_0 * pp->sig_0;
     489          QQ_sav(iPar,iPar)   = pp->sig0 * pp->sig0;
    480490        }
    481491      }
  • trunk/BNC/combination/bnccomb.h

    r3438 r3439  
    2323  QString prn;
    2424  double  xx;
    25   double  sig_0;
    26   double  sig_P;
     25  double  sig0;
     26  double  sigP;
     27  bool    epoSpec;
    2728};
    2829
Note: See TracChangeset for help on using the changeset viewer.