Initializing SAMURAI 3D
iMin iMax iIncr jMin jMax jIncr kMin kMax kIncr
10 35 0.5 10 40 0.5 0 10 0.5
DEBUG: center date = 2018-01-07 00:00:00.000000000, time = 0
Found matching reference time 00:42:01 at 41.6016655, 1.40277696
Physical (mish) State size = 3971968
Nodal State size = 537579
Grid dimensions: (51, 61, 21)
Found zero C level at 1836.31949 based on reference sounding
Found 34 data files to read...
Processing 20180107.cen of type cen
Processing cfrad.20180107_004225.662_to_20180107_004240.662_CREU_DEL_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 30699, Radius problem: 0
obVector size: 949
Processed 949 observations from 31772 entries (0.497817785%) file: cfrad.20180107_004225.662_to_20180107_004240.662_CREU_DEL_SUR.nc
949 total observations. ( 2 of 34 files processed )
Processing cfrad.20180107_004230.125_to_20180107_004245.125_PUIG_BER_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 24831, Radius problem: 0
obVector size: 1668
Processed 719 observations from 25568 entries (0.468684814%) file: cfrad.20180107_004230.125_to_20180107_004245.125_PUIG_BER_SUR.nc
1668 total observations. ( 3 of 34 files processed )
Processing cfrad.20180107_004244.092_to_20180107_004259.092_CREU_DEL_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 32227, Radius problem: 0
obVector size: 2667
Processed 999 observations from 33332 entries (0.499519981%) file: cfrad.20180107_004244.092_to_20180107_004259.092_CREU_DEL_SUR.nc
2667 total observations. ( 4 of 34 files processed )
Processing cfrad.20180107_004247.823_to_20180107_004301.823_PUIG_BER_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 27700, Radius problem: 0
obVector size: 3397
Processed 730 observations from 28446 entries (0.427710985%) file: cfrad.20180107_004247.823_to_20180107_004301.823_PUIG_BER_SUR.nc
3397 total observations. ( 5 of 34 files processed )
Processing cfrad.20180107_004301.347_to_20180107_004316.347_CREU_DEL_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 32671, Radius problem: 0
obVector size: 4395
Processed 998 observations from 33786 entries (0.492314371%) file: cfrad.20180107_004301.347_to_20180107_004316.347_CREU_DEL_SUR.nc
4395 total observations. ( 6 of 34 files processed )
Processing cfrad.20180107_004304.358_to_20180107_004319.358_PUIG_BER_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 30682, Radius problem: 0
obVector size: 5163
Processed 768 observations from 31456 entries (0.406917599%) file: cfrad.20180107_004304.358_to_20180107_004319.358_PUIG_BER_SUR.nc
5163 total observations. ( 7 of 34 files processed )
Processing cfrad.20180107_004318.235_to_20180107_004332.235_CREU_DEL_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 31713, Radius problem: 0
obVector size: 6195
Processed 1032 observations from 32856 entries (0.523496469%) file: cfrad.20180107_004318.235_to_20180107_004332.235_CREU_DEL_SUR.nc
6195 total observations. ( 8 of 34 files processed )
Processing cfrad.20180107_004322.342_to_20180107_004337.342_PUIG_BER_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 33360, Radius problem: 0
obVector size: 6954
Processed 759 observations from 34130 entries (0.370641664%) file: cfrad.20180107_004322.342_to_20180107_004337.342_PUIG_BER_SUR.nc
6954 total observations. ( 9 of 34 files processed )
Processing cfrad.20180107_004335.672_to_20180107_004350.672_CREU_DEL_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 29558, Radius problem: 0
obVector size: 7992
Processed 1038 observations from 30725 entries (0.563059398%) file: cfrad.20180107_004335.672_to_20180107_004350.672_CREU_DEL_SUR.nc
7992 total observations. ( 10 of 34 files processed )
Processing cfrad.20180107_004340.760_to_20180107_004355.760_PUIG_BER_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 33634, Radius problem: 0
obVector size: 8749
Processed 757 observations from 34400 entries (0.366763566%) file: cfrad.20180107_004340.760_to_20180107_004355.760_PUIG_BER_SUR.nc
8749 total observations. ( 11 of 34 files processed )
Processing cfrad.20180107_004351.684_to_20180107_004406.684_CREU_DEL_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 27510, Radius problem: 0
obVector size: 9762
Processed 1013 observations from 28657 entries (0.589152156%) file: cfrad.20180107_004351.684_to_20180107_004406.684_CREU_DEL_SUR.nc
9762 total observations. ( 12 of 34 files processed )
Processing cfrad.20180107_004359.360_to_20180107_004414.360_PUIG_BER_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 31107, Radius problem: 0
obVector size: 10463
Processed 701 observations from 31837 entries (0.366973438%) file: cfrad.20180107_004359.360_to_20180107_004414.360_PUIG_BER_SUR.nc
10463 total observations. ( 13 of 34 files processed )
Processing cfrad.20180107_004409.010_to_20180107_004423.010_CREU_DEL_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 21448, Radius problem: 0
obVector size: 11439
Processed 976 observations from 22545 entries (0.721519923%) file: cfrad.20180107_004409.010_to_20180107_004423.010_CREU_DEL_SUR.nc
11439 total observations. ( 14 of 34 files processed )
Processing cfrad.20180107_004416.403_to_20180107_004431.403_PUIG_BER_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 27802, Radius problem: 0
obVector size: 12075
Processed 636 observations from 28455 entries (0.372518011%) file: cfrad.20180107_004416.403_to_20180107_004431.403_PUIG_BER_SUR.nc
12075 total observations. ( 15 of 34 files processed )
Processing cfrad.20180107_004426.910_to_20180107_004441.910_CREU_DEL_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 18078, Radius problem: 0
obVector size: 12968
Processed 893 observations from 19115 entries (0.77862063%) file: cfrad.20180107_004426.910_to_20180107_004441.910_CREU_DEL_SUR.nc
12968 total observations. ( 16 of 34 files processed )
Processing cfrad.20180107_004435.419_to_20180107_004450.419_PUIG_BER_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 25612, Radius problem: 0
obVector size: 13595
Processed 627 observations from 26255 entries (0.398019425%) file: cfrad.20180107_004435.419_to_20180107_004450.419_PUIG_BER_SUR.nc
13595 total observations. ( 17 of 34 files processed )
Processing cfrad.20180107_004444.218_to_20180107_004459.218_CREU_DEL_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 14849, Radius problem: 0
obVector size: 14473
Processed 878 observations from 15823 entries (0.924814089%) file: cfrad.20180107_004444.218_to_20180107_004459.218_CREU_DEL_SUR.nc
14473 total observations. ( 18 of 34 files processed )
Processing cfrad.20180107_004453.420_to_20180107_004508.420_PUIG_BER_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 21862, Radius problem: 0
obVector size: 15054
Processed 581 observations from 22475 entries (0.430849092%) file: cfrad.20180107_004453.420_to_20180107_004508.420_PUIG_BER_SUR.nc
15054 total observations. ( 19 of 34 files processed )
Processing cfrad.20180107_004501.640_to_20180107_004516.640_CREU_DEL_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 11660, Radius problem: 0
obVector size: 15902
Processed 848 observations from 12604 entries (1.12133714%) file: cfrad.20180107_004501.640_to_20180107_004516.640_CREU_DEL_SUR.nc
15902 total observations. ( 20 of 34 files processed )
Processing cfrad.20180107_004510.725_to_20180107_004525.725_PUIG_BER_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 19200, Radius problem: 0
obVector size: 16431
Processed 529 observations from 19775 entries (0.445849136%) file: cfrad.20180107_004510.725_to_20180107_004525.725_PUIG_BER_SUR.nc
16431 total observations. ( 21 of 34 files processed )
Processing cfrad.20180107_004519.306_to_20180107_004534.306_CREU_DEL_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 7960, Radius problem: 0
obVector size: 17173
Processed 742 observations from 8742 entries (1.41462671%) file: cfrad.20180107_004519.306_to_20180107_004534.306_CREU_DEL_SUR.nc
17173 total observations. ( 22 of 34 files processed )
Processing cfrad.20180107_004529.202_to_20180107_004544.202_PUIG_BER_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 16330, Radius problem: 0
obVector size: 17482
Processed 309 observations from 16683 entries (0.308697476%) file: cfrad.20180107_004529.202_to_20180107_004544.202_PUIG_BER_SUR.nc
17482 total observations. ( 23 of 34 files processed )
Processing cfrad.20180107_004535.929_to_20180107_004550.929_CREU_DEL_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 6141, Radius problem: 0
obVector size: 18107
Processed 625 observations from 6861 entries (1.51824321%) file: cfrad.20180107_004535.929_to_20180107_004550.929_CREU_DEL_SUR.nc
18107 total observations. ( 24 of 34 files processed )
Processing cfrad.20180107_004546.178_to_20180107_004601.178_PUIG_BER_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 14452, Radius problem: 0
obVector size: 18136
Processed 29 observations from 14487 entries (0.0333632452%) file: cfrad.20180107_004546.178_to_20180107_004601.178_PUIG_BER_SUR.nc
18136 total observations. ( 25 of 34 files processed )
Processing cfrad.20180107_004554.649_to_20180107_004608.649_CREU_DEL_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 4552, Radius problem: 0
obVector size: 18249
Processed 113 observations from 4776 entries (0.394332775%) file: cfrad.20180107_004554.649_to_20180107_004608.649_CREU_DEL_SUR.nc
18249 total observations. ( 26 of 34 files processed )
Processing cfrad.20180107_004605.722_to_20180107_004620.722_PUIG_BER_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 12388, Radius problem: 0
obVector size: 18249
No valid observations in file
18249 total observations. ( 27 of 34 files processed )
Processing cfrad.20180107_004612.565_to_20180107_004626.565_CREU_DEL_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 3852, Radius problem: 0
obVector size: 18259
Processed 10 observations from 3889 entries (0.0428559184%) file: cfrad.20180107_004612.565_to_20180107_004626.565_CREU_DEL_SUR.nc
18259 total observations. ( 28 of 34 files processed )
Processing cfrad.20180107_004625.394_to_20180107_004640.394_PUIG_BER_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 10608, Radius problem: 0
obVector size: 18259
No valid observations in file
18259 total observations. ( 29 of 34 files processed )
Processing cfrad.20180107_004630.909_to_20180107_004645.909_CREU_DEL_SUR.nc of type cfrad
i = 29
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 3077, Radius problem: 0
obVector size: 18259
No valid observations in file
18259 total observations. ( 30 of 34 files processed )
Processing cfrad.20180107_004645.119_to_20180107_004700.119_PUIG_BER_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 9293, Radius problem: 0
obVector size: 18259
No valid observations in file
18259 total observations. ( 31 of 34 files processed )
Processing cfrad.20180107_004649.801_to_20180107_004704.801_CREU_DEL_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 2914, Radius problem: 0
obVector size: 18259
No valid observations in file
18259 total observations. ( 32 of 34 files processed )
Processing cfrad.20180107_004704.310_to_20180107_004719.310_PUIG_BER_SUR.nc of type cfrad
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 8402, Radius problem: 0
obVector size: 18259
No valid observations in file
18259 total observations. ( 33 of 34 files processed )
Processing terrain.hgt of type hgt
i = 33
Observation problem: 0, Time problem: 0, Coordinate problem: 0, Domain problem: 104, Radius problem: 0
obVector size: 27280
Processed 9021 observations from 3111 entries (48.3285117%) file: terrain.hgt
27280 total observations. ( 34 of 34 files processed )
92611 total observations including pseudo-obs for W and mass continuity
Finished preprocessing 33 files.
Number of New Observations: 92611
runTime: loadMetObs: 14.363112
done with stof conversion...
kRankMax: 23
runTime: initObCost3D: 0.01922
runTime: gridDependentInit: 14.394375
Outer Loop Iteration: 1
Initializing state vector...
4th Order Recursive Filter coefficients computed for lengthscale 4 Delta X
4th Order Recursive Filter coefficients computed for lengthscale 4 Delta X
4th Order Recursive Filter coefficients computed for lengthscale 4 Delta X
i Spline cutoff set to 2
j Spline cutoff set to 2
k Spline cutoff set to 2
Mass continuity weight set to 1
Initializing background...
Variable 0 RMS = 0 BG Error = 5.05923221e-321 ( Infinite! %)
Variable 1 RMS = 0 BG Error = 0 ( Infinite! %)
Variable 2 RMS = 0 BG Error = 0 ( Infinite! %)
Variable 3 RMS = 0 BG Error = 4.66019637e-310 ( Infinite! %)
Variable 4 RMS = 0 BG Error = 1.38338381e-322 ( Infinite! %)
Variable 5 RMS = 0 BG Error = 6.90503772e-310 ( Infinite! %)
Variable 6 RMS = 4.06861428 BG Error = 0 ( 0 %)
Build H transform matrix...
calcHmatrix: Grid dimensions: (53, 63, 23)
sizeof(integer): 8
Non-zero entries in sparse H matrix: 7560194 = 0.0151854641 %
Memory usage for [H] (Mbytes): 57.6797028
Memory usage for [obsVector] (Mbytes): 24.729805
Memory usage for [obsData] (Mbytes): 0.706565857
Memory usage for [HCq] (Mbytes): 1.29248047
Memory usage for [mPtr,mVal,I2H] (Mbytes): 119.460815
Memory usage for [IH,JH] (Mbytes): 58.3862762
Memory usage for [state] (Mbytes): 4.10140228
Initializing innovation vector...
Innovation RMS : 4.14597356
Outputting background...
runTime: Cost3D Init: 4.130128
Beginning minimize()
Solver conv. tolerance = 0.0001
SOLVER: Samurai Truncated Newton
Newton Iteration: 0 J = 162295.75 Residual = 8824.86594
CG iteration 0: r_norm = 8824.86594150 rel_resid = 1.0000000000
CG iteration 1: r_norm = 5615.89841651 rel_resid = 0.6363720938
CG iteration 2: r_norm = 7253.96362984 rel_resid = 0.8219913682
CG iteration 3: r_norm = 3719.27506046 rel_resid = 0.4214540011
CG iteration 4: r_norm = 3431.15302209 rel_resid = 0.3888051155
CG iteration 5: r_norm = 3431.65020968 rel_resid = 0.3888614549
CG iteration 6: r_norm = 2570.45376292 rel_resid = 0.2912739729
CG iteration 7: r_norm = 2538.87481416 rel_resid = 0.2876955674
CG iteration 8: r_norm = 1992.34761062 rel_resid = 0.2257651985
CG iteration 9: r_norm = 1689.11150344 rel_resid = 0.1914036445
CG iteration 10: r_norm = 1761.64803318 rel_resid = 0.1996232062
CG iteration 11: r_norm = 1538.37004739 rel_resid = 0.1743222002
CG iteration 12: r_norm = 1605.00467853 rel_resid = 0.1818729813
CG iteration 13: r_norm = 1443.10645291 rel_resid = 0.1635272946
CG iteration 14: r_norm = 1228.59333689 rel_resid = 0.1392194901
CG iteration 15: r_norm = 1390.72749523 rel_resid = 0.1575919118
CG iteration 16: r_norm = 1218.73565071 rel_resid = 0.1381024549
CG iteration 17: r_norm = 981.90111401 rel_resid = 0.1112652725
CG iteration 18: r_norm = 1014.00705665 rel_resid = 0.1149033949
CG iteration 19: r_norm = 1058.33085927 rel_resid = 0.1199259985
CG iteration 20: r_norm = 881.77321618 rel_resid = 0.0999191627
CG iteration 21: r_norm = 939.83333144 rel_resid = 0.1064983126
CG iteration 22: r_norm = 718.47392112 rel_resid = 0.0814147122
CG iteration 23: r_norm = 960.82647771 rel_resid = 0.1088771755
CG iteration 24: r_norm = 703.36095568 rel_resid = 0.0797021689
CG iteration 25: r_norm = 678.01820828 rel_resid = 0.0768304258
CG iteration 26: r_norm = 626.67179983 rel_resid = 0.0710120475
CG iteration 27: r_norm = 552.19382315 rel_resid = 0.0625724886
CG iteration 28: r_norm = 707.75468638 rel_resid = 0.0802000496
CG iteration 29: r_norm = 541.74635610 rel_resid = 0.0613886216
CG iteration 30: r_norm = 574.86443956 rel_resid = 0.0651414360
CG iteration 31: r_norm = 473.85288570 rel_resid = 0.0536951937
CG iteration 32: r_norm = 545.11389650 rel_resid = 0.0617702184
CG iteration 33: r_norm = 452.88419724 rel_resid = 0.0513191022
CG iteration 34: r_norm = 668.67522952 rel_resid = 0.0757717153
CG iteration 35: r_norm = 519.55850739 rel_resid = 0.0588743796
CG iteration 36: r_norm = 429.82448774 rel_resid = 0.0487060643
CG iteration 37: r_norm = 586.88989789 rel_resid = 0.0665041148
CG iteration 38: r_norm = 440.25855214 rel_resid = 0.0498884125
CG iteration 39: r_norm = 331.76330402 rel_resid = 0.0375941466
CG iteration 40: r_norm = 324.40613048 rel_resid = 0.0367604599
CG iteration 41: r_norm = 300.03256556 rel_resid = 0.0339985409
CG iteration 42: r_norm = 305.22813044 rel_resid = 0.0345872824
CG iteration 43: r_norm = 406.78665347 rel_resid = 0.0460955052
CG iteration 44: r_norm = 234.18100511 rel_resid = 0.0265364943
CG iteration 45: r_norm = 247.22923589 rel_resid = 0.0280150699
CG iteration 46: r_norm = 389.17699582 rel_resid = 0.0441000462
CG iteration 47: r_norm = 261.47904260 rel_resid = 0.0296298034
CG iteration 48: r_norm = 236.49154438 rel_resid = 0.0267983158
CG iteration 49: r_norm = 376.21524727 rel_resid = 0.0426312705
CG iteration 50: r_norm = 225.42893415 rel_resid = 0.0255447432
CG iteration 51: r_norm = 194.86760693 rel_resid = 0.0220816507
CG iteration 52: r_norm = 215.19874239 rel_resid = 0.0243854971
CG iteration 53: r_norm = 287.76081745 rel_resid = 0.0326079534
CG iteration 54: r_norm = 197.79760506 rel_resid = 0.0224136668
CG iteration 55: r_norm = 185.83492832 rel_resid = 0.0210581021
CG iteration 56: r_norm = 243.86758187 rel_resid = 0.0276341401
CG iteration 57: r_norm = 173.24874471 rel_resid = 0.0196318840
CG iteration 58: r_norm = 253.36085600 rel_resid = 0.0287098816
CG iteration 59: r_norm = 169.77387937 rel_resid = 0.0192381256
CG iteration 60: r_norm = 158.23201608 rel_resid = 0.0179302459
CG iteration 61: r_norm = 208.56118994 rel_resid = 0.0236333550
CG iteration 62: r_norm = 137.76569088 rel_resid = 0.0156110803
CG iteration 63: r_norm = 236.83214040 rel_resid = 0.0268369108
CG iteration 64: r_norm = 127.30258176 rel_resid = 0.0144254409
CG iteration 65: r_norm = 164.14864235 rel_resid = 0.0186006953
CG iteration 66: r_norm = 131.90117197 rel_resid = 0.0149465355
CG iteration 67: r_norm = 110.91050109 rel_resid = 0.0125679531
CG iteration 68: r_norm = 144.39577105 rel_resid = 0.0163623756
CG iteration 69: r_norm = 96.62031752 rel_resid = 0.0109486442
CG iteration 70: r_norm = 98.49916777 rel_resid = 0.0111615483
CG iteration 71: r_norm = 138.49123157 rel_resid = 0.0156932958
CG iteration 72: r_norm = 120.86417679 rel_resid = 0.0136958655
CG iteration 73: r_norm = 124.52467712 rel_resid = 0.0141106594
CG iteration 74: r_norm = 81.25911684 rel_resid = 0.0092079718
CG iteration 75: r_norm = 90.75868943 rel_resid = 0.0102844270
CG iteration 76: r_norm = 104.07402960 rel_resid = 0.0117932703
CG iteration 77: r_norm = 70.12596942 rel_resid = 0.0079464062
CG iteration 78: r_norm = 118.99382582 rel_resid = 0.0134839245
CG iteration 79: r_norm = 69.45493503 rel_resid = 0.0078703672
CG iteration 80: r_norm = 105.94122819 rel_resid = 0.0120048541
CG iteration 81: r_norm = 88.93113561 rel_resid = 0.0100773356
CG iteration 82: r_norm = 76.71855697 rel_resid = 0.0086934530
CG iteration 83: r_norm = 87.97351094 rel_resid = 0.0099688212
CG iteration 84: r_norm = 58.34458877 rel_resid = 0.0066113853
CG iteration 85: r_norm = 52.67166495 rel_resid = 0.0059685513
CG iteration 86: r_norm = 62.50258947 rel_resid = 0.0070825540
CG iteration 87: r_norm = 62.02462742 rel_resid = 0.0070283932
CG iteration 88: r_norm = 88.73856375 rel_resid = 0.0100555141
CG iteration 89: r_norm = 50.35501554 rel_resid = 0.0057060374
CG iteration 90: r_norm = 66.36236497 rel_resid = 0.0075199290
CG iteration 91: r_norm = 61.48398756 rel_resid = 0.0069671299
CG iteration 92: r_norm = 58.47102767 rel_resid = 0.0066257129
CG iteration 93: r_norm = 42.62680415 rel_resid = 0.0048303061
CG iteration 94: r_norm = 56.07672008 rel_resid = 0.0063543991
CG iteration 95: r_norm = 36.06532961 rel_resid = 0.0040867850
CG iteration 96: r_norm = 35.37625167 rel_resid = 0.0040087013
CG iteration 97: r_norm = 51.57801829 rel_resid = 0.0058446234
CG iteration 98: r_norm = 41.83336626 rel_resid = 0.0047403968
CG iteration 99: r_norm = 41.20600784 rel_resid = 0.0046693069
CG iteration 100: r_norm = 45.76146183 rel_resid = 0.0051855135
CG iteration 101: r_norm = 35.76087598 rel_resid = 0.0040522855
CG iteration 102: r_norm = 32.15313366 rel_resid = 0.0036434699
CG iteration 103: r_norm = 45.98787615 rel_resid = 0.0052111699
CG iteration 104: r_norm = 46.99801377 rel_resid = 0.0053256349
CG iteration 105: r_norm = 28.82702720 rel_resid = 0.0032665683
CG iteration 106: r_norm = 25.43263163 rel_resid = 0.0028819284
CG iteration 107: r_norm = 28.16531310 rel_resid = 0.0031915854
CG iteration 108: r_norm = 29.95134475 rel_resid = 0.0033939716
CG iteration 109: r_norm = 30.77700134 rel_resid = 0.0034875319
CG iteration 110: r_norm = 21.37121955 rel_resid = 0.0024217047
CG iteration 111: r_norm = 24.46411719 rel_resid = 0.0027721800
CG iteration 112: r_norm = 29.66593552 rel_resid = 0.0033616302
CG iteration 113: r_norm = 25.20063967 rel_resid = 0.0028556399
CG iteration 114: r_norm = 23.26443675 rel_resid = 0.0026362368
CG iteration 115: r_norm = 17.47873943 rel_resid = 0.0019806238
CG iteration 116: r_norm = 24.67364332 rel_resid = 0.0027959227
CG iteration 117: r_norm = 18.27643370 rel_resid = 0.0020710154
CG iteration 118: r_norm = 18.24765147 rel_resid = 0.0020677540
CG iteration 119: r_norm = 19.95634112 rel_resid = 0.0022613761
CG iteration 120: r_norm = 17.93586214 rel_resid = 0.0020324232
CG iteration 121: r_norm = 14.08204592 rel_resid = 0.0015957235
CG iteration 122: r_norm = 14.44470256 rel_resid = 0.0016368184
CG iteration 123: r_norm = 17.02149068 rel_resid = 0.0019288101
CG iteration 124: r_norm = 19.77628569 rel_resid = 0.0022409729
CG iteration 125: r_norm = 13.14771323 rel_resid = 0.0014898485
CG iteration 126: r_norm = 19.71696249 rel_resid = 0.0022342507
CG iteration 127: r_norm = 11.08626557 rel_resid = 0.0012562531
CG iteration 128: r_norm = 15.00003568 rel_resid = 0.0016997466
CG iteration 129: r_norm = 14.63790738 rel_resid = 0.0016587116
CG iteration 130: r_norm = 14.11250775 rel_resid = 0.0015991753
CG iteration 131: r_norm = 11.54584858 rel_resid = 0.0013083313
CG iteration 132: r_norm = 9.55908846 rel_resid = 0.0010831993
CG iteration 133: r_norm = 10.06208086 rel_resid = 0.0011401965
CG iteration 134: r_norm = 8.55682997 rel_resid = 0.0009696272
CG iteration 135: r_norm = 7.42669296 rel_resid = 0.0008415644
CG iteration 136: r_norm = 11.01115193 rel_resid = 0.0012477416
CG iteration 137: r_norm = 10.91470095 rel_resid = 0.0012368121
CG iteration 138: r_norm = 9.10785455 rel_resid = 0.0010320672
CG iteration 139: r_norm = 10.34154941 rel_resid = 0.0011718648
CG iteration 140: r_norm = 8.51499874 rel_resid = 0.0009648870
CG iteration 141: r_norm = 7.65966156 rel_resid = 0.0008679635
CG iteration 142: r_norm = 6.29202168 rel_resid = 0.0007129878
CG iteration 143: r_norm = 8.27504196 rel_resid = 0.0009376961
CG iteration 144: r_norm = 8.05174652 rel_resid = 0.0009123931
CG iteration 145: r_norm = 7.13579119 rel_resid = 0.0008086005
CG iteration 146: r_norm = 8.22453095 rel_resid = 0.0009319723
CG iteration 147: r_norm = 7.86619343 rel_resid = 0.0008913669
CG iteration 148: r_norm = 7.20886998 rel_resid = 0.0008168815
CG iteration 149: r_norm = 5.87314888 rel_resid = 0.0006655227
CG iteration 150: r_norm = 5.55234329 rel_resid = 0.0006291703
CG iteration 151: r_norm = 7.28056957 rel_resid = 0.0008250063
CG iteration 152: r_norm = 6.41635268 rel_resid = 0.0007270765
CG iteration 153: r_norm = 4.59359688 rel_resid = 0.0005205288
CG iteration 154: r_norm = 4.04930476 rel_resid = 0.0004588517
CG iteration 155: r_norm = 5.03025264 rel_resid = 0.0005700090
CG iteration 156: r_norm = 4.35901882 rel_resid = 0.0004939473
CG iteration 157: r_norm = 5.03632665 rel_resid = 0.0005706972
CG iteration 158: r_norm = 5.64617856 rel_resid = 0.0006398033
CG iteration 159: r_norm = 3.22142933 rel_resid = 0.0003650400
CG iteration 160: r_norm = 3.05328987 rel_resid = 0.0003459871
CG iteration 161: r_norm = 4.12198049 rel_resid = 0.0004670870
CG iteration 162: r_norm = 3.05290605 rel_resid = 0.0003459436
CG iteration 163: r_norm = 5.03626813 rel_resid = 0.0005706906
CG iteration 164: r_norm = 2.75024701 rel_resid = 0.0003116475
CG iteration 165: r_norm = 5.76258314 rel_resid = 0.0006529938
CG iteration 166: r_norm = 2.68324107 rel_resid = 0.0003040546
CG iteration 167: r_norm = 2.41854946 rel_resid = 0.0002740608
CG iteration 168: r_norm = 4.23095019 rel_resid = 0.0004794351
CG iteration 169: r_norm = 3.55951775 rel_resid = 0.0004033509
CG iteration 170: r_norm = 2.81041974 rel_resid = 0.0003184660
CG iteration 171: r_norm = 2.64979449 rel_resid = 0.0003002646
CG iteration 172: r_norm = 2.67075879 rel_resid = 0.0003026402
CG iteration 173: r_norm = 2.00937033 rel_resid = 0.0002276941
CG iteration 174: r_norm = 2.09722403 rel_resid = 0.0002376494
CG iteration 175: r_norm = 2.69225255 rel_resid = 0.0003050757
CG iteration 176: r_norm = 2.34224956 rel_resid = 0.0002654147
CG iteration 177: r_norm = 2.79358790 rel_resid = 0.0003165587
CG iteration 178: r_norm = 1.84253016 rel_resid = 0.0002087885
CG iteration 179: r_norm = 2.18936334 rel_resid = 0.0002480903
CG iteration 180: r_norm = 2.03963223 rel_resid = 0.0002311233
CG iteration 181: r_norm = 2.27059905 rel_resid = 0.0002572956
CG iteration 182: r_norm = 1.70054072 rel_resid = 0.0001926988
CG iteration 183: r_norm = 2.17775824 rel_resid = 0.0002467752
CG iteration 184: r_norm = 1.57575370 rel_resid = 0.0001785584
CG iteration 185: r_norm = 1.74949890 rel_resid = 0.0001982465
CG iteration 186: r_norm = 1.49043286 rel_resid = 0.0001688901
CG iteration 187: r_norm = 1.25994887 rel_resid = 0.0001427726
CG iteration 188: r_norm = 1.16790644 rel_resid = 0.0001323427
CG iteration 189: r_norm = 1.30090705 rel_resid = 0.0001474138
CG iteration 190: r_norm = 1.62162029 rel_resid = 0.0001837558
CG iteration 191: r_norm = 1.11138509 rel_resid = 0.0001259379
CG iteration 192: r_norm = 1.38256998 rel_resid = 0.0001566675
CG iteration 193: r_norm = 1.63187337 rel_resid = 0.0001849176
CG iteration 194: r_norm = 0.86796514 rel_resid = 0.0000983545
MT LS: Number of function evals = 1
MT LS: step = 1.0000000000
Newton Iteration: 1 J = 20529.5642703829 Residual = 0.8679651393
Minimum J: 20529.5642703829
Found minimum in 1 outer Newton iterations.
(and a total of 194 inner CG iterations.)
FINAL ||g(X)||/||g(X0)|| = 0.0000983545
(Linesearch iterations = 1 )
runTime: Cost3D minimize: 109.7625850000
runTime: Cost3d update: 6.4667280000
Analysis successful!
PropList::~PropList - H5Pclose failed
PropList::~PropList - H5Pclose failed
PropList::~PropList - H5Pclose failed
PropList::~PropList - H5Pclose failed
PropList::~PropList - H5Pclose failed
PropList::~PropList - H5Pclose failed
PropList::~PropList - H5Pclose failed
PropList::~PropList - H5Pclose failed