Changeset 8903 in ntrip for trunk/BNC/src/bncutils.cpp
- Timestamp:
- Mar 18, 2020, 11:08:50 AM (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNC/src/bncutils.cpp
r8801 r8903 317 317 const ColumnVector& rsw, ColumnVector& xyz) { 318 318 319 ColumnVector along = vv / vv. norm_Frobenius();320 ColumnVector cross = crossproduct(rr, vv); cross /= cross. norm_Frobenius();319 ColumnVector along = vv / vv.NormFrobenius(); 320 ColumnVector cross = crossproduct(rr, vv); cross /= cross.NormFrobenius(); 321 321 ColumnVector radial = crossproduct(along, cross); 322 322 … … 334 334 const ColumnVector& xyz, ColumnVector& rsw) { 335 335 336 ColumnVector along = vv / vv. norm_Frobenius();337 ColumnVector cross = crossproduct(rr, vv); cross /= cross. norm_Frobenius();336 ColumnVector along = vv / vv.NormFrobenius(); 337 ColumnVector cross = crossproduct(rr, vv); cross /= cross.NormFrobenius(); 338 338 ColumnVector radial = crossproduct(along, cross); 339 339 … … 808 808 } 809 809 810 // 811 //////////////////////////////////////////////////////////////////////////// 810 812 double accuracyFromIndex(int index, t_eph::e_type type) { 811 813 double accuracy = -1.0; … … 869 871 } 870 872 873 // 874 //////////////////////////////////////////////////////////////////////////// 871 875 int indexFromAccuracy(double accuracy, t_eph::e_type type) { 872 876 … … 944 948 945 949 return (type == t_eph::Galileo) ? 255 : 15; 950 } 951 952 // Returns fit interval in hours from flag 953 //////////////////////////////////////////////////////////////////////////// 954 double fitIntervalFromFlag(int flag, double iodc, t_eph::e_type type) { 955 double fitInterval = 0.0; 956 957 switch (flag) { 958 case 0: 959 if (type == t_eph::GPS) { 960 fitInterval = 4.0; 961 } 962 else if (type == t_eph::QZSS) { 963 fitInterval = 2.0; 964 } 965 break; 966 case 1: 967 if (type == t_eph::GPS) { 968 if (iodc >= 240 && iodc <= 247) { 969 fitInterval = 8.0; 970 } 971 else if ((iodc >= 248 && iodc <= 255) || 972 (iodc == 496) ) { 973 fitInterval = 14.0; 974 } 975 else if ((iodc >= 497 && iodc <= 503) || 976 (iodc >= 2021 && iodc <= 1023) ) { 977 fitInterval = 26.0; 978 } 979 else { 980 fitInterval = 6.0; 981 } 982 } 983 break; 984 } 985 return fitInterval; 946 986 } 947 987
Note:
See TracChangeset
for help on using the changeset viewer.