Changeset 2246 in ntrip


Ignore:
Timestamp:
Jan 12, 2010, 5:22:33 PM (14 years ago)
Author:
mervart
Message:

* empty log message *

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/bncmodel.cpp

    r2245 r2246  
    103103  // ---------------
    104104  else if (type == RECCLK_GPS) {
    105     return 1.0;
     105    if (satData->prn[0] == 'G') {
     106      return 1.0;
     107    }
     108    else {
     109      return 0.0;
     110    }
    106111  }
    107112  else if (type == RECCLK_GLO) {
     
    662667    vv    = ll - AA * dx;
    663668
    664 //    //// beg test
    665 //    {
    666 //      cout.setf(ios::fixed);
    667 //      ColumnVector vv_code(epoData->sizeGPS());
    668 //      ColumnVector vv_phase(epoData->sizeGPS());
    669 //      ColumnVector vv_glo(epoData->sizeGlo());
    670 //
    671 //      for (int iobs = 1; iobs <= epoData->sizeGPS(); ++iobs) {
    672 //        vv_code(iobs)  = vv(2*iobs-1);
    673 //        vv_phase(iobs) = vv(2*iobs);
    674 //      }
    675 //      for (int iobs = 1; iobs <= epoData->sizeGlo(); ++iobs) {
    676 //        vv_glo(iobs)  = vv(2*epoData->sizeGPS()+iobs);
    677 //      }
    678 //
    679 //      cout << "residuals code  " << setprecision(3) << vv_code.t();
    680 //      cout << "residuals phase " << setprecision(3) << vv_phase.t();
    681 //      cout << "residuals glo   " << setprecision(3) << vv_glo.t();
    682 //    }
    683 //    //// end test
     669    //// beg test
     670    {
     671      cout.setf(ios::fixed);
     672      ColumnVector vv_code(epoData->sizeGPS());
     673      ColumnVector vv_phase(epoData->sizeGPS());
     674      ColumnVector vv_glo(epoData->sizeGlo());
     675
     676      for (unsigned iobs = 1; iobs <= epoData->sizeGPS(); ++iobs) {
     677        vv_code(iobs)  = vv(2*iobs-1);
     678        vv_phase(iobs) = vv(2*iobs);
     679      }
     680      for (unsigned iobs = 1; iobs <= epoData->sizeGlo(); ++iobs) {
     681        vv_glo(iobs)  = vv(2*epoData->sizeGPS()+iobs);
     682      }
     683
     684      cout << "residuals code  " << setprecision(3) << vv_code.t();
     685      cout << "residuals phase " << setprecision(3) << vv_phase.t();
     686      cout << "residuals glo   " << setprecision(3) << vv_glo.t();
     687    }
     688    //// end test
    684689
    685690  } while (outlierDetection(QQsav, vv, epoData->satDataGPS,
Note: See TracChangeset for help on using the changeset viewer.