Stats_T.cpp
Go to the documentation of this file.
1 //==============================================================================
2 //
3 // This file is part of GNSSTk, the ARL:UT GNSS Toolkit.
4 //
5 // The GNSSTk is free software; you can redistribute it and/or modify
6 // it under the terms of the GNU Lesser General Public License as published
7 // by the Free Software Foundation; either version 3.0 of the License, or
8 // any later version.
9 //
10 // The GNSSTk is distributed in the hope that it will be useful,
11 // but WITHOUT ANY WARRANTY; without even the implied warranty of
12 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
13 // GNU Lesser General Public License for more details.
14 //
15 // You should have received a copy of the GNU Lesser General Public
16 // License along with GNSSTk; if not, write to the Free Software Foundation,
17 // Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110, USA
18 //
19 // This software was developed by Applied Research Laboratories at the
20 // University of Texas at Austin.
21 // Copyright 2004-2022, The Board of Regents of The University of Texas System
22 //
23 //==============================================================================
24 
25 //==============================================================================
26 //
27 // This software was developed by Applied Research Laboratories at the
28 // University of Texas at Austin, under contract to an agency or agencies
29 // within the U.S. Department of Defense. The U.S. Government retains all
30 // rights to use, duplicate, distribute, disclose, or release this software.
31 //
32 // Pursuant to DoD Directive 523024
33 //
34 // DISTRIBUTION STATEMENT A: This software has been approved for public
35 // release, distribution is unlimited.
36 //
37 //==============================================================================
38 
40 
41 #include <vector>
42 #include "Vector.hpp"
43 #include "Stats.hpp"
44 #include "logstream.hpp"
45 
46 //------------------------------------------------------------------------------------
47 using namespace std;
48 using namespace gnsstk;
49 
50 //------------------------------------------------------------------------------------
51 int main(int argc, char **argv)
52 {
53 try {
54  // This include file defines
55  // const unsigned int M=734; double data[5*M] = { };
56  // which is 5 vectors of length M:
57  //for(i=0; i<M; i++) {
58  // sow[i] = data[5*i];
59  // L1[i] = data[5*i+1];
60  // L2[i] = data[5*i+2];
61  // P1[i] = data[5*i+3];
62  // P2[i] = data[5*i+4];
63  // L1-L2 = wl1*data[5*i+1] - wl2*data[5*i+2];
64  // L1-P1 = wl1*data[5*i+1] - data[5*i+3];
65  // L2-P2 = wl2*data[5*i+2] - data[5*i+4];
66  //}
67 
68  // This data comes from MSN data at StLouis, file s011072a.07o,
69  // satellite G01, first complete pass.
70  // The data is seconds of week, L1, L2, P1, P2, L in cycles, P in meters.
71  const unsigned int M=734;
72  double data[5*M] = {
73  86400.000, 123335818.734, 96105816.038, 23470039.004, 23470042.870,
74  86430.000, 123265156.115, 96050754.313, 23456591.718, 23456595.395,
75  86460.000, 123194673.332, 95995832.714, 23443178.907, 23443183.384,
76  86490.000, 123124370.538, 95941051.365, 23429800.786, 23429805.463,
77  86520.000, 123054247.657, 95886410.203, 23416457.467, 23416461.670,
78  86550.000, 122984304.730, 95831909.269, 23403147.839, 23403151.894,
79  86580.000, 122914541.707, 95777548.509, 23389871.985, 23389875.676,
80  86610.000, 122844958.703, 95723328.035, 23376630.318, 23376634.588,
81  86640.000, 122775555.760, 95669247.865, 23363423.840, 23363427.962,
82  86670.000, 122706332.734, 95615307.878, 23350251.545, 23350255.081,
83  86700.000, 122637289.983, 95561508.370, 23337112.954, 23337116.936,
84  86730.000, 122568427.349, 95507849.223, 23324008.454, 23324012.788,
85  86760.000, 122499744.875, 95454330.466, 23310938.555, 23310942.689,
86  86790.000, 122431242.589, 95400952.098, 23297903.060, 23297907.173,
87  86820.000, 122362920.352, 95347714.021, 23284901.758, 23284904.901,
88  86850.000, 122294778.318, 95294616.362, 23271934.552, 23271938.345,
89  86880.000, 122226816.352, 95241659.037, 23259001.663, 23259005.725,
90  86910.000, 122159034.652, 95188842.186, 23246103.080, 23246107.502,
91  86940.000, 122091433.256, 95136165.812, 23233239.621, 23233243.320,
92  86970.000, 122024012.029, 95083629.828, 23220409.796, 23220413.464,
93  87000.000, 121956770.904, 95031234.188, 23207613.948, 23207617.509,
94  87030.000, 121889709.873, 94978978.878, 23194852.197, 23194856.350,
95  87060.000, 121822829.087, 94926864.024, 23182125.608, 23182129.376,
96  87090.000, 121756128.227, 94874889.374, 23169433.515, 23169436.665,
97  87120.000, 121689607.354, 94823054.965, 23156774.799, 23156778.296,
98  87150.000, 121623266.648, 94771360.946, 23144150.226, 23144153.495,
99  87180.000, 121557105.814, 94719807.103, 23131559.984, 23131563.345,
100  87210.000, 121491124.876, 94668393.421, 23119004.323, 23119007.681,
101  87240.000, 121425323.835, 94617119.925, 23106482.905, 23106486.377,
102  87270.000, 121359702.442, 94565986.421, 23093995.156, 23093999.520,
103  87300.000, 121294260.899, 94514993.062, 23081541.646, 23081546.020,
104  87330.000, 121228999.036, 94464139.709, 23069122.632, 23069126.706,
105  87360.000, 121163916.542, 94413426.110, 23056738.239, 23056742.189,
106  87390.000, 121099013.729, 94362852.518, 23044387.541, 23044390.974,
107  87420.000, 121034290.339, 94312418.734, 23032071.001, 23032075.119,
108  87450.000, 120969746.257, 94262124.684, 23019788.634, 23019792.820,
109  87480.000, 120905381.503, 94211970.365, 23007540.387, 23007544.888,
110  87510.000, 120841196.134, 94161955.832, 22995326.573, 22995330.493,
111  87540.000, 120777190.031, 94112080.977, 22983146.449, 22983150.086,
112  87570.000, 120713362.959, 94062345.623, 22971000.419, 22971004.006,
113  87600.000, 120649714.905, 94012749.762, 22958888.432, 22958892.127,
114  87630.000, 120586245.703, 93963293.279, 22946810.686, 22946814.778,
115  87660.000, 120522955.365, 93913976.159, 22934767.203, 22934770.787,
116  87690.000, 120459843.656, 93864798.241, 22922757.114, 22922761.467,
117  87720.000, 120396910.616, 93815759.547, 22910781.490, 22910785.295,
118  87750.000, 120334156.093, 93766859.959, 22898839.704, 22898843.437,
119  87780.000, 120271580.124, 93718099.492, 22886931.935, 22886935.211,
120  87810.000, 120209182.662, 93669478.128, 22875058.489, 22875061.290,
121  87840.000, 120146963.278, 93620995.521, 22863218.290, 22863221.703,
122  87870.000, 120084921.911, 93572651.643, 22851411.776, 22851415.677,
123  87900.000, 120023058.595, 93524446.497, 22839639.431, 22839643.587,
124  87930.000, 119961373.257, 93476380.044, 22827901.338, 22827904.862,
125  87960.000, 119899865.631, 93428452.048, 22816196.982, 22816200.200,
126  87990.000, 119838535.745, 93380662.554, 22804526.360, 22804529.370,
127  88020.000, 119777383.399, 93333011.417, 22792889.143, 22792892.833,
128  88050.000, 119716408.460, 93285498.526, 22781286.333, 22781289.757,
129  88080.000, 119655610.682, 93238123.671, 22769716.416, 22769720.343,
130  88110.000, 119594990.234, 93190887.001, 22758181.238, 22758184.757,
131  88140.000, 119534546.886, 93143788.323, 22746678.996, 22746682.795,
132  88170.000, 119474280.531, 93096827.552, 22735210.590, 22735214.167,
133  88200.000, 119414190.899, 93050004.499, 22723775.772, 22723779.293,
134  88230.000, 119354277.815, 93003319.014, 22712374.780, 22712378.172,
135  88260.000, 119294541.355, 92956771.165, 22701007.058, 22701010.784,
136  88290.000, 119234981.220, 92910360.699, 22689673.555, 22689677.123,
137  88320.000, 119175597.336, 92864087.585, 22678372.676, 22678376.603,
138  88350.000, 119116389.460, 92817951.602, 22667106.105, 22667109.588,
139  88380.000, 119057357.399, 92771952.634, 22655872.539, 22655875.823,
140  88410.000, 118998501.181, 92726090.674, 22644672.486, 22644676.166,
141  88440.000, 118939820.527, 92680365.526, 22633505.646, 22633509.033,
142  88470.000, 118881315.355, 92634777.103, 22622372.721, 22622376.284,
143  88500.000, 118822985.597, 92589325.377, 22611273.083, 22611276.575,
144  88530.000, 118764831.063, 92544010.182, 22600206.665, 22600209.565,
145  88560.000, 118706851.712, 92498831.501, 22589173.549, 22589176.797,
146  88590.000, 118649047.336, 92453789.160, 22578173.788, 22578176.679,
147  88620.000, 118591417.649, 92408882.943, 22567206.913, 22567210.031,
148  88650.000, 118533962.330, 92364112.600, 22556273.765, 22556276.867,
149  88680.000, 118476681.373, 92319478.119, 22545373.341, 22545376.451,
150  88710.000, 118419574.763, 92274979.481, 22534506.552, 22534509.433,
151  88740.000, 118362642.308, 92230616.560, 22523672.169, 22523675.939,
152  88770.000, 118305883.852, 92186389.224, 22512871.659, 22512874.860,
153  88800.000, 118249299.289, 92142297.390, 22502103.823, 22502107.151,
154  88830.000, 118192888.331, 92098340.831, 22491369.020, 22491372.276,
155  88860.000, 118136650.959, 92054519.528, 22480667.861, 22480670.649,
156  88890.000, 118080586.919, 92010833.288, 22469999.027, 22470002.063,
157  88920.000, 118024695.902, 91967281.870, 22459363.320, 22459366.664,
158  88950.000, 117968977.764, 91923865.164, 22448760.360, 22448763.646,
159  88980.000, 117913432.555, 91880583.208, 22438190.625, 22438193.703,
160  89010.000, 117858059.902, 91837435.712, 22427653.501, 22427656.684,
161  89040.000, 117802859.793, 91794422.665, 22417149.023, 22417152.115,
162  89070.000, 117747832.328, 91751544.139, 22406677.906, 22406680.719,
163  89100.000, 117692977.034, 91708799.776, 22396239.215, 22396242.478,
164  89130.000, 117638293.727, 91666189.437, 22385833.397, 22385836.769,
165  89160.000, 117583782.204, 91623712.952, 22375460.089, 22375463.619,
166  89190.000, 117529442.407, 91581370.278, 22365119.606, 22365122.975,
167  89220.000, 117475274.092, 91539161.237, 22354811.578, 22354814.976,
168  89250.000, 117421277.069, 91497085.658, 22344536.302, 22344539.454,
169  89280.000, 117367451.297, 91455143.526, 22334293.348, 22334296.510,
170  89310.000, 117313796.330, 91413334.480, 22324083.516, 22324086.403,
171  89340.000, 117260312.171, 91371658.539, 22313905.526, 22313908.782,
172  89370.000, 117206998.564, 91330115.491, 22303759.901, 22303763.755,
173  89400.000, 117153855.499, 91288705.340, 22293646.996, 22293651.011,
174  89430.000, 117100882.857, 91247427.985, 22283566.769, 22283570.561,
175  89460.000, 117048080.380, 91206283.216, 22273519.045, 22273522.548,
176  89490.000, 116995447.967, 91165270.967, 22263503.337, 22263506.589,
177  89520.000, 116942985.367, 91124391.040, 22253519.956, 22253523.413,
178  89550.000, 116890692.454, 91083643.340, 22243568.938, 22243572.526,
179  89580.000, 116838569.213, 91043027.857, 22233650.186, 22233653.652,
180  89610.000, 116786615.470, 91002544.446, 22223763.697, 22223767.365,
181  89640.000, 116734830.933, 90962192.875, 22213909.685, 22213913.191,
182  89670.000, 116683215.455, 90921973.049, 22204087.323, 22204090.880,
183  89700.000, 116631768.764, 90881884.741, 22194297.346, 22194300.742,
184  89730.000, 116580490.905, 90841927.986, 22184539.644, 22184542.920,
185  89760.000, 116529381.585, 90802102.566, 22174813.728, 22174816.898,
186  89790.000, 116478440.787, 90762408.458, 22165120.144, 22165123.114,
187  89820.000, 116427668.346, 90722845.537, 22155458.327, 22155461.507,
188  89850.000, 116377063.984, 90683413.591, 22145828.756, 22145831.625,
189  89880.000, 116326627.740, 90644112.641, 22136230.852, 22136233.871,
190  89910.000, 116276359.421, 90604942.548, 22126664.840, 22126668.118,
191  89940.000, 116226258.710, 90565903.056, 22117131.047, 22117134.337,
192  89970.000, 116176325.612, 90526994.181, 22107629.005, 22107632.158,
193  90000.000, 116126559.703, 90488215.569, 22098159.081, 22098161.944,
194  90030.000, 116076961.130, 90449567.347, 22088721.018, 22088723.636,
195  90060.000, 116027529.539, 90411049.247, 22079314.361, 22079317.152,
196  90090.000, 115978264.948, 90372661.272, 22069939.570, 22069942.414,
197  90120.000, 115929167.074, 90334403.207, 22060596.535, 22060599.465,
198  90150.000, 115880235.913, 90296275.059, 22051285.035, 22051288.363,
199  90180.000, 115831471.382, 90258276.739, 22042005.624, 22042009.002,
200  90210.000, 115782873.275, 90220408.106, 22032757.680, 22032760.827,
201  90240.000, 115734441.495, 90182669.071, 22023541.499, 22023544.517,
202  90270.000, 115686175.860, 90145059.509, 22014356.789, 22014359.952,
203  90300.000, 115638076.152, 90107579.229, 22005203.773, 22005206.640,
204  90330.000, 115590142.372, 90070228.252, 21996082.151, 21996085.266,
205  90360.000, 115542374.377, 90033006.457, 21986992.266, 21986995.408,
206  90390.000, 115494771.924, 89995913.654, 21977933.654, 21977936.998,
207  90420.000, 115447334.855, 89958949.723, 21968906.941, 21968910.027,
208  90450.000, 115400063.397, 89922114.838, 21959911.622, 21959914.762,
209  90480.000, 115352957.128, 89885408.676, 21950947.336, 21950950.509,
210  90510.000, 115306015.881, 89848831.094, 21942014.807, 21942017.584,
211  90540.000, 115259239.632, 89812382.094, 21933113.360, 21933116.356,
212  90570.000, 115212628.360, 89776061.636, 21924243.421, 21924246.305,
213  90600.000, 115166182.098, 89739869.764, 21915405.064, 21915407.928,
214  90630.000, 115119900.599, 89703806.276, 21906598.186, 21906601.128,
215  90660.000, 115073783.617, 89667870.983, 21897822.452, 21897825.236,
216  90690.000, 115027831.142, 89632063.874, 21889077.875, 21889080.716,
217  90720.000, 114982043.284, 89596385.044, 21880364.567, 21880367.605,
218  90750.000, 114936419.750, 89560834.262, 21871682.428, 21871685.810,
219  90780.000, 114890960.493, 89525411.482, 21863031.802, 21863034.990,
220  90810.000, 114845665.328, 89490116.570, 21854412.422, 21854415.429,
221  90840.000, 114800534.219, 89454949.485, 21845824.782, 21845827.478,
222  90870.000, 114755567.009, 89419910.113, 21837267.709, 21837270.385,
223  90900.000, 114710763.660, 89384998.427, 21828741.806, 21828744.762,
224  90930.000, 114666124.204, 89350214.447, 21820247.349, 21820250.086,
225  90960.000, 114621648.379, 89315557.977, 21811783.443, 21811786.796,
226  90990.000, 114577336.296, 89281029.101, 21803351.056, 21803354.448,
227  91020.000, 114533187.590, 89246627.538, 21794949.784, 21794953.083,
228  91050.000, 114489202.406, 89212353.374, 21786580.042, 21786582.846,
229  91080.000, 114445380.721, 89178206.625, 21778240.979, 21778243.657,
230  91110.000, 114401722.420, 89144187.186, 21769933.053, 21769935.933,
231  91140.000, 114358227.329, 89110294.916, 21761656.330, 21761659.035,
232  91170.000, 114314895.665, 89076529.993, 21753410.448, 21753413.220,
233  91200.000, 114271727.259, 89042892.292, 21745195.653, 21745198.880,
234  91230.000, 114228722.132, 89009381.818, 21737012.128, 21737015.316,
235  91260.000, 114185880.381, 88975998.657, 21728859.652, 21728862.910,
236  91290.000, 114143201.388, 88942742.314, 21720738.111, 21720741.319,
237  91320.000, 114100685.559, 88909613.106, 21712647.693, 21712650.603,
238  91350.000, 114058332.731, 88876610.915, 21704588.394, 21704591.127,
239  91380.000, 114016142.875, 88843735.713, 21696559.676, 21696562.574,
240  91410.000, 113974115.991, 88810987.505, 21688562.169, 21688565.015,
241  91440.000, 113932251.998, 88778366.222, 21680595.537, 21680598.509,
242  91470.000, 113890550.914, 88745871.875, 21672660.160, 21672663.141,
243  91500.000, 113849012.879, 88713504.598, 21664755.591, 21664758.858,
244  91530.000, 113807637.629, 88681264.160, 21656882.382, 21656885.470,
245  91560.000, 113766425.370, 88649150.724, 21649039.760, 21649042.989,
246  91590.000, 113725375.943, 88617164.172, 21641228.456, 21641231.313,
247  91620.000, 113684489.279, 88585304.452, 21633447.719, 21633450.854,
248  91650.000, 113643765.457, 88553571.610, 21625698.306, 21625701.409,
249  91680.000, 113603204.642, 88521965.795, 21617979.741, 21617982.845,
250  91710.000, 113562806.768, 88490486.937, 21610292.230, 21610295.281,
251  91740.000, 113522571.835, 88459135.053, 21602635.749, 21602638.925,
252  91770.000, 113482499.780, 88427910.089, 21595010.399, 21595013.612,
253  91800.000, 113442590.607, 88396812.037, 21587416.178, 21587419.171,
254  91830.000, 113402844.256, 88365840.868, 21579852.670, 21579855.703,
255  91860.000, 113363260.970, 88334996.762, 21572320.123, 21572323.144,
256  91890.000, 113323840.757, 88304279.726, 21564818.853, 21564821.691,
257  91920.000, 113284583.684, 88273689.816, 21557348.381, 21557351.237,
258  91950.000, 113245489.756, 88243227.023, 21549909.179, 21549911.861,
259  91980.000, 113206559.122, 88212891.478, 21542500.885, 21542503.582,
260  92010.000, 113167791.752, 88182683.146, 21535123.710, 21535126.276,
261  92040.000, 113129187.566, 88152601.975, 21527777.460, 21527780.140,
262  92070.000, 113090746.772, 88122648.119, 21520462.462, 21520465.203,
263  92100.000, 113052469.256, 88092821.492, 21513178.580, 21513181.303,
264  92130.000, 113014355.302, 88063122.312, 21505925.819, 21505928.587,
265  92160.000, 112976404.863, 88033550.550, 21498704.157, 21498706.834,
266  92190.000, 112938618.105, 88004106.333, 21491513.565, 21491516.157,
267  92220.000, 112900995.094, 87974789.711, 21484354.039, 21484356.666,
268  92250.000, 112863535.983, 87945600.800, 21477225.865, 21477228.455,
269  92280.000, 112826240.741, 87916539.580, 21470128.869, 21470131.307,
270  92310.000, 112789109.657, 87887606.275, 21463062.865, 21463065.249,
271  92340.000, 112752142.579, 87858800.764, 21456028.307, 21456030.867,
272  92370.000, 112715339.676, 87830123.182, 21449024.954, 21449027.516,
273  92400.000, 112678701.148, 87801573.685, 21442052.798, 21442055.506,
274  92430.000, 112642227.038, 87773152.310, 21435112.030, 21435114.730,
275  92460.000, 112605917.613, 87744859.255, 21428202.542, 21428205.274,
276  92490.000, 112569772.858, 87716694.518, 21421324.270, 21421327.162,
277  92520.000, 112533792.998, 87688658.276, 21414477.615, 21414480.316,
278  92550.000, 112497978.007, 87660750.497, 21407662.306, 21407664.857,
279  92580.000, 112462328.167, 87632971.407, 21400878.480, 21400880.871,
280  92610.000, 112426843.584, 87605321.095, 21394125.998, 21394128.499,
281  92640.000, 112391524.369, 87577799.634, 21387405.038, 21387407.470,
282  92670.000, 112356370.648, 87550407.134, 21380715.347, 21380717.906,
283  92700.000, 112321382.768, 87523143.866, 21374057.156, 21374059.959,
284  92730.000, 112286560.812, 87496009.881, 21367430.753, 21367433.701,
285  92760.000, 112251904.919, 87469005.297, 21360836.038, 21360838.734,
286  92790.000, 112217415.264, 87442130.254, 21354272.889, 21354275.654,
287  92820.000, 112183092.073, 87415384.912, 21347741.502, 21347744.082,
288  92850.000, 112148935.454, 87388769.371, 21341241.672, 21341244.205,
289  92880.000, 112114945.554, 87362283.735, 21334773.670, 21334775.995,
290  92910.000, 112081122.699, 87335928.274, 21328337.098, 21328339.912,
291  92940.000, 112047466.905, 87309702.988, 21321932.617, 21321935.429,
292  92970.000, 112013978.401, 87283608.051, 21315559.916, 21315562.776,
293  93000.000, 111980657.644, 87257643.835, 21309219.272, 21309222.100,
294  93030.000, 111947504.680, 87231810.357, 21302910.481, 21302913.319,
295  93060.000, 111914519.827, 87206107.882, 21296633.944, 21296636.434,
296  93090.000, 111881703.555, 87180536.766, 21290389.222, 21290391.554,
297  93120.000, 111849055.788, 87155096.951, 21284176.385, 21284179.135,
298  93150.000, 111816576.637, 87129788.526, 21277995.726, 21277998.368,
299  93180.000, 111784266.399, 87104611.717, 21271847.128, 21271850.023,
300  93210.000, 111752125.405, 87079566.789, 21265730.923, 21265733.884,
301  93240.000, 111720153.985, 87054654.002, 21259646.780, 21259649.701,
302  93270.000, 111688352.340, 87029873.499, 21253595.336, 21253598.047,
303  93300.000, 111656720.822, 87005225.566, 21247576.148, 21247578.888,
304  93330.000, 111625259.618, 86980710.347, 21241589.269, 21241592.014,
305  93360.000, 111593968.846, 86956327.929, 21235635.009, 21235637.435,
306  93390.000, 111562848.820, 86932078.559, 21229713.071, 21229715.418,
307  93420.000, 111531899.872, 86907962.498, 21223823.605, 21223826.111,
308  93450.000, 111501122.326, 86883979.998, 21217966.793, 21217969.318,
309  93480.000, 111470516.614, 86860131.396, 21212142.695, 21212145.240,
310  93510.000, 111440082.876, 86836416.799, 21206351.266, 21206353.976,
311  93540.000, 111409821.498, 86812836.512, 21200592.554, 21200595.478,
312  93570.000, 111379732.815, 86789390.789, 21194866.878, 21194869.848,
313  93600.000, 111349817.127, 86766079.863, 21189174.260, 21189176.956,
314  93630.000, 111320074.694, 86742903.952, 21183514.598, 21183517.105,
315  93660.000, 111290505.928, 86719863.355, 21177887.923, 21177890.342,
316  93690.000, 111261111.175, 86696958.359, 21172294.326, 21172296.677,
317  93720.000, 111231890.741, 86674189.190, 21166733.882, 21166736.273,
318  93750.000, 111202844.958, 86651556.117, 21161206.575, 21161209.165,
319  93780.000, 111173974.150, 86629059.388, 21155712.560, 21155715.346,
320  93810.000, 111145278.704, 86606699.305, 21150251.847, 21150254.706,
321  93840.000, 111116758.881, 86584476.068, 21144824.743, 21144827.599,
322  93870.000, 111088415.024, 86562389.953, 21139430.856, 21139433.950,
323  93900.000, 111060247.684, 86540441.387, 21134070.828, 21134073.782,
324  93930.000, 111032257.118, 86518630.556, 21128744.490, 21128747.223,
325  93960.000, 111004443.894, 86496957.922, 21123451.869, 21123454.554,
326  93990.000, 110976808.374, 86475423.754, 21118193.172, 21118195.560,
327  94020.000, 110949350.920, 86454028.337, 21112968.114, 21112970.533,
328  94050.000, 110922071.800, 86432771.883, 21107777.203, 21107779.570,
329  94080.000, 110894971.468, 86411654.746, 21102620.086, 21102622.612,
330  94110.000, 110868050.302, 86390677.216, 21097497.093, 21097499.741,
331  94140.000, 110841308.818, 86369839.698, 21092408.201, 21092410.983,
332  94170.000, 110814747.654, 86349142.694, 21087353.537, 21087356.444,
333  94200.000, 110788366.897, 86328586.264, 21082333.532, 21082336.543,
334  94230.000, 110762167.227, 86308170.939, 21077347.675, 21077350.904,
335  94260.000, 110736149.021, 86287897.012, 21072396.649, 21072399.812,
336  94290.000, 110710312.709, 86267764.823, 21067480.169, 21067483.530,
337  94320.000, 110684658.533, 86247774.558, 21062598.415, 21062601.583,
338  94350.000, 110659187.213, 86227926.778, 21057751.553, 21057754.404,
339  94380.000, 110633899.085, 86208221.742, 21052939.411, 21052942.317,
340  94410.000, 110608794.607, 86188659.809, 21048162.374, 21048165.139,
341  94440.000, 110583874.246, 86169241.344, 21043420.079, 21043422.978,
342  94470.000, 110559138.511, 86149966.749, 21038712.965, 21038715.931,
343  94500.000, 110534587.906, 86130836.411, 21034041.107, 21034044.104,
344  94530.000, 110510222.902, 86111850.689, 21029404.568, 21029407.465,
345  94560.000, 110486044.007, 86093009.999, 21024803.425, 21024806.266,
346  94590.000, 110462051.674, 86074314.674, 21020237.762, 21020240.679,
347  94620.000, 110438246.434, 86055765.144, 21015707.721, 21015710.540,
348  94650.000, 110414628.715, 86037361.730, 21011213.533, 21011216.221,
349  94680.000, 110391199.074, 86019104.866, 21006754.987, 21006757.699,
350  94710.000, 110367958.275, 86000995.153, 21002332.504, 21002335.372,
351  94740.000, 110344906.678, 85983032.867, 20997945.916, 20997948.755,
352  94770.000, 110322044.698, 85965218.344, 20993595.588, 20993598.238,
353  94800.000, 110299373.009, 85947552.093, 20989281.314, 20989284.080,
354  94830.000, 110276892.012, 85930034.437, 20985003.286, 20985006.131,
355  94860.000, 110254602.384, 85912665.897, 20980761.668, 20980764.490,
356  94890.000, 110232504.479, 85895446.751, 20976556.568, 20976559.329,
357  94920.000, 110210598.918, 85878377.490, 20972387.930, 20972390.935,
358  94950.000, 110188886.299, 85861458.561, 20968256.116, 20968259.150,
359  94980.000, 110167367.282, 85844690.506, 20964161.102, 20964164.220,
360  95010.000, 110146042.525, 85828073.809, 20960103.087, 20960106.188,
361  95040.000, 110124912.532, 85811608.883, 20956082.227, 20956085.324,
362  95070.000, 110103977.737, 85795296.053, 20952098.601, 20952101.562,
363  95100.000, 110083238.805, 85779135.839, 20948152.098, 20948154.980,
364  95130.000, 110062696.312, 85763128.700, 20944243.094, 20944245.866,
365  95160.000, 110042351.085, 85747275.280, 20940371.610, 20940374.331,
366  95190.000, 110022203.438, 85731575.818, 20936537.604, 20936540.446,
367  95220.000, 110002254.134, 85716030.904, 20932741.458, 20932744.261,
368  95250.000, 109982503.789, 85700641.032, 20928983.104, 20928985.850,
369  95280.000, 109962952.881, 85685406.559, 20925262.807, 20925265.444,
370  95310.000, 109943602.086, 85670328.024, 20921580.397, 20921583.054,
371  95340.000, 109924452.214, 85655406.050, 20917936.293, 20917938.909,
372  95370.000, 109905503.590, 85640640.892, 20914330.503, 20914333.147,
373  95400.000, 109886757.130, 85626033.261, 20910763.191, 20910765.780,
374  95430.000, 109868213.467, 85611583.657, 20907234.386, 20907236.940,
375  95460.000, 109849873.279, 85597292.607, 20903744.235, 20903747.044,
376  95490.000, 109831737.163, 85583160.580, 20900292.918, 20900295.778,
377  95520.000, 109813805.664, 85569187.988, 20896880.770, 20896883.547,
378  95550.000, 109796079.547, 85555375.430, 20893507.528, 20893510.376,
379  95580.000, 109778559.377, 85541723.350, 20890173.816, 20890176.496,
380  95610.000, 109761245.858, 85528232.298, 20886879.127, 20886881.669,
381  95640.000, 109744139.612, 85514902.753, 20883623.942, 20883626.499,
382  95670.000, 109727241.373, 85501735.296, 20880408.370, 20880410.797,
383  95700.000, 109710551.738, 85488730.387, 20877232.487, 20877234.994,
384  95730.000, 109694071.508, 85475888.647, 20874096.303, 20874098.644,
385  95760.000, 109677801.582, 85463210.774, 20871000.327, 20871002.763,
386  95790.000, 109661742.440, 85450697.160, 20867944.328, 20867946.910,
387  95820.000, 109645894.678, 85438348.251, 20864928.775, 20864931.177,
388  95850.000, 109630259.299, 85426164.834, 20861953.509, 20861955.872,
389  95880.000, 109614836.759, 85414147.272, 20859018.631, 20859021.207,
390  95910.000, 109599627.823, 85402296.151, 20856124.336, 20856127.045,
391  95940.000, 109584633.185, 85390612.011, 20853271.106, 20853273.621,
392  95970.000, 109569853.706, 85379095.526, 20850458.382, 20850461.188,
393  96000.000, 109555289.997, 85367747.165, 20847686.904, 20847689.876,
394  96030.000, 109540942.761, 85356567.497, 20844956.655, 20844959.536,
395  96060.000, 109526812.654, 85345557.021, 20842267.720, 20842270.751,
396  96090.000, 109512900.582, 85334716.438, 20839620.370, 20839623.366,
397  96120.000, 109499207.330, 85324046.365, 20837014.626, 20837017.768,
398  96150.000, 109485733.578, 85313547.338, 20834450.633, 20834453.734,
399  96180.000, 109472480.034, 85303219.900, 20831928.441, 20831931.542,
400  96210.000, 109459447.535, 85293064.698, 20829448.559, 20829451.644,
401  96240.000, 109446636.731, 85283082.248, 20827010.745, 20827013.872,
402  96270.000, 109434048.379, 85273273.137, 20824615.327, 20824618.283,
403  96300.000, 109421683.332, 85263638.035, 20822262.347, 20822265.321,
404  96330.000, 109409542.421, 85254177.575, 20819952.145, 20819954.942,
405  96360.000, 109397626.292, 85244892.273, 20817684.551, 20817687.331,
406  96390.000, 109385935.793, 85235782.796, 20815459.928, 20815462.687,
407  96420.000, 109374471.793, 85226849.809, 20813278.447, 20813281.191,
408  96450.000, 109363234.881, 85218093.777, 20811140.273, 20811142.957,
409  96480.000, 109352225.954, 85209515.395, 20809045.329, 20809048.128,
410  96510.000, 109341445.659, 85201115.161, 20806993.942, 20806996.523,
411  96540.000, 109330894.979, 85192893.845, 20804986.146, 20804988.912,
412  96570.000, 109320574.521, 85184851.933, 20803022.077, 20803024.893,
413  96600.000, 109310485.241, 85176990.153, 20801102.421, 20801104.998,
414  96630.000, 109300627.847, 85169309.069, 20799226.634, 20799229.009,
415  96660.000, 109291003.254, 85161809.379, 20797395.107, 20797397.544,
416  96690.000, 109281612.216, 85154491.676, 20795608.107, 20795610.568,
417  96720.000, 109272455.567, 85147356.619, 20793865.832, 20793868.097,
418  96750.000, 109263534.001, 85140404.735, 20792168.113, 20792170.461,
419  96780.000, 109254848.443, 85133636.760, 20790515.338, 20790517.633,
420  96810.000, 109246399.711, 85127053.322, 20788907.399, 20788909.924,
421  96840.000, 109238188.729, 85120655.145, 20787344.912, 20787347.545,
422  96870.000, 109230216.267, 85114442.828, 20785828.009, 20785830.448,
423  96900.000, 109222483.008, 85108416.907, 20784356.325, 20784358.773,
424  96930.000, 109214990.074, 85102578.254, 20782930.392, 20782932.857,
425  96960.000, 109207738.149, 85096927.398, 20781550.488, 20781552.965,
426  96990.000, 109200727.842, 85091464.818, 20780216.302, 20780218.921,
427  97020.000, 109193960.322, 85086191.425, 20778928.553, 20778931.155,
428  97050.000, 109187436.004, 85081107.531, 20777687.000, 20777689.618,
429  97080.000, 109181156.050, 85076214.052, 20776491.938, 20776494.562,
430  97110.000, 109175121.429, 85071511.743, 20775343.657, 20775346.148,
431  97140.000, 109169332.785, 85067001.111, 20774242.176, 20774244.645,
432  97170.000, 109163790.927, 85062682.771, 20773187.497, 20773190.109,
433  97200.000, 109158496.855, 85058557.519, 20772180.083, 20772182.651,
434  97230.000, 109153451.498, 85054626.072, 20771220.026, 20771222.547,
435  97260.000, 109148655.350, 85050888.810, 20770307.261, 20770309.894,
436  97290.000, 109144109.624, 85047346.686, 20769442.209, 20769444.919,
437  97320.000, 109139814.924, 85044000.168, 20768624.917, 20768627.623,
438  97350.000, 109135772.143, 85040849.941, 20767855.654, 20767858.365,
439  97380.000, 109131982.196, 85037896.750, 20767134.486, 20767137.139,
440  97410.000, 109128446.088, 85035141.370, 20766461.435, 20766464.168,
441  97440.000, 109125164.668, 85032584.428, 20765836.873, 20765839.679,
442  97470.000, 109122138.479, 85030226.370, 20765261.039, 20765263.701,
443  97500.000, 109119368.766, 85028068.153, 20764733.979, 20764736.519,
444  97530.000, 109116856.119, 85026110.254, 20764256.036, 20764258.537,
445  97560.000, 109114601.741, 85024353.612, 20763826.936, 20763829.579,
446  97590.000, 109112606.160, 85022798.624, 20763447.185, 20763449.809,
447  97620.000, 109110870.370, 85021446.066, 20763116.717, 20763119.675,
448  97650.000, 109109395.087, 85020296.479, 20762836.008, 20762838.817,
449  97680.000, 109108181.513, 85019350.830, 20762605.274, 20762607.800,
450  97710.000, 109107230.269, 85018609.581, 20762424.195, 20762426.804,
451  97740.000, 109106542.425, 85018073.581, 20762293.342, 20762295.928,
452  97770.000, 109106118.778, 85017743.477, 20762212.570, 20762215.370,
453  97800.000, 109105960.000, 85017619.746, 20762182.311, 20762185.200,
454  97830.000, 109106067.160, 85017703.245, 20762202.946, 20762205.544,
455  97860.000, 109106441.228, 85017994.735, 20762274.075, 20762276.761,
456  97890.000, 109107082.817, 85018494.678, 20762396.233, 20762398.896,
457  97920.000, 109107992.921, 85019203.858, 20762569.458, 20762572.042,
458  97950.000, 109109172.573, 85020123.062, 20762793.866, 20762796.475,
459  97980.000, 109110622.537, 85021252.881, 20763069.746, 20763072.442,
460  98010.000, 109112343.831, 85022594.140, 20763397.178, 20763399.829,
461  98040.000, 109114337.141, 85024147.368, 20763776.530, 20763779.212,
462  98070.000, 109116603.456, 85025913.316, 20764207.780, 20764210.544,
463  98100.000, 109119143.524, 85027892.587, 20764691.263, 20764694.042,
464  98130.000, 109121958.449, 85030086.038, 20765226.956, 20765229.636,
465  98160.000, 109125048.762, 85032494.079, 20765815.021, 20765817.758,
466  98190.000, 109128415.760, 85035117.726, 20766455.660, 20766458.457,
467  98220.000, 109132059.905, 85037957.330, 20767149.096, 20767151.986,
468  98250.000, 109135982.331, 85041013.774, 20767895.593, 20767898.398,
469  98280.000, 109140183.880, 85044287.707, 20768695.049, 20768697.889,
470  98310.000, 109144665.348, 85047779.761, 20769547.901, 20769550.627,
471  98340.000, 109149427.542, 85051490.567, 20770454.124, 20770456.779,
472  98370.000, 109154471.525, 85055420.940, 20771413.976, 20771416.727,
473  98400.000, 109159798.139, 85059571.545, 20772427.484, 20772430.360,
474  98430.000, 109165408.043, 85063942.899, 20773495.093, 20773497.852,
475  98460.000, 109171302.206, 85068535.756, 20774616.626, 20774619.365,
476  98490.000, 109177481.554, 85073350.831, 20775792.640, 20775795.357,
477  98520.000, 109183946.896, 85078388.752, 20777022.939, 20777025.476,
478  98550.000, 109190699.303, 85083650.368, 20778307.806, 20778310.563,
479  98580.000, 109197739.412, 85089136.163, 20779647.522, 20779650.339,
480  98610.000, 109205068.150, 85094846.849, 20781042.071, 20781044.889,
481  98640.000, 109212686.398, 85100783.135, 20782491.778, 20782494.474,
482  98670.000, 109220595.060, 85106945.727, 20783996.704, 20783999.485,
483  98700.000, 109228795.016, 85113335.305, 20785557.178, 20785559.919,
484  98730.000, 109237286.999, 85119952.426, 20787173.156, 20787175.977,
485  98760.000, 109246072.078, 85126797.936, 20788844.840, 20788847.775,
486  98790.000, 109255150.911, 85133872.345, 20790572.563, 20790575.378,
487  98820.000, 109264524.499, 85141176.433, 20792356.264, 20792359.102,
488  98850.000, 109274193.292, 85148710.554, 20794196.331, 20794199.058,
489  98880.000, 109284158.577, 85156475.701, 20796092.602, 20796095.354,
490  98910.000, 109294420.973, 85164472.372, 20798045.480, 20798048.228,
491  98940.000, 109304981.478, 85172701.330, 20800055.026, 20800057.787,
492  98970.000, 109315840.742, 85181163.085, 20802121.464, 20802124.213,
493  99000.000, 109326999.618, 85189858.306, 20804244.964, 20804247.806,
494  99030.000, 109338458.963, 85198787.661, 20806425.533, 20806428.522,
495  99060.000, 109350219.664, 85207951.842, 20808663.663, 20808666.488,
496  99090.000, 109362282.488, 85217351.439, 20810959.197, 20810961.947,
497  99120.000, 109374648.216, 85226987.072, 20813312.295, 20813314.999,
498  99150.000, 109387317.866, 85236859.531, 20815723.103, 20815725.986,
499  99180.000, 109400292.229, 85246969.423, 20818191.983, 20818194.901,
500  99210.000, 109413571.982, 85257317.281, 20820719.070, 20820722.090,
501  99240.000, 109427157.924, 85267903.725, 20823304.453, 20823307.274,
502  99270.000, 109441050.894, 85278729.409, 20825948.160, 20825951.036,
503  99300.000, 109455251.852, 85289795.083, 20828650.531, 20828653.496,
504  99330.000, 109469761.440, 85301101.246, 20831411.650, 20831414.598,
505  99360.000, 109484580.478, 85312648.537, 20834231.672, 20834234.554,
506  99390.000, 109499709.760, 85324437.579, 20837110.685, 20837113.482,
507  99420.000, 109515150.224, 85336469.106, 20840048.833, 20840051.689,
508  99450.000, 109530902.517, 85348743.611, 20843046.412, 20843049.357,
509  99480.000, 109546967.282, 85361261.600, 20846103.532, 20846106.398,
510  99510.000, 109563345.534, 85374023.866, 20849220.037, 20849223.022,
511  99540.000, 109580038.021, 85387030.991, 20852396.677, 20852399.551,
512  99570.000, 109597045.352, 85400283.442, 20855633.025, 20855636.021,
513  99600.000, 109614368.450, 85413781.951, 20858929.493, 20858932.310,
514  99630.000, 109632007.869, 85427526.946, 20862286.158, 20862289.142,
515  99660.000, 109649964.550, 85441519.157, 20865703.301, 20865706.233,
516  99690.000, 109668239.039, 85455759.010, 20869180.989, 20869183.587,
517  99720.000, 109686832.394, 85470247.319, 20872719.074, 20872721.906,
518  99750.000, 109705745.186, 85484984.544, 20876317.965, 20876320.907,
519  99780.000, 109724978.056, 85499971.191, 20879977.708, 20879980.793,
520  99810.000, 109744532.007, 85515208.049, 20883698.898, 20883701.810,
521  99840.000, 109764407.486, 85530695.419, 20887481.060, 20887484.120,
522  99870.000, 109784605.522, 85546434.135, 20891324.603, 20891327.514,
523  99900.000, 109805126.766, 85562424.709, 20895229.527, 20895232.680,
524  99930.000, 109825971.656, 85578667.470, 20899196.224, 20899199.214,
525  99960.000, 109847141.017, 85595163.072, 20903224.729, 20903227.704,
526  99990.000, 109868635.892, 85611912.339, 20907315.031, 20907318.067,
527  100020.000, 109890456.614, 85628915.486, 20911467.590, 20911470.379,
528  100050.000, 109912603.869, 85646173.073, 20915681.958, 20915684.735,
529  100080.000, 109935078.478, 85663685.749, 20919958.823, 20919961.534,
530  100110.000, 109957880.998, 85681453.944, 20924297.722, 20924300.794,
531  100140.000, 109981012.168, 85699478.224, 20928699.519, 20928702.398,
532  100170.000, 110004472.883, 85717759.304, 20933163.999, 20933166.777,
533  100200.000, 110028263.493, 85736297.439, 20937691.185, 20937694.098,
534  100230.000, 110052384.783, 85755093.246, 20942281.289, 20942284.144,
535  100260.000, 110076837.257, 85774147.117, 20946934.417, 20946937.322,
536  100290.000, 110101621.623, 85793459.608, 20951650.852, 20951653.705,
537  100320.000, 110126738.613, 85813031.285, 20956430.475, 20956433.242,
538  100350.000, 110152188.852, 85832862.635, 20961273.462, 20961276.377,
539  100380.000, 110177972.883, 85852954.093, 20966179.937, 20966182.833,
540  100410.000, 110204091.348, 85873306.141, 20971150.116, 20971152.976,
541  100440.000, 110230544.843, 85893919.246, 20976184.173, 20976186.876,
542  100470.000, 110257333.974, 85914793.891, 20981281.848, 20981284.686,
543  100500.000, 110284459.261, 85935930.464, 20986443.651, 20986446.476,
544  100530.000, 110311921.449, 85957329.572, 20991669.461, 20991672.516,
545  100560.000, 110339721.184, 85978991.709, 20996959.593, 20996962.583,
546  100590.000, 110367858.770, 86000917.100, 21002314.012, 21002317.088,
547  100620.000, 110396334.895, 86023106.277, 21007732.906, 21007735.932,
548  100650.000, 110425150.149, 86045559.719, 21013216.345, 21013219.328,
549  100680.000, 110454305.070, 86068277.834, 21018764.239, 21018767.347,
550  100710.000, 110483800.298, 86091261.121, 21024377.053, 21024380.176,
551  100740.000, 110513636.175, 86114509.847, 21030054.630, 21030057.789,
552  100770.000, 110543813.328, 86138024.511, 21035797.268, 21035800.296,
553  100800.000, 110574332.241, 86161805.476, 21041604.724, 21041607.820,
554  100830.000, 110605193.504, 86185853.204, 21047477.436, 21047480.564,
555  100860.000, 110636397.414, 86210167.931, 21053415.327, 21053418.446,
556  100890.000, 110667944.783, 86234750.292, 21059418.614, 21059421.660,
557  100920.000, 110699836.022, 86259600.589, 21065487.227, 21065490.445,
558  100950.000, 110732071.387, 86284719.039, 21071621.445, 21071624.552,
559  100980.000, 110764651.717, 86310106.298, 21077821.208, 21077824.337,
560  101010.000, 110797577.200, 86335762.519, 21084086.707, 21084089.793,
561  101040.000, 110830848.736, 86361688.369, 21090418.155, 21090421.324,
562  101070.000, 110864466.467, 86387883.997, 21096815.457, 21096818.626,
563  101100.000, 110898430.860, 86414349.740, 21103278.663, 21103281.745,
564  101130.000, 110932742.345, 86441085.956, 21109807.883, 21109810.948,
565  101160.000, 110967401.229, 86468092.869, 21116403.263, 21116406.294,
566  101190.000, 111002408.098, 86495370.947, 21123064.795, 21123067.839,
567  101220.000, 111037763.272, 86522920.428, 21129792.750, 21129795.734,
568  101250.000, 111073467.398, 86550741.812, 21136586.912, 21136590.171,
569  101280.000, 111109520.655, 86578835.256, 21143447.580, 21143450.757,
570  101310.000, 111145923.508, 86607201.118, 21150374.798, 21150377.863,
571  101340.000, 111182676.251, 86635839.609, 21157368.650, 21157371.744,
572  101370.000, 111219779.249, 86664751.025, 21164429.336, 21164432.288,
573  101400.000, 111257232.918, 86693935.687, 21171556.630, 21171559.413,
574  101430.000, 111295037.606, 86723393.883, 21178750.469, 21178753.482,
575  101460.000, 111333193.600, 86753125.819, 21186011.330, 21186014.417,
576  101490.000, 111371701.296, 86783131.813, 21193339.203, 21193342.181,
577  101520.000, 111410560.993, 86813412.083, 21200733.985, 21200736.952,
578  101550.000, 111449772.931, 86843966.818, 21208195.866, 21208198.734,
579  101580.000, 111489337.665, 86874796.478, 21215724.648, 21215727.649,
580  101610.000, 111529255.256, 86905901.084, 21223320.847, 21223323.743,
581  101640.000, 111569525.970, 86937280.865, 21230984.088, 21230987.098,
582  101670.000, 111610150.317, 86968936.193, 21238714.620, 21238717.398,
583  101700.000, 111651128.234, 87000867.029, 21246512.315, 21246515.401,
584  101730.000, 111692460.194, 87033073.743, 21254377.514, 21254380.402,
585  101760.000, 111734146.493, 87065556.573, 21262310.171, 21262313.067,
586  101790.000, 111776187.381, 87098315.698, 21270310.435, 21270313.200,
587  101820.000, 111818582.744, 87131351.037, 21278377.820, 21278380.961,
588  101850.000, 111861333.123, 87164663.009, 21286513.028, 21286515.976,
589  101880.000, 111904438.833, 87198251.863, 21294715.793, 21294718.866,
590  101910.000, 111947899.881, 87232117.613, 21302986.142, 21302989.146,
591  101940.000, 111991716.368, 87266260.320, 21311324.228, 21311327.130,
592  101970.000, 112035888.650, 87300680.269, 21319729.934, 21319732.996,
593  102000.000, 112080416.600, 87335377.353, 21328203.263, 21328206.454,
594  102030.000, 112125300.597, 87370351.876, 21336744.360, 21336747.598,
595  102060.000, 112170540.874, 87405604.027, 21345353.484, 21345356.654,
596  102090.000, 112216137.508, 87441133.849, 21354030.141, 21354033.341,
597  102120.000, 112262090.693, 87476941.507, 21362774.946, 21362778.128,
598  102150.000, 112308400.334, 87513026.931, 21371587.354, 21371590.456,
599  102180.000, 112355066.744, 87549390.347, 21380467.535, 21380470.683,
600  102210.000, 112402090.035, 87586031.855, 21389415.858, 21389419.045,
601  102240.000, 112449470.149, 87622951.407, 21398431.982, 21398435.003,
602  102270.000, 112497207.369, 87660149.211, 21407516.078, 21407519.167,
603  102300.000, 112545301.413, 87697625.077, 21416668.162, 21416670.956,
604  102330.000, 112593752.666, 87735379.287, 21425888.119, 21425890.923,
605  102360.000, 112642561.200, 87773411.911, 21435176.042, 21435179.110,
606  102390.000, 112691726.794, 87811722.741, 21444531.854, 21444534.961,
607  102420.000, 112741249.673, 87850311.982, 21453955.780, 21453958.955,
608  102450.000, 112791129.810, 87889179.607, 21463447.886, 21463450.732,
609  102480.000, 112841367.130, 87928325.558, 21473007.657, 21473010.463,
610  102510.000, 112891961.636, 87967749.830, 21482635.479, 21482638.306,
611  102540.000, 112942913.318, 88007452.424, 21492331.170, 21492334.218,
612  102570.000, 112994222.068, 88047433.239, 21502094.960, 21502098.082,
613  102600.000, 113045887.887, 88087692.291, 21511926.576, 21511929.740,
614  102630.000, 113097910.738, 88128229.560, 21521826.217, 21521829.432,
615  102660.000, 113150290.545, 88169044.972, 21531793.947, 21531797.203,
616  102690.000, 113203027.442, 88210138.617, 21541829.418, 21541832.548,
617  102720.000, 113256121.099, 88251510.287, 21551932.819, 21551935.839,
618  102750.000, 113309571.677, 88293160.073, 21562103.963, 21562107.179,
619  102780.000, 113363378.689, 88335087.601, 21572343.278, 21572346.438,
620  102810.000, 113417542.229, 88377292.935, 21582650.286, 21582653.359,
621  102840.000, 113472062.016, 88419775.855, 21593024.954, 21593028.296,
622  102870.000, 113526938.032, 88462536.373, 21603467.429, 21603470.934,
623  102900.000, 113582170.201, 88505574.402, 21613977.732, 21613981.396,
624  102930.000, 113637758.298, 88548889.784, 21624555.873, 21624559.295,
625  102960.000, 113693702.069, 88592482.311, 21635201.933, 21635205.150,
626  102990.000, 113750001.400, 88636351.895, 21645915.341, 21645918.490,
627  103020.000, 113806655.950, 88680498.258, 21656696.238, 21656699.563,
628  103050.000, 113863665.995, 88724921.646, 21667545.066, 21667548.167,
629  103080.000, 113921030.918, 88769621.555, 21678461.070, 21678464.428,
630  103110.000, 113978750.804, 88814598.071, 21689444.816, 21689448.058,
631  103140.000, 114036825.146, 88859850.784, 21700495.942, 21700499.197,
632  103170.000, 114095253.894, 88905379.661, 21711614.747, 21711617.745,
633  103200.000, 114154036.576, 88951184.316, 21722800.727, 21722803.816,
634  103230.000, 114213172.964, 88997264.623, 21734053.935, 21734057.185,
635  103260.000, 114272662.684, 89043620.231, 21745374.526, 21745377.675,
636  103290.000, 114332505.441, 89090250.931, 21756762.180, 21756765.541,
637  103320.000, 114392701.028, 89137156.561, 21768217.345, 21768220.445,
638  103350.000, 114453249.015, 89184336.802, 21779739.256, 21779742.436,
639  103380.000, 114514149.046, 89231791.341, 21791328.043, 21791331.273,
640  103410.000, 114575400.955, 89279520.065, 21802983.830, 21802987.169,
641  103440.000, 114637004.264, 89327522.617, 21814706.246, 21814709.785,
642  103470.000, 114698958.673, 89375798.742, 21826495.967, 21826499.333,
643  103500.000, 114761263.729, 89424348.101, 21838352.302, 21838355.652,
644  103530.000, 114823919.190, 89473170.521, 21850275.463, 21850278.669,
645  103560.000, 114886924.511, 89522265.529, 21862265.035, 21862267.988,
646  103590.000, 114950279.438, 89571632.987, 21874321.156, 21874324.059,
647  103620.000, 115013983.531, 89621272.493, 21886443.434, 21886446.521,
648  103650.000, 115078036.106, 89671183.565, 21898632.199, 21898635.476,
649  103680.000, 115142436.994, 89721366.034, 21910887.387, 21910890.602,
650  103710.000, 115207185.511, 89771819.377, 21923208.735, 21923211.928,
651  103740.000, 115272281.374, 89822543.403, 21935596.118, 21935599.372,
652  103770.000, 115337723.884, 89873537.537, 21948049.593, 21948052.401,
653  103800.000, 115403512.669, 89924801.518, 21960568.478, 21960571.605,
654  103830.000, 115469647.184, 89976334.874, 21973153.439, 21973156.850,
655  103860.000, 115536126.950, 90028137.247, 21985804.160, 21985807.526,
656  103890.000, 115602951.375, 90080208.188, 21998520.642, 21998523.808,
657  103920.000, 115670119.893, 90132547.271, 22011302.423, 22011305.668,
658  103950.000, 115737631.984, 90185154.070, 22024149.551, 22024152.728,
659  103980.000, 115805487.185, 90238028.223, 22037061.847, 22037065.188,
660  104010.000, 115873684.886, 90291169.263, 22050039.400, 22050042.784,
661  104040.000, 115942224.581, 90344576.774, 22063082.187, 22063085.762,
662  104070.000, 116011105.588, 90398250.254, 22076189.833, 22076193.378,
663  104100.000, 116080327.358, 90452189.257, 22089362.246, 22089365.925,
664  104130.000, 116149889.042, 90506393.123, 22102599.595, 22102602.869,
665  104160.000, 116219790.221, 90560861.535, 22115901.048, 22115904.681,
666  104190.000, 116290030.303, 90615594.020, 22129267.395, 22129270.808,
667  104220.000, 116360608.494, 90670589.955, 22142697.965, 22142701.238,
668  104250.000, 116431524.208, 90725848.919, 22156193.254, 22156196.197,
669  104280.000, 116502776.622, 90781370.232, 22169752.265, 22169755.370,
670  104310.000, 116574365.188, 90837153.487, 22183375.185, 22183378.306,
671  104340.000, 116646289.121, 90893198.077, 22197062.028, 22197064.997,
672  104370.000, 116718547.899, 90949503.596, 22210812.219, 22210815.542,
673  104400.000, 116791140.542, 91006069.256, 22224625.817, 22224629.535,
674  104430.000, 116864066.500, 91062894.640, 22238503.046, 22238506.523,
675  104460.000, 116937325.086, 91119979.217, 22252443.796, 22252447.254,
676  104490.000, 117010915.216, 91177322.147, 22266447.284, 22266451.134,
677  104520.000, 117084836.404, 91234923.044, 22280514.263, 22280517.835,
678  104550.000, 117159087.702, 91292781.186, 22294643.904, 22294647.425,
679  104580.000, 117233668.268, 91350895.888, 22308836.429, 22308839.759,
680  104610.000, 117308577.262, 91409266.514, 22323090.559, 22323094.280,
681  104640.000, 117383813.945, 91467892.475, 22337407.770, 22337411.575,
682  104670.000, 117459377.518, 91526773.148, 22351787.185, 22351791.053,
683  104700.000, 117535267.116, 91585907.861, 22366228.878, 22366232.398,
684  104730.000, 117611482.006, 91645296.075, 22380732.025, 22380735.972,
685  104760.000, 117688021.230, 91704937.000, 22395297.095, 22395300.830,
686  104790.000, 117764883.892, 91764829.962, 22409923.596, 22409927.214,
687  104820.000, 117842069.205, 91824974.330, 22424611.664, 22424615.207,
688  104850.000, 117919576.096, 91885369.290, 22439360.623, 22439363.996,
689  104880.000, 117997403.895, 91946014.304, 22454170.995, 22454174.302,
690  104910.000, 118075551.573, 92006908.570, 22469041.932, 22469045.101,
691  104940.000, 118154018.188, 92068051.341, 22483973.178, 22483976.482,
692  104970.000, 118232802.690, 92129441.832, 22498965.414, 22498968.868,
693  105000.000, 118311904.196, 92191079.332, 22514018.299, 22514021.454,
694  105030.000, 118391321.710, 92252963.090, 22529130.550, 22529133.785,
695  105060.000, 118471054.437, 92315092.441, 22544303.273, 22544306.380,
696  105090.000, 118551101.233, 92377466.530, 22559535.894, 22559538.826,
697  105120.000, 118631461.104, 92440084.581, 22574828.002, 22574831.235,
698  105150.000, 118712133.284, 92502945.990, 22590179.491, 22590182.564,
699  105180.000, 118793116.470, 92566049.749, 22605589.917, 22605593.139,
700  105210.000, 118874409.752, 92629395.149, 22621059.709, 22621063.120,
701  105240.000, 118956012.046, 92692981.324, 22636587.987, 22636591.490,
702  105270.000, 119037922.364, 92756807.497, 22652174.888, 22652178.791,
703  105300.000, 119120139.791, 92820872.992, 22667820.286, 22667824.137,
704  105330.000, 119202663.242, 92885176.950, 22683524.026, 22683527.968,
705  105360.000, 119285491.432, 92949718.350, 22699285.935, 22699289.588,
706  105390.000, 119368623.385, 93014496.460, 22715105.128, 22715109.193,
707  105420.000, 119452057.851, 93079510.283, 22730982.011, 22730986.228,
708  105450.000, 119535793.880, 93144759.089, 22746916.382, 22746920.655,
709  105480.000, 119619830.233, 93210241.942, 22762908.533, 22762912.078,
710  105510.000, 119704165.797, 93275957.927, 22778956.702, 22778960.730,
711  105540.000, 119788799.503, 93341906.216, 22795062.394, 22795066.010,
712  105570.000, 119873730.237, 93408085.957, 22811224.160, 22811228.091,
713  105600.000, 119958956.609, 93474496.081, 22827442.312, 22827446.041,
714  105630.000, 120044477.590, 93541135.765, 22843716.498, 22843720.231,
715  105660.000, 120130292.064, 93608004.154, 22860046.378, 22860050.214,
716  105690.000, 120216398.939, 93675100.397, 22876432.237, 22876435.798,
717  105720.000, 120302796.756, 93742423.349, 22892872.846, 22892876.875,
718  105750.000, 120389484.356, 93809972.099, 22909369.159, 22909372.825,
719  105780.000, 120476460.647, 93877745.788, 22925919.629, 22925924.059,
720  105810.000, 120563724.317, 93945743.411, 22942525.983, 22942529.615,
721  105840.000, 120651274.005, 94013963.926, 22959185.512, 22959189.907,
722  105870.000, 120739108.732, 94082406.536, 22975900.455, 22975904.428,
723  105900.000, 120827226.873, 94151069.982, 22992668.411, 22992673.112,
724  105930.000, 120915627.454, 94219953.520, 23009490.589, 23009495.055,
725  105960.000, 121004309.103, 94289056.073, 23026366.407, 23026370.254,
726  105990.000, 121093270.555, 94358376.645, 23043295.421, 23043299.091,
727  106020.000, 121182510.300, 94427914.079, 23060277.338, 23060281.034,
728  106050.000, 121272027.133, 94497667.431, 23077311.680, 23077315.484,
729  106080.000, 121361819.829, 94567635.716, 23094398.459, 23094402.478,
730  106110.000, 121451886.960, 94637817.860, 23111537.689, 23111541.905,
731  106140.000, 121542227.259, 94708212.856, 23128729.487, 23128733.002,
732  106170.000, 121632839.332, 94778819.647, 23145972.052, 23145975.879,
733  106200.000, 121723721.880, 94849637.179, 23163266.221, 23163270.247,
734  106230.000, 121814873.529, 94920664.385, 23180611.413, 23180616.129,
735  106260.000, 121906292.932, 94991900.248, 23198008.391, 23198012.546,
736  106290.000, 121997978.681, 95063343.652, 23215455.758, 23215459.829,
737  106320.000, 122089929.310, 95134993.462, 23232953.263, 23232957.531,
738  106350.000, 122182143.299, 95206848.462, 23250501.280, 23250505.506,
739  106380.000, 122274619.605, 95278907.893, 23268099.429, 23268102.842,
740  106410.000, 122367356.599, 95351170.466, 23285746.244, 23285750.139,
741  106440.000, 122460352.870, 95423635.059, 23303442.603, 23303447.066,
742  106470.000, 122553607.025, 95496300.602, 23321188.435, 23321192.764,
743  106500.000, 122647117.549, 95569165.915, 23338983.207, 23338987.285,
744  106530.000, 122740882.974, 95642229.847, 23356826.164, 23356830.194,
745  106560.000, 122834901.784, 95715491.220, 23374716.887, 23374721.279,
746  106590.000, 122929172.576, 95788948.955, 23392656.152, 23392660.459,
747  106620.000, 123023693.984, 95862601.983, 23410643.358, 23410647.006,
748  106650.000, 123118464.274, 95936448.951, 23428677.668, 23428681.580,
749  106680.000, 123213481.945, 96010488.686, 23446758.750, 23446763.010,
750  106710.000, 123308745.531, 96084720.028, 23464886.324, 23464890.807,
751  106740.000, 123404253.569, 96159141.844, 23483061.364, 23483065.326,
752  106770.000, 123500004.455, 96233752.911, 23501282.329, 23501286.471,
753  106800.000, 123595996.808, 96308552.127, 23519549.200, 23519553.281,
754  106830.000, 123692228.995, 96383538.220, 23537861.589, 23537865.817,
755  106860.000, 123788699.389, 96458709.945, 23556218.545, 23556223.738,
756  106890.000, 123885406.533, 96534066.122, 23574621.453, 23574626.334,
757  106920.000, 123982348.826, 96609605.536, 23593069.261, 23593073.500,
758  106950.000, 124079524.821, 96685327.089, 23611561.091, 23611565.473,
759  106980.000, 124176932.689, 96761229.287, 23630097.713, 23630102.034,
760  107010.000, 124274570.939, 96837311.001, 23648677.853, 23648681.542,
761  107040.000, 124372438.055, 96913571.070, 23667301.103, 23667304.938,
762  107070.000, 124470532.653, 96990008.405, 23685967.866, 23685972.133,
763  107100.000, 124568852.655, 97066621.349, 23704676.985, 23704681.264,
764  107130.000, 124667396.621, 97143408.830, 23723429.739, 23723433.618,
765  107160.000, 124766162.852, 97220369.504, 23742224.731, 23742228.386,
766  107190.000, 124865149.824, 97297502.178, 23761061.641, 23761065.073,
767  107220.000, 124964355.816, 97374805.519, 23779939.186, 23779943.183,
768  107250.000, 125063779.360, 97452278.367, 23798859.193, 23798862.912,
769  107280.000, 125163418.792, 97529919.455, 23817819.711, 23817823.943,
770  107310.000, 125263272.466, 97607727.507, 23836820.750, 23836825.416,
771  107340.000, 125363338.626, 97685701.086, 23855863.303, 23855868.027,
772  107370.000, 125463615.473, 97763838.852, 23874944.589, 23874950.004,
773  107400.000, 125564101.347, 97842139.497, 23894067.164, 23894071.664,
774  107430.000, 125664794.765, 97920601.845, 23913228.466, 23913232.434,
775  107460.000, 125765693.724, 97999224.356, 23932429.411, 23932432.795,
776  107490.000, 125866796.838, 98078005.980, 23951667.779, 23951671.945,
777  107520.000, 125968102.099, 98156945.101, 23970946.045, 23970949.778,
778  107550.000, 126069608.027, 98236040.592, 23990261.658, 23990266.709,
779  107580.000, 126171312.865, 98315291.095, 24009616.133, 24009620.745,
780  107610.000, 126273214.778, 98394695.165, 24029007.641, 24029011.505,
781  107670.000, 126477602.735, 98553958.498, 24067901.128, 24067904.830,
782  107700.000, 126580085.328, 98633815.056, 24087402.625, 24087407.188,
783  107730.000, 126682758.084, 98713819.776, 24106940.049, 24106944.927,
784  107760.000, 126785619.162, 98793971.231, 24126515.207, 24126518.730,
785  107790.000, 126888666.862, 98874268.154, 24146124.132, 24146127.789,
786  107820.000, 126991899.422, 98954709.105, 24165768.415, 24165773.258,
787  107880.000, 127198911.969, 99116017.602, 24205160.934, 24205166.891,
788  107910.000, 127302688.288, 99196882.263, 24224908.769, 24224913.950,
789  107940.000, 127406642.293, 99277885.403, 24244691.128, 24244696.492,
790  107970.000, 127510772.168, 99359025.596, 24264506.707, 24264511.340,
791  108000.000, 127615075.765, 99440301.138, 24284354.303, 24284358.972,
792  108060.000, 127824197.777, 99603253.364, 24324149.748, 24324154.054,
793  108090.000, 127929012.405, 99684927.098, 24344093.917, 24344099.457,
794  108180.000, 128244449.308, 99930721.989, 24404121.328, 24404125.347,
795  108210.000, 128349919.628, 100012906.642, 24424190.437, 24424195.906,
796  108240.000, 128455549.461, 100095215.550, 24444291.273, 24444296.146,
797  108270.000, 128561336.767, 100177647.197, 24464422.282, 24464426.916,
798  108300.000, 128667279.658, 100260200.104, 24484583.446, 24484587.304,
799  108330.000, 128773376.298, 100342872.808, 24504772.689, 24504776.554,
800  108360.000, 128879624.939, 100425663.935, 24524990.946, 24524994.743,
801  108420.000, 129092570.222, 100591595.281, 24565512.452, 24565517.667,
802  108450.000, 129199263.176, 100674732.620, 24585816.386, 24585819.277,
803  108480.000, 129306100.499, 100757982.489, 24606145.979, 24606149.902,
804  108540.000, 129520200.790, 100924813.902, 24646889.215, 24646892.419,
805  108570.000, 129627459.842, 101008392.409, 24667298.694, 24667303.302,
806  108600.000, 129734855.706, 101092077.499, 24687735.887, 24687740.707
807  };
808 
809  size_t i;
810  const double wl1=0.190293672798; // GPS L1 wl
811  const double wl2=0.244210213425; // GPS L2 wl
812  Stats<double> S,S1,S2,S3;
813  WtdStats<double> WS,WS1,WS2,WS3;
814  SeqStats<double> SS,SS1,SS2,SS3;
815  TwoSampleStats<double> TSS,TSS1,TSS2,TSS3;
816  vector<unsigned int> vi;
817  vector<double> vd;
818  Vector<double> Vd(M);
819 
820  for(i=0; i<M; i++) { // all the data
821  S.Add(wl1*data[5*i+1] - wl2*data[5*i+2]); // Stats L1-L2
822  WS.Add(wl1*data[5*i+1] - wl2*data[5*i+2],
823  (wl1*data[5*i+1] - data[5*i+3])/9.94); // WtdStats L1-L2 wt L1-P1/9.94
824  //SS.Add(wl1*data[5*i+1] - data[5*i+3]); // L1-P1
825  SS.Add(wl1*data[5*i+1] - wl2*data[5*i+2]); // SeqStats L1-L2
826  TSS.Add(data[5*i],wl2*data[5*i+2] - data[5*i+4]); // TwoSampleStats sow,L2-P2
827  vd.push_back(wl2*data[5*i+2] - data[5*i+4]); // L2-P2
828  Vd[i] = wl2*data[5*i+2] - data[5*i+4]; // L2-P2
829 
830  if(i < M/2) { // first half of data
831  S1.Add(wl1*data[5*i+1] - wl2*data[5*i+2]);
832  WS1.Add(wl1*data[5*i+1] - wl2*data[5*i+2],
833  (wl1*data[5*i+1] - data[5*i+3])/9.94);
834  //SS1.Add(wl1*data[5*i+1] - data[5*i+3]);
835  SS1.Add(wl1*data[5*i+1] - wl2*data[5*i+2]);
836  TSS1.Add(data[5*i],wl2*data[5*i+2] - data[5*i+4]);
837  }
838  else { // second half of data
839  S2.Add(wl1*data[5*i+1] - wl2*data[5*i+2]);
840  WS2.Add(wl1*data[5*i+1] - wl2*data[5*i+2],
841  (wl1*data[5*i+1] - data[5*i+3])/9.94);
842  //SS2.Add(wl1*data[5*i+1] - data[5*i+3]);
843  SS2.Add(wl1*data[5*i+1] - wl2*data[5*i+2]);
844  TSS2.Add(data[5*i],wl2*data[5*i+2] - data[5*i+4]);
845  }
846  }
847  LOG(INFO) << fixed << setprecision(8)
848  << "\n Stats " << S.asString("",0,8) << endl
849  << " WtdStats " << WS.asString("",0,8) << endl
850  << " SeqStats " << SS.asString("",0,8) << endl
851  << " TwoSampleStats\n" << TSS.asString("",0,8) << endl
852  << "\n median and mad (vector) " << median(vd) << " " << mad(vd) << endl
853  << " median and mad (Vector) " << median(Vd) << " " << mad(Vd);
854 
855  // -------------------
856  LOG(INFO) << "\nDump and Load ----------------------";
857  S.Dump(vi,vd); S.Load(vi,vd);
858  WS.Dump(vi,vd); WS.Load(vi,vd);
859  SS.Dump(vi,vd); SS.Load(vi,vd);
860  TSS.Dump(vi,vd); TSS.Load(vi,vd);
861  LOG(INFO) << fixed << setprecision(8)
862  << "\n Stats " << S.asString("",0,8) << endl
863  << " WtdStats " << WS.asString("",0,8) << endl
864  << " SeqStats " << SS.asString("",0,8) << endl
865  << " TwoSampleStats\n" << TSS.asString("",0,8);
866 
867  // ----------------------
868  LOG(INFO) << "\noperator =,+= --------------------------";
869  S = S1; S += S2;
870  WS = WS1; WS += WS2;
871  SS = SS1; SS += SS2;
872  TSS = TSS1; TSS += TSS2;
873  LOG(INFO) << fixed << setprecision(8)
874  << "\n Stats " << S.asString("",0,8) << endl
875  << " WtdStats " << WS.asString("",0,8) << endl
876  << " SeqStats " << SS.asString("",0,8) << endl
877  << " TwoSampleStats\n" << TSS.asString("",0,8);
878 
879  // -------------------
880  LOG(INFO) << "\noperator-=() ( NB min and max lost ) --------";
881  // subtract second half of data, first via operator-=() ...
882  S3 = S; WS3 = WS; SS3 = SS; TSS3 = TSS; // save total in *3
883  S3 -= S2; WS3 -= WS2; SS3 -= SS2; TSS3 -= TSS2; // subtract *2
884  LOG(INFO) << fixed << setprecision(8)
885  << "\n Stats " << S3.asShortString("",0,8) << endl
886  << " Compare: " << S1.asShortString("",0,8) << endl
887  << "\n WtdStats " << WS3.asShortString("",0,8) << endl
888  << " Compare: " << WS1.asShortString("",0,8) << endl
889  << "\n SeqStats " << SS3.asShortString("",0,8) << endl
890  << " Compare: " << SS1.asShortString("",0,8) << endl
891  << "\n TwoSampleStats\n" << TSS3.asShortString("",0,8) << endl
892  << TSS1.asShortString("",0,8);
893 
894  // -------------------
895  LOG(INFO) << "\nSubtract() ( NB min and max lost ) --------";
896  // ... then by Subtract()
897  for(i=M/2; i<M; i++) {
898  S.Subtract(wl1*data[5*i+1] - wl2*data[5*i+2]); // L1-L2
899  WS.Subtract(wl1*data[5*i+1] - wl2*data[5*i+2],
900  (wl1*data[5*i+1] - data[5*i+3])/9.94); // L1-L2 wtd by L1-P1
901  //SS.Subtract(wl1*data[5*i+1] - data[5*i+3]); // L1-P1
902  SS.Subtract(wl1*data[5*i+1] - wl2*data[5*i+2]); // L1-L2
903  TSS.Subtract(data[5*i],wl2*data[5*i+2] - data[5*i+4]); // sow,L2-P2
904  }
905  LOG(INFO) << fixed << setprecision(8)
906  << "\n Stats " << S.asShortString("",0,8) << endl
907  << " Compare: " << S1.asShortString("",0,8) << endl
908  << "\n WtdStats " << WS.asShortString("",0,8) << endl
909  << " Compare: " << WS1.asShortString("",0,8) << endl
910  << "\n SeqStats " << SS.asShortString("",0,8) << endl
911  << " Compare: " << SS1.asShortString("",0,8) << endl
912  << "\n TwoSampleStats:\n" << TSS.asShortString("",0,8) << endl
913  << TSS1.asShortString("",0,8);
914 
915  return 0;
916 }
917 catch(Exception& e) { cout << "Exception: " << e << endl; }
918 }
919 
920 /* the result
921  Stats stats(con): N 734 Ave 3.20392238 Std 0.48508513 Var 0.23530758 Min 2.74167747 Max 4.37292143 P2P 1.63124395
922  WtdStats stats(wtd): N 734 Ave 3.27797344 Std 0.51425401 Var 0.26445718 Min 2.74167747 Max 4.37292143 P2P 1.63124395 Wts 731.38968840
923  SeqStats stats(seq): N 734 Ave 3.20392238 Std 0.48508513 Var 0.23530758 Min 2.74167747 Max 4.37292143 P2P 1.63124395
924  TwoSampleStats
925 stats(con): N 734 Ave 97398.92370572 Std 6367.63431582 Var 40546766.77999742 Min 86400.00000000 Max 108600.00000000 P2P 22200.00000000 (X)
926 stats(con): N 734 Ave -16.31245032 Std 2.51395081 Var 6.31994869 Min -23.07857952 Max -13.63482460 P2P 9.44375493 (Y)
927 stats(tss): N 734 Int -3.31282812 Slp -0.00013347 +- 0.00001373 CSig 2.36755348 Corr -0.33806321
928 
929  median and mad (vector) -15.17627371 1.13523974
930  median and mad (Vector) -15.17627371 1.13523974
931 
932 Dump and Load ----------------------
933 
934  Stats stats(con): N 734 Ave 3.20392238 Std 0.48508513 Var 0.23530758 Min 2.74167747 Max 4.37292143 P2P 1.63124395
935  WtdStats stats(wtd): N 734 Ave 3.27797344 Std 0.51425401 Var 0.26445718 Min 2.74167747 Max 4.37292143 P2P 1.63124395 Wts 731.38968840
936  SeqStats stats(seq): N 734 Ave 3.20392238 Std 0.48508513 Var 0.23530758 Min 2.74167747 Max 4.37292143 P2P 1.63124395
937  TwoSampleStats
938 stats(con): N 734 Ave 97398.92370572 Std 6367.63431582 Var 40546766.77999742 Min 86400.00000000 Max 108600.00000000 P2P 22200.00000000 (X)
939 stats(con): N 734 Ave -16.31245032 Std 2.51395081 Var 6.31994869 Min -23.07857952 Max -13.63482460 P2P 9.44375493 (Y)
940 stats(tss): N 734 Int -3.31282812 Slp -0.00013347 +- 0.00001373 CSig 2.36755348 Corr -0.33806321
941 
942 operator =,+= --------------------------
943 
944  Stats stats(con): N 734 Ave 3.20392238 Std 0.48508513 Var 0.23530758 Min 2.74167747 Max 4.37292143 P2P 1.63124395
945  WtdStats stats(wtd): N 734 Ave 3.27797344 Std 0.51425401 Var 0.26445718 Min 2.74167747 Max 4.37292143 P2P 1.63124395 Wts 731.38968840
946  SeqStats stats(seq): N 734 Ave 3.20392238 Std 0.48508513 Var 0.23530758 Min 2.74167747 Max 4.37292143 P2P 1.63124395
947  TwoSampleStats
948 stats(con): N 734 Ave 97398.92370572 Std 6367.63431582 Var 40546766.77998585 Min 86400.00000000 Max 108600.00000000 P2P 22200.00000000 (X)
949 stats(con): N 734 Ave -16.31245032 Std 2.51395081 Var 6.31994869 Min -23.07857952 Max -13.63482460 P2P 9.44375493 (Y)
950 stats(tss): N 734 Int -3.31282812 Slp -0.00013347 +- 0.00001373 CSig 2.36755348 Corr -0.33806321
951 
952 operator-=() ( NB min and max lost ) --------
953 
954  Stats N 367 Ave 3.06723777 Std 0.37991196
955  Compare: N 367 Ave 3.06723777 Std 0.37991196
956 
957  WtdStats N 367 Ave 3.11435451 Std 0.40327050
958  Compare: N 367 Ave 3.11435451 Std 0.40327050
959 
960  SeqStats N 367 Ave 3.06723777 Std 0.37991196
961  Compare: N 367 Ave 3.06723777 Std 0.37991196
962 
963  TwoSampleStats
964  N 367 Ave 91890.00000000 Std 3182.64041324 (X); N 367 Ave -15.59581385 Std 1.96793508 (Y); Int -66.28330763 Slp 0.00055161 +- 0.00001462 CSig 0.89043841 Corr 0.89209144
965  N 367 Ave 91890.00000000 Std 3182.64041324 (X); N 367 Ave -15.59581385 Std 1.96793508 (Y); Int -66.28330763 Slp 0.00055161 +- 0.00001462 CSig 0.89043841 Corr 0.89209144
966 
967 Subtract() ( NB min and max lost ) --------
968 
969  Stats N 367 Ave 3.06723777 Std 0.37991196
970  Compare: N 367 Ave 3.06723777 Std 0.37991196
971 
972  WtdStats N 367 Ave 3.11435451 Std 0.40327050
973  Compare: N 367 Ave 3.11435451 Std 0.40327050
974 
975  SeqStats N 367 Ave 3.06723777 Std 0.37991196
976  Compare: N 367 Ave 3.06723777 Std 0.37991196
977 
978  TwoSampleStats:
979  N 367 Ave 91890.00000000 Std 3182.64041324 (X); N 367 Ave -15.59581385 Std 1.96793508 (Y); Int -66.28330763 Slp 0.00055161 +- 0.00001462 CSig 0.89043841 Corr 0.89209144
980  N 367 Ave 91890.00000000 Std 3182.64041324 (X); N 367 Ave -15.59581385 Std 1.96793508 (Y); Int -66.28330763 Slp 0.00055161 +- 0.00001462 CSig 0.89043841 Corr 0.89209144
981  */
gnsstk::Stats::asShortString
std::string asShortString(std::string msg=std::string(), int w=0, int p=3) const
Write Stats N,ave,sig to a short single-line string.
Definition: Stats.hpp:306
gnsstk::SeqStats::Add
void Add(const T &x)
Definition: Stats.hpp:404
gnsstk::WtdStats::Add
void Add(const T &x, const T &wt_in)
Definition: Stats.hpp:638
gnsstk::TwoSampleStats::Add
void Add(const T &x, const T &y)
Definition: Stats.hpp:867
gnsstk::Stats::Subtract
void Subtract(const T &x)
Definition: Stats.hpp:178
gnsstk::WtdStats::Dump
void Dump(std::vector< unsigned int > &vuint, std::vector< T > &vT)
Definition: Stats.hpp:767
gnsstk::SeqStats::asShortString
std::string asShortString(std::string msg=std::string(), int w=0, int p=3) const
Write SeqStats N,ave,sig to a short single-line string.
Definition: Stats.hpp:556
gnsstk::TwoSampleStats::Load
void Load(std::vector< unsigned int > &vuint, std::vector< T > &vT)
Definition: Stats.hpp:977
logstream.hpp
gnsstk::SeqStats::Subtract
void Subtract(const T x)
Definition: Stats.hpp:425
gnsstk::TwoSampleStats::Subtract
void Subtract(const T &x, const T &y)
Definition: Stats.hpp:879
gnsstk
For Sinex::InputHistory.
Definition: BasicFramework.cpp:50
gnsstk::Stats< double >
gnsstk::WtdStats::asShortString
std::string asShortString(std::string msg=std::string(), int w=0, int p=3) const
Write WtdStats N,ave,sig to a short single-line string.
Definition: Stats.hpp:812
Stats.hpp
gnsstk::WtdStats
Definition: Stats.hpp:620
gnsstk::Exception
Definition: Exception.hpp:151
gnsstk::TwoSampleStats::asShortString
std::string asShortString(std::string msg=std::string(), int w=0, int p=3) const
Write TwoSampleStats as a short 1-line string.
Definition: Stats.hpp:1018
gnsstk::WtdStats::Load
void Load(std::vector< unsigned int > &vuint, std::vector< T > &vT)
Definition: Stats.hpp:779
gnsstk::SeqStats::Load
void Load(const std::vector< unsigned int > &vuint, const std::vector< T > &vT)
Definition: Stats.hpp:523
gnsstk::WtdStats::Subtract
void Subtract(const T &x, const T &wt_in)
Definition: Stats.hpp:667
wl1
double wl1
Definition: DiscCorr.cpp:656
gnsstk::SeqStats::asString
std::string asString(std::string msg=std::string(), int w=7, int p=4) const
Write SeqStats to a single-line string.
Definition: Stats.hpp:540
gnsstk::Stats::asString
std::string asString(std::string msg=std::string(), int w=7, int p=4) const
Write Stats to a single-line string.
Definition: Stats.hpp:290
gnsstk::Stats::Load
void Load(const std::vector< unsigned int > &vuint, const std::vector< T > &vT)
Definition: Stats.hpp:269
gnsstk::WtdStats::asString
std::string asString(std::string msg=std::string(), int w=7, int p=4) const
Write WtdStats to a single-line string.
Definition: Stats.hpp:794
wl2
double wl2
Definition: DiscCorr.cpp:656
gnsstk::Vector< double >
example3.data
data
Definition: example3.py:22
LOG
#define LOG(level)
define the macro that is used to write to the log stream
Definition: logstream.hpp:315
gnsstk::TwoSampleStats::asString
std::string asString(std::string msg=std::string(), int w=7, int p=4) const
Write TwoSampleStats to a three-line string.
Definition: Stats.hpp:1001
std
Definition: Angle.hpp:142
gnsstk::SeqStats::Dump
void Dump(std::vector< unsigned int > &vuint, std::vector< T > &vT)
Definition: Stats.hpp:507
gnsstk::INFO
@ INFO
Definition: logstream.hpp:57
gnsstk::Stats::Add
void Add(const T &x)
Definition: Stats.hpp:158
gnsstk::TwoSampleStats::Dump
void Dump(std::vector< unsigned int > &vuint, std::vector< T > &vT)
Definition: Stats.hpp:951
gnsstk::SeqStats
Definition: Stats.hpp:388
main
int main(int argc, char **argv)
Definition: Stats_T.cpp:51
gnsstk::median
T median(const Vector< T > &v)
Compute the median of a gnsstk::Vector.
Definition: Stats.hpp:59
gnsstk::TwoSampleStats< double >
gnsstk::mad
T mad(const gnsstk::Vector< T > &v)
median absolute deviation of a gnsstk::Vector
Definition: Stats.hpp:85
Vector.hpp
gnsstk::Stats::Dump
void Dump(std::vector< unsigned int > &vuint, std::vector< T > &vT)
Definition: Stats.hpp:252

gnsstk
Author(s):
autogenerated on Wed Oct 25 2023 02:40:41