Changeset 8542 in ntrip for trunk/BNC/src/ephemeris.cpp
- Timestamp:
- Nov 15, 2018, 11:16:46 AM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNC/src/ephemeris.cpp
r8489 r8542 65 65 << 900 << 1800 << 3600 << 7200 66 66 << 10800; 67 xc.ReSize( 7);67 xc.ReSize(6); 68 68 vv.ReSize(3); 69 69 if (position(tt.gpsw(), tt.gpssec(), xc.data(), vv.data()) != success) { … … 320 320 static const double gmGRS = 398.6005e12; 321 321 322 memset(xc, 0, 7*sizeof(double));322 memset(xc, 0, 6*sizeof(double)); 323 323 memset(vv, 0, 3*sizeof(double)); 324 324 … … 368 368 double tc = tt - _TOC; 369 369 xc[3] = _clock_bias + _clock_drift*tc + _clock_driftrate*tc*tc; 370 371 xc[4] = _clock_bias;372 xc[5] = _clock_drift;373 xc[6] = _clock_driftrate;374 370 375 371 // Velocity … … 402 398 xc[3] -= 2.0 * (xc[0]*vv[0] + xc[1]*vv[1] + xc[2]*vv[2]) / t_CST::c / t_CST::c; 403 399 400 xc[4] = _clock_drift + _clock_driftrate*tc; 401 xc[5] = _clock_driftrate; 402 404 403 return success; 405 404 } … … 618 617 static const double nominalStep = 10.0; 619 618 620 memset(xc, 0, 7*sizeof(double));619 memset(xc, 0, 6*sizeof(double)); 621 620 memset(vv, 0, 3*sizeof(double)); 622 621 … … 655 654 xc[3] = -_tau + _gamma * dtClk; 656 655 657 xc[4] = -_tau; 658 xc[5] = _gamma; 659 xc[6] = 0.0; 656 xc[4] = _gamma; 657 xc[5] = 0.0; 660 658 661 659 return success; … … 923 921 static const double gmWGS = 398.6004418e12; 924 922 925 memset(xc, 0, 7*sizeof(double));923 memset(xc, 0, 6*sizeof(double)); 926 924 memset(vv, 0, 3*sizeof(double)); 927 925 … … 971 969 double tc = tt - _TOC; 972 970 xc[3] = _clock_bias + _clock_drift*tc + _clock_driftrate*tc*tc; 973 974 xc[4] = _clock_bias;975 xc[5] = _clock_drift;976 xc[6] = _clock_driftrate;977 971 978 972 // Velocity … … 1007 1001 //xc[3] -= 2.0 * (xc[0]*vv[0] + xc[1]*vv[1] + xc[2]*vv[2]) / t_CST::c / t_CST::c; 1008 1002 1003 xc[4] = _clock_drift + _clock_driftrate*tc; 1004 xc[5] = _clock_driftrate; 1005 1009 1006 return success; 1010 1007 } … … 1319 1316 xc[3] = _agf0 + _agf1 * dt; 1320 1317 1321 xc[4] = _agf0; 1322 xc[5] = _agf1; 1323 xc[6] = 0.0; 1318 xc[4] = _agf1; 1319 xc[5] = 0.0; 1324 1320 1325 1321 return success; … … 1673 1669 xc[3] = _clock_bias + _clock_drift*tc + _clock_driftrate*tc*tc; 1674 1670 1675 xc[4] = _clock_bias;1676 xc[5] = _clock_drift;1677 xc[6] = _clock_driftrate;1678 1679 1671 // dotC = _clock_drift + _clock_driftrate*tc 1680 1672 // - 4.442807633e-10*_e*sqrt(a0)*cos(E) * dEdM * n; … … 1687 1679 // xc[3] -= 2.0 * (xc[0]*vv[0] + xc[1]*vv[1] + xc[2]*vv[2]) / t_CST::c / t_CST::c; 1688 1680 1681 xc[4] = _clock_drift + _clock_driftrate*tc; 1682 xc[5] = _clock_driftrate; 1689 1683 return success; 1690 1684 }
Note:
See TracChangeset
for help on using the changeset viewer.