Generate new starting conditions of numbers and f-at-age.
monteCarloStock.RdMethod to sample from the variance-co-variance matrix and based on a multi-variante normal distribution generate new starting conditions of numbers-at-age and f-at-age.
monteCarloStock(stck,sam,realisations,run.dir=tempdir())Arguments
- stck
An object of class
FLStock- sam
An object of class
FLSAM- realisations
Number of newly generated starting conditions.
- run.dir
Path to which the temporary files are to be written to or read from. Default is the R temporary directory provided by tempdir()
Details
Large number of realisations might take longer to create.
Value
An FLStock with iter=number of realisations is returned
See also
Examples
data(NSH)
data(NSH.sam)
#Stock with one iter
summary(NSH)#iter = 1
#> An object of class "FLStock"
#>
#> Name: North Sea Herring
#> Description: Imported from a VPA file. ( ./data/index.txt ). Fri May 4 12:16:13 2012 [...]
#> Quant: age
#> Dims: age year unit season area iter
#> 9 20 1 1 1 1
#>
#> Range: min max pgroup minyear maxyear minfbar maxfbar
#> 0 8 8 1992 2011 2 6
#>
#> Metrics:
#> rec: 2.1e+07 - 7.8e+07 (NA)
#> ssb: 748114 - 2281590 (NA)
#> catch: 168443 - 716799 (tonnes)
#> fbar: 0.082 - 0.471 (f)
plot(NSH)
#Generate new stock with 100 iters
newStock <- monteCarloStock(NSH,NSH.tun, NSH.sam,100, ncores=2L)
#> Loading required package: doParallel
#> Loading required package: foreach
#> Loading required package: parallel
#> Order of parameters:
#> [1] "logFpar" "logQpow" "logSdLogFsta"
#> [4] "logSdLogN" "logSdLogP" "logSdLogObs"
#> [7] "logSdLogTotalObs" "transfIRARdist" "sigmaObsParUS"
#> [10] "rec_loga" "rec_logb" "itrans_rho"
#> [13] "rhop" "logScale" "logitReleaseSurvival"
#> [16] "logitRecapturePhi" "logAlphaSCB" "logF"
#> [19] "logN" "logPS" "missing"
#> Not matching template order:
#> [1] "logFpar" "logQpow" "logSdLogFsta"
#> [4] "logSdLogN" "logSdLogP" "logSdLogObs"
#> [7] "logSdLogTotalObs" "transfIRARdist" "sigmaObsParUS"
#> [10] "rec_loga" "rec_logb" "itrans_rho"
#> [13] "rhop" "logScale" "logitReleaseSurvival"
#> [16] "logitRecapturePhi" "logAlphaSCB" "logF"
#> [19] "logN" "missing"
#> Your parameter list has been re-ordered.
#> (Disable this warning with checkParameterOrder=FALSE)
#> iter: 1 value: 29858.94 mgc: 102.6853 ustep: 0.1843228
#> iter: 2 value: 6427.944 mgc: 41.19916 ustep: 0.4293854
#> iter: 3 value: 3643.672 mgc: 11.96224 ustep: 0.6553096
#> iter: 4 value: 2383.509 mgc: 9.832191 ustep: 0.809531
#> iter: 5 value: 2252.316 mgc: 1.584616 ustep: 0.8997494
#> iter: 6 value: 2250.408 mgc: 0.2357221 ustep: 0.9485564
#> iter: 7 value: 2250.401 mgc: 0.01217538 ustep: 0.9739412
#> iter: 8 value: 2250.401 mgc: 0.0007627327 ustep: 0.9868859
#> iter: 9 value: 2250.401 mgc: 2.411668e-05 ustep: 0.993422
#> iter: 10 value: 2250.401 mgc: 3.787673e-07 ustep: 0.9967059
#> iter: 11 mgc: 2.962014e-09
#> iter: 1 mgc: 2.962014e-09
#> Matching hessian patterns... Done
#> outer mgc: 937.0649
#> 0: 2288.8818: -5.00000 -5.00000 -5.00000 -5.00000 -5.00000 -5.00000 -5.00000 -5.00000 -5.00000 -5.00000 -0.700000 -0.350000 -0.350000 -0.350000 -0.350000 -0.350000 -0.350000 -0.350000 0.500000
#> iter: 1 value: 1295.734 mgc: 3.191557 ustep: 1
#> iter: 2 value: 1295.689 mgc: 0.1784069 ustep: 1
#> iter: 3 value: 1295.689 mgc: 0.003145601 ustep: 1
#> iter: 4 value: 1295.689 mgc: 1.676632e-06 ustep: 1
#> iter: 5 mgc: 3.715361e-13
#> iter: 1 mgc: 3.715361e-13
#> outer mgc: 344.246
#> 1: 1270.1871: -5.08863 -4.91055 -4.94563 -4.96799 -4.97773 -4.98310 -4.98710 -4.98254 -5.08762 -5.11268 -0.718198 -0.351086 0.0769365 -0.342072 -0.0408757 0.230886 -0.0324865 0.141923 0.508719
#> iter: 1 value: 916.0222 mgc: 6.34312 ustep: 1
#> iter: 2 value: 910.0898 mgc: 1.137912 ustep: 1
#> iter: 3 value: 909.8531 mgc: 0.2259272 ustep: 1
#> iter: 4 value: 909.8526 mgc: 0.01513482 ustep: 1
#> iter: 5 value: 909.8526 mgc: 6.401704e-05 ustep: 1
#> iter: 6 mgc: 1.148572e-09
#> iter: 1 mgc: 1.148572e-09
#> outer mgc: 163.7715
#> 2: 942.62357: -5.26901 -4.85738 -4.88709 -4.91877 -4.93595 -4.94694 -4.95573 -4.92846 -5.23754 -5.24126 -0.845336 -0.380952 -0.224249 -0.649964 0.305213 0.687157 0.150118 0.730410 0.559968
#> iter: 1 value: 711.6633 mgc: 2.375197 ustep: 1
#> iter: 2 value: 710.0988 mgc: 0.8324954 ustep: 1
#> iter: 3 value: 710.0852 mgc: 0.09314666 ustep: 1
#> iter: 4 value: 710.0852 mgc: 0.001066194 ustep: 1
#> iter: 5 value: 710.0852 mgc: 1.415094e-07 ustep: 1
#> iter: 6 mgc: 3.896883e-14
#> iter: 1 mgc: 3.896883e-14
#> outer mgc: 57.50863
#> 3: 786.26128: -5.33686 -4.80028 -4.81374 -4.85811 -4.88974 -4.91120 -4.92796 -4.88979 -5.45167 -5.33225 -0.965713 -0.411656 -0.177317 -0.901893 0.299597 1.02425 0.733917 1.35011 0.620378
#> iter: 1 value: 492.9034 mgc: 2.889296 ustep: 1
#> iter: 2 value: 491.1143 mgc: 0.9347992 ustep: 1
#> iter: 3 value: 491.1137 mgc: 0.03301362 ustep: 1
#> iter: 4 value: 491.1137 mgc: 1.489686e-05 ustep: 1
#> iter: 5 mgc: 2.430611e-12
#> iter: 1 value: 43.9683 mgc: 29.49475 ustep: 1
#> iter: 2 value: 43.83872 mgc: 0.4673514 ustep: 1
#> iter: 3 value: 43.83871 mgc: 0.0003426394 ustep: 1
#> iter: 4 mgc: 6.38518e-10
#> iter: 1 mgc: 6.38518e-10
#> outer mgc: 126.6576
#> 4: 576.80251: -6.06882 -4.51133 -4.50742 -4.59489 -4.66793 -4.72204 -4.76780 -4.63555 -6.11355 -5.70707 -1.67300 -0.604378 -2.03271 -2.03541 1.15094 0.338267 2.63337 3.40849 1.01895
#> iter: 1 value: 10.90599 mgc: 287.6597 ustep: 1
#> iter: 2 value: 10.57287 mgc: 2.012651 ustep: 1
#> iter: 3 value: 10.57287 mgc: 0.001464435 ustep: 1
#> iter: 4 mgc: 3.178302e-11
#> iter: 1 mgc: 3.178302e-11
#> outer mgc: 93.30105
#> 5: 498.37476: -6.65570 -4.31143 -4.25414 -4.34687 -4.46008 -4.55356 -4.64757 -4.63072 -6.56274 -5.99569 0.137165 -0.711751 -1.27442 -1.74676 1.71751 -2.56440 3.31441 4.42017 1.98770
#> iter: 1 value: 484.7178 mgc: 113.5134 ustep: 1
#> iter: 2 value: 484.5256 mgc: 0.5652743 ustep: 1
#> iter: 3 value: 484.5256 mgc: 0.0001485126 ustep: 1
#> iter: 4 mgc: 1.70769e-12
#> iter: 1 value: 106.9895 mgc: 55.12887 ustep: 1
#> iter: 2 value: 106.9844 mgc: 0.2052389 ustep: 1
#> iter: 3 value: 106.9844 mgc: 7.489573e-06 ustep: 1
#> iter: 4 mgc: 4.911627e-13
#> iter: 1 mgc: 4.911627e-13
#> outer mgc: 76.85225
#> 6: 491.00656: -6.66999 -4.26083 -4.34105 -4.09262 -4.42427 -4.45011 -4.37306 -5.23851 -6.56441 -5.99600 -0.128228 -0.670972 -0.390569 -1.18392 1.56334 -2.22251 3.12617 4.22461 2.07530
#> iter: 1 value: -0.5745773 mgc: 29.00865 ustep: 1
#> iter: 2 value: -0.5746843 mgc: 0.02304346 ustep: 1
#> iter: 3 value: -0.5746843 mgc: 2.547398e-07 ustep: 1
#> mgc: 3.312906e-13
#> iter: 1 value: 49.36313 mgc: 32.37045 ustep: 1
#> iter: 2 value: 49.36204 mgc: 0.0462509 ustep: 1
#> iter: 3 value: 49.36204 mgc: 1.097814e-06 ustep: 1
#> iter: 4 mgc: 4.707346e-13
#> iter: 1 mgc: 3.312906e-13
#> outer mgc: 27.28687
#> 7: 434.75037: -6.68411 -4.26905 -4.24485 -4.20237 -4.34531 -4.41311 -4.61056 -5.07841 -6.56616 -5.99633 -0.315710 -0.675389 -0.912721 -1.27809 1.46269 -2.22922 2.99468 4.08550 2.15432
#> iter: 1 value: -31.74598 mgc: 29.47695 ustep: 1
#> iter: 2 value: -31.74606 mgc: 0.01982511 ustep: 1
#> iter: 3 value: -31.74606 mgc: 1.125765e-07 ustep: 1
#> mgc: 2.922107e-13
#> iter: 1 mgc: 2.922107e-13
#> outer mgc: 74.91487
#> 8: 428.91412: -6.71488 -4.22055 -4.13926 -4.23411 -4.27124 -4.53460 -4.89226 -4.83272 -6.57027 -5.99712 -0.505705 -0.670007 -1.24980 -1.27787 1.29098 -2.13682 2.76149 3.83463 2.17510
#> iter: 1 value: -44.24781 mgc: 20.36353 ustep: 1
#> iter: 2 value: -44.24787 mgc: 0.03235738 ustep: 1
#> iter: 3 value: -44.24787 mgc: 9.909799e-07 ustep: 1
#> iter: 4 mgc: 2.950973e-13
#> iter: 1 mgc: 2.950973e-13
#> outer mgc: 45.7465
#> 9: 408.29444: -6.76537 -4.08991 -4.07067 -4.12277 -4.54335 -4.64501 -4.59675 -4.93873 -6.61104 -6.02762 -0.472540 -0.671571 -1.21432 -1.47061 1.08332 -1.90196 2.52168 3.59189 2.22633
#> iter: 1 value: -6.717147 mgc: 8.573656 ustep: 1
#> iter: 2 value: -6.717177 mgc: 0.0284937 ustep: 1
#> iter: 3 value: -6.717177 mgc: 1.096477e-06 ustep: 1
#> iter: 4 mgc: 1.815215e-13
#> iter: 1 mgc: 1.815215e-13
#> outer mgc: 32.87295
#> 10: 394.31589: -6.81731 -4.02510 -4.10146 -4.30707 -4.27650 -4.40279 -4.80134 -5.00047 -6.63440 -6.04003 -0.430123 -0.659987 -1.01709 -1.47696 0.951682 -1.55221 2.34211 3.37177 2.18546
#> iter: 1 value: -50.1307 mgc: 11.00374 ustep: 1
#> iter: 2 value: -50.13088 mgc: 0.06387695 ustep: 1
#> iter: 3 value: -50.13088 mgc: 4.803896e-06 ustep: 1
#> iter: 4 mgc: 2.082778e-13
#> iter: 1 mgc: 2.082778e-13
#> outer mgc: 24.99611
#> 11: 375.69112: -6.95112 -4.09117 -4.16527 -3.97784 -4.23833 -4.50069 -4.71779 -4.99852 -6.69564 -6.06941 -0.375082 -0.640590 -1.33727 -1.66785 0.839080 -1.35530 2.13643 3.07604 2.18809
#> iter: 1 value: -89.61545 mgc: 7.139076 ustep: 1
#> iter: 2 value: -89.61546 mgc: 0.005349777 ustep: 1
#> iter: 3 value: -89.61546 mgc: 3.020959e-08 ustep: 1
#> mgc: 2.930989e-13
#> iter: 1 value: -75.98503 mgc: 3.155036 ustep: 1
#> iter: 2 value: -75.98503 mgc: 0.001101492 ustep: 1
#> iter: 3 mgc: 1.265883e-09
#> iter: 1 mgc: 1.265883e-09
#> outer mgc: 14.81626
#> 12: 368.86281: -7.00322 -4.09849 -3.98769 -4.13110 -4.26981 -4.44866 -4.70371 -4.97953 -6.70813 -6.07275 -0.483879 -0.632218 -1.35571 -1.70698 0.807091 -1.35541 2.06263 2.96205 2.19733
#> iter: 1 value: -89.45378 mgc: 2.242948 ustep: 1
#> iter: 2 value: -89.45378 mgc: 0.0004036388 ustep: 1
#> iter: 3 mgc: 2.206004e-10
#> iter: 1 mgc: 2.206004e-10
#> outer mgc: 17.48676
#> 13: 364.45262: -7.08605 -4.02318 -4.10072 -3.99898 -4.27613 -4.48015 -4.67286 -4.95153 -6.73024 -6.07918 -0.576026 -0.616547 -1.37192 -1.72953 0.758923 -1.34833 1.96556 2.80012 2.15309
#> iter: 1 value: -93.32494 mgc: 2.971631 ustep: 1
#> iter: 2 value: -93.32494 mgc: 0.001642562 ustep: 1
#> iter: 3 mgc: 4.984663e-09
#> iter: 1 mgc: 4.984663e-09
#> outer mgc: 14.02576
#> 14: 359.47725: -7.20614 -4.03265 -3.95273 -4.10361 -4.22457 -4.42130 -4.65876 -4.92430 -6.77471 -6.09713 -0.617694 -0.598782 -1.35516 -1.76957 0.708356 -1.32967 1.90315 2.65320 2.07248
#> iter: 1 value: -91.91236 mgc: 5.9552 ustep: 1
#> iter: 2 value: -91.91236 mgc: 0.004405361 ustep: 1
#> iter: 3 mgc: 7.01302e-09
#> iter: 1 mgc: 7.01302e-09
#> outer mgc: 19.38089
#> 15: 352.76546: -7.56094 -3.97994 -3.93135 -3.86215 -4.22535 -4.30388 -4.58847 -4.85509 -6.91783 -6.16650 -0.654940 -0.562017 -1.37058 -1.83444 0.552728 -1.34904 1.83075 2.41184 1.84740
#> iter: 1 value: -104.1327 mgc: 4.331121 ustep: 1
#> iter: 2 value: -104.1327 mgc: 0.002005606 ustep: 1
#> iter: 3 mgc: 1.427948e-09
#> iter: 1 mgc: 1.427948e-09
#> outer mgc: 31.57602
#> 16: 347.80653: -7.91191 -3.82429 -3.90400 -3.89727 -3.95320 -4.39887 -4.42622 -4.76992 -7.05715 -6.24440 -0.775335 -0.542026 -1.37067 -1.86604 0.347147 -1.35639 1.75840 2.28037 1.74785
#> iter: 1 value: -119.9602 mgc: 7.89804 ustep: 1
#> iter: 2 value: -119.9602 mgc: 0.002476018 ustep: 1
#> iter: 3 mgc: 7.03952e-09
#> iter: 1 mgc: 7.03952e-09
#> outer mgc: 15.24421
#> 17: 337.06396: -8.23055 -3.69979 -3.76931 -3.89042 -4.04133 -4.15305 -4.47697 -4.61065 -7.17794 -6.33048 -0.917952 -0.546071 -1.31631 -1.81735 0.0563173 -1.26940 1.67922 2.32910 1.79631
#> iter: 1 value: -128.0141 mgc: 3.705085 ustep: 1
#> iter: 2 value: -128.0141 mgc: 0.005910159 ustep: 1
#> iter: 3 value: -128.0141 mgc: 5.766538e-08 ustep: 1
#> mgc: 2.593481e-13
#> iter: 1 mgc: 2.593481e-13
#> outer mgc: 21.27486
#> 18: 330.75662: -8.52589 -3.73179 -3.66388 -3.63825 -3.94609 -4.17030 -4.26442 -4.70957 -7.30689 -6.41707 -0.922718 -0.546740 -1.36373 -1.82273 -0.287121 -1.31255 1.60092 2.34005 1.74733
#> iter: 1 value: -113.7405 mgc: 4.681769 ustep: 1
#> iter: 2 value: -113.7405 mgc: 0.007929936 ustep: 1
#> iter: 3 value: -113.7405 mgc: 3.70615e-08 ustep: 1
#> mgc: 2.939871e-13
#> iter: 1 mgc: 2.939871e-13
#> outer mgc: 20.35934
#> 19: 328.72342: -8.75901 -3.55010 -3.63257 -3.68571 -3.82468 -4.02991 -4.36348 -4.52223 -7.47698 -6.52248 -0.926637 -0.551111 -1.40081 -1.76677 -0.674361 -1.31909 1.55323 2.32432 1.62436
#> iter: 1 value: -105.1224 mgc: 2.615154 ustep: 1
#> iter: 2 value: -105.1224 mgc: 0.001958143 ustep: 1
#> iter: 3 mgc: 2.446879e-09
#> iter: 1 mgc: 2.446879e-09
#> outer mgc: 20.73529
#> 20: 326.49498: -8.26251 -3.59899 -3.53822 -3.70491 -3.86662 -4.07014 -4.32019 -4.57849 -7.74471 -6.66499 -0.913209 -0.576408 -1.34422 -1.78385 -0.771078 -1.29652 1.54353 2.25617 1.55594
#> iter: 1 value: -146.2882 mgc: 6.933896 ustep: 1
#> iter: 2 value: -146.2882 mgc: 0.004359399 ustep: 1
#> iter: 3 value: -146.2882 mgc: 1.198158e-08 ustep: 1
#> mgc: 3.501643e-13
#> iter: 1 mgc: 3.501643e-13
#> outer mgc: 28.24945
#> 21: 325.76336: -8.47716 -3.55525 -3.59217 -3.45450 -3.80316 -3.99310 -4.15402 -4.45759 -8.03177 -6.79974 -1.06966 -0.593885 -1.33108 -1.83845 -0.543479 -1.44367 1.51393 2.19427 1.65453
#> iter: 1 value: -121.088 mgc: 3.901788 ustep: 1
#> iter: 2 value: -121.088 mgc: 0.002751137 ustep: 1
#> iter: 3 mgc: 9.98716e-09
#> iter: 1 mgc: 9.98716e-09
#> outer mgc: 14.28193
#> 22: 321.61859: -8.44559 -3.53345 -3.47948 -3.64057 -3.73643 -3.94716 -4.21083 -4.44575 -8.05942 -6.81560 -0.988772 -0.582858 -1.32856 -1.77818 -0.562312 -1.29295 1.49091 2.20282 1.61179
#> iter: 1 value: -120.9591 mgc: 3.337906 ustep: 1
#> iter: 2 value: -120.9591 mgc: 0.006136845 ustep: 1
#> iter: 3 value: -120.9591 mgc: 6.944051e-08 ustep: 1
#> mgc: 2.268186e-13
#> iter: 1 mgc: 2.268186e-13
#> outer mgc: 8.657306
#> 23: 320.37471: -8.41722 -3.41488 -3.51204 -3.53411 -3.78226 -3.93260 -4.12064 -4.46144 -8.21725 -6.89580 -0.929805 -0.560096 -1.43551 -1.77095 -0.592733 -1.28817 1.39448 2.20387 1.58906
#> iter: 1 value: -132.4649 mgc: 2.857943 ustep: 1
#> iter: 2 value: -132.465 mgc: 0.004922483 ustep: 1
#> iter: 3 value: -132.465 mgc: 2.533174e-08 ustep: 1
#> mgc: 2.222666e-13
#> iter: 1 mgc: 2.222666e-13
#> outer mgc: 8.729763
#> 24: 317.89702: -8.34163 -3.46818 -3.41859 -3.46740 -3.69224 -3.92286 -4.18331 -4.40322 -8.38674 -6.98170 -0.976545 -0.556986 -1.36650 -1.79024 -0.627434 -1.35619 1.30963 2.22639 1.62533
#> iter: 1 value: -128.5957 mgc: 2.105914 ustep: 1
#> iter: 2 value: -128.5957 mgc: 0.001561296 ustep: 1
#> iter: 3 mgc: 1.704339e-09
#> iter: 1 mgc: 1.704339e-09
#> outer mgc: 10.15385
#> 25: 316.18660: -8.29247 -3.39915 -3.41452 -3.48458 -3.66107 -3.84682 -4.05005 -4.40763 -8.61061 -7.07517 -0.973352 -0.571162 -1.33072 -1.76766 -0.645933 -1.34054 1.27949 2.19424 1.63696
#> iter: 1 value: -127.2518 mgc: 2.48407 ustep: 1
#> iter: 2 value: -127.2518 mgc: 0.003728131 ustep: 1
#> iter: 3 mgc: 7.041147e-09
#> iter: 1 mgc: 7.041147e-09
#> outer mgc: 5.502011
#> 26: 313.87951: -8.26513 -3.36100 -3.36476 -3.40136 -3.64077 -3.83444 -4.04281 -4.29303 -8.84358 -7.16603 -0.956284 -0.578390 -1.37253 -1.74037 -0.621510 -1.30635 1.23812 2.15774 1.64033
#> iter: 1 value: -142.1829 mgc: 1.918994 ustep: 1
#> iter: 2 value: -142.1829 mgc: 0.001291332 ustep: 1
#> iter: 3 mgc: 2.208718e-09
#> iter: 1 mgc: 2.208718e-09
#> outer mgc: 9.203617
#> 27: 312.18646: -8.21269 -3.30876 -3.30120 -3.39999 -3.55819 -3.75326 -4.02829 -4.31859 -9.05277 -7.25355 -0.958675 -0.581249 -1.40120 -1.78751 -0.619837 -1.39367 1.14042 2.16750 1.64104
#> iter: 1 value: -134.5646 mgc: 1.845753 ustep: 1
#> iter: 2 value: -134.5646 mgc: 0.0009035202 ustep: 1
#> iter: 3 mgc: 1.396582e-09
#> iter: 1 mgc: 1.396582e-09
#> outer mgc: 6.880071
#> 28: 309.26288: -8.17991 -3.27052 -3.27249 -3.29969 -3.53794 -3.74857 -3.96316 -4.26603 -9.27656 -7.34173 -0.948206 -0.577773 -1.35801 -1.77187 -0.620578 -1.33921 1.03339 2.17279 1.62875
#> iter: 1 value: -139.6183 mgc: 1.661759 ustep: 1
#> iter: 2 value: -139.6183 mgc: 0.002179619 ustep: 1
#> iter: 3 mgc: 4.229629e-09
#> iter: 1 mgc: 4.229629e-09
#> outer mgc: 11.59808
#> 29: 307.24905: -8.14062 -3.21573 -3.19751 -3.33850 -3.48886 -3.68563 -3.96909 -4.17091 -9.50627 -7.42008 -0.952029 -0.573050 -1.38139 -1.77313 -0.630285 -1.34481 0.943999 2.14461 1.63274
#> iter: 1 value: -143.4753 mgc: 1.44656 ustep: 1
#> iter: 2 value: -143.4753 mgc: 0.0007161764 ustep: 1
#> iter: 3 mgc: 7.584497e-10
#> iter: 1 value: -156.1174 mgc: 1.711774 ustep: 1
#> iter: 2 value: -156.1174 mgc: 0.004236756 ustep: 1
#> iter: 3 value: -156.1174 mgc: 4.181891e-08 ustep: 1
#> mgc: 2.202682e-13
#> iter: 1 value: -183.0237 mgc: 7.851342 ustep: 1
#> iter: 2 value: -183.0254 mgc: 0.09407709 ustep: 1
#> iter: 3 value: -183.0254 mgc: 9.809408e-06 ustep: 1
#> iter: 4 mgc: 3.663736e-13
#> iter: 1 mgc: 2.202682e-13
#> outer mgc: 8.834843
#> 30: 297.04518: -8.00270 -3.08167 -3.07672 -3.17299 -3.36456 -3.55267 -3.79083 -4.06906 -10.2615 -7.67695 -0.959654 -0.567855 -1.39435 -1.77839 -0.618232 -1.34191 0.660549 2.11339 1.67168
#> iter: 1 value: -93.99458 mgc: 7.944217 ustep: 1
#> iter: 2 value: -93.99476 mgc: 0.02458971 ustep: 1
#> iter: 3 value: -93.99476 mgc: 3.664464e-07 ustep: 1
#> iter: 4 mgc: 1.617595e-13
#> iter: 1 mgc: 1.617595e-13
#> outer mgc: 45.88069
#> 31: 292.48806: -7.79262 -2.86688 -3.06734 -2.60530 -3.11122 -3.47281 -3.54477 -3.81831 -11.5469 -8.08964 -0.791124 -0.566227 -1.11294 -1.70546 -0.675987 -1.22095 -0.282994 1.94694 1.53924
#> iter: 1 value: -66.84359 mgc: 8.067476 ustep: 1
#> iter: 2 value: -66.84394 mgc: 0.03355627 ustep: 1
#> iter: 3 value: -66.84394 mgc: 2.198293e-06 ustep: 1
#> iter: 4 mgc: 2.4869e-13
#> iter: 1 value: -134.938 mgc: 3.032645 ustep: 1
#> iter: 2 value: -134.938 mgc: 0.002094367 ustep: 1
#> iter: 3 value: -134.938 mgc: 1.468709e-08 ustep: 1
#> mgc: 1.833256e-13
#> iter: 1 mgc: 1.833256e-13
#> outer mgc: 21.33665
#> 32: 276.49998: -7.78031 -2.88558 -2.93767 -2.78351 -3.10280 -3.39107 -3.56612 -3.82270 -11.5415 -8.10172 -0.855847 -0.568832 -1.21542 -1.74531 -0.671257 -1.25207 -0.342219 1.95379 1.57979
#> iter: 1 value: -174.5127 mgc: 3.076463 ustep: 1
#> iter: 2 value: -174.5127 mgc: 0.002701163 ustep: 1
#> iter: 3 value: -174.5127 mgc: 1.183342e-08 ustep: 1
#> mgc: 2.005063e-13
#> iter: 1 mgc: 2.005063e-13
#> outer mgc: 13.7304
#> 33: 269.39586: -7.76584 -2.89292 -2.85086 -2.92392 -3.11023 -3.31931 -3.56610 -3.82656 -11.5299 -8.12187 -0.923375 -0.571123 -1.35268 -1.79142 -0.664593 -1.34607 -0.439068 1.96270 1.62234
#> iter: 1 value: -186.7597 mgc: 1.962868 ustep: 1
#> iter: 2 value: -186.7597 mgc: 0.001294899 ustep: 1
#> iter: 3 mgc: 2.603285e-09
#> iter: 1 mgc: 2.603285e-09
#> outer mgc: 10.2687
#> 34: 266.55096: -7.74789 -2.87672 -2.92308 -2.89635 -3.14220 -3.31681 -3.50760 -3.80687 -11.4968 -8.17464 -0.956717 -0.571551 -1.41078 -1.78895 -0.649106 -1.33926 -0.680699 1.97723 1.63956
#> iter: 1 value: -182.7085 mgc: 1.74895 ustep: 1
#> iter: 2 value: -182.7085 mgc: 0.0007275734 ustep: 1
#> iter: 3 mgc: 3.263985e-09
#> iter: 1 value: -182.4933 mgc: 1.143937 ustep: 1
#> iter: 2 value: -182.4933 mgc: 0.0003562506 ustep: 1
#> iter: 3 mgc: 2.585114e-10
#> iter: 1 mgc: 2.585114e-10
#> outer mgc: 8.770688
#> 35: 265.42975: -7.74772 -2.90215 -2.85540 -2.95156 -3.10323 -3.30080 -3.54058 -3.78459 -11.4971 -8.19639 -0.950847 -0.573025 -1.37928 -1.77952 -0.646706 -1.32155 -0.757919 1.97652 1.63545
#> iter: 1 value: -185.887 mgc: 1.089121 ustep: 1
#> iter: 2 value: -185.887 mgc: 0.0002185108 ustep: 1
#> iter: 3 mgc: 1.096838e-10
#> iter: 1 mgc: 1.096838e-10
#> outer mgc: 6.784882
#> 36: 264.66675: -7.73746 -2.87941 -2.89431 -2.88388 -3.12054 -3.29784 -3.49821 -3.79696 -11.5228 -8.23329 -0.951006 -0.575159 -1.38738 -1.78515 -0.644191 -1.33798 -0.848158 1.97429 1.63733
#> iter: 1 value: -184.3164 mgc: 0.8553228 ustep: 1
#> iter: 2 value: -184.3164 mgc: 0.0004203745 ustep: 1
#> iter: 3 mgc: 1.159748e-10
#> iter: 1 mgc: 1.159748e-10
#> outer mgc: 7.92419
#> 37: 264.41075: -7.74976 -2.90672 -2.88048 -2.93322 -3.10848 -3.30232 -3.53206 -3.78594 -11.4147 -8.24317 -0.947332 -0.578652 -1.38388 -1.78099 -0.645801 -1.33242 -0.899837 1.97161 1.63215
#> iter: 1 value: -185.1135 mgc: 0.9089372 ustep: 1
#> iter: 2 value: -185.1135 mgc: 0.0002229392 ustep: 1
#> iter: 3 mgc: 1.59357e-11
#> iter: 1 mgc: 1.59357e-11
#> outer mgc: 4.863725
#> 38: 263.79077: -7.74078 -2.84363 -2.85794 -2.89984 -3.11812 -3.29836 -3.50603 -3.79742 -11.5050 -8.27609 -0.946351 -0.580842 -1.38704 -1.77967 -0.643925 -1.33388 -0.960188 1.96717 1.63149
#> iter: 1 value: -184.7091 mgc: 1.295468 ustep: 1
#> iter: 2 value: -184.7091 mgc: 0.000565141 ustep: 1
#> iter: 3 mgc: 1.894276e-10
#> iter: 1 value: -185.5498 mgc: 0.5064377 ustep: 1
#> iter: 2 value: -185.5498 mgc: 7.617522e-05 ustep: 1
#> iter: 3 mgc: 3.038236e-12
#> iter: 1 mgc: 3.038236e-12
#> outer mgc: 3.613651
#> 39: 263.41970: -7.73622 -2.87452 -2.85726 -2.90983 -3.09581 -3.29174 -3.52225 -3.77941 -11.4727 -8.29533 -0.949094 -0.581810 -1.38448 -1.77688 -0.642654 -1.33502 -0.987490 1.96426 1.63333
#> iter: 1 value: -185.698 mgc: 0.5333973 ustep: 1
#> iter: 2 value: -185.698 mgc: 6.395254e-05 ustep: 1
#> iter: 3 mgc: 1.025358e-11
#> iter: 1 mgc: 1.025358e-11
#> outer mgc: 2.920778
#> 40: 263.18124: -7.72778 -2.86676 -2.86101 -2.89510 -3.10191 -3.28461 -3.49483 -3.78523 -11.4725 -8.33035 -0.948210 -0.583783 -1.38653 -1.77426 -0.641367 -1.33601 -1.03101 1.95919 1.63303
#> iter: 1 value: -185.4519 mgc: 0.4610834 ustep: 1
#> iter: 2 value: -185.4519 mgc: 0.0001580008 ustep: 1
#> iter: 3 mgc: 1.487521e-11
#> iter: 1 mgc: 1.487521e-11
#> outer mgc: 3.543853
#> 41: 263.01404: -7.72137 -2.86752 -2.84274 -2.90779 -3.07969 -3.27493 -3.51126 -3.75769 -11.4762 -8.36585 -0.948203 -0.585562 -1.38778 -1.76939 -0.640359 -1.33389 -1.06130 1.95446 1.63416
#> iter: 1 value: -185.5589 mgc: 0.5889439 ustep: 1
#> iter: 2 value: -185.5589 mgc: 0.0001078573 ustep: 1
#> iter: 3 mgc: 2.94067e-11
#> iter: 1 mgc: 2.94067e-11
#> outer mgc: 4.185684
#> 42: 262.89354: -7.71099 -2.86175 -2.84294 -2.89001 -3.08478 -3.26583 -3.47137 -3.76930 -11.4918 -8.40128 -0.947171 -0.587053 -1.38649 -1.76938 -0.639329 -1.33207 -1.08543 1.95070 1.63595
#> iter: 1 value: -186.2026 mgc: 0.4665431 ustep: 1
#> iter: 2 value: -186.2026 mgc: 0.0001728082 ustep: 1
#> iter: 3 mgc: 2.399414e-11
#> iter: 1 mgc: 2.399414e-11
#> outer mgc: 3.213317
#> 43: 262.76925: -7.71020 -2.87136 -2.85822 -2.88326 -3.07655 -3.26521 -3.49788 -3.74462 -11.4471 -8.42469 -0.948928 -0.588488 -1.39062 -1.76717 -0.639548 -1.33713 -1.09043 1.94602 1.63627
#> iter: 1 value: -185.5791 mgc: 0.5099855 ustep: 1
#> iter: 2 value: -185.5791 mgc: 8.986419e-05 ustep: 1
#> iter: 3 mgc: 1.503242e-11
#> iter: 1 mgc: 1.503242e-11
#> outer mgc: 4.187853
#> 44: 262.65670: -7.70264 -2.85569 -2.82521 -2.89967 -3.06385 -3.25601 -3.48366 -3.74217 -11.4761 -8.45918 -0.946908 -0.589705 -1.38927 -1.76498 -0.638986 -1.33243 -1.10462 1.94163 1.63662
#> iter: 1 value: -185.9749 mgc: 0.5602958 ustep: 1
#> iter: 2 value: -185.9749 mgc: 6.939509e-05 ustep: 1
#> iter: 3 mgc: 7.931877e-12
#> iter: 1 mgc: 7.931877e-12
#> outer mgc: 3.478932
#> 45: 262.54946: -7.69962 -2.85372 -2.84978 -2.86328 -3.07874 -3.25753 -3.46203 -3.75506 -11.4466 -8.48538 -0.947823 -0.591136 -1.38818 -1.76591 -0.639264 -1.33587 -1.10653 1.93635 1.63631
#> iter: 1 value: -186.3045 mgc: 0.5819925 ustep: 1
#> iter: 2 value: -186.3045 mgc: 0.0002454965 ustep: 1
#> iter: 3 mgc: 3.587797e-11
#> iter: 1 mgc: 3.587797e-11
#> outer mgc: 5.752161
#> 46: 262.43988: -7.69354 -2.85696 -2.83766 -2.87375 -3.05613 -3.24528 -3.49151 -3.71981 -11.4425 -8.51979 -0.948475 -0.591910 -1.39217 -1.76297 -0.638643 -1.33614 -1.11494 1.93168 1.63810
#> iter: 1 value: -185.6816 mgc: 0.6182255 ustep: 1
#> iter: 2 value: -185.6816 mgc: 0.0002255229 ustep: 1
#> iter: 3 mgc: 6.522916e-11
#> iter: 1 mgc: 6.522916e-11
#> outer mgc: 3.245662
#> 47: 262.28620: -7.68961 -2.83524 -2.81555 -2.87683 -3.05902 -3.24938 -3.45097 -3.74815 -11.4503 -8.54830 -0.947125 -0.593162 -1.38748 -1.76333 -0.639018 -1.33385 -1.11491 1.92590 1.63594
#> iter: 1 value: -186.3732 mgc: 0.4556932 ustep: 1
#> iter: 2 value: -186.3732 mgc: 8.493742e-05 ustep: 1
#> iter: 3 mgc: 5.296541e-12
#> iter: 1 mgc: 5.296541e-12
#> outer mgc: 3.129478
#> 48: 262.08187: -7.68151 -2.83627 -2.83057 -2.84861 -3.05767 -3.23673 -3.46245 -3.72688 -11.4293 -8.59337 -0.948625 -0.594040 -1.39042 -1.76181 -0.638335 -1.33703 -1.11909 1.91885 1.63718
#> iter: 1 value: -185.9522 mgc: 0.5217167 ustep: 1
#> iter: 2 value: -185.9522 mgc: 0.0001093261 ustep: 1
#> iter: 3 mgc: 3.597878e-11
#> iter: 1 mgc: 3.597878e-11
#> outer mgc: 4.34768
#> 49: 261.94336: -7.67691 -2.82156 -2.79563 -2.87122 -3.03617 -3.23492 -3.45904 -3.72142 -11.4305 -8.63634 -0.948210 -0.595094 -1.38979 -1.75975 -0.638707 -1.33615 -1.11900 1.91063 1.63558
#> iter: 1 value: -186.2604 mgc: 0.4669873 ustep: 1
#> iter: 2 value: -186.2604 mgc: 7.266764e-05 ustep: 1
#> iter: 3 mgc: 6.816159e-12
#> iter: 1 mgc: 6.816159e-12
#> outer mgc: 3.081353
#> 50: 261.69175: -7.66757 -2.81543 -2.80666 -2.84007 -3.04174 -3.22529 -3.44108 -3.72037 -11.4171 -8.68615 -0.948873 -0.595631 -1.38931 -1.75953 -0.638625 -1.33709 -1.12068 1.90215 1.63551
#> iter: 1 value: -186.3536 mgc: 0.3614256 ustep: 1
#> iter: 2 value: -186.3536 mgc: 9.26844e-05 ustep: 1
#> iter: 3 mgc: 1.531841e-11
#> iter: 1 mgc: 1.531841e-11
#> outer mgc: 3.663214
#> 51: 261.55294: -7.66257 -2.81019 -2.78512 -2.85047 -3.02005 -3.21647 -3.45178 -3.69778 -11.4093 -8.73567 -0.948727 -0.596239 -1.39121 -1.75803 -0.639002 -1.33772 -1.12131 1.89301 1.63536
#> iter: 1 value: -186.5698 mgc: 0.4138182 ustep: 1
#> iter: 2 value: -186.5698 mgc: 5.14084e-05 ustep: 1
#> iter: 3 mgc: 7.23821e-12
#> iter: 1 mgc: 7.23821e-12
#> outer mgc: 3.137757
#> 52: 261.28681: -7.65178 -2.79918 -2.78799 -2.82739 -3.02316 -3.20920 -3.42460 -3.70338 -11.4025 -8.78724 -0.949023 -0.596381 -1.38963 -1.75863 -0.638977 -1.33743 -1.12203 1.88411 1.63506
#> iter: 1 value: -187.0582 mgc: 0.5957186 ustep: 1
#> iter: 2 value: -187.0582 mgc: 0.0002517181 ustep: 1
#> iter: 3 mgc: 1.074403e-10
#> iter: 1 mgc: 1.074403e-10
#> outer mgc: 5.461532
#> 53: 260.92060: -7.63478 -2.78410 -2.76309 -2.81937 -2.99699 -3.18697 -3.43287 -3.66619 -11.3809 -8.90076 -0.949893 -0.596801 -1.39270 -1.75662 -0.638998 -1.33891 -1.12239 1.86413 1.63485
#> iter: 1 value: -187.9182 mgc: 1.253023 ustep: 1
#> iter: 2 value: -187.9182 mgc: 0.0008833191 ustep: 1
#> iter: 3 mgc: 1.050337e-09
#> iter: 1 mgc: 1.050337e-09
#> outer mgc: 8.35194
#> 54: 260.23746: -7.59212 -2.76255 -2.71291 -2.81232 -2.93961 -3.16003 -3.34980 -3.64121 -11.3657 -9.12922 -0.949310 -0.595275 -1.39073 -1.75878 -0.639922 -1.33596 -1.13180 1.82568 1.63702
#> iter: 1 value: -189.386 mgc: 1.361969 ustep: 1
#> iter: 2 value: -189.386 mgc: 0.0009211608 ustep: 1
#> iter: 3 mgc: 1.624163e-09
#> iter: 1 mgc: 1.624163e-09
#> outer mgc: 4.128509
#> 55: 259.04994: -7.55141 -2.70617 -2.71296 -2.72312 -2.94548 -3.10986 -3.32617 -3.61292 -11.3231 -9.35347 -0.950933 -0.595942 -1.39255 -1.75946 -0.638044 -1.33331 -1.13264 1.78576 1.63701
#> iter: 1 value: -190.9283 mgc: 1.067596 ustep: 1
#> iter: 2 value: -190.9283 mgc: 0.001143683 ustep: 1
#> iter: 3 mgc: 2.542368e-09
#> iter: 1 value: -195.4124 mgc: 3.680597 ustep: 1
#> iter: 2 value: -195.4124 mgc: 0.008395325 ustep: 1
#> iter: 3 value: -195.4124 mgc: 1.948353e-07 ustep: 1
#> iter: 4 mgc: 2.629008e-13
#> iter: 1 value: -208.1192 mgc: 15.42634 ustep: 1
#> iter: 2 value: -208.1273 mgc: 0.1776816 ustep: 1
#> iter: 3 value: -208.1273 mgc: 9.688305e-05 ustep: 1
#> iter: 4 mgc: 3.668321e-11
#> iter: 1 mgc: 3.668321e-11
#> outer mgc: 42.38944
#> 56: 248.99510: -6.90357 -1.88032 -2.07584 -2.08418 -2.41015 -2.64114 -2.62991 -3.25452 -10.6804 -13.1098 -0.988683 -0.621804 -1.39628 -1.78499 -0.616259 -1.33503 -1.22981 1.11153 1.63978
#> iter: 1 value: -186.8762 mgc: 15.67179 ustep: 1
#> iter: 2 value: -186.9324 mgc: 0.5685008 ustep: 1
#> iter: 3 value: -186.9324 mgc: 0.0009063712 ustep: 1
#> iter: 4 mgc: 5.113241e-09
#> iter: 1 value: -220.211 mgc: 3.803081 ustep: 1
#> iter: 2 value: -220.2112 mgc: 0.03003581 ustep: 1
#> iter: 3 value: -220.2112 mgc: 2.46931e-06 ustep: 1
#> iter: 4 mgc: 2.526868e-13
#> iter: 1 mgc: 2.526868e-13
#> outer mgc: 24.27286
#> 57: 234.41811: -6.77584 -1.80260 -1.93020 -1.96311 -2.24747 -2.46113 -2.55938 -3.02892 -10.5187 -13.9901 -0.975339 -0.625746 -1.38220 -1.79204 -0.617257 -1.28382 -1.23007 0.878594 1.65396
#> iter: 1 value: -226.2729 mgc: 5.436528 ustep: 1
#> iter: 2 value: -226.2732 mgc: 0.02661691 ustep: 1
#> iter: 3 value: -226.2732 mgc: 2.189867e-06 ustep: 1
#> iter: 4 mgc: 2.344513e-13
#> iter: 1 mgc: 2.344513e-13
#> outer mgc: 26.01495
#> 58: 224.46534: -6.68330 -1.84058 -1.79356 -1.88635 -2.04668 -2.23740 -2.59143 -2.67674 -10.3821 -14.7754 -0.950251 -0.619929 -1.38607 -1.76996 -0.622582 -1.30650 -1.19639 0.499982 1.65860
#> iter: 1 value: -225.0864 mgc: 5.852736 ustep: 1
#> iter: 2 value: -225.0867 mgc: 0.02807385 ustep: 1
#> iter: 3 value: -225.0867 mgc: 2.151257e-06 ustep: 1
#> iter: 4 mgc: 2.271516e-13
#> iter: 1 value: -231.4924 mgc: 3.476877 ustep: 1
#> iter: 2 value: -231.4924 mgc: 0.004774861 ustep: 1
#> iter: 3 value: -231.4924 mgc: 4.960289e-08 ustep: 1
#> mgc: 2.78888e-13
#> iter: 1 mgc: 2.78888e-13
#> outer mgc: 18.59683
#> 59: 220.94253: -6.65476 -1.79641 -1.81703 -1.82277 -2.04551 -2.24553 -2.34611 -2.76157 -10.4600 -14.8777 -0.939591 -0.617503 -1.40110 -1.78147 -0.629024 -1.31836 -1.17624 0.324951 1.65491
#> iter: 1 value: -234.5052 mgc: 2.300451 ustep: 1
#> iter: 2 value: -234.5052 mgc: 0.002512284 ustep: 1
#> iter: 3 value: -234.5052 mgc: 2.006456e-08 ustep: 1
#> mgc: 2.194911e-13
#> iter: 1 mgc: 2.194911e-13
#> outer mgc: 17.1191
#> 60: 216.65087: -6.69249 -1.89226 -1.85352 -1.89150 -2.05780 -2.22835 -2.48986 -2.68386 -10.4927 -14.6104 -0.939162 -0.613562 -1.41568 -1.76536 -0.632505 -1.32408 -1.15896 0.248280 1.65026
#> iter: 1 value: -236.2125 mgc: 2.900023 ustep: 1
#> iter: 2 value: -236.2125 mgc: 0.01408111 ustep: 1
#> iter: 3 value: -236.2125 mgc: 7.412505e-07 ustep: 1
#> iter: 4 mgc: 2.382261e-13
#> iter: 1 value: -238.93 mgc: 2.171608 ustep: 1
#> iter: 2 value: -238.93 mgc: 0.001644609 ustep: 1
#> iter: 3 value: -238.93 mgc: 1.550874e-08 ustep: 1
#> mgc: 2.629008e-13
#> iter: 1 mgc: 2.629008e-13
#> outer mgc: 21.38179
#> 61: 213.69580: -6.69553 -1.91440 -1.89954 -1.89274 -2.08345 -2.27300 -2.35419 -2.76422 -10.5349 -14.4594 -0.937384 -0.613653 -1.41846 -1.76890 -0.637775 -1.34492 -1.14503 0.0122615 1.64873
#> iter: 1 value: -240.6288 mgc: 2.846181 ustep: 1
#> iter: 2 value: -240.6289 mgc: 0.0138144 ustep: 1
#> iter: 3 value: -240.6289 mgc: 6.799435e-07 ustep: 1
#> iter: 4 mgc: 2.179923e-13
#> iter: 1 value: -241.4081 mgc: 1.974102 ustep: 1
#> iter: 2 value: -241.4081 mgc: 0.002982863 ustep: 1
#> iter: 3 value: -241.4081 mgc: 3.6016e-08 ustep: 1
#> mgc: 2.44027e-13
#> iter: 1 mgc: 2.44027e-13
#> outer mgc: 12.31044
#> 62: 210.30148: -6.65756 -1.88965 -1.85381 -1.88032 -2.04189 -2.19581 -2.43381 -2.64624 -10.4892 -14.6492 -0.933875 -0.617809 -1.41896 -1.76482 -0.638512 -1.31893 -1.14605 -0.169081 1.65855
#> iter: 1 value: -244.5177 mgc: 2.765156 ustep: 1
#> iter: 2 value: -244.5177 mgc: 0.007969463 ustep: 1
#> iter: 3 value: -244.5177 mgc: 2.355039e-07 ustep: 1
#> iter: 4 mgc: 2.460254e-13
#> iter: 1 mgc: 2.460254e-13
#> outer mgc: 12.87486
#> 63: 209.68898: -6.72882 -1.96844 -1.92097 -1.92959 -2.08904 -2.28139 -2.45479 -2.73823 -10.6074 -14.0725 -0.926288 -0.626170 -1.42560 -1.76929 -0.636753 -1.35033 -1.13773 -0.347285 1.66279
#> iter: 1 value: -243.1138 mgc: 4.110487 ustep: 1
#> iter: 2 value: -243.1138 mgc: 0.006563811 ustep: 1
#> iter: 3 value: -243.1138 mgc: 4.820297e-07 ustep: 1
#> iter: 4 mgc: 2.442491e-13
#> iter: 1 mgc: 2.442491e-13
#> outer mgc: 17.42471
#> 64: 205.48480: -6.66686 -1.85024 -1.87260 -1.88507 -2.09422 -2.20744 -2.37456 -2.70831 -10.3738 -14.4631 -0.914862 -0.632004 -1.39591 -1.75553 -0.637807 -1.32398 -1.13279 -0.762326 1.66161
#> iter: 1 value: -242.6526 mgc: 4.414207 ustep: 1
#> iter: 2 value: -242.6533 mgc: 0.06714703 ustep: 1
#> iter: 3 value: -242.6533 mgc: 1.661153e-05 ustep: 1
#> iter: 4 mgc: 1.163514e-12
#> iter: 1 value: -247.1272 mgc: 1.227091 ustep: 1
#> iter: 2 value: -247.1272 mgc: 0.0007084528 ustep: 1
#> iter: 3 mgc: 3.667273e-09
#> iter: 1 mgc: 3.667273e-09
#> outer mgc: 4.819222
#> 65: 203.36656: -6.66953 -1.84948 -1.84185 -1.88435 -2.04402 -2.21797 -2.43017 -2.65893 -10.4748 -14.4027 -0.919640 -0.633357 -1.41125 -1.75368 -0.637653 -1.33888 -1.12650 -0.782312 1.65851
#> iter: 1 value: -247.4729 mgc: 1.285416 ustep: 1
#> iter: 2 value: -247.4729 mgc: 0.000761605 ustep: 1
#> iter: 3 mgc: 1.264633e-09
#> iter: 1 value: -247.6612 mgc: 0.6021967 ustep: 1
#> iter: 2 value: -247.6612 mgc: 0.00016919 ustep: 1
#> iter: 3 mgc: 6.257128e-11
#> iter: 1 mgc: 6.257128e-11
#> outer mgc: 3.451336
#> 66: 203.02678: -6.66501 -1.87746 -1.85461 -1.86453 -2.04600 -2.20890 -2.40027 -2.67702 -10.4815 -14.3636 -0.922722 -0.634835 -1.41290 -1.75257 -0.637194 -1.34639 -1.13110 -0.814205 1.65223
#> iter: 1 value: -247.0176 mgc: 0.6215022 ustep: 1
#> iter: 2 value: -247.0176 mgc: 0.0002910204 ustep: 1
#> iter: 3 mgc: 1.195392e-10
#> iter: 1 mgc: 1.195392e-10
#> outer mgc: 2.979586
#> 67: 202.85608: -6.65351 -1.86097 -1.83774 -1.87375 -2.02041 -2.20145 -2.40920 -2.64442 -10.4911 -14.3496 -0.923336 -0.640982 -1.41761 -1.74591 -0.635357 -1.34883 -1.13532 -0.862078 1.64489
#> iter: 1 value: -245.7694 mgc: 0.6136477 ustep: 1
#> iter: 2 value: -245.7694 mgc: 0.0004699427 ustep: 1
#> iter: 3 mgc: 6.852101e-10
#> iter: 1 mgc: 6.852101e-10
#> outer mgc: 2.029581
#> 68: 202.65536: -6.63524 -1.85073 -1.83013 -1.84594 -2.02257 -2.17804 -2.38249 -2.63839 -10.4559 -14.3199 -0.926761 -0.650602 -1.41205 -1.74080 -0.632496 -1.34796 -1.13772 -0.883352 1.63663
#> iter: 1 value: -247.1845 mgc: 0.6252408 ustep: 1
#> iter: 2 value: -247.1845 mgc: 0.0003957892 ustep: 1
#> iter: 3 mgc: 6.133949e-10
#> iter: 1 mgc: 6.133949e-10
#> outer mgc: 4.032042
#> 69: 202.61987: -6.62334 -1.82844 -1.80857 -1.83644 -1.98722 -2.17947 -2.35381 -2.62646 -10.4626 -14.3210 -0.930703 -0.651603 -1.41745 -1.74202 -0.634020 -1.35583 -1.14144 -0.925838 1.63471
#> iter: 1 value: -244.8033 mgc: 0.6371171 ustep: 1
#> iter: 2 value: -244.8033 mgc: 0.000546944 ustep: 1
#> iter: 3 mgc: 1.033094e-09
#> iter: 1 mgc: 1.033094e-09
#> outer mgc: 4.064391
#> 70: 202.46878: -6.61132 -1.82078 -1.79402 -1.82183 -1.97390 -2.14119 -2.36210 -2.59391 -10.4344 -14.2942 -0.933864 -0.644140 -1.41644 -1.73106 -0.635993 -1.35054 -1.14511 -0.935862 1.62294
#> iter: 1 value: -245.4264 mgc: 0.8387518 ustep: 1
#> iter: 2 value: -245.4264 mgc: 0.0003983801 ustep: 1
#> iter: 3 mgc: 8.306641e-10
#> iter: 1 mgc: 8.306641e-10
#> outer mgc: 3.939819
#> 71: 202.45378: -6.58717 -1.80923 -1.78991 -1.78694 -1.98193 -2.12818 -2.31131 -2.59551 -10.4149 -14.2853 -0.934896 -0.646533 -1.41082 -1.73742 -0.635926 -1.35216 -1.14391 -0.945157 1.62284
#> iter: 1 value: -245.4419 mgc: 0.3303159 ustep: 1
#> iter: 2 value: -245.4419 mgc: 9.126391e-05 ustep: 1
#> iter: 3 mgc: 3.469086e-11
#> iter: 1 mgc: 3.469086e-11
#> outer mgc: 0.7741606
#> 72: 202.31246: -6.58415 -1.80263 -1.77892 -1.79884 -1.96189 -2.12546 -2.32557 -2.57661 -10.4153 -14.2791 -0.935234 -0.647438 -1.41361 -1.73590 -0.635964 -1.35159 -1.14479 -0.948367 1.62109
#> iter: 1 value: -245.3795 mgc: 0.3478868 ustep: 1
#> iter: 2 value: -245.3795 mgc: 0.000127421 ustep: 1
#> iter: 3 mgc: 6.151271e-11
#> iter: 1 mgc: 6.151271e-11
#> outer mgc: 1.473726
#> 73: 202.26611: -6.57602 -1.79191 -1.77369 -1.78416 -1.95525 -2.11413 -2.30933 -2.57207 -10.4063 -14.2651 -0.936403 -0.649765 -1.41610 -1.73431 -0.635485 -1.35401 -1.14782 -0.960267 1.61688
#> iter: 1 value: -245.3241 mgc: 0.5158818 ustep: 1
#> iter: 2 value: -245.3241 mgc: 0.0008983711 ustep: 1
#> iter: 3 mgc: 3.498579e-09
#> iter: 1 value: -246.2358 mgc: 1.501823 ustep: 1
#> iter: 2 value: -246.2358 mgc: 0.007590922 ustep: 1
#> iter: 3 value: -246.2358 mgc: 2.477182e-07 ustep: 1
#> iter: 4 mgc: 2.700062e-13
#> iter: 1 value: -249.042 mgc: 6.051362 ustep: 1
#> iter: 2 value: -249.0499 mgc: 0.1439669 ustep: 1
#> iter: 3 value: -249.0499 mgc: 0.0001052933 ustep: 1
#> iter: 4 mgc: 6.158096e-11
#> iter: 1 mgc: 6.158096e-11
#> outer mgc: 11.31223
#> 74: 200.80403: -6.20276 -1.49308 -1.40953 -1.38330 -1.61773 -1.70385 -1.89753 -2.12465 -10.0806 -13.9582 -0.967360 -0.663580 -1.40844 -1.74558 -0.658881 -1.36626 -1.15451 -1.01481 1.58777
#> iter: 1 value: -254.5399 mgc: 23.5967 ustep: 1
#> iter: 2 value: -259.4766 mgc: 3.85769 ustep: 1
#> iter: 3 value: -259.4808 mgc: 0.121422 ustep: 1
#> iter: 4 value: -259.4808 mgc: 0.0003793811 ustep: 1
#> iter: 5 mgc: 1.369375e-09
#> iter: 1 value: 4880.884 mgc: 78.50725 ustep: 0.1336527
#> iter: 2 value: 1830.565 mgc: 1392.328 ustep: 0.3656489
#> iter: 3 value: 1291.88 mgc: 397.0308 ustep: 0.6047286
#> iter: 4 value: 1196.644 mgc: 152.598 ustep: 0.7776652
#> iter: 5 value: 1186.624 mgc: 46.74395 ustep: 0.8818651
#> iter: 6 value: 1186.373 mgc: 7.791282 ustep: 0.9390828
#> iter: 7 value: 1186.373 mgc: 0.2917598 ustep: 0.969066
#> iter: 8 value: 1186.373 mgc: 0.0004788979 ustep: 0.984413
#> iter: 9 value: 1186.373 mgc: 3.302259e-07 ustep: 0.9921767
#> iter: 10 mgc: 1.988492e-09
#> iter: 1 mgc: 1.369375e-09
#> outer mgc: 23.62966
#> 75: 192.70173: -5.15198 -0.614955 -0.421697 -0.312042 -0.593685 -0.587578 -0.739963 -0.916545 -9.16627 -13.0766 -1.02820 -0.664201 -1.42007 -1.75127 -0.678929 -1.38054 -1.17023 -1.14634 1.51117
#> iter: 1 value: 3160.221 mgc: 49.58408 ustep: 0.4367355
#> iter: 2 value: 665.3024 mgc: 658.7087 ustep: 0.6608936
#> iter: 3 value: 158.7354 mgc: 201.166 ustep: 0.8129723
#> iter: 4 value: 60.54585 mgc: 73.82437 ustep: 0.9016596
#> iter: 5 value: 52.00882 mgc: 20.98796 ustep: 0.9495626
#> iter: 6 value: 51.88574 mgc: 2.807289 ustep: 0.9744576
#> iter: 7 value: 51.88569 mgc: 0.06211556 ustep: 0.9871475
#> iter: 8 value: 51.88569 mgc: 3.786837e-05 ustep: 0.9935536
#> iter: 9 mgc: 2.664112e-09
#> iter: 1 value: -263.1113 mgc: 5.026129 ustep: 1
#> iter: 2 value: -263.1599 mgc: 0.198795 ustep: 1
#> iter: 3 value: -263.1599 mgc: 0.0008600085 ustep: 1
#> iter: 4 value: -263.1599 mgc: 1.666584e-08 ustep: 1
#> mgc: 2.469136e-13
#> iter: 1 mgc: 2.469136e-13
#> outer mgc: 24.08432
#> 76: 189.68030: -4.93752 -0.429718 -0.233694 -0.119206 -0.373481 -0.368606 -0.505469 -0.667493 -8.98887 -12.9021 -1.03542 -0.664497 -1.43735 -1.75032 -0.681436 -1.37941 -1.17206 -1.16960 1.49480
#> iter: 1 value: -202.4227 mgc: 22.00633 ustep: 1
#> iter: 2 value: -202.697 mgc: 1.603614 ustep: 1
#> iter: 3 value: -202.6971 mgc: 0.03297971 ustep: 1
#> iter: 4 value: -202.6971 mgc: 2.078969e-05 ustep: 1
#> iter: 5 mgc: 6.104783e-12
#> iter: 1 value: -271.8149 mgc: 3.646308 ustep: 1
#> iter: 2 value: -271.8152 mgc: 0.04016352 ustep: 1
#> iter: 3 value: -271.8152 mgc: 1.85058e-05 ustep: 1
#> iter: 4 mgc: 1.943334e-12
#> iter: 1 mgc: 1.943334e-12
#> outer mgc: 10.27746
#> 77: 184.06221: -4.90588 -0.355902 -0.298775 -0.263685 -0.274527 -0.392124 -0.479609 -0.592420 -9.02054 -12.9086 -0.999952 -0.662780 -1.55300 -1.73815 -0.669112 -1.36054 -1.16735 -1.15033 1.49427
#> iter: 1 value: -285.3616 mgc: 3.022192 ustep: 1
#> iter: 2 value: -285.3624 mgc: 0.05556201 ustep: 1
#> iter: 3 value: -285.3624 mgc: 4.149847e-05 ustep: 1
#> iter: 4 mgc: 1.833644e-11
#> iter: 1 mgc: 1.833644e-11
#> outer mgc: 23.01971
#> 78: 182.47200: -4.87216 -0.417702 -0.337493 -0.209873 -0.381811 -0.303108 -0.382840 -0.501560 -9.01893 -12.9236 -0.980712 -0.663803 -1.67675 -1.73775 -0.652937 -1.38358 -1.16122 -1.11488 1.51742
#> iter: 1 value: -285.8851 mgc: 3.052018 ustep: 1
#> iter: 2 value: -285.8874 mgc: 0.0681079 ustep: 1
#> iter: 3 value: -285.8874 mgc: 0.0001221851 ustep: 1
#> iter: 4 mgc: 1.747833e-10
#> iter: 1 mgc: 1.747833e-10
#> outer mgc: 7.566694
#> 79: 179.22617: -4.80134 -0.371712 -0.278265 -0.217381 -0.208416 -0.265025 -0.294466 -0.388664 -9.00705 -12.9023 -0.978911 -0.663358 -1.68774 -1.71825 -0.648606 -1.35569 -1.15606 -1.09909 1.51821
#> iter: 1 value: -295.7796 mgc: 4.47093 ustep: 1
#> iter: 2 value: -295.8584 mgc: 0.4761009 ustep: 1
#> iter: 3 value: -295.8584 mgc: 0.003849439 ustep: 1
#> iter: 4 value: -295.8584 mgc: 1.850165e-07 ustep: 1
#> iter: 5 mgc: 4.294343e-13
#> iter: 1 value: -295.927 mgc: 7.06292 ustep: 1
#> iter: 2 value: -297.2038 mgc: 2.485095 ustep: 1
#> iter: 3 value: -297.2049 mgc: 0.1110587 ustep: 1
#> iter: 4 value: -297.2049 mgc: 0.0001676536 ustep: 1
#> iter: 5 mgc: 3.764749e-10
#> iter: 1 mgc: 3.764749e-10
#> outer mgc: 26.43212
#> 80: 172.41069: -4.29941 -0.0916132 0.0318505 0.223366 0.200516 0.290355 0.317261 0.225136 -8.70292 -12.6244 -1.00596 -0.653796 -1.78845 -1.66383 -0.633573 -1.37122 -1.12140 -1.06862 1.48455
#> iter: 1 value: -20.22358 mgc: 26.67408 ustep: 1
#> iter: 2 value: -250.9567 mgc: 154.8237 ustep: 1
#> iter: 3 value: -286.456 mgc: 45.5921 ustep: 1
#> iter: 4 value: -288.498 mgc: 10.39774 ustep: 1
#> iter: 5 value: -288.5086 mgc: 0.7434939 ustep: 1
#> iter: 6 value: -288.5086 mgc: 0.00426641 ustep: 1
#> iter: 7 value: -288.5086 mgc: 1.535395e-07 ustep: 1
#> mgc: 4.249934e-13
#> iter: 1 value: -310.2553 mgc: 3.00397 ustep: 1
#> iter: 2 value: -310.2594 mgc: 0.1814172 ustep: 1
#> iter: 3 value: -310.2594 mgc: 0.0005821772 ustep: 1
#> iter: 4 mgc: 5.564128e-09
#> iter: 1 mgc: 5.564128e-09
#> outer mgc: 12.47406
#> 81: 169.20175: -4.31189 -0.0763548 0.0385915 0.138012 0.189090 0.223543 0.256227 0.379597 -8.73622 -12.6446 -1.02173 -0.660633 -1.83183 -1.67658 -0.631375 -1.37077 -1.12805 -1.06558 1.51061
#> iter: 1 value: -325.9441 mgc: 2.893596 ustep: 1
#> iter: 2 value: -325.9455 mgc: 0.09384114 ustep: 1
#> iter: 3 value: -325.9455 mgc: 0.000207897 ustep: 1
#> iter: 4 mgc: 8.545255e-10
#> iter: 1 mgc: 8.545255e-10
#> outer mgc: 8.15124
#> 82: 168.23377: -4.31705 -0.117685 -0.0202771 0.169137 0.156543 0.294264 0.380995 0.426249 -8.77138 -12.6952 -1.02283 -0.696769 -1.89982 -1.68643 -0.639930 -1.36723 -1.13283 -1.05589 1.57765
#> iter: 1 value: -328.6732 mgc: 2.890587 ustep: 1
#> iter: 2 value: -328.6746 mgc: 0.1460109 ustep: 1
#> iter: 3 value: -328.6746 mgc: 0.0002640566 ustep: 1
#> iter: 4 mgc: 1.094781e-09
#> iter: 1 mgc: 1.094781e-09
#> outer mgc: 6.135686
#> 83: 167.98732: -4.39421 -0.176300 -0.00568997 0.114120 0.176806 0.248974 0.301797 0.371181 -8.81890 -12.7540 -0.981509 -0.731559 -1.96254 -1.64842 -0.633270 -1.32186 -1.14752 -1.04804 1.63844
#> iter: 1 value: -328.5402 mgc: 2.354379 ustep: 1
#> iter: 2 value: -328.5452 mgc: 0.2517745 ustep: 1
#> iter: 3 value: -328.5452 mgc: 0.001015325 ustep: 1
#> iter: 4 value: -328.5452 mgc: 1.349816e-08 ustep: 1
#> mgc: 4.884981e-13
#> iter: 1 value: -331.6944 mgc: 0.966842 ustep: 1
#> iter: 2 value: -331.6944 mgc: 0.02603713 ustep: 1
#> iter: 3 value: -331.6944 mgc: 1.001434e-05 ustep: 1
#> iter: 4 mgc: 1.254108e-12
#> iter: 1 mgc: 1.254108e-12
#> outer mgc: 3.131887
#> 84: 167.55442: -4.36746 -0.159371 -0.0350965 0.122230 0.159027 0.235076 0.327904 0.404200 -8.79457 -12.7143 -0.979952 -0.724482 -1.97680 -1.65437 -0.640721 -1.35231 -1.15750 -1.03502 1.63222
#> iter: 1 value: -331.5399 mgc: 0.9287581 ustep: 1
#> iter: 2 value: -331.5399 mgc: 0.007704394 ustep: 1
#> iter: 3 value: -331.5399 mgc: 1.482021e-06 ustep: 1
#> iter: 4 mgc: 3.756995e-13
#> iter: 1 value: -331.7277 mgc: 0.4308506 ustep: 1
#> iter: 2 value: -331.7277 mgc: 0.0008838052 ustep: 1
#> iter: 3 value: -331.7277 mgc: 1.287325e-08 ustep: 1
#> mgc: 3.410605e-13
#> iter: 1 mgc: 3.410605e-13
#> outer mgc: 2.449023
#> 85: 167.51857: -4.35694 -0.154598 -0.0158827 0.112231 0.160880 0.259013 0.323796 0.397167 -8.80149 -12.7150 -0.979986 -0.720924 -1.98019 -1.65483 -0.641567 -1.34635 -1.16308 -1.02849 1.62997
#> iter: 1 value: -332.4861 mgc: 0.4074901 ustep: 1
#> iter: 2 value: -332.4861 mgc: 0.004308336 ustep: 1
#> iter: 3 value: -332.4861 mgc: 3.551683e-07 ustep: 1
#> iter: 4 mgc: 5.155876e-13
#> iter: 1 mgc: 5.155876e-13
#> outer mgc: 1.060942
#> 86: 167.46408: -4.34366 -0.149871 -0.0168779 0.125677 0.173692 0.253410 0.335692 0.415569 -8.79171 -12.7091 -0.980970 -0.709324 -1.98531 -1.66002 -0.644185 -1.34343 -1.16641 -1.02516 1.62550
#> iter: 1 value: -331.2944 mgc: 0.9280322 ustep: 1
#> iter: 2 value: -331.2944 mgc: 0.001192809 ustep: 1
#> iter: 3 mgc: 5.731144e-09
#> iter: 1 value: -332.9273 mgc: 0.1477072 ustep: 1
#> iter: 2 value: -332.9273 mgc: 5.119212e-05 ustep: 1
#> iter: 3 mgc: 3.836337e-11
#> iter: 1 mgc: 3.836337e-11
#> outer mgc: 0.9358055
#> 87: 167.45689: -4.34076 -0.148790 -0.0177489 0.127596 0.166855 0.260846 0.338005 0.413392 -8.79167 -12.7084 -0.983246 -0.709197 -1.98693 -1.65901 -0.644195 -1.34643 -1.16971 -1.02126 1.62569
#> iter: 1 value: -332.3881 mgc: 0.2338707 ustep: 1
#> iter: 2 value: -332.3881 mgc: 0.0003353921 ustep: 1
#> iter: 3 mgc: 1.594557e-09
#> iter: 1 mgc: 1.594557e-09
#> outer mgc: 0.3073223
#> 88: 167.45039: -4.33880 -0.148311 -0.0152766 0.126990 0.172666 0.258721 0.339382 0.418209 -8.79094 -12.7077 -0.985971 -0.710713 -1.98811 -1.65164 -0.643764 -1.34790 -1.17174 -1.01690 1.62351
#> iter: 1 value: -333.4255 mgc: 0.2302157 ustep: 1
#> iter: 2 value: -333.4255 mgc: 0.0001004639 ustep: 1
#> iter: 3 mgc: 8.174794e-11
#> iter: 1 mgc: 8.174794e-11
#> outer mgc: 0.3000414
#> 89: 167.44786: -4.33718 -0.147548 -0.0148759 0.129265 0.173037 0.261164 0.339980 0.418680 -8.78907 -12.7061 -0.988185 -0.714443 -1.98983 -1.65915 -0.644791 -1.34304 -1.17664 -1.01220 1.62306
#> iter: 1 value: -333.0047 mgc: 0.2014773 ustep: 1
#> iter: 2 value: -333.0047 mgc: 0.0001890199 ustep: 1
#> iter: 3 mgc: 4.915446e-10
#> iter: 1 mgc: 4.915446e-10
#> outer mgc: 0.425712
#> 90: 167.44657: -4.33487 -0.146907 -0.0162807 0.129334 0.171222 0.260938 0.343243 0.420618 -8.79133 -12.7080 -0.988570 -0.708078 -1.99229 -1.65428 -0.644182 -1.34756 -1.17962 -1.00736 1.61949
#> iter: 1 value: -333.1656 mgc: 0.2896099 ustep: 1
#> iter: 2 value: -333.1656 mgc: 8.274036e-05 ustep: 1
#> iter: 3 mgc: 1.336377e-10
#> iter: 1 value: -333.0199 mgc: 0.06214588 ustep: 1
#> iter: 2 value: -333.0199 mgc: 1.520976e-05 ustep: 1
#> iter: 3 mgc: 4.209522e-12
#> iter: 1 mgc: 4.209522e-12
#> outer mgc: 0.2243897
#> 91: 167.44524: -4.33496 -0.148448 -0.0150324 0.128363 0.173835 0.261369 0.340988 0.421004 -8.79010 -12.7076 -0.989228 -0.708412 -1.99208 -1.65418 -0.644213 -1.34529 -1.17995 -1.00709 1.62011
#> iter: 1 value: -333.1738 mgc: 0.1188129 ustep: 1
#> iter: 2 value: -333.1738 mgc: 2.410588e-05 ustep: 1
#> iter: 3 mgc: 6.7919e-12
#> iter: 1 mgc: 6.7919e-12
#> outer mgc: 0.1673713
#> 92: 167.44521: -4.33516 -0.147578 -0.0143027 0.129770 0.172892 0.261378 0.342565 0.420140 -8.78939 -12.7062 -0.989562 -0.710490 -1.99260 -1.65245 -0.644451 -1.34490 -1.18078 -1.00693 1.62272
#> iter: 1 value: -333.2904 mgc: 0.0348062 ustep: 1
#> iter: 2 value: -333.2904 mgc: 1.934389e-05 ustep: 1
#> iter: 3 mgc: 5.823786e-12
#> iter: 1 mgc: 5.823786e-12
#> outer mgc: 0.1659818
#> 93: 167.44491: -4.33478 -0.147156 -0.0151749 0.129248 0.173678 0.261679 0.341552 0.421277 -8.78971 -12.7066 -0.989321 -0.710614 -1.99288 -1.65338 -0.644464 -1.34560 -1.18095 -1.00670 1.62255
#> iter: 1 value: -333.2805 mgc: 0.07839555 ustep: 1
#> iter: 2 value: -333.2805 mgc: 2.819603e-06 ustep: 1
#> iter: 3 mgc: 4.442557e-13
#> iter: 1 mgc: 4.442557e-13
#> outer mgc: 0.07631744
#> 94: 167.44480: -4.33469 -0.147619 -0.0146718 0.129266 0.173379 0.261622 0.342528 0.421083 -8.79000 -12.7068 -0.989125 -0.710300 -1.99304 -1.65477 -0.644554 -1.34527 -1.18116 -1.00638 1.62105
#> iter: 1 value: -333.3213 mgc: 0.03199395 ustep: 1
#> iter: 2 value: -333.3213 mgc: 3.665691e-06 ustep: 1
#> iter: 3 mgc: 4.218847e-13
#> iter: 1 value: -333.2625 mgc: 0.010314 ustep: 1
#> iter: 2 value: -333.2625 mgc: 1.549262e-07 ustep: 1
#> iter: 3 mgc: 3.259615e-13
#> iter: 1 mgc: 3.259615e-13
#> outer mgc: 0.02922481
#> 95: 167.44475: -4.33460 -0.147586 -0.0150537 0.129408 0.173442 0.261748 0.342168 0.421267 -8.78984 -12.7069 -0.989216 -0.710260 -1.99296 -1.65457 -0.644544 -1.34507 -1.18123 -1.00628 1.62115
#> iter: 1 value: -333.2852 mgc: 0.01771874 ustep: 1
#> iter: 2 value: -333.2852 mgc: 1.691873e-07 ustep: 1
#> iter: 3 mgc: 5.537792e-13
#> iter: 1 mgc: 5.537792e-13
#> outer mgc: 0.03942717
#> 96: 167.44474: -4.33443 -0.147607 -0.0149097 0.129312 0.173444 0.261736 0.342465 0.421179 -8.78991 -12.7070 -0.989313 -0.709961 -1.99290 -1.65440 -0.644594 -1.34506 -1.18138 -1.00611 1.62157
#> iter: 1 value: -333.2778 mgc: 0.01081539 ustep: 1
#> iter: 2 value: -333.2778 mgc: 1.71786e-06 ustep: 1
#> iter: 3 mgc: 4.045653e-13
#> iter: 1 mgc: 4.045653e-13
#> outer mgc: 0.02944055
#> 97: 167.44474: -4.33442 -0.147592 -0.0147930 0.129450 0.173488 0.261901 0.342430 0.421511 -8.78983 -12.7069 -0.989119 -0.710119 -1.99298 -1.65415 -0.644463 -1.34516 -1.18156 -1.00584 1.62184
#> iter: 1 value: -333.2908 mgc: 0.007363089 ustep: 1
#> iter: 2 value: -333.2908 mgc: 1.532354e-06 ustep: 1
#> iter: 3 mgc: 4.63185e-13
#> iter: 1 value: -333.287 mgc: 0.003520262 ustep: 1
#> iter: 2 value: -333.287 mgc: 1.46816e-07 ustep: 1
#> mgc: 5.391243e-13
#> iter: 1 mgc: 5.391243e-13
#> outer mgc: 0.0158172
#> 98: 167.44473: -4.33437 -0.147544 -0.0149830 0.129454 0.173585 0.261836 0.342490 0.421425 -8.78982 -12.7069 -0.989149 -0.710118 -1.99301 -1.65419 -0.644474 -1.34510 -1.18158 -1.00580 1.62184
#> iter: 1 value: -333.2926 mgc: 0.002802459 ustep: 1
#> iter: 2 value: -333.2926 mgc: 1.218213e-07 ustep: 1
#> iter: 3 mgc: 6.212808e-13
#> iter: 1 mgc: 6.212808e-13
#> outer mgc: 0.004498458
#> 99: 167.44473: -4.33439 -0.147616 -0.0149535 0.129414 0.173500 0.261794 0.342363 0.421362 -8.78987 -12.7070 -0.989196 -0.710041 -1.99303 -1.65424 -0.644563 -1.34509 -1.18162 -1.00571 1.62178
#> iter: 1 value: -333.2868 mgc: 0.003193784 ustep: 1
#> iter: 2 mgc: 7.788166e-09
#> iter: 1 value: -333.293 mgc: 0.001260446 ustep: 1
#> iter: 2 mgc: 7.980207e-09
#> iter: 1 value: -333.2928 mgc: 0.0005981356 ustep: 1
#> iter: 2 mgc: 1.797115e-09
#> iter: 1 mgc: 1.797115e-09
#> outer mgc: 0.002671636
#> 100: 167.44473: -4.33439 -0.147608 -0.0149547 0.129396 0.173492 0.261805 0.342389 0.421350 -8.78987 -12.7070 -0.989192 -0.710051 -1.99304 -1.65425 -0.644555 -1.34508 -1.18163 -1.00569 1.62178
#> iter: 1 value: -333.2943 mgc: 0.0005902048 ustep: 1
#> iter: 2 mgc: 4.147604e-09
#> iter: 1 mgc: 4.147604e-09
#> outer mgc: 0.003926587
#> 101: 167.44473: -4.33438 -0.147614 -0.0149589 0.129406 0.173497 0.261800 0.342371 0.421372 -8.78987 -12.7070 -0.989191 -0.710065 -1.99305 -1.65426 -0.644544 -1.34508 -1.18163 -1.00566 1.62177
#> iter: 1 value: -333.2922 mgc: 0.001459419 ustep: 1
#> iter: 2 mgc: 2.843849e-09
#> iter: 1 mgc: 2.843849e-09
#> outer mgc: 0.001822165
#> 102: 167.44473: -4.33438 -0.147629 -0.0149539 0.129404 0.173491 0.261802 0.342384 0.421352 -8.78988 -12.7070 -0.989206 -0.710138 -1.99305 -1.65424 -0.644526 -1.34508 -1.18163 -1.00564 1.62174
#> iter: 1 value: -333.2951 mgc: 0.001413151 ustep: 1
#> iter: 2 mgc: 5.301692e-09
#> iter: 1 value: -333.2939 mgc: 0.0005192793 ustep: 1
#> iter: 2 mgc: 2.06576e-09
#> iter: 1 mgc: 2.06576e-09
#> outer mgc: 0.002001365
#> 103: 167.44473: -4.33438 -0.147611 -0.0149649 0.129390 0.173499 0.261804 0.342382 0.421366 -8.78988 -12.7070 -0.989210 -0.710123 -1.99306 -1.65425 -0.644524 -1.34507 -1.18164 -1.00562 1.62174
#> iter: 1 value: -333.2958 mgc: 0.000453828 ustep: 1
#> iter: 2 mgc: 1.090661e-09
#> iter: 1 mgc: 1.090661e-09
#> outer mgc: 0.001402065
#> 104: 167.44473: -4.33438 -0.147616 -0.0149616 0.129406 0.173494 0.261812 0.342393 0.421367 -8.78988 -12.7070 -0.989205 -0.710095 -1.99307 -1.65426 -0.644522 -1.34507 -1.18164 -1.00561 1.62174
#> iter: 1 value: -333.296 mgc: 0.0003083809 ustep: 1
#> iter: 2 mgc: 1.315131e-09
#> iter: 1 mgc: 1.315131e-09
#> outer mgc: 0.001619388
#> 105: 167.44473: -4.33437 -0.147618 -0.0149474 0.129404 0.173502 0.261812 0.342399 0.421375 -8.78987 -12.7070 -0.989198 -0.710109 -1.99308 -1.65425 -0.644547 -1.34507 -1.18164 -1.00560 1.62175
#> iter: 1 value: -333.2968 mgc: 0.0004732299 ustep: 1
#> iter: 2 mgc: 1.200474e-09
#> iter: 1 value: -333.2964 mgc: 0.0001998044 ustep: 1
#> iter: 2 mgc: 2.035009e-10
#> iter: 1 mgc: 2.035009e-10
#> outer mgc: 0.0009009657
#> 106: 167.44473: -4.33437 -0.147613 -0.0149578 0.129405 0.173502 0.261815 0.342393 0.421379 -8.78987 -12.7070 -0.989203 -0.710107 -1.99307 -1.65425 -0.644538 -1.34507 -1.18165 -1.00560 1.62176
#> iter: 1 value: -333.2963 mgc: 0.0001168454 ustep: 1
#> iter: 2 mgc: 1.902016e-10
#> iter: 1 mgc: 1.902016e-10
#> outer mgc: 0.0005125852
#> 107: 167.44473: -4.33437 -0.147614 -0.0149537 0.129403 0.173502 0.261810 0.342396 0.421375 -8.78987 -12.7070 -0.989204 -0.710102 -1.99307 -1.65425 -0.644523 -1.34507 -1.18165 -1.00560 1.62176
#> iter: 1 value: -333.296 mgc: 0.0001092468 ustep: 1
#> iter: 2 mgc: 1.238676e-10
#> 108: 167.44473: -4.33437 -0.147614 -0.0149537 0.129403 0.173502 0.261810 0.342396 0.421375 -8.78987 -12.7070 -0.989204 -0.710102 -1.99307 -1.65425 -0.644523 -1.34507 -1.18165 -1.00560 1.62176
#> iter: 1 mgc: 1.902016e-10
#> iter: 1 mgc: 1.902016e-10
#> outer mgc: 0.0005125852
#> iter: 1 value: -333.2968 mgc: 0.001706019 ustep: 1
#> iter: 2 value: -333.2968 mgc: 2.131715e-08 ustep: 1
#> mgc: 3.614886e-13
#> outer mgc: 0.07017049
#> iter: 1 value: -333.2958 mgc: 0.001706019 ustep: 1
#> iter: 2 value: -333.2958 mgc: 2.131584e-08 ustep: 1
#> mgc: 3.166356e-13
#> outer mgc: 0.07036467
#> iter: 1 value: -333.2964 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2964 mgc: 1.521443e-07 ustep: 1
#> iter: 3 mgc: 3.121947e-13
#> outer mgc: 0.1618723
#> iter: 1 value: -333.296 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.296 mgc: 1.521288e-07 ustep: 1
#> iter: 3 mgc: 4.551914e-13
#> outer mgc: 0.1617207
#> iter: 1 value: -333.2967 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2967 mgc: 1.229146e-06 ustep: 1
#> iter: 3 mgc: 4.565237e-13
#> outer mgc: 0.2199124
#> iter: 1 value: -333.2957 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2957 mgc: 1.228506e-06 ustep: 1
#> iter: 3 mgc: 5.582201e-13
#> outer mgc: 0.2195578
#> iter: 1 value: -333.2968 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2968 mgc: 1.178571e-06 ustep: 1
#> iter: 3 mgc: 3.159695e-13
#> outer mgc: 0.2211451
#> iter: 1 value: -333.2956 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2956 mgc: 1.177963e-06 ustep: 1
#> iter: 3 mgc: 5.839773e-13
#> outer mgc: 0.2218815
#> iter: 1 value: -333.2968 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2968 mgc: 1.331866e-06 ustep: 1
#> iter: 3 mgc: 3.224088e-13
#> outer mgc: 0.2214605
#> iter: 1 value: -333.2956 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2956 mgc: 1.331132e-06 ustep: 1
#> iter: 3 mgc: 3.717027e-13
#> outer mgc: 0.2208204
#> iter: 1 value: -333.2968 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2968 mgc: 1.447861e-06 ustep: 1
#> iter: 3 mgc: 4.289902e-13
#> outer mgc: 0.2178739
#> iter: 1 value: -333.2956 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2956 mgc: 1.446994e-06 ustep: 1
#> iter: 3 mgc: 3.743672e-13
#> outer mgc: 0.2188981
#> iter: 1 value: -333.2966 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2966 mgc: 1.80015e-06 ustep: 1
#> iter: 3 mgc: 3.193001e-13
#> outer mgc: 0.2104465
#> iter: 1 value: -333.2958 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2958 mgc: 1.798867e-06 ustep: 1
#> iter: 3 mgc: 4.116707e-13
#> outer mgc: 0.2101655
#> iter: 1 value: -333.2966 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2966 mgc: 1.184153e-05 ustep: 1
#> iter: 3 mgc: 2.19913e-12
#> outer mgc: 0.147267
#> iter: 1 value: -333.2959 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2959 mgc: 1.182156e-05 ustep: 1
#> iter: 3 mgc: 2.19913e-12
#> outer mgc: 0.147355
#> iter: 1 value: -333.2966 mgc: 0.01062593 ustep: 1
#> iter: 2 value: -333.2966 mgc: 3.380456e-07 ustep: 1
#> iter: 3 mgc: 4.249934e-13
#> outer mgc: 0.1767063
#> iter: 1 value: -333.2958 mgc: 0.01062593 ustep: 1
#> iter: 2 value: -333.2958 mgc: 3.379826e-07 ustep: 1
#> iter: 3 mgc: 3.63265e-13
#> outer mgc: 0.1765883
#> iter: 1 value: -333.2966 mgc: 0.007472262 ustep: 1
#> iter: 2 value: -333.2966 mgc: 2.317905e-07 ustep: 1
#> iter: 3 mgc: 4.110046e-13
#> outer mgc: 0.1187934
#> iter: 1 value: -333.2959 mgc: 0.007472262 ustep: 1
#> iter: 2 value: -333.2959 mgc: 2.3175e-07 ustep: 1
#> iter: 3 mgc: 4.343192e-13
#> outer mgc: 0.1184172
#> iter: 1 value: -333.1929 mgc: 0.02626242 ustep: 1
#> iter: 2 value: -333.1929 mgc: 5.145017e-07 ustep: 1
#> iter: 3 mgc: 3.517187e-13
#> outer mgc: 0.07859371
#> iter: 1 value: -333.3997 mgc: 0.026315 ustep: 1
#> iter: 2 value: -333.3997 mgc: 5.154156e-07 ustep: 1
#> iter: 3 mgc: 3.287925e-13
#> outer mgc: 0.07873186
#> iter: 1 value: -333.2936 mgc: 0.01464721 ustep: 1
#> iter: 2 value: -333.2936 mgc: 2.077272e-08 ustep: 1
#> mgc: 4.112266e-13
#> outer mgc: 0.02861136
#> iter: 1 value: -333.299 mgc: 0.01467653 ustep: 1
#> iter: 2 value: -333.299 mgc: 2.083e-08 ustep: 1
#> mgc: 3.606004e-13
#> outer mgc: 0.02856361
#> iter: 1 value: -333.1836 mgc: 0.02122111 ustep: 1
#> iter: 2 value: -333.1836 mgc: 1.554966e-06 ustep: 1
#> iter: 3 mgc: 4.574119e-13
#> outer mgc: 0.05249009
#> iter: 1 value: -333.409 mgc: 0.0212636 ustep: 1
#> iter: 2 value: -333.409 mgc: 1.558626e-06 ustep: 1
#> iter: 3 mgc: 5.402345e-13
#> outer mgc: 0.05241999
#> iter: 1 value: -333.1873 mgc: 0.02569488 ustep: 1
#> iter: 2 value: -333.1873 mgc: 1.144447e-06 ustep: 1
#> iter: 3 mgc: 4.054534e-13
#> outer mgc: 0.08226291
#> iter: 1 value: -333.4053 mgc: 0.02574632 ustep: 1
#> iter: 2 value: -333.4053 mgc: 1.145373e-06 ustep: 1
#> iter: 3 mgc: 4.092282e-13
#> outer mgc: 0.08233307
#> iter: 1 value: -333.296 mgc: 0.002474377 ustep: 1
#> iter: 2 value: -333.296 mgc: 3.148051e-08 ustep: 1
#> mgc: 3.872458e-13
#> outer mgc: 0.037938
#> iter: 1 value: -333.2966 mgc: 0.00247933 ustep: 1
#> iter: 2 value: -333.2966 mgc: 3.160221e-08 ustep: 1
#> mgc: 4.427569e-13
#> outer mgc: 0.03744565
#> iter: 1 value: -333.2597 mgc: 0.01703845 ustep: 1
#> iter: 2 value: -333.2597 mgc: 2.323142e-06 ustep: 1
#> iter: 3 mgc: 4.209966e-13
#> outer mgc: 0.1710651
#> iter: 1 value: -333.3328 mgc: 0.01707256 ustep: 1
#> iter: 2 value: -333.3328 mgc: 2.326511e-06 ustep: 1
#> iter: 3 mgc: 5.089262e-13
#> outer mgc: 0.1710632
#> iter: 1 value: -333.2932 mgc: 0.0123519 ustep: 1
#> iter: 2 value: -333.2932 mgc: 1.660562e-08 ustep: 1
#> mgc: 3.175238e-13
#> outer mgc: 0.03179032
#> iter: 1 value: -333.2994 mgc: 0.01237663 ustep: 1
#> iter: 2 value: -333.2994 mgc: 1.665806e-08 ustep: 1
#> mgc: 4.405365e-13
#> outer mgc: 0.03184436
#> iter: 1 value: -333.293 mgc: 0.01440225 ustep: 1
#> iter: 2 value: -333.293 mgc: 8.081473e-08 ustep: 1
#> mgc: 4.414247e-13
#> outer mgc: 0.02941761
#> iter: 1 value: -333.2996 mgc: 0.01443108 ustep: 1
#> iter: 2 value: -333.2996 mgc: 8.104812e-08 ustep: 1
#> mgc: 4.82725e-13
#> outer mgc: 0.02961914
#> iter: 1 value: -333.3936 mgc: 0.02594613 ustep: 1
#> iter: 2 value: -333.3936 mgc: 2.434221e-07 ustep: 1
#> iter: 3 mgc: 4.396483e-13
#> outer mgc: 0.04180977
#> iter: 1 value: -333.199 mgc: 0.02589355 ustep: 1
#> iter: 2 value: -333.199 mgc: 2.431137e-07 ustep: 1
#> iter: 3 mgc: 3.765876e-13
#> outer mgc: 0.04165187
#> iter: 1 value: -333.2962 mgc: 7.811658e-05 ustep: 1
#> iter: 2 mgc: 5.983836e-11
#> iter: 1 mgc: 5.983836e-11
#> outer mgc: 2.406246e-09
#> iter: 1 value: -333.2967 mgc: 0.001706046 ustep: 1
#> iter: 2 value: -333.2967 mgc: 2.131766e-08 ustep: 1
#> mgc: 4.962697e-13
#> outer mgc: 0.07026863
#> iter: 1 value: -333.2957 mgc: 0.001706046 ustep: 1
#> iter: 2 value: -333.2957 mgc: 2.131643e-08 ustep: 1
#> mgc: 3.671508e-13
#> outer mgc: 0.07026861
#> iter: 1 value: -333.2964 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2964 mgc: 1.521444e-07 ustep: 1
#> iter: 3 mgc: 4.600764e-13
#> outer mgc: 0.1617967
#> iter: 1 value: -333.296 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.296 mgc: 1.52129e-07 ustep: 1
#> iter: 3 mgc: 4.596323e-13
#> outer mgc: 0.1617963
#> iter: 1 value: -333.2967 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2967 mgc: 1.229144e-06 ustep: 1
#> iter: 3 mgc: 4.192202e-13
#> outer mgc: 0.2197355
#> iter: 1 value: -333.2956 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2956 mgc: 1.228504e-06 ustep: 1
#> iter: 3 mgc: 4.085621e-13
#> outer mgc: 0.2197349
#> iter: 1 value: -333.2967 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2967 mgc: 1.178569e-06 ustep: 1
#> iter: 3 mgc: 4.476419e-13
#> outer mgc: 0.2215138
#> iter: 1 value: -333.2956 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2956 mgc: 1.177962e-06 ustep: 1
#> iter: 3 mgc: 4.697354e-13
#> outer mgc: 0.2215132
#> iter: 1 value: -333.2967 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2967 mgc: 1.331864e-06 ustep: 1
#> iter: 3 mgc: 4.201084e-13
#> outer mgc: 0.2211409
#> iter: 1 value: -333.2955 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2955 mgc: 1.33113e-06 ustep: 1
#> iter: 3 mgc: 4.236611e-13
#> outer mgc: 0.2211402
#> iter: 1 value: -333.2967 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2967 mgc: 1.447861e-06 ustep: 1
#> iter: 3 mgc: 4.058975e-13
#> outer mgc: 0.2183867
#> iter: 1 value: -333.2956 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2956 mgc: 1.446995e-06 ustep: 1
#> iter: 3 mgc: 3.53495e-13
#> outer mgc: 0.2183857
#> iter: 1 value: -333.2965 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2965 mgc: 1.800152e-06 ustep: 1
#> iter: 3 mgc: 4.933831e-13
#> outer mgc: 0.2103066
#> iter: 1 value: -333.2957 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2957 mgc: 1.798869e-06 ustep: 1
#> iter: 3 mgc: 3.813616e-13
#> outer mgc: 0.2103056
#> iter: 1 value: -333.2965 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2965 mgc: 1.184155e-05 ustep: 1
#> iter: 3 mgc: 2.149392e-12
#> outer mgc: 0.1473114
#> iter: 1 value: -333.2959 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2959 mgc: 1.182159e-05 ustep: 1
#> iter: 3 mgc: 2.151168e-12
#> outer mgc: 0.1473108
#> iter: 1 value: -333.2966 mgc: 0.01062593 ustep: 1
#> iter: 2 value: -333.2966 mgc: 3.380448e-07 ustep: 1
#> iter: 3 mgc: 5.608847e-13
#> outer mgc: 0.1766477
#> iter: 1 value: -333.2957 mgc: 0.01062593 ustep: 1
#> iter: 2 value: -333.2957 mgc: 3.379817e-07 ustep: 1
#> iter: 3 mgc: 4.813927e-13
#> outer mgc: 0.1766469
#> iter: 1 value: -333.2965 mgc: 0.007472208 ustep: 1
#> iter: 2 value: -333.2965 mgc: 2.31787e-07 ustep: 1
#> iter: 3 mgc: 5.400125e-13
#> outer mgc: 0.1186051
#> iter: 1 value: -333.2959 mgc: 0.007472208 ustep: 1
#> iter: 2 value: -333.2959 mgc: 2.317466e-07 ustep: 1
#> iter: 3 mgc: 4.192202e-13
#> outer mgc: 0.1186042
#> iter: 1 value: -333.1928 mgc: 0.02626238 ustep: 1
#> iter: 2 value: -333.1928 mgc: 5.14503e-07 ustep: 1
#> iter: 3 mgc: 3.828049e-13
#> outer mgc: 0.0786511
#> iter: 1 value: -333.3996 mgc: 0.02631496 ustep: 1
#> iter: 2 value: -333.3996 mgc: 5.154168e-07 ustep: 1
#> iter: 3 mgc: 4.178879e-13
#> outer mgc: 0.0786745
#> iter: 1 value: -333.2936 mgc: 0.01464723 ustep: 1
#> iter: 2 value: -333.2936 mgc: 2.077278e-08 ustep: 1
#> mgc: 4.840572e-13
#> outer mgc: 0.02857573
#> iter: 1 value: -333.2989 mgc: 0.01467656 ustep: 1
#> iter: 2 value: -333.2989 mgc: 2.083004e-08 ustep: 1
#> mgc: 3.863576e-13
#> outer mgc: 0.02859926
#> iter: 1 value: -333.1836 mgc: 0.02122111 ustep: 1
#> iter: 2 value: -333.1836 mgc: 1.554973e-06 ustep: 1
#> iter: 3 mgc: 4.072298e-13
#> outer mgc: 0.05249249
#> iter: 1 value: -333.4089 mgc: 0.02126359 ustep: 1
#> iter: 2 value: -333.4089 mgc: 1.558633e-06 ustep: 1
#> iter: 3 mgc: 3.73479e-13
#> outer mgc: 0.05241714
#> iter: 1 value: -333.1873 mgc: 0.02569484 ustep: 1
#> iter: 2 value: -333.1873 mgc: 1.144451e-06 ustep: 1
#> iter: 3 mgc: 3.941292e-13
#> outer mgc: 0.0823399
#> iter: 1 value: -333.4052 mgc: 0.02574628 ustep: 1
#> iter: 2 value: -333.4052 mgc: 1.145377e-06 ustep: 1
#> iter: 3 mgc: 4.116707e-13
#> outer mgc: 0.08225642
#> iter: 1 value: -333.2959 mgc: 0.002474418 ustep: 1
#> iter: 2 value: -333.2959 mgc: 3.148145e-08 ustep: 1
#> mgc: 4.451994e-13
#> outer mgc: 0.03765903
#> iter: 1 value: -333.2966 mgc: 0.002479371 ustep: 1
#> iter: 2 value: -333.2966 mgc: 3.160319e-08 ustep: 1
#> mgc: 4.485301e-13
#> outer mgc: 0.03772561
#> iter: 1 value: -333.2596 mgc: 0.01703843 ustep: 1
#> iter: 2 value: -333.2596 mgc: 2.323146e-06 ustep: 1
#> iter: 3 mgc: 4.112266e-13
#> outer mgc: 0.1710082
#> iter: 1 value: -333.3327 mgc: 0.01707255 ustep: 1
#> iter: 2 value: -333.3327 mgc: 2.326515e-06 ustep: 1
#> iter: 3 mgc: 4.576339e-13
#> outer mgc: 0.1711201
#> iter: 1 value: -333.2931 mgc: 0.01235191 ustep: 1
#> iter: 2 value: -333.2931 mgc: 1.660574e-08 ustep: 1
#> mgc: 3.472778e-13
#> outer mgc: 0.03180161
#> iter: 1 value: -333.2994 mgc: 0.01237664 ustep: 1
#> iter: 2 value: -333.2994 mgc: 1.665824e-08 ustep: 1
#> mgc: 5.071499e-13
#> outer mgc: 0.03183304
#> iter: 1 value: -333.2929 mgc: 0.01440214 ustep: 1
#> iter: 2 value: -333.2929 mgc: 8.081477e-08 ustep: 1
#> mgc: 4.805045e-13
#> outer mgc: 0.02950621
#> iter: 1 value: -333.2995 mgc: 0.01443098 ustep: 1
#> iter: 2 value: -333.2995 mgc: 8.104804e-08 ustep: 1
#> mgc: 4.178879e-13
#> outer mgc: 0.02953045
#> iter: 1 value: -333.3935 mgc: 0.02594609 ustep: 1
#> iter: 2 value: -333.3935 mgc: 2.434224e-07 ustep: 1
#> iter: 3 mgc: 4.449774e-13
#> outer mgc: 0.0417333
#> iter: 1 value: -333.199 mgc: 0.02589351 ustep: 1
#> iter: 2 value: -333.199 mgc: 2.431139e-07 ustep: 1
#> iter: 3 mgc: 3.641532e-13
#> outer mgc: 0.04172827
#> iter: 1 mgc: 5.090741e-10
#> iter: 1 mgc: 5.090741e-10
#> outer mgc: 2.803283e-11
#> iter: 1 value: -333.2967 mgc: 0.001706046 ustep: 1
#> iter: 2 value: -333.2967 mgc: 2.131772e-08 ustep: 1
#> mgc: 5.207501e-13
#> outer mgc: 0.07026863
#> iter: 1 value: -333.2957 mgc: 0.001706046 ustep: 1
#> iter: 2 value: -333.2957 mgc: 2.131642e-08 ustep: 1
#> mgc: 5.604406e-13
#> outer mgc: 0.07026861
#> iter: 1 value: -333.2964 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2964 mgc: 1.521443e-07 ustep: 1
#> iter: 3 mgc: 4.28324e-13
#> outer mgc: 0.1617967
#> iter: 1 value: -333.296 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.296 mgc: 1.52129e-07 ustep: 1
#> iter: 3 mgc: 4.040102e-13
#> outer mgc: 0.1617963
#> iter: 1 value: -333.2967 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2967 mgc: 1.229144e-06 ustep: 1
#> iter: 3 mgc: 3.903544e-13
#> outer mgc: 0.2197355
#> iter: 1 value: -333.2956 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2956 mgc: 1.228504e-06 ustep: 1
#> iter: 3 mgc: 5.160317e-13
#> outer mgc: 0.2197349
#> iter: 1 value: -333.2967 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2967 mgc: 1.178569e-06 ustep: 1
#> iter: 3 mgc: 5.12923e-13
#> outer mgc: 0.2215138
#> iter: 1 value: -333.2956 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2956 mgc: 1.177962e-06 ustep: 1
#> iter: 3 mgc: 4.951595e-13
#> outer mgc: 0.2215132
#> iter: 1 value: -333.2967 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2967 mgc: 1.331864e-06 ustep: 1
#> iter: 3 mgc: 5.133671e-13
#> outer mgc: 0.2211409
#> iter: 1 value: -333.2955 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2955 mgc: 1.33113e-06 ustep: 1
#> iter: 3 mgc: 3.717027e-13
#> outer mgc: 0.2211402
#> iter: 1 value: -333.2967 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2967 mgc: 1.447861e-06 ustep: 1
#> iter: 3 mgc: 5.391243e-13
#> outer mgc: 0.2183867
#> iter: 1 value: -333.2956 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2956 mgc: 1.446995e-06 ustep: 1
#> iter: 3 mgc: 4.094503e-13
#> outer mgc: 0.2183857
#> iter: 1 value: -333.2965 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2965 mgc: 1.800152e-06 ustep: 1
#> iter: 3 mgc: 4.116707e-13
#> outer mgc: 0.2103066
#> iter: 1 value: -333.2957 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2957 mgc: 1.798869e-06 ustep: 1
#> iter: 3 mgc: 3.508305e-13
#> outer mgc: 0.2103056
#> iter: 1 value: -333.2965 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2965 mgc: 1.184155e-05 ustep: 1
#> iter: 3 mgc: 2.184919e-12
#> outer mgc: 0.1473114
#> iter: 1 value: -333.2959 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2959 mgc: 1.182159e-05 ustep: 1
#> iter: 3 mgc: 2.122746e-12
#> outer mgc: 0.1473108
#> iter: 1 value: -333.2966 mgc: 0.01062593 ustep: 1
#> iter: 2 value: -333.2966 mgc: 3.380447e-07 ustep: 1
#> iter: 3 mgc: 4.565237e-13
#> outer mgc: 0.1766477
#> iter: 1 value: -333.2957 mgc: 0.01062593 ustep: 1
#> iter: 2 value: -333.2957 mgc: 3.379817e-07 ustep: 1
#> iter: 3 mgc: 4.662937e-13
#> outer mgc: 0.1766469
#> iter: 1 value: -333.2965 mgc: 0.007472208 ustep: 1
#> iter: 2 value: -333.2965 mgc: 2.317871e-07 ustep: 1
#> iter: 3 mgc: 4.498624e-13
#> outer mgc: 0.1186051
#> iter: 1 value: -333.2959 mgc: 0.007472208 ustep: 1
#> iter: 2 value: -333.2959 mgc: 2.317465e-07 ustep: 1
#> iter: 3 mgc: 4.13003e-13
#> outer mgc: 0.1186042
#> iter: 1 value: -333.1928 mgc: 0.02626238 ustep: 1
#> iter: 2 value: -333.1928 mgc: 5.14503e-07 ustep: 1
#> iter: 3 mgc: 3.579359e-13
#> outer mgc: 0.0786511
#> iter: 1 value: -333.3996 mgc: 0.02631496 ustep: 1
#> iter: 2 value: -333.3996 mgc: 5.154169e-07 ustep: 1
#> iter: 3 mgc: 3.372858e-13
#> outer mgc: 0.0786745
#> iter: 1 value: -333.2936 mgc: 0.01464723 ustep: 1
#> iter: 2 value: -333.2936 mgc: 2.077279e-08 ustep: 1
#> mgc: 3.979039e-13
#> outer mgc: 0.02857573
#> iter: 1 value: -333.2989 mgc: 0.01467656 ustep: 1
#> iter: 2 value: -333.2989 mgc: 2.083002e-08 ustep: 1
#> mgc: 4.464207e-13
#> outer mgc: 0.02859926
#> iter: 1 value: -333.1836 mgc: 0.02122111 ustep: 1
#> iter: 2 value: -333.1836 mgc: 1.554973e-06 ustep: 1
#> iter: 3 mgc: 5.124234e-13
#> outer mgc: 0.05249249
#> iter: 1 value: -333.4089 mgc: 0.02126359 ustep: 1
#> iter: 2 value: -333.4089 mgc: 1.558633e-06 ustep: 1
#> iter: 3 mgc: 3.654854e-13
#> outer mgc: 0.05241714
#> iter: 1 value: -333.1873 mgc: 0.02569484 ustep: 1
#> iter: 2 value: -333.1873 mgc: 1.144451e-06 ustep: 1
#> iter: 3 mgc: 4.152234e-13
#> outer mgc: 0.08233989
#> iter: 1 value: -333.4052 mgc: 0.02574628 ustep: 1
#> iter: 2 value: -333.4052 mgc: 1.145377e-06 ustep: 1
#> iter: 3 mgc: 4.245493e-13
#> outer mgc: 0.08225642
#> iter: 1 value: -333.2959 mgc: 0.002474418 ustep: 1
#> iter: 2 value: -333.2959 mgc: 3.148148e-08 ustep: 1
#> mgc: 5.273559e-13
#> outer mgc: 0.03765903
#> iter: 1 value: -333.2966 mgc: 0.002479371 ustep: 1
#> iter: 2 value: -333.2966 mgc: 3.160319e-08 ustep: 1
#> mgc: 5.764278e-13
#> outer mgc: 0.03772561
#> iter: 1 value: -333.2596 mgc: 0.01703843 ustep: 1
#> iter: 2 value: -333.2596 mgc: 2.323146e-06 ustep: 1
#> iter: 3 mgc: 4.654055e-13
#> outer mgc: 0.1710082
#> iter: 1 value: -333.3327 mgc: 0.01707254 ustep: 1
#> iter: 2 value: -333.3327 mgc: 2.326515e-06 ustep: 1
#> iter: 3 mgc: 4.52971e-13
#> outer mgc: 0.1711201
#> iter: 1 value: -333.2931 mgc: 0.01235191 ustep: 1
#> iter: 2 value: -333.2931 mgc: 1.660574e-08 ustep: 1
#> mgc: 3.366196e-13
#> outer mgc: 0.03180161
#> iter: 1 value: -333.2994 mgc: 0.01237664 ustep: 1
#> iter: 2 value: -333.2994 mgc: 1.665816e-08 ustep: 1
#> mgc: 4.760636e-13
#> outer mgc: 0.03183304
#> iter: 1 value: -333.2929 mgc: 0.01440214 ustep: 1
#> iter: 2 value: -333.2929 mgc: 8.081471e-08 ustep: 1
#> mgc: 5.018208e-13
#> outer mgc: 0.02950621
#> iter: 1 value: -333.2995 mgc: 0.01443098 ustep: 1
#> iter: 2 value: -333.2995 mgc: 8.104805e-08 ustep: 1
#> mgc: 5.084821e-13
#> outer mgc: 0.02953045
#> iter: 1 value: -333.3935 mgc: 0.02594609 ustep: 1
#> iter: 2 value: -333.3935 mgc: 2.434224e-07 ustep: 1
#> iter: 3 mgc: 4.618528e-13
#> outer mgc: 0.0417333
#> iter: 1 value: -333.199 mgc: 0.02589351 ustep: 1
#> iter: 2 value: -333.199 mgc: 2.431141e-07 ustep: 1
#> iter: 3 mgc: 5.133671e-13
#> outer mgc: 0.04172827
#> iter: 1 mgc: 5.177463e-10
#> iter: 1 mgc: 5.177463e-10
#> outer mgc: 1.008388e-12
#> iter: 1 value: -333.2967 mgc: 0.001706046 ustep: 1
#> iter: 2 value: -333.2967 mgc: 2.131773e-08 ustep: 1
#> mgc: 4.662937e-13
#> outer mgc: 0.07026863
#> iter: 1 value: -333.2957 mgc: 0.001706046 ustep: 1
#> iter: 2 value: -333.2957 mgc: 2.131641e-08 ustep: 1
#> mgc: 4.618528e-13
#> outer mgc: 0.07026861
#> iter: 1 value: -333.2964 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2964 mgc: 1.521443e-07 ustep: 1
#> iter: 3 mgc: 5.105361e-13
#> outer mgc: 0.1617967
#> iter: 1 value: -333.296 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.296 mgc: 1.52129e-07 ustep: 1
#> iter: 3 mgc: 3.708145e-13
#> outer mgc: 0.1617963
#> iter: 1 value: -333.2967 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2967 mgc: 1.229144e-06 ustep: 1
#> iter: 3 mgc: 4.209966e-13
#> outer mgc: 0.2197355
#> iter: 1 value: -333.2956 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2956 mgc: 1.228504e-06 ustep: 1
#> iter: 3 mgc: 3.739231e-13
#> outer mgc: 0.2197349
#> iter: 1 value: -333.2967 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2967 mgc: 1.178569e-06 ustep: 1
#> iter: 3 mgc: 3.93019e-13
#> outer mgc: 0.2215138
#> iter: 1 value: -333.2956 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2956 mgc: 1.177962e-06 ustep: 1
#> iter: 3 mgc: 4.634071e-13
#> outer mgc: 0.2215132
#> iter: 1 value: -333.2967 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2967 mgc: 1.331864e-06 ustep: 1
#> iter: 3 mgc: 4.121148e-13
#> outer mgc: 0.2211409
#> iter: 1 value: -333.2955 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2955 mgc: 1.33113e-06 ustep: 1
#> iter: 3 mgc: 3.552714e-13
#> outer mgc: 0.2211402
#> iter: 1 value: -333.2967 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2967 mgc: 1.447861e-06 ustep: 1
#> iter: 3 mgc: 4.08118e-13
#> outer mgc: 0.2183867
#> iter: 1 value: -333.2956 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2956 mgc: 1.446995e-06 ustep: 1
#> iter: 3 mgc: 5.506706e-13
#> outer mgc: 0.2183857
#> iter: 1 value: -333.2965 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2965 mgc: 1.800152e-06 ustep: 1
#> iter: 3 mgc: 4.142242e-13
#> outer mgc: 0.2103066
#> iter: 1 value: -333.2957 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2957 mgc: 1.798869e-06 ustep: 1
#> iter: 3 mgc: 4.567458e-13
#> outer mgc: 0.2103056
#> iter: 1 value: -333.2965 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2965 mgc: 1.184155e-05 ustep: 1
#> iter: 3 mgc: 2.17959e-12
#> outer mgc: 0.1473114
#> iter: 1 value: -333.2959 mgc: 0.01473375 ustep: 1
#> iter: 2 value: -333.2959 mgc: 1.182159e-05 ustep: 1
#> iter: 3 mgc: 2.067679e-12
#> outer mgc: 0.1473108
#> iter: 1 value: -333.2966 mgc: 0.01062593 ustep: 1
#> iter: 2 value: -333.2966 mgc: 3.380447e-07 ustep: 1
#> iter: 3 mgc: 3.739231e-13
#> outer mgc: 0.1766477
#> iter: 1 value: -333.2957 mgc: 0.01062593 ustep: 1
#> iter: 2 value: -333.2957 mgc: 3.379818e-07 ustep: 1
#> iter: 3 mgc: 3.810285e-13
#> outer mgc: 0.1766469
#> iter: 1 value: -333.2965 mgc: 0.007472208 ustep: 1
#> iter: 2 value: -333.2965 mgc: 2.317871e-07 ustep: 1
#> iter: 3 mgc: 4.005685e-13
#> outer mgc: 0.1186051
#> iter: 1 value: -333.2959 mgc: 0.007472208 ustep: 1
#> iter: 2 value: -333.2959 mgc: 2.317466e-07 ustep: 1
#> iter: 3 mgc: 3.206324e-13
#> outer mgc: 0.1186042
#> iter: 1 value: -333.1928 mgc: 0.02626238 ustep: 1
#> iter: 2 value: -333.1928 mgc: 5.145029e-07 ustep: 1
#> iter: 3 mgc: 4.840572e-13
#> outer mgc: 0.0786511
#> iter: 1 value: -333.3996 mgc: 0.02631496 ustep: 1
#> iter: 2 value: -333.3996 mgc: 5.154171e-07 ustep: 1
#> iter: 3 mgc: 4.636291e-13
#> outer mgc: 0.0786745
#> iter: 1 value: -333.2936 mgc: 0.01464723 ustep: 1
#> iter: 2 value: -333.2936 mgc: 2.077281e-08 ustep: 1
#> mgc: 5.684342e-13
#> outer mgc: 0.02857573
#> iter: 1 value: -333.2989 mgc: 0.01467656 ustep: 1
#> iter: 2 value: -333.2989 mgc: 2.083003e-08 ustep: 1
#> mgc: 4.849454e-13
#> outer mgc: 0.02859926
#> iter: 1 value: -333.1836 mgc: 0.02122111 ustep: 1
#> iter: 2 value: -333.1836 mgc: 1.554973e-06 ustep: 1
#> iter: 3 mgc: 4.356515e-13
#> outer mgc: 0.05249249
#> iter: 1 value: -333.4089 mgc: 0.02126359 ustep: 1
#> iter: 2 value: -333.4089 mgc: 1.558633e-06 ustep: 1
#> iter: 3 mgc: 4.086731e-13
#> outer mgc: 0.05241714
#> iter: 1 value: -333.1873 mgc: 0.02569484 ustep: 1
#> iter: 2 value: -333.1873 mgc: 1.144451e-06 ustep: 1
#> iter: 3 mgc: 3.925749e-13
#> outer mgc: 0.08233989
#> iter: 1 value: -333.4052 mgc: 0.02574628 ustep: 1
#> iter: 2 value: -333.4052 mgc: 1.145377e-06 ustep: 1
#> iter: 3 mgc: 4.074519e-13
#> outer mgc: 0.08225642
#> iter: 1 value: -333.2959 mgc: 0.002474418 ustep: 1
#> iter: 2 value: -333.2959 mgc: 3.148148e-08 ustep: 1
#> mgc: 4.813927e-13
#> outer mgc: 0.03765903
#> iter: 1 value: -333.2966 mgc: 0.002479371 ustep: 1
#> iter: 2 value: -333.2966 mgc: 3.160312e-08 ustep: 1
#> mgc: 3.63487e-13
#> outer mgc: 0.03772561
#> iter: 1 value: -333.2596 mgc: 0.01703843 ustep: 1
#> iter: 2 value: -333.2596 mgc: 2.323146e-06 ustep: 1
#> iter: 3 mgc: 3.810285e-13
#> outer mgc: 0.1710082
#> iter: 1 value: -333.3327 mgc: 0.01707254 ustep: 1
#> iter: 2 value: -333.3327 mgc: 2.326515e-06 ustep: 1
#> iter: 3 mgc: 4.063416e-13
#> outer mgc: 0.1711201
#> iter: 1 value: -333.2931 mgc: 0.01235191 ustep: 1
#> iter: 2 value: -333.2931 mgc: 1.660586e-08 ustep: 1
#> mgc: 4.266587e-13
#> outer mgc: 0.03180161
#> iter: 1 value: -333.2994 mgc: 0.01237664 ustep: 1
#> iter: 2 value: -333.2994 mgc: 1.665817e-08 ustep: 1
#> mgc: 4.640732e-13
#> outer mgc: 0.03183304
#> iter: 1 value: -333.2929 mgc: 0.01440214 ustep: 1
#> iter: 2 value: -333.2929 mgc: 8.081478e-08 ustep: 1
#> mgc: 5.795364e-13
#> outer mgc: 0.02950621
#> iter: 1 value: -333.2995 mgc: 0.01443098 ustep: 1
#> iter: 2 value: -333.2995 mgc: 8.104804e-08 ustep: 1
#> mgc: 4.476419e-13
#> outer mgc: 0.02953045
#> iter: 1 value: -333.3935 mgc: 0.02594609 ustep: 1
#> iter: 2 value: -333.3935 mgc: 2.434224e-07 ustep: 1
#> iter: 3 mgc: 5.37792e-13
#> outer mgc: 0.0417333
#> iter: 1 value: -333.199 mgc: 0.02589351 ustep: 1
#> iter: 2 value: -333.199 mgc: 2.431141e-07 ustep: 1
#> iter: 3 mgc: 4.276579e-13
#> outer mgc: 0.04172827
#> outer mgc: 1
#Now does it look?
summary(newStock)#iter = 100
#> An object of class "FLStock"
#>
#> Name: North Sea Herring
#> Description: Imported from a VPA file. ( ./data/index.txt ). Fri May 4 12:16:13 2012 [...]
#> Quant: age
#> Dims: age year unit season area iter
#> 9 21 1 1 1 100
#>
#> Range: min max pgroup minyear maxyear minfbar maxfbar
#> 0 8 8 1992 2012 2 6
#>
#> Metrics:
#> rec: 2.1e+07 - 8.1e+07 (NA)
#> ssb: NA - NA (NA)
#> catch: NA - NA (tonnes)
#> fbar: 0.084 - 0.565 (f)
plot(newStock) #Voila!
