diff --git a/unit_test/test_nse_net/ci-benchmarks/nse_net_unit_test.out b/unit_test/test_nse_net/ci-benchmarks/nse_net_unit_test.out index 18e17a9981..2af0ee8961 100644 --- a/unit_test/test_nse_net/ci-benchmarks/nse_net_unit_test.out +++ b/unit_test/test_nse_net/ci-benchmarks/nse_net_unit_test.out @@ -1,5 +1,5 @@ -Initializing AMReX (23.07-395-ge6c93bf22695)... -AMReX (23.07-395-ge6c93bf22695) initialized +Initializing AMReX (23.07-445-gc09c45c31007)... +AMReX (23.07-445-gc09c45c31007) initialized starting the single zone burn... reading in network electron-capture / beta-decay tables... chemical potential of proton is -3 @@ -8,28 +8,28 @@ State Density (g/cm^3): 10000000 State Temperature (K): 6000000000 electron fraction is 0.5 NSE state: -n : 0.000607569351 -H1 : 0.001992432526 -He4 : 0.5190078752 -C12 : 1.326101988e-05 -N13 : 9.477376579e-11 -N14 : 2.56548683e-09 -O16 : 3.225455214e-05 -F18 : 5.422064122e-11 -Ne20 : 7.427509352e-07 -Ne21 : 4.99891257e-08 -Na22 : 2.608008157e-09 -Na23 : 1.008842941e-06 -Mg24 : 0.0001022222498 -Al27 : 0.0005548826564 -Si28 : 0.03680966808 -P31 : 0.0422922376 -S32 : 0.0386577195 -Ar36 : 0.02464742425 -Ca40 : 0.0288515261 -Ti44 : 0.001661231048 -Cr48 : 0.009967017622 -Fe52 : 0.05968005276 -Ni56 : 0.2351208186 +n : 0.0006075779147 +H1 : 0.001992190328 +He4 : 0.519153618 +C12 : 1.326242057e-05 +N13 : 9.476378074e-11 +N14 : 2.565252696e-09 +O16 : 3.225118799e-05 +F18 : 5.421427915e-11 +Ne20 : 7.426541003e-07 +Ne21 : 4.998331295e-08 +Na22 : 2.607634144e-09 +Na23 : 1.008712481e-06 +Mg24 : 0.0001022062575 +Al27 : 0.0005547964335 +Si28 : 0.03680294966 +P31 : 0.0422845632 +S32 : 0.03864965592 +Ar36 : 0.02464164049 +Ca40 : 0.02884400365 +Ti44 : 0.001660754609 +Cr48 : 0.009963899266 +Fe52 : 0.05965982508 +Ni56 : 0.2350349989 We're in NSE. -AMReX (23.07-395-ge6c93bf22695) finalized +AMReX (23.07-445-gc09c45c31007) finalized diff --git a/unit_test/test_nse_net/make_table/burn_cell.H b/unit_test/test_nse_net/make_table/burn_cell.H index 886338bf59..17f056ef06 100644 --- a/unit_test/test_nse_net/make_table/burn_cell.H +++ b/unit_test/test_nse_net/make_table/burn_cell.H @@ -37,9 +37,13 @@ void burn_cell_c() state.rho = rho; state.y_e = Ye; - if (state.y_e > 0.52_rt){ - state.mu_p = -4.0_rt; - state.mu_n = -14.0_rt; + if (state.y_e >= 0.63_rt) { + state.mu_p = -3.0_rt; + state.mu_n = -14.9_rt; + } + else if (state.y_e > 0.52_rt){ + state.mu_p = -6.0_rt; + state.mu_n = -11.0_rt; } else if (state.y_e > 0.48_rt){ state.mu_p = -8.0_rt; diff --git a/unit_test/test_nse_net/make_table/ci-benchmarks/nse_net_make_table_unit_test.out b/unit_test/test_nse_net/make_table/ci-benchmarks/nse_net_make_table_unit_test.out index 7ec0293319..b6c5f387b3 100644 --- a/unit_test/test_nse_net/make_table/ci-benchmarks/nse_net_make_table_unit_test.out +++ b/unit_test/test_nse_net/make_table/ci-benchmarks/nse_net_make_table_unit_test.out @@ -1,355 +1,355 @@ -Initializing AMReX (23.07-395-ge6c93bf22695)... -AMReX (23.07-395-ge6c93bf22695) initialized +Initializing AMReX (23.07-445-gc09c45c31007)... +AMReX (23.07-445-gc09c45c31007) initialized starting the single zone burn... reading in network electron-capture / beta-decay tables... - 1.0000000000e+06 4.0000000000e+09 0.4000000000 -12.9555025204 -4.7372062842 0 - 1.0000000000e+06 4.4286927156e+09 0.4000000000 -12.4096892935 -5.3017297034 1 - 1.0000000000e+06 4.9033297922e+09 0.4000000000 -11.8124863350 -5.9247157509 1 - 1.0000000000e+06 5.4288352332e+09 0.4000000000 -11.2186398049 -6.5894751704 1 - 1.0000000000e+06 6.0106607628e+09 0.4000000000 -10.8221627699 -7.3603511111 1 - 1.0000000000e+06 6.6548423840e+09 0.4000000000 -10.4363763722 -8.2252373664 1 - 1.0000000000e+06 7.3680629973e+09 0.4000000000 -10.1445056150 -9.0849662772 1 - 1.0000000000e+06 8.1577217310e+09 0.4000000000 -10.2582950656 -9.7983309243 1 - 1.0000000000e+06 9.0320107014e+09 0.4000000000 -11.0251962339 -10.7011623998 1 - 1.0000000000e+06 1.0000000000e+10 0.4000000000 -12.3061952746 -11.9583139126 1 - 1.0000000000e+07 4.0000000000e+09 0.4000000000 -13.7539410035 -3.9436985166 0 - 1.0000000000e+07 4.4286927156e+09 0.4000000000 -13.2910063988 -4.4245064973 1 - 1.0000000000e+07 4.9033297922e+09 0.4000000000 -12.7738571048 -4.9622430425 1 - 1.0000000000e+07 5.4288352332e+09 0.4000000000 -12.2022519657 -5.5597835137 1 - 1.0000000000e+07 6.0106607628e+09 0.4000000000 -11.5968451667 -6.2091582807 1 - 1.0000000000e+07 6.6548423840e+09 0.4000000000 -11.0838653206 -6.9172956472 1 - 1.0000000000e+07 7.3680629973e+09 0.4000000000 -10.7117195040 -7.7493017684 1 - 1.0000000000e+07 8.1577217310e+09 0.4000000000 -10.3608224687 -8.6311809609 1 - 1.0000000000e+07 9.0320107014e+09 0.4000000000 -10.2512028838 -9.4182840030 1 - 1.0000000000e+07 1.0000000000e+10 0.4000000000 -10.6008153728 -10.1495273727 1 - 1.0000000000e+08 4.0000000000e+09 0.4000000000 -14.5607973487 -3.1500296846 0 - 1.0000000000e+08 4.4286927156e+09 0.4000000000 -14.1824332333 -3.5459083093 1 - 1.0000000000e+08 4.9033297922e+09 0.4000000000 -13.7577662180 -3.9903922193 1 - 1.0000000000e+08 5.4288352332e+09 0.4000000000 -13.2819583494 -4.4889815352 1 - 1.0000000000e+08 6.0106607628e+09 0.4000000000 -12.7514980414 -5.0464198084 1 - 1.0000000000e+08 6.6548423840e+09 0.4000000000 -12.1712468695 -5.6630196541 1 - 1.0000000000e+08 7.3680629973e+09 0.4000000000 -11.5722427319 -6.3289690263 1 - 1.0000000000e+08 8.1577217310e+09 0.4000000000 -11.0954919155 -7.0732199085 1 - 1.0000000000e+08 9.0320107014e+09 0.4000000000 -10.7313790718 -7.9307043937 1 - 1.0000000000e+08 1.0000000000e+10 0.4000000000 -10.4459234786 -8.7985272846 1 + 1.0000000000e+06 4.0000000000e+09 0.4000000000 -12.9555025985 -4.7372062322 0 + 1.0000000000e+06 4.4286927156e+09 0.4000000000 -12.4096886951 -5.3017300984 1 + 1.0000000000e+06 4.9033297922e+09 0.4000000000 -11.8124827277 -5.9247179903 1 + 1.0000000000e+06 5.4288352332e+09 0.4000000000 -11.2186195287 -6.5894809028 1 + 1.0000000000e+06 6.0106607628e+09 0.4000000000 -10.8221238417 -7.3603510805 1 + 1.0000000000e+06 6.6548423840e+09 0.4000000000 -10.4363399360 -8.2252368745 1 + 1.0000000000e+06 7.3680629973e+09 0.4000000000 -10.1444749829 -9.0849628652 1 + 1.0000000000e+06 8.1577217310e+09 0.4000000000 -10.2582729101 -9.7983255913 1 + 1.0000000000e+06 9.0320107014e+09 0.4000000000 -11.0251835300 -10.7011614674 1 + 1.0000000000e+06 1.0000000000e+10 0.4000000000 -12.3061845155 -11.9583138919 1 + 1.0000000000e+07 4.0000000000e+09 0.4000000000 -13.7539410130 -3.9436985103 0 + 1.0000000000e+07 4.4286927156e+09 0.4000000000 -13.2910064791 -4.4245064438 1 + 1.0000000000e+07 4.9033297922e+09 0.4000000000 -12.7738576654 -4.9622426711 1 + 1.0000000000e+07 5.4288352332e+09 0.4000000000 -12.2022552335 -5.5597814123 1 + 1.0000000000e+07 6.0106607628e+09 0.4000000000 -11.5968613147 -6.2091497038 1 + 1.0000000000e+07 6.6548423840e+09 0.4000000000 -11.0839445473 -6.9172943556 1 + 1.0000000000e+07 7.3680629973e+09 0.4000000000 -10.7118071972 -7.7493029627 1 + 1.0000000000e+07 8.1577217310e+09 0.4000000000 -10.3609046402 -8.6311911606 1 + 1.0000000000e+07 9.0320107014e+09 0.4000000000 -10.2511246531 -9.4182694259 1 + 1.0000000000e+07 1.0000000000e+10 0.4000000000 -10.6007602543 -10.1495148201 1 + 1.0000000000e+08 4.0000000000e+09 0.4000000000 -14.5607973496 -3.1500296840 0 + 1.0000000000e+08 4.4286927156e+09 0.4000000000 -14.1824332412 -3.5459083040 1 + 1.0000000000e+08 4.9033297922e+09 0.4000000000 -13.7577662756 -3.9903921809 1 + 1.0000000000e+08 5.4288352332e+09 0.4000000000 -13.2819587028 -4.4889813003 1 + 1.0000000000e+08 6.0106607628e+09 0.4000000000 -12.7514998816 -5.0464186008 1 + 1.0000000000e+08 6.6548423840e+09 0.4000000000 -12.1712550423 -5.6630145848 1 + 1.0000000000e+08 7.3680629973e+09 0.4000000000 -11.5722750148 -6.3289545463 1 + 1.0000000000e+08 8.1577217310e+09 0.4000000000 -11.0956062224 -7.0732197715 1 + 1.0000000000e+08 9.0320107014e+09 0.4000000000 -10.7314991991 -7.9307107472 1 + 1.0000000000e+08 1.0000000000e+10 0.4000000000 -10.4460221082 -8.7985576148 1 1.0000000000e+09 4.0000000000e+09 0.4000000000 -15.3871516246 -2.3563444183 0 - 1.0000000000e+09 4.4286927156e+09 0.4000000000 -15.0933118585 -2.6671723129 1 - 1.0000000000e+09 4.9033297922e+09 0.4000000000 -14.7621490253 -3.0175641915 1 - 1.0000000000e+09 5.4288352332e+09 0.4000000000 -14.3892566411 -3.4124048451 1 - 1.0000000000e+09 6.0106607628e+09 0.4000000000 -13.9696328335 -3.8570453479 1 - 1.0000000000e+09 6.6548423840e+09 0.4000000000 -13.4990892461 -4.3569663734 1 - 1.0000000000e+09 7.3680629973e+09 0.4000000000 -12.9751070035 -4.9165801598 1 - 1.0000000000e+09 8.1577217310e+09 0.4000000000 -12.4029615922 -5.5362921583 1 - 1.0000000000e+09 9.0320107014e+09 0.4000000000 -11.8075492548 -6.2097916688 1 - 1.0000000000e+09 1.0000000000e+10 0.4000000000 -11.2831770772 -6.9503077459 1 - 1.0000000000e+10 4.0000000000e+09 0.4000000000 -16.2582598176 -1.5626567041 0 - 1.0000000000e+10 4.4286927156e+09 0.4000000000 -16.0486738826 -1.7884208125 1 - 1.0000000000e+10 4.9033297922e+09 0.4000000000 -15.8108280448 -2.0446360334 1 - 1.0000000000e+10 5.4288352332e+09 0.4000000000 -15.5412395186 -2.3352407233 1 - 1.0000000000e+10 6.0106607628e+09 0.4000000000 -15.2357440992 -2.6646600961 1 - 1.0000000000e+10 6.6548423840e+09 0.4000000000 -14.8903553277 -3.0378300141 1 - 1.0000000000e+10 7.3680629973e+09 0.4000000000 -14.5003230778 -3.4601446829 1 - 1.0000000000e+10 8.1577217310e+09 0.4000000000 -14.0614570153 -3.9371817304 1 - 1.0000000000e+10 9.0320107014e+09 0.4000000000 -13.5711724217 -4.4739432697 1 - 1.0000000000e+10 1.0000000000e+10 0.4000000000 -13.0312175162 -5.0733655763 1 - 1.0000000000e+06 4.0000000000e+09 0.4500000000 -12.6199618053 -5.0493825987 0 - 1.0000000000e+06 4.4286927156e+09 0.4500000000 -12.0411581961 -5.6445822151 1 - 1.0000000000e+06 4.9033297922e+09 0.4500000000 -11.4232312199 -6.2868479094 1 - 1.0000000000e+06 5.4288352332e+09 0.4500000000 -10.8518428398 -6.9298657584 1 - 1.0000000000e+06 6.0106607628e+09 0.4500000000 -10.4348708730 -7.7172921982 1 - 1.0000000000e+06 6.6548423840e+09 0.4500000000 -10.0424621871 -8.5877955647 1 - 1.0000000000e+06 7.3680629973e+09 0.4500000000 -9.8849946892 -9.3214605955 1 - 1.0000000000e+06 8.1577217310e+09 0.4500000000 -10.1350747780 -9.9053651295 1 - 1.0000000000e+06 9.0320107014e+09 0.4500000000 -10.9317730375 -10.7720016856 1 - 1.0000000000e+06 1.0000000000e+10 0.4500000000 -12.2046629248 -12.0333675970 1 - 1.0000000000e+07 4.0000000000e+09 0.4500000000 -13.4182964409 -4.2562148091 0 - 1.0000000000e+07 4.4286927156e+09 0.4500000000 -12.9196141769 -4.7702516152 1 - 1.0000000000e+07 4.9033297922e+09 0.4500000000 -12.3646700494 -5.3431239685 1 - 1.0000000000e+07 5.4288352332e+09 0.4500000000 -11.7607919864 -5.9706982948 1 - 1.0000000000e+07 6.0106607628e+09 0.4500000000 -11.1529755129 -6.6222947524 1 - 1.0000000000e+07 6.6548423840e+09 0.4500000000 -10.6533363165 -7.3144541390 1 - 1.0000000000e+07 7.3680629973e+09 0.4500000000 -10.2538115130 -8.1710937198 1 - 1.0000000000e+07 8.1577217310e+09 0.4500000000 -9.9753897162 -8.9845986143 1 - 1.0000000000e+07 9.0320107014e+09 0.4500000000 -10.0361542054 -9.6115697958 1 - 1.0000000000e+07 1.0000000000e+10 0.4500000000 -10.4770937073 -10.2529776199 1 - 1.0000000000e+08 4.0000000000e+09 0.4500000000 -14.2257599760 -3.4625787945 0 - 1.0000000000e+08 4.4286927156e+09 0.4500000000 -13.8113628847 -3.8919431610 1 - 1.0000000000e+08 4.9033297922e+09 0.4500000000 -13.3469605437 -4.3733521393 1 - 1.0000000000e+08 5.4288352332e+09 0.4500000000 -12.8281601561 -4.9118963225 1 - 1.0000000000e+08 6.0106607628e+09 0.4500000000 -12.2549827771 -5.5090675693 1 - 1.0000000000e+08 6.6548423840e+09 0.4500000000 -11.6452553291 -6.1532183846 1 - 1.0000000000e+08 7.3680629973e+09 0.4500000000 -11.0519169747 -6.8138046108 1 - 1.0000000000e+08 8.1577217310e+09 0.4500000000 -10.5744119695 -7.5537001211 1 - 1.0000000000e+08 9.0320107014e+09 0.4500000000 -10.2097691952 -8.4106502242 1 - 1.0000000000e+08 1.0000000000e+10 0.4500000000 -10.0527955898 -9.1577337848 1 - 1.0000000000e+09 4.0000000000e+09 0.4500000000 -15.0536228396 -2.6688948472 0 - 1.0000000000e+09 4.4286927156e+09 0.4500000000 -14.7237134072 -3.0132309283 1 - 1.0000000000e+09 4.9033297922e+09 0.4500000000 -14.3526160867 -3.4007197563 1 - 1.0000000000e+09 5.4288352332e+09 0.4500000000 -13.9355877638 -3.8365577467 1 - 1.0000000000e+09 6.0106607628e+09 0.4500000000 -13.4675881563 -4.3261360138 1 - 1.0000000000e+09 6.6548423840e+09 0.4500000000 -12.9456617325 -4.8738185708 1 - 1.0000000000e+09 7.3680629973e+09 0.4500000000 -12.3724764785 -5.4792958152 1 - 1.0000000000e+09 8.1577217310e+09 0.4500000000 -11.7661588560 -6.1311892901 1 - 1.0000000000e+09 9.0320107014e+09 0.4500000000 -11.1668612132 -6.8085657276 1 - 1.0000000000e+09 1.0000000000e+10 0.4500000000 -10.6512855625 -7.5351949184 1 - 1.0000000000e+10 4.0000000000e+09 0.4500000000 -15.9281485305 -1.8752022445 0 - 1.0000000000e+10 4.4286927156e+09 0.4500000000 -15.6824856145 -2.1344670415 1 - 1.0000000000e+10 4.9033297922e+09 0.4500000000 -15.4047001813 -2.4277705489 1 - 1.0000000000e+10 5.4288352332e+09 0.4500000000 -15.0909188396 -2.7594143274 1 - 1.0000000000e+10 6.0106607628e+09 0.4500000000 -14.7366091103 -3.1341694060 1 - 1.0000000000e+10 6.6548423840e+09 0.4500000000 -14.3375648682 -3.5571890304 1 - 1.0000000000e+10 7.3680629973e+09 0.4500000000 -13.8893261475 -4.0336347952 1 - 1.0000000000e+10 8.1577217310e+09 0.4500000000 -13.3893312974 -4.5676563215 1 - 1.0000000000e+10 9.0320107014e+09 0.4500000000 -12.8393727587 -5.1603436508 1 - 1.0000000000e+10 1.0000000000e+10 0.4500000000 -12.2492192359 -5.8073345674 1 - 1.0000000000e+06 4.0000000000e+09 0.5000000000 -8.0664287063 -9.4890300885 0 - 1.0000000000e+06 4.4286927156e+09 0.5000000000 -8.2111123981 -9.3757297125 1 - 1.0000000000e+06 4.9033297922e+09 0.5000000000 -8.3918399895 -9.2344583793 1 - 1.0000000000e+06 5.4288352332e+09 0.5000000000 -8.7773708532 -8.9447667662 1 - 1.0000000000e+06 6.0106607628e+09 0.5000000000 -9.0637265807 -9.0647441208 1 - 1.0000000000e+06 6.6548423840e+09 0.5000000000 -9.3058731093 -9.3070059375 1 - 1.0000000000e+06 7.3680629973e+09 0.5000000000 -9.5983915599 -9.5996554213 1 - 1.0000000000e+06 8.1577217310e+09 0.5000000000 -10.0168541391 -10.0182626012 1 - 1.0000000000e+06 9.0320107014e+09 0.5000000000 -10.8474197344 -10.8489878418 1 - 1.0000000000e+06 1.0000000000e+10 0.5000000000 -12.1138199092 -12.1155643414 1 - 1.0000000000e+07 4.0000000000e+09 0.5000000000 -7.8839123533 -9.6536792188 0 - 1.0000000000e+07 4.4286927156e+09 0.5000000000 -8.0019180002 -9.5630334034 1 - 1.0000000000e+07 4.9033297922e+09 0.5000000000 -8.1354121681 -9.4609849615 1 - 1.0000000000e+07 5.4288352332e+09 0.5000000000 -8.2961986080 -9.3384887105 1 - 1.0000000000e+07 6.0106607628e+09 0.5000000000 -8.5443326081 -9.1484326209 1 - 1.0000000000e+07 6.6548423840e+09 0.5000000000 -8.9669523913 -8.9714790676 1 - 1.0000000000e+07 7.3680629973e+09 0.5000000000 -9.2004322394 -9.2015703723 1 - 1.0000000000e+07 8.1577217310e+09 0.5000000000 -9.4726257222 -9.4739144064 1 - 1.0000000000e+07 9.0320107014e+09 0.5000000000 -9.8194771269 -9.8209313702 1 - 1.0000000000e+07 1.0000000000e+10 0.5000000000 -10.3609081088 -10.3625443016 1 - 1.0000000000e+08 4.0000000000e+09 0.5000000000 -7.7564252605 -9.7786841281 0 - 1.0000000000e+08 4.4286927156e+09 0.5000000000 -7.8492019112 -9.7096717498 1 - 1.0000000000e+08 4.9033297922e+09 0.5000000000 -7.9527961073 -9.6330305105 1 - 1.0000000000e+08 5.4288352332e+09 0.5000000000 -8.0692464391 -9.5476093419 1 - 1.0000000000e+08 6.0106607628e+09 0.5000000000 -8.2063900317 -9.4473474639 1 - 1.0000000000e+08 6.6548423840e+09 0.5000000000 -8.3860213927 -9.3159021644 1 - 1.0000000000e+08 7.3680629973e+09 0.5000000000 -8.6892088477 -9.0970021780 1 - 1.0000000000e+08 8.1577217310e+09 0.5000000000 -9.0475207044 -9.0501554594 1 - 1.0000000000e+08 9.0320107014e+09 0.5000000000 -9.2986740372 -9.2997699354 1 - 1.0000000000e+08 1.0000000000e+10 0.5000000000 -9.5981877903 -9.5994826015 1 - 1.0000000000e+09 4.0000000000e+09 0.5000000000 -7.7470638015 -9.8201507833 0 - 1.0000000000e+09 4.4286927156e+09 0.5000000000 -7.8147769492 -9.7725305169 1 - 1.0000000000e+09 4.9033297922e+09 0.5000000000 -7.8899734067 -9.7200977838 1 - 1.0000000000e+09 5.4288352332e+09 0.5000000000 -7.9729788613 -9.6630648524 1 - 1.0000000000e+09 6.0106607628e+09 0.5000000000 -8.0659663127 -9.5999364916 1 - 1.0000000000e+09 6.6548423840e+09 0.5000000000 -8.1710732556 -9.5300539621 1 - 1.0000000000e+09 7.3680629973e+09 0.5000000000 -8.2980125161 -9.4462720608 1 - 1.0000000000e+09 8.1577217310e+09 0.5000000000 -8.4687660372 -9.3330282318 1 - 1.0000000000e+09 9.0320107014e+09 0.5000000000 -8.7320500861 -9.1630033338 1 - 1.0000000000e+09 1.0000000000e+10 0.5000000000 -9.0632943617 -9.0828318763 1 - 1.0000000000e+10 4.0000000000e+09 0.5000000000 -7.9879953396 -9.6879468798 1 - 1.0000000000e+10 4.4286927156e+09 0.5000000000 -8.0302314490 -9.6619337577 1 - 1.0000000000e+10 4.9033297922e+09 0.5000000000 -8.0763284304 -9.6342740546 1 - 1.0000000000e+10 5.4288352332e+09 0.5000000000 -8.1258459897 -9.6058175046 1 - 1.0000000000e+10 6.0106607628e+09 0.5000000000 -8.1795088599 -9.5762491105 1 - 1.0000000000e+10 6.6548423840e+09 0.5000000000 -8.2365078557 -9.5472461611 1 - 1.0000000000e+10 7.3680629973e+09 0.5000000000 -8.2988571338 -9.5177368349 1 - 1.0000000000e+10 8.1577217310e+09 0.5000000000 -8.3695493416 -9.4863947153 1 - 1.0000000000e+10 9.0320107014e+09 0.5000000000 -8.4559431120 -9.4487939586 1 - 1.0000000000e+10 1.0000000000e+10 0.5000000000 -8.5723779233 -9.3960835543 1 - 1.0000000000e+06 4.0000000000e+09 0.5500000000 -4.9263567583 -12.6306847771 0 - 1.0000000000e+06 4.4286927156e+09 0.5500000000 -5.5125982821 -12.0759967435 1 - 1.0000000000e+06 4.9033297922e+09 0.5500000000 -6.1679059616 -11.4605293767 1 - 1.0000000000e+06 5.4288352332e+09 0.5500000000 -6.9002510955 -10.8361899606 1 - 1.0000000000e+06 6.0106607628e+09 0.5500000000 -7.7162795920 -10.4358987387 1 - 1.0000000000e+06 6.6548423840e+09 0.5500000000 -8.5866665410 -10.0436084115 1 - 1.0000000000e+06 7.3680629973e+09 0.5500000000 -9.3202024829 -9.8862686067 1 - 1.0000000000e+06 8.1577217310e+09 0.5500000000 -9.9039624201 -10.1364903176 1 - 1.0000000000e+06 9.0320107014e+09 0.5500000000 -10.7704366546 -10.9333444696 1 - 1.0000000000e+06 1.0000000000e+10 0.5500000000 -12.0316256213 -12.2064099419 1 - 1.0000000000e+07 4.0000000000e+09 0.5500000000 -4.1328499283 -13.4067273031 0 - 1.0000000000e+07 4.4286927156e+09 0.5500000000 -4.6340211015 -12.9330486665 1 - 1.0000000000e+07 4.9033297922e+09 0.5500000000 -5.1951522570 -12.4035178568 1 - 1.0000000000e+07 5.4288352332e+09 0.5500000000 -5.8233286811 -11.8138168049 1 - 1.0000000000e+07 6.0106607628e+09 0.5500000000 -6.5263096328 -11.1696374790 1 - 1.0000000000e+07 6.6548423840e+09 0.5500000000 -7.3118205884 -10.6533164642 1 - 1.0000000000e+07 7.3680629973e+09 0.5500000000 -8.1699653523 -10.2549768923 1 - 1.0000000000e+07 8.1577217310e+09 0.5500000000 -8.9833218276 -9.9767011458 1 - 1.0000000000e+07 9.0320107014e+09 0.5500000000 -9.6101289590 -10.0376261258 1 - 1.0000000000e+07 1.0000000000e+10 0.5500000000 -10.2513560243 -10.4787469547 1 - 1.0000000000e+08 4.0000000000e+09 0.5500000000 -3.3394222565 -14.1986157657 0 - 1.0000000000e+08 4.4286927156e+09 0.5500000000 -3.7555058127 -13.8064185762 1 - 1.0000000000e+08 4.9033297922e+09 0.5500000000 -4.2224338348 -13.3665796932 1 - 1.0000000000e+08 5.4288352332e+09 0.5500000000 -4.7463244462 -12.8738625707 1 - 1.0000000000e+08 6.0106607628e+09 0.5500000000 -5.3340170370 -12.3231580847 1 - 1.0000000000e+08 6.6548423840e+09 0.5500000000 -5.9930630460 -11.7121630303 1 - 1.0000000000e+08 7.3680629973e+09 0.5500000000 -6.7314189359 -11.0586117888 1 - 1.0000000000e+08 8.1577217310e+09 0.5500000000 -7.5518998378 -10.5748765826 1 - 1.0000000000e+08 9.0320107014e+09 0.5500000000 -8.4095851780 -10.2109427835 1 - 1.0000000000e+08 1.0000000000e+10 0.5500000000 -9.1564767533 -10.0541540547 1 - 1.0000000000e+09 4.0000000000e+09 0.5500000000 -2.5470603946 -15.0251960133 0 - 1.0000000000e+09 4.4286927156e+09 0.5500000000 -2.8780137697 -14.7144476551 1 - 1.0000000000e+09 4.9033297922e+09 0.5500000000 -3.2506978413 -14.3646520162 1 - 1.0000000000e+09 5.4288352332e+09 0.5500000000 -3.6702538871 -13.9712053620 1 - 1.0000000000e+09 6.0106607628e+09 0.5500000000 -4.1424458459 -13.5290133725 1 - 1.0000000000e+09 6.6548423840e+09 0.5500000000 -4.6737326992 -13.0330537776 1 - 1.0000000000e+09 7.3680629973e+09 0.5500000000 -5.2713212052 -12.4785205588 1 - 1.0000000000e+09 8.1577217310e+09 0.5500000000 -5.9429385089 -11.8635582114 1 - 1.0000000000e+09 9.0320107014e+09 0.5500000000 -6.6949519374 -11.2026555889 1 - 1.0000000000e+09 1.0000000000e+10 0.5500000000 -7.5247904700 -10.6481506638 1 - 1.0000000000e+10 4.0000000000e+09 0.5500000000 -1.7575285290 -15.9281119586 1 - 1.0000000000e+10 4.4286927156e+09 0.5500000000 -2.0032273072 -15.6987367857 1 - 1.0000000000e+10 4.9033297922e+09 0.5500000000 -2.2815466531 -15.4389665570 1 - 1.0000000000e+10 5.4288352332e+09 0.5500000000 -2.5966508213 -15.1450469031 1 - 1.0000000000e+10 6.0106607628e+09 0.5500000000 -2.9532220861 -14.8127045193 1 - 1.0000000000e+10 6.6548423840e+09 0.5500000000 -3.3565239474 -14.4375393899 1 - 1.0000000000e+10 7.3680629973e+09 0.5500000000 -3.8124716671 -14.0145582613 1 - 1.0000000000e+10 8.1577217310e+09 0.5500000000 -4.3277071808 -13.5387148150 1 - 1.0000000000e+10 9.0320107014e+09 0.5500000000 -4.9096377922 -13.0053181043 1 - 1.0000000000e+10 1.0000000000e+10 0.5500000000 -5.5661097741 -12.4114337141 1 - 1.0000000000e+06 4.0000000000e+09 0.6000000000 -4.6874372498 -12.8713326831 0 - 1.0000000000e+06 4.4286927156e+09 0.6000000000 -5.2480726287 -12.3424870400 1 - 1.0000000000e+06 4.9033297922e+09 0.6000000000 -5.8750376557 -11.7561606661 1 - 1.0000000000e+06 5.4288352332e+09 0.6000000000 -6.5760898652 -11.1834318133 1 - 1.0000000000e+06 6.0106607628e+09 0.6000000000 -7.3593425290 -10.8232051193 1 - 1.0000000000e+06 6.6548423840e+09 0.6000000000 -8.2241108560 -10.4375374976 1 - 1.0000000000e+06 7.3680629973e+09 0.6000000000 -9.0837111804 -10.1457926302 1 - 1.0000000000e+06 8.1577217310e+09 0.6000000000 -9.7969327525 -10.2597190165 1 - 1.0000000000e+06 9.0320107014e+09 0.6000000000 -10.6996001970 -11.0267712343 1 - 1.0000000000e+06 1.0000000000e+10 0.6000000000 -11.9565742814 -12.3079450256 1 - 1.0000000000e+07 4.0000000000e+09 0.6000000000 -3.8939383372 -13.6477417081 0 - 1.0000000000e+07 4.4286927156e+09 0.6000000000 -4.3695027901 -13.1998278106 1 - 1.0000000000e+07 4.9033297922e+09 0.6000000000 -4.9022832206 -12.6988733383 1 - 1.0000000000e+07 5.4288352332e+09 0.6000000000 -5.4990786383 -12.1410754239 1 - 1.0000000000e+07 6.0106607628e+09 0.6000000000 -6.1674541533 -11.5347440371 1 - 1.0000000000e+07 6.6548423840e+09 0.6000000000 -6.9156298388 -11.0830753295 1 - 1.0000000000e+07 7.3680629973e+09 0.6000000000 -7.7481817497 -10.7129141315 1 - 1.0000000000e+07 8.1577217310e+09 0.6000000000 -8.6299116983 -10.3621615417 1 - 1.0000000000e+07 9.0320107014e+09 0.6000000000 -9.4168075327 -10.2525633805 1 - 1.0000000000e+07 1.0000000000e+10 0.6000000000 -10.1478709457 -10.6024085069 1 - 1.0000000000e+08 4.0000000000e+09 0.6000000000 -3.1005259273 -14.4405028547 0 - 1.0000000000e+08 4.4286927156e+09 0.6000000000 -3.4910018429 -14.0740545856 1 - 1.0000000000e+08 4.9033297922e+09 0.6000000000 -3.9295779936 -13.6627315310 1 - 1.0000000000e+08 5.4288352332e+09 0.6000000000 -4.4220785916 -13.2016151666 1 - 1.0000000000e+08 6.0106607628e+09 0.6000000000 -4.9750223848 -12.6860143020 1 - 1.0000000000e+08 6.6548423840e+09 0.6000000000 -5.5956833129 -12.1143859321 1 - 1.0000000000e+08 7.3680629973e+09 0.6000000000 -6.2923470974 -11.5100241165 1 - 1.0000000000e+08 8.1577217310e+09 0.6000000000 -7.0719798685 -11.0955902794 1 - 1.0000000000e+08 9.0320107014e+09 0.6000000000 -7.9296619962 -10.7326405101 1 - 1.0000000000e+08 1.0000000000e+10 0.6000000000 -8.7972934615 -10.4473618809 1 - 1.0000000000e+09 4.0000000000e+09 0.6000000000 -2.3082235668 -15.2690142179 0 - 1.0000000000e+09 4.4286927156e+09 0.6000000000 -2.6135665219 -14.9840071714 1 - 1.0000000000e+09 4.9033297922e+09 0.6000000000 -2.9578960346 -14.6627124682 1 - 1.0000000000e+09 5.4288352332e+09 0.6000000000 -3.3460592682 -14.3008241528 1 - 1.0000000000e+09 6.0106607628e+09 0.6000000000 -3.7834946071 -13.8935836540 1 - 1.0000000000e+09 6.6548423840e+09 0.6000000000 -4.2763022554 -13.4363612027 1 - 1.0000000000e+09 7.3680629973e+09 0.6000000000 -4.8313168013 -12.9248402164 1 - 1.0000000000e+09 8.1577217310e+09 0.6000000000 -5.4561228019 -12.3577854824 1 - 1.0000000000e+09 9.0320107014e+09 0.6000000000 -6.1587039833 -11.7516690135 1 - 1.0000000000e+09 1.0000000000e+10 0.6000000000 -6.9454524603 -11.2745753097 1 - 1.0000000000e+10 4.0000000000e+09 0.6000000000 -1.5188716438 -16.1761351108 0 - 1.0000000000e+10 4.4286927156e+09 0.6000000000 -1.7389525113 -15.9724978714 1 - 1.0000000000e+10 4.9033297922e+09 0.6000000000 -1.9889099833 -15.7412229679 1 - 1.0000000000e+10 5.4288352332e+09 0.6000000000 -2.2726142074 -15.4788512292 1 - 1.0000000000e+10 6.0106607628e+09 0.6000000000 -2.5944218278 -15.1814362147 1 - 1.0000000000e+10 6.6548423840e+09 0.6000000000 -2.9592358252 -14.8449412444 1 - 1.0000000000e+10 7.3680629973e+09 0.6000000000 -3.3725725553 -14.4647794265 1 - 1.0000000000e+10 8.1577217310e+09 0.6000000000 -3.8406360372 -14.0363659318 1 - 1.0000000000e+10 9.0320107014e+09 0.6000000000 -4.3703910764 -13.5555230831 1 - 1.0000000000e+10 1.0000000000e+10 0.6000000000 -4.9695596343 -13.0196507108 1 - 1.0000000000e+06 4.0000000000e+09 0.6500000000 -4.5476796742 -13.0130039947 0 - 1.0000000000e+06 4.4286927156e+09 0.6500000000 -5.0933362712 -12.4994174293 1 - 1.0000000000e+06 4.9033297922e+09 0.6500000000 -5.7037185995 -11.9307011924 1 - 1.0000000000e+06 5.4288352332e+09 0.6500000000 -6.3864250027 -11.4022341176 1 - 1.0000000000e+06 6.0106607628e+09 0.6500000000 -7.1497090923 -11.0673819184 1 - 1.0000000000e+06 6.6548423840e+09 0.6500000000 -7.9986396847 -10.7004385555 1 - 1.0000000000e+06 7.3680629973e+09 0.6500000000 -8.8932281298 -10.3705502418 1 - 1.0000000000e+06 8.1577217310e+09 0.6500000000 -9.6962470719 -10.3881605331 1 - 1.0000000000e+06 9.0320107014e+09 0.6500000000 -10.6340099963 -11.1317007071 1 - 1.0000000000e+06 1.0000000000e+10 0.6500000000 -11.8875232376 -12.4230320085 1 - 1.0000000000e+07 4.0000000000e+09 0.6500000000 -3.7541883557 -13.7897569631 0 - 1.0000000000e+07 4.4286927156e+09 0.6500000000 -4.2147736231 -13.3570050888 1 - 1.0000000000e+07 4.9033297922e+09 0.6500000000 -4.7309696318 -12.8729015431 1 - 1.0000000000e+07 5.4288352332e+09 0.6500000000 -5.3094047388 -12.3341136606 1 - 1.0000000000e+07 6.0106607628e+09 0.6500000000 -5.9574765589 -11.7528091119 1 - 1.0000000000e+07 6.6548423840e+09 0.6500000000 -6.6833511744 -11.3535641629 1 - 1.0000000000e+07 7.3680629973e+09 0.6500000000 -7.4942271332 -11.0088492097 1 - 1.0000000000e+07 8.1577217310e+09 0.6500000000 -8.3793161375 -10.6554902076 1 - 1.0000000000e+07 9.0320107014e+09 0.6500000000 -9.2440396955 -10.4605434295 1 - 1.0000000000e+07 1.0000000000e+10 0.6500000000 -10.0502569485 -10.7352236557 1 - 1.0000000000e+08 4.0000000000e+09 0.6500000000 -2.9607906232 -14.5833440181 0 - 1.0000000000e+08 4.4286927156e+09 0.6500000000 -3.3362864793 -14.2320397904 1 - 1.0000000000e+08 4.9033297922e+09 0.6500000000 -3.7582773168 -13.8374930298 1 - 1.0000000000e+08 5.4288352332e+09 0.6500000000 -4.2324153362 -13.3949878626 1 - 1.0000000000e+08 6.0106607628e+09 0.6500000000 -4.7650296834 -12.9001483764 1 - 1.0000000000e+08 6.6548423840e+09 0.6500000000 -5.3631981018 -12.3522858919 1 - 1.0000000000e+08 7.3680629973e+09 0.6500000000 -6.0351054434 -11.7844816604 1 - 1.0000000000e+08 8.1577217310e+09 0.6500000000 -6.7883154852 -11.4259937475 1 - 1.0000000000e+08 9.0320107014e+09 0.6500000000 -7.6264687221 -11.0863271605 1 - 1.0000000000e+08 1.0000000000e+10 0.6500000000 -8.5199483965 -10.7735064153 1 - 1.0000000000e+09 4.0000000000e+09 0.6500000000 -2.1685452785 -15.4136828685 0 - 1.0000000000e+09 4.4286927156e+09 0.6500000000 -2.4589055197 -15.1438122255 1 - 1.0000000000e+09 4.9033297922e+09 0.6500000000 -2.7866471569 -14.8392773913 1 - 1.0000000000e+09 5.4288352332e+09 0.6500000000 -3.1564453048 -14.4959492528 1 - 1.0000000000e+09 6.0106607628e+09 0.6500000000 -3.5735478746 -14.1092731402 1 - 1.0000000000e+09 6.6548423840e+09 0.6500000000 -4.0438444754 -13.6748739586 1 - 1.0000000000e+09 7.3680629973e+09 0.6500000000 -4.5739409491 -13.1888195133 1 - 1.0000000000e+09 8.1577217310e+09 0.6500000000 -5.1712163628 -12.6507561248 1 - 1.0000000000e+09 9.0320107014e+09 0.6500000000 -5.8437444564 -12.0829813034 1 - 1.0000000000e+09 1.0000000000e+10 0.6500000000 -6.5994405224 -11.6765113113 1 - 1.0000000000e+10 4.0000000000e+09 0.6500000000 -1.3793655235 -16.3247814756 0 - 1.0000000000e+10 4.4286927156e+09 0.6500000000 -1.5844565484 -16.1362774695 1 - 1.0000000000e+10 4.9033297922e+09 0.6500000000 -1.8178191791 -15.9217568381 1 - 1.0000000000e+10 5.4288352332e+09 0.6500000000 -2.0831515252 -15.6779340065 1 - 1.0000000000e+10 6.0106607628e+09 0.6500000000 -2.3846197513 -15.4010550224 1 - 1.0000000000e+10 6.6548423840e+09 0.6500000000 -2.7269159436 -15.0872983686 1 - 1.0000000000e+10 7.3680629973e+09 0.6500000000 -3.1153229999 -14.7323215524 1 - 1.0000000000e+10 8.1577217310e+09 0.6500000000 -3.5557871014 -14.3318302438 1 - 1.0000000000e+10 9.0320107014e+09 0.6500000000 -4.0549949963 -13.8820225244 1 - 1.0000000000e+10 1.0000000000e+10 0.6500000000 -4.6204256792 -13.3808512952 1 - 1.0000000000e+06 4.0000000000e+09 0.7000000000 -4.4485209891 -13.1143251183 0 - 1.0000000000e+06 4.4286927156e+09 0.7000000000 -4.9835499142 -12.6117081124 1 - 1.0000000000e+06 4.9033297922e+09 0.7000000000 -5.5821660601 -12.0561394505 1 - 1.0000000000e+06 5.4288352332e+09 0.7000000000 -6.2518496554 -11.5720368075 1 - 1.0000000000e+06 6.0106607628e+09 0.7000000000 -7.0008334984 -11.2561651451 1 - 1.0000000000e+06 6.6548423840e+09 0.7000000000 -7.8361333603 -10.9067344717 1 - 1.0000000000e+06 7.3680629973e+09 0.7000000000 -8.7392243787 -10.5687792383 1 - 1.0000000000e+06 8.1577217310e+09 0.7000000000 -9.6021916703 -10.5229827098 1 - 1.0000000000e+06 9.0320107014e+09 0.7000000000 -10.5729938190 -11.2517560750 1 - 1.0000000000e+06 1.0000000000e+10 0.7000000000 -11.8235850796 -12.5558647574 1 - 1.0000000000e+07 4.0000000000e+09 0.7000000000 -3.6550369761 -13.8914003282 0 - 1.0000000000e+07 4.4286927156e+09 0.7000000000 -4.1049942046 -13.4694919130 1 - 1.0000000000e+07 4.9033297922e+09 0.7000000000 -4.6094232304 -12.9974776854 1 - 1.0000000000e+07 5.4288352332e+09 0.7000000000 -5.1748308102 -12.4725550651 1 - 1.0000000000e+07 6.0106607628e+09 0.7000000000 -5.8084870128 -11.9131933897 1 - 1.0000000000e+07 6.6548423840e+09 0.7000000000 -6.5184625267 -11.5626253210 1 - 1.0000000000e+07 7.3680629973e+09 0.7000000000 -7.3127774909 -11.2390485970 1 - 1.0000000000e+07 8.1577217310e+09 0.7000000000 -8.1901397748 -10.8966710259 1 - 1.0000000000e+07 9.0320107014e+09 0.7000000000 -9.0920672122 -10.6613615921 1 - 1.0000000000e+07 1.0000000000e+10 0.7000000000 -9.9582875767 -10.8797496537 1 - 1.0000000000e+08 4.0000000000e+09 0.7000000000 -2.8616533998 -14.6857719715 0 - 1.0000000000e+08 4.4286927156e+09 0.7000000000 -3.2265203806 -14.3452898486 1 - 1.0000000000e+08 4.9033297922e+09 0.7000000000 -3.6367434104 -13.9627369387 1 - 1.0000000000e+08 5.4288352332e+09 0.7000000000 -4.0978526310 -13.5335594957 1 - 1.0000000000e+08 6.0106607628e+09 0.7000000000 -4.6160423318 -13.0536842248 1 - 1.0000000000e+08 6.6548423840e+09 0.7000000000 -5.1982455454 -12.5235029605 1 - 1.0000000000e+08 7.3680629973e+09 0.7000000000 -5.8525289991 -11.9932157006 1 - 1.0000000000e+08 8.1577217310e+09 0.7000000000 -6.5865650606 -11.6818564067 1 - 1.0000000000e+08 9.0320107014e+09 0.7000000000 -7.4069463043 -11.3651455549 1 - 1.0000000000e+08 1.0000000000e+10 0.7000000000 -8.3014810042 -11.0535468230 1 - 1.0000000000e+09 4.0000000000e+09 0.7000000000 -2.0694628674 -15.5178474652 0 - 1.0000000000e+09 4.4286927156e+09 0.7000000000 -2.3491916894 -15.2587909978 1 - 1.0000000000e+09 4.9033297922e+09 0.7000000000 -2.6651630630 -14.9662311049 1 - 1.0000000000e+09 5.4288352332e+09 0.7000000000 -3.0219300320 -14.6361670520 1 - 1.0000000000e+09 6.0106607628e+09 0.7000000000 -3.4246053577 -14.2641972416 1 - 1.0000000000e+09 6.6548423840e+09 0.7000000000 -3.8789291126 -13.8461550124 1 - 1.0000000000e+09 7.3680629973e+09 0.7000000000 -4.3913430475 -13.3784836554 1 - 1.0000000000e+09 8.1577217310e+09 0.7000000000 -4.9690605512 -12.8619701183 1 - 1.0000000000e+09 9.0320107014e+09 0.7000000000 -5.6200778085 -12.3291976596 1 - 1.0000000000e+09 1.0000000000e+10 0.7000000000 -6.3527292745 -11.9889184561 1 - 1.0000000000e+10 4.0000000000e+09 0.7000000000 -1.2804483941 -16.4327248844 0 - 1.0000000000e+10 4.4286927156e+09 0.7000000000 -1.4749011906 -16.2550319759 1 - 1.0000000000e+10 4.9033297922e+09 0.7000000000 -1.6964869006 -16.0524804719 1 - 1.0000000000e+10 5.4288352332e+09 0.7000000000 -1.9487815730 -15.8219089869 1 - 1.0000000000e+10 6.0106607628e+09 0.7000000000 -2.2358162273 -15.5597016802 1 - 1.0000000000e+10 6.6548423840e+09 0.7000000000 -2.5621333174 -15.2621928113 1 - 1.0000000000e+10 7.3680629973e+09 0.7000000000 -2.9328500820 -14.9252211738 1 - 1.0000000000e+10 8.1577217310e+09 0.7000000000 -3.3537294917 -14.5447208945 1 - 1.0000000000e+10 9.0320107014e+09 0.7000000000 -3.8312578887 -14.1172251767 1 - 1.0000000000e+10 1.0000000000e+10 0.7000000000 -4.3727140958 -13.6413288908 1 -AMReX (23.07-395-ge6c93bf22695) finalized + 1.0000000000e+09 4.4286927156e+09 0.4000000000 -15.0933118588 -2.6671723127 1 + 1.0000000000e+09 4.9033297922e+09 0.4000000000 -14.7621490285 -3.0175641894 1 + 1.0000000000e+09 5.4288352332e+09 0.4000000000 -14.3892566633 -3.4124048303 1 + 1.0000000000e+09 6.0106607628e+09 0.4000000000 -13.9696329645 -3.8570452607 1 + 1.0000000000e+09 6.6548423840e+09 0.4000000000 -13.4990899024 -4.3569659385 1 + 1.0000000000e+09 7.3680629973e+09 0.4000000000 -12.9751097903 -4.9165783433 1 + 1.0000000000e+09 8.1577217310e+09 0.4000000000 -12.4029718543 -5.5362858494 1 + 1.0000000000e+09 9.0320107014e+09 0.4000000000 -11.8075831340 -6.2097753499 1 + 1.0000000000e+09 1.0000000000e+10 0.4000000000 -11.2833005757 -6.9502991890 1 + 1.0000000000e+10 4.0000000000e+09 0.4000000000 -16.2582552424 -1.5626567049 0 + 1.0000000000e+10 4.4286927156e+09 0.4000000000 -16.0486695196 -1.7884208146 1 + 1.0000000000e+10 4.9033297922e+09 0.4000000000 -15.8108280439 -2.0446360339 1 + 1.0000000000e+10 5.4288352332e+09 0.4000000000 -15.5412395159 -2.3352407251 1 + 1.0000000000e+10 6.0106607628e+09 0.4000000000 -15.2357440912 -2.6646601015 1 + 1.0000000000e+10 6.6548423840e+09 0.4000000000 -14.8903553070 -3.0378300279 1 + 1.0000000000e+10 7.3680629973e+09 0.4000000000 -14.5003230284 -3.4601447158 1 + 1.0000000000e+10 8.1577217310e+09 0.4000000000 -14.0614573455 -3.9371815118 1 + 1.0000000000e+10 9.0320107014e+09 0.4000000000 -13.5711738956 -4.4739423061 1 + 1.0000000000e+10 1.0000000000e+10 0.4000000000 -13.0312228621 -5.0733622012 1 + 1.0000000000e+06 4.0000000000e+09 0.4500000000 -12.6199619882 -5.0493824492 0 + 1.0000000000e+06 4.4286927156e+09 0.4500000000 -12.0411567544 -5.6445833881 1 + 1.0000000000e+06 4.9033297922e+09 0.4500000000 -11.4232232546 -6.2868541638 1 + 1.0000000000e+06 5.4288352332e+09 0.4500000000 -10.8518135789 -6.9298795402 1 + 1.0000000000e+06 6.0106607628e+09 0.4500000000 -10.4348297049 -7.7172920610 1 + 1.0000000000e+06 6.6548423840e+09 0.4500000000 -10.0424249574 -8.5877937253 1 + 1.0000000000e+06 7.3680629973e+09 0.4500000000 -9.8849653906 -9.3214540023 1 + 1.0000000000e+06 8.1577217310e+09 0.4500000000 -10.1350525691 -9.9053582092 1 + 1.0000000000e+06 9.0320107014e+09 0.4500000000 -10.9317596586 -10.7720005416 1 + 1.0000000000e+06 1.0000000000e+10 0.4500000000 -12.2046515153 -12.0333675716 1 + 1.0000000000e+07 4.0000000000e+09 0.4500000000 -13.4182964630 -4.2562147910 0 + 1.0000000000e+07 4.4286927156e+09 0.4500000000 -12.9196143636 -4.7702514625 1 + 1.0000000000e+07 4.9033297922e+09 0.4500000000 -12.3646713416 -5.3431229151 1 + 1.0000000000e+07 5.4288352332e+09 0.4500000000 -11.7607991622 -5.9706925457 1 + 1.0000000000e+07 6.0106607628e+09 0.4500000000 -11.1530047519 -6.6222738557 1 + 1.0000000000e+07 6.6548423840e+09 0.4500000000 -10.6534176006 -7.3144511945 1 + 1.0000000000e+07 7.3680629973e+09 0.4500000000 -10.2538965760 -8.1710984842 1 + 1.0000000000e+07 8.1577217310e+09 0.4500000000 -9.9754563098 -8.9846245478 1 + 1.0000000000e+07 9.0320107014e+09 0.4500000000 -10.0361948008 -9.6116151464 1 + 1.0000000000e+07 1.0000000000e+10 0.4500000000 -10.4770374370 -10.2529619738 1 + 1.0000000000e+08 4.0000000000e+09 0.4500000000 -14.2257599781 -3.4625787929 0 + 1.0000000000e+08 4.4286927156e+09 0.4500000000 -13.8113629028 -3.8919431463 1 + 1.0000000000e+08 4.9033297922e+09 0.4500000000 -13.3469606757 -4.3733520313 1 + 1.0000000000e+08 5.4288352332e+09 0.4500000000 -12.8281609635 -4.9118956632 1 + 1.0000000000e+08 6.0106607628e+09 0.4500000000 -12.2549868876 -5.5090642377 1 + 1.0000000000e+08 6.6548423840e+09 0.4500000000 -11.6452720540 -6.1532052811 1 + 1.0000000000e+08 7.3680629973e+09 0.4500000000 -11.0519682942 -6.8137717934 1 + 1.0000000000e+08 8.1577217310e+09 0.4500000000 -10.5745250527 -7.5537012214 1 + 1.0000000000e+08 9.0320107014e+09 0.4500000000 -10.2098734359 -8.4106723339 1 + 1.0000000000e+08 1.0000000000e+10 0.4500000000 -10.0528603240 -9.1577967123 1 + 1.0000000000e+09 4.0000000000e+09 0.4500000000 -15.0536228397 -2.6688948470 0 + 1.0000000000e+09 4.4286927156e+09 0.4500000000 -14.7237134075 -3.0132309281 1 + 1.0000000000e+09 4.9033297922e+09 0.4500000000 -14.3526160907 -3.4007197531 1 + 1.0000000000e+09 5.4288352332e+09 0.4500000000 -13.9355878021 -3.8365577153 1 + 1.0000000000e+09 6.0106607628e+09 0.4500000000 -13.4675884154 -4.3261358021 1 + 1.0000000000e+09 6.6548423840e+09 0.4500000000 -12.9456631769 -4.8738173937 1 + 1.0000000000e+09 7.3680629973e+09 0.4500000000 -12.3724824027 -5.4792910380 1 + 1.0000000000e+09 8.1577217310e+09 0.4500000000 -11.7661786775 -6.1311738817 1 + 1.0000000000e+09 9.0320107014e+09 0.4500000000 -11.1669148875 -6.8085299247 1 + 1.0000000000e+09 1.0000000000e+10 0.4500000000 -10.6514166226 -7.5351746507 1 + 1.0000000000e+10 4.0000000000e+09 0.4500000000 -15.9281437071 -1.8752022471 0 + 1.0000000000e+10 4.4286927156e+09 0.4500000000 -15.6824809806 -2.1344670468 1 + 1.0000000000e+10 4.9033297922e+09 0.4500000000 -15.4047001663 -2.4277705612 1 + 1.0000000000e+10 5.4288352332e+09 0.4500000000 -15.0909187931 -2.7594143654 1 + 1.0000000000e+10 6.0106607628e+09 0.4500000000 -14.7366089802 -3.1341695123 1 + 1.0000000000e+10 6.6548423840e+09 0.4500000000 -14.3375645420 -3.5571892970 1 + 1.0000000000e+10 7.3680629973e+09 0.4500000000 -13.8893254086 -4.0336353982 1 + 1.0000000000e+10 8.1577217310e+09 0.4500000000 -13.3893298655 -4.5676574853 1 + 1.0000000000e+10 9.0320107014e+09 0.4500000000 -12.8393736234 -5.1603429531 1 + 1.0000000000e+10 1.0000000000e+10 0.4500000000 -12.2492265137 -5.8073288375 1 + 1.0000000000e+06 4.0000000000e+09 0.5000000000 -8.0664194354 -9.4890393815 0 + 1.0000000000e+06 4.4286927156e+09 0.5000000000 -8.2111046670 -9.3757376097 1 + 1.0000000000e+06 4.9033297922e+09 0.5000000000 -8.3918234083 -9.2344738095 1 + 1.0000000000e+06 5.4288352332e+09 0.5000000000 -8.7772909062 -8.9448205348 1 + 1.0000000000e+06 6.0106607628e+09 0.5000000000 -9.0636971271 -9.0647301302 1 + 1.0000000000e+06 6.6548423840e+09 0.5000000000 -9.3058452726 -9.3069927949 1 + 1.0000000000e+06 7.3680629973e+09 0.5000000000 -9.5983657200 -9.5996435483 1 + 1.0000000000e+06 8.1577217310e+09 0.5000000000 -10.0168321560 -10.0182538936 1 + 1.0000000000e+06 9.0320107014e+09 0.5000000000 -10.8474057609 -10.8489864867 1 + 1.0000000000e+06 1.0000000000e+10 0.5000000000 -12.1138078856 -12.1155643114 1 + 1.0000000000e+07 4.0000000000e+09 0.5000000000 -7.8838825215 -9.6537090536 0 + 1.0000000000e+07 4.4286927156e+09 0.5000000000 -8.0018897889 -9.5630616368 1 + 1.0000000000e+07 4.9033297922e+09 0.5000000000 -8.1353862644 -9.4610110058 1 + 1.0000000000e+07 5.4288352332e+09 0.5000000000 -8.2961790243 -9.3385091392 1 + 1.0000000000e+07 6.0106607628e+09 0.5000000000 -8.5443481247 -9.1484241687 1 + 1.0000000000e+07 6.6548423840e+09 0.5000000000 -8.9669746383 -8.9715401000 1 + 1.0000000000e+07 7.3680629973e+09 0.5000000000 -9.2004551976 -9.2016374220 1 + 1.0000000000e+07 8.1577217310e+09 0.5000000000 -9.4726511197 -9.4739817205 1 + 1.0000000000e+07 9.0320107014e+09 0.5000000000 -9.8195004568 -9.8209945474 1 + 1.0000000000e+07 1.0000000000e+10 0.5000000000 -10.3608512486 -10.3625253205 1 + 1.0000000000e+08 4.0000000000e+09 0.5000000000 -7.7564783284 -9.7786310606 0 + 1.0000000000e+08 4.4286927156e+09 0.5000000000 -7.8492582957 -9.7096153678 1 + 1.0000000000e+08 4.9033297922e+09 0.5000000000 -7.9528574806 -9.6329691536 1 + 1.0000000000e+08 5.4288352332e+09 0.5000000000 -8.0693126926 -9.5475431816 1 + 1.0000000000e+08 6.0106607628e+09 0.5000000000 -8.2064635264 -9.4472744455 1 + 1.0000000000e+08 6.6548423840e+09 0.5000000000 -8.3859624449 -9.3159632798 1 + 1.0000000000e+08 7.3680629973e+09 0.5000000000 -8.6892137197 -9.0970157579 1 + 1.0000000000e+08 8.1577217310e+09 0.5000000000 -9.0475146736 -9.0502773705 1 + 1.0000000000e+08 9.0320107014e+09 0.5000000000 -9.2986738307 -9.2998953589 1 + 1.0000000000e+08 1.0000000000e+10 0.5000000000 -9.5981915843 -9.5996058564 1 + 1.0000000000e+09 4.0000000000e+09 0.5000000000 -7.7470944343 -9.8201166535 0 + 1.0000000000e+09 4.4286927156e+09 0.5000000000 -7.8148151276 -9.7724887691 1 + 1.0000000000e+09 4.9033297922e+09 0.5000000000 -7.8900133105 -9.7200543142 1 + 1.0000000000e+09 5.4288352332e+09 0.5000000000 -7.9730314841 -9.6630087712 1 + 1.0000000000e+09 6.0106607628e+09 0.5000000000 -8.0660308619 -9.5998690277 1 + 1.0000000000e+09 6.6548423840e+09 0.5000000000 -8.1711502365 -9.5299746078 1 + 1.0000000000e+09 7.3680629973e+09 0.5000000000 -8.2981193470 -9.4461659460 1 + 1.0000000000e+09 8.1577217310e+09 0.5000000000 -8.4688942226 -9.3329031919 1 + 1.0000000000e+09 9.0320107014e+09 0.5000000000 -8.7322319942 -9.1628388393 1 + 1.0000000000e+09 1.0000000000e+10 0.5000000000 -9.0635018458 -9.0827468205 1 + 1.0000000000e+10 4.0000000000e+09 0.5000000000 -7.9880361426 -9.6879018121 1 + 1.0000000000e+10 4.4286927156e+09 0.5000000000 -8.0302837106 -9.6618767633 1 + 1.0000000000e+10 4.9033297922e+09 0.5000000000 -8.0763904321 -9.6342068332 1 + 1.0000000000e+10 5.4288352332e+09 0.5000000000 -8.1258757774 -9.6057820019 1 + 1.0000000000e+10 6.0106607628e+09 0.5000000000 -8.1795463024 -9.5762054669 1 + 1.0000000000e+10 6.6548423840e+09 0.5000000000 -8.2365564026 -9.5471909647 1 + 1.0000000000e+10 7.3680629973e+09 0.5000000000 -8.2989177606 -9.5176692137 1 + 1.0000000000e+10 8.1577217310e+09 0.5000000000 -8.3696249312 -9.4863119797 1 + 1.0000000000e+10 9.0320107014e+09 0.5000000000 -8.4560400376 -9.4486905350 1 + 1.0000000000e+10 1.0000000000e+10 0.5000000000 -8.5725042842 -9.3959527205 1 + 1.0000000000e+06 4.0000000000e+09 0.5500000000 -4.9263369035 -12.6307046567 0 + 1.0000000000e+06 4.4286927156e+09 0.5500000000 -5.5125794078 -12.0760158044 1 + 1.0000000000e+06 4.9033297922e+09 0.5500000000 -6.1678880195 -11.4605459733 1 + 1.0000000000e+06 5.4288352332e+09 0.5500000000 -6.9002340024 -10.8361713231 1 + 1.0000000000e+06 6.0106607628e+09 0.5500000000 -7.7162632259 -10.4358694538 1 + 1.0000000000e+06 6.6548423840e+09 0.5500000000 -8.5866490983 -10.0435826768 1 + 1.0000000000e+06 7.3680629973e+09 0.5500000000 -9.3201805492 -9.8862508729 1 + 1.0000000000e+06 8.1577217310e+09 0.5500000000 -9.9039408495 -10.1364796930 1 + 1.0000000000e+06 9.0320107014e+09 0.5500000000 -10.7704221565 -10.9333429140 1 + 1.0000000000e+06 1.0000000000e+10 0.5500000000 -12.0316130148 -12.2064099076 1 + 1.0000000000e+07 4.0000000000e+09 0.5500000000 -4.1327873451 -13.4067898896 0 + 1.0000000000e+07 4.4286927156e+09 0.5500000000 -4.6339615895 -12.9331082034 1 + 1.0000000000e+07 4.9033297922e+09 0.5500000000 -5.1950956687 -12.4035746041 1 + 1.0000000000e+07 5.4288352332e+09 0.5500000000 -5.8232748773 -11.8138715673 1 + 1.0000000000e+07 6.0106607628e+09 0.5500000000 -6.5262585311 -11.1696968936 1 + 1.0000000000e+07 6.6548423840e+09 0.5500000000 -7.3117723591 -10.6534502323 1 + 1.0000000000e+07 7.3680629973e+09 0.5500000000 -8.1699241087 -10.2551099463 1 + 1.0000000000e+07 8.1577217310e+09 0.5500000000 -8.9833050294 -9.9768125421 1 + 1.0000000000e+07 9.0320107014e+09 0.5500000000 -9.6101346160 -10.0377081495 1 + 1.0000000000e+07 1.0000000000e+10 0.5500000000 -10.2512990136 -10.4787244918 1 + 1.0000000000e+08 4.0000000000e+09 0.5500000000 -3.3395295007 -14.1985085219 0 + 1.0000000000e+08 4.4286927156e+09 0.5500000000 -3.7556211435 -13.8063032478 1 + 1.0000000000e+08 4.9033297922e+09 0.5500000000 -4.2225575794 -13.3664559653 1 + 1.0000000000e+08 5.4288352332e+09 0.5500000000 -4.7464569494 -12.8737301627 1 + 1.0000000000e+08 6.0106607628e+09 0.5500000000 -5.3341586613 -12.3230169469 1 + 1.0000000000e+08 6.6548423840e+09 0.5500000000 -5.9932141535 -11.7120144662 1 + 1.0000000000e+08 7.3680629973e+09 0.5500000000 -6.7312739156 -11.0587794004 1 + 1.0000000000e+08 8.1577217310e+09 0.5500000000 -7.5517647875 -10.5751303759 1 + 1.0000000000e+08 9.0320107014e+09 0.5500000000 -8.4094768040 -10.2111789283 1 + 1.0000000000e+08 1.0000000000e+10 0.5500000000 -9.1564178754 -10.0543419536 1 + 1.0000000000e+09 4.0000000000e+09 0.5500000000 -2.5471797995 -15.0250730353 0 + 1.0000000000e+09 4.4286927156e+09 0.5500000000 -2.8781461534 -14.7143116038 1 + 1.0000000000e+09 4.9033297922e+09 0.5500000000 -3.2508437389 -14.3645024308 1 + 1.0000000000e+09 5.4288352332e+09 0.5500000000 -3.6704138417 -13.9710418050 1 + 1.0000000000e+09 6.0106607628e+09 0.5500000000 -4.1426204134 -13.5288354140 1 + 1.0000000000e+09 6.6548423840e+09 0.5500000000 -4.6739224526 -13.0328615008 1 + 1.0000000000e+09 7.3680629973e+09 0.5500000000 -5.2715267295 -12.4783156933 1 + 1.0000000000e+09 8.1577217310e+09 0.5500000000 -5.9431603152 -11.8633393484 1 + 1.0000000000e+09 9.0320107014e+09 0.5500000000 -6.6951903623 -11.2024350566 1 + 1.0000000000e+09 1.0000000000e+10 0.5500000000 -7.5250446792 -10.6480296449 1 + 1.0000000000e+10 4.0000000000e+09 0.5500000000 -1.7576144396 -15.9280217909 1 + 1.0000000000e+10 4.4286927156e+09 0.5500000000 -2.0033314082 -15.6986279517 1 + 1.0000000000e+10 4.9033297922e+09 0.5500000000 -2.2816705223 -15.4388374567 1 + 1.0000000000e+10 5.4288352332e+09 0.5500000000 -2.5967959027 -15.1448960799 1 + 1.0000000000e+10 6.0106607628e+09 0.5500000000 -2.9533897196 -14.8125306367 1 + 1.0000000000e+10 6.6548423840e+09 0.5500000000 -3.3567153951 -14.4373412142 1 + 1.0000000000e+10 7.3680629973e+09 0.5500000000 -3.8126881350 -14.0143346733 1 + 1.0000000000e+10 8.1577217310e+09 0.5500000000 -4.3279498369 -13.5384648003 1 + 1.0000000000e+10 9.0320107014e+09 0.5500000000 -4.9099077591 -13.0050412453 1 + 1.0000000000e+10 1.0000000000e+10 0.5500000000 -5.5664079579 -12.4111302237 1 + 1.0000000000e+06 4.0000000000e+09 0.6000000000 -4.6874165138 -12.8713534471 0 + 1.0000000000e+06 4.4286927156e+09 0.6000000000 -5.2480529163 -12.3425069636 1 + 1.0000000000e+06 4.9033297922e+09 0.6000000000 -5.8750189170 -11.7561777822 1 + 1.0000000000e+06 5.4288352332e+09 0.6000000000 -6.5760720464 -11.1834054668 1 + 1.0000000000e+06 6.0106607628e+09 0.6000000000 -7.3593255585 -10.8231744084 1 + 1.0000000000e+06 6.6548423840e+09 0.6000000000 -8.2240941596 -10.4375089980 1 + 1.0000000000e+06 7.3680629973e+09 0.6000000000 -9.0836917067 -10.1457704344 1 + 1.0000000000e+06 8.1577217310e+09 0.6000000000 -9.7969116846 -10.2597064355 1 + 1.0000000000e+06 9.0320107014e+09 0.6000000000 -10.6995852342 -11.0267695004 1 + 1.0000000000e+06 1.0000000000e+10 0.6000000000 -11.9565611188 -12.3079449877 1 + 1.0000000000e+07 4.0000000000e+09 0.6000000000 -3.8938729829 -13.6478070662 0 + 1.0000000000e+07 4.4286927156e+09 0.6000000000 -4.3694406419 -13.1998899866 1 + 1.0000000000e+07 4.9033297922e+09 0.6000000000 -4.9022241247 -12.6989326129 1 + 1.0000000000e+07 5.4288352332e+09 0.6000000000 -5.4990224481 -12.1411327078 1 + 1.0000000000e+07 6.0106607628e+09 0.6000000000 -6.1674007417 -11.5348079253 1 + 1.0000000000e+07 6.6548423840e+09 0.6000000000 -6.9155791417 -11.0832126364 1 + 1.0000000000e+07 7.3680629973e+09 0.6000000000 -7.7481347730 -10.7130540280 1 + 1.0000000000e+07 8.1577217310e+09 0.6000000000 -8.6298769859 -10.3622929168 1 + 1.0000000000e+07 9.0320107014e+09 0.6000000000 -9.4167978029 -10.2526624488 1 + 1.0000000000e+07 1.0000000000e+10 0.6000000000 -10.1478611795 -10.6024621642 1 + 1.0000000000e+08 4.0000000000e+09 0.6000000000 -3.1006340152 -14.4403947671 0 + 1.0000000000e+08 4.4286927156e+09 0.6000000000 -3.4911181626 -14.0739382686 1 + 1.0000000000e+08 4.9033297922e+09 0.6000000000 -3.9297028770 -13.6626066662 1 + 1.0000000000e+08 5.4288352332e+09 0.6000000000 -4.4222123883 -13.2014814753 1 + 1.0000000000e+08 6.0106607628e+09 0.6000000000 -4.9751654635 -12.6858717662 1 + 1.0000000000e+08 6.6548423840e+09 0.6000000000 -5.5958360589 -12.1142361002 1 + 1.0000000000e+08 7.3680629973e+09 0.6000000000 -6.2921951408 -11.5102059086 1 + 1.0000000000e+08 8.1577217310e+09 0.6000000000 -7.0718360681 -11.0958538127 1 + 1.0000000000e+08 9.0320107014e+09 0.6000000000 -7.9295315244 -10.7329004795 1 + 1.0000000000e+08 1.0000000000e+10 0.6000000000 -8.7971964340 -10.4475908920 1 + 1.0000000000e+09 4.0000000000e+09 0.6000000000 -2.3083427797 -15.2688913657 0 + 1.0000000000e+09 4.4286927156e+09 0.6000000000 -2.6136989393 -14.9838710003 1 + 1.0000000000e+09 4.9033297922e+09 0.6000000000 -2.9580422061 -14.6625625008 1 + 1.0000000000e+09 5.4288352332e+09 0.6000000000 -3.3462197500 -14.3006599290 1 + 1.0000000000e+09 6.0106607628e+09 0.6000000000 -3.7836699672 -13.8934047622 1 + 1.0000000000e+09 6.6548423840e+09 0.6000000000 -4.2764930786 -13.4361677469 1 + 1.0000000000e+09 7.3680629973e+09 0.6000000000 -4.8315236918 -12.9246340459 1 + 1.0000000000e+09 8.1577217310e+09 0.6000000000 -5.4563463669 -12.3575652127 1 + 1.0000000000e+09 9.0320107014e+09 0.6000000000 -6.1589448250 -11.7514500907 1 + 1.0000000000e+09 1.0000000000e+10 0.6000000000 -6.9457107700 -11.2744549299 1 + 1.0000000000e+10 4.0000000000e+09 0.6000000000 -1.5189550687 -16.1760474376 0 + 1.0000000000e+10 4.4286927156e+09 0.6000000000 -1.7390541688 -15.9723914831 1 + 1.0000000000e+10 4.9033297922e+09 0.6000000000 -1.9890315435 -15.7410961693 1 + 1.0000000000e+10 5.4288352332e+09 0.6000000000 -2.2727571929 -15.4787024815 1 + 1.0000000000e+10 6.0106607628e+09 0.6000000000 -2.5945876455 -15.1812641094 1 + 1.0000000000e+10 6.6548423840e+09 0.6000000000 -2.9594257944 -14.8447444844 1 + 1.0000000000e+10 7.3680629973e+09 0.6000000000 -3.3727879321 -14.4645568357 1 + 1.0000000000e+10 8.1577217310e+09 0.6000000000 -3.8408780345 -14.0361164439 1 + 1.0000000000e+10 9.0320107014e+09 0.6000000000 -4.3706608771 -13.5552459151 1 + 1.0000000000e+10 1.0000000000e+10 0.6000000000 -4.9698583638 -13.0193464360 1 + 1.0000000000e+06 4.0000000000e+09 0.6500000000 -4.5476580932 -13.0130256076 0 + 1.0000000000e+06 4.4286927156e+09 0.6500000000 -5.0933157554 -12.4994381879 1 + 1.0000000000e+06 4.9033297922e+09 0.6500000000 -5.7036990969 -11.9307186936 1 + 1.0000000000e+06 5.4288352332e+09 0.6500000000 -6.3864064622 -11.4022027916 1 + 1.0000000000e+06 6.0106607628e+09 0.6500000000 -7.1496914530 -11.0673499344 1 + 1.0000000000e+06 6.6548423840e+09 0.6500000000 -7.9986226626 -10.7004084995 1 + 1.0000000000e+06 7.3680629973e+09 0.6500000000 -8.8932097425 -10.3705250415 1 + 1.0000000000e+06 8.1577217310e+09 0.6500000000 -9.6962265139 -10.3881460503 1 + 1.0000000000e+06 9.0320107014e+09 0.6500000000 -10.6339946188 -11.1316988308 1 + 1.0000000000e+06 1.0000000000e+10 0.6500000000 -11.8875095425 -12.4230319679 1 + 1.0000000000e+07 4.0000000000e+09 0.6500000000 -3.7541203442 -13.7898249789 0 + 1.0000000000e+07 4.4286927156e+09 0.6500000000 -4.2147089469 -13.3570697965 1 + 1.0000000000e+07 4.9033297922e+09 0.6500000000 -4.7309081311 -12.8729632474 1 + 1.0000000000e+07 5.4288352332e+09 0.6500000000 -5.3093462612 -12.3341734101 1 + 1.0000000000e+07 6.0106607628e+09 0.6500000000 -5.9574209643 -11.7528786941 1 + 1.0000000000e+07 6.6548423840e+09 0.6500000000 -6.6832983524 -11.3537044311 1 + 1.0000000000e+07 7.3680629973e+09 0.6500000000 -7.4941774420 -11.0089926944 1 + 1.0000000000e+07 8.1577217310e+09 0.6500000000 -8.3792738800 -10.6556306767 1 + 1.0000000000e+07 9.0320107014e+09 0.6500000000 -9.2440179618 -10.4606562951 1 + 1.0000000000e+07 1.0000000000e+10 0.6500000000 -10.0502406180 -10.7352829937 1 + 1.0000000000e+08 4.0000000000e+09 0.6500000000 -2.9608994723 -14.5832351693 0 + 1.0000000000e+08 4.4286927156e+09 0.6500000000 -3.3364036964 -14.2319225764 1 + 1.0000000000e+08 4.9033297922e+09 0.6500000000 -3.7584032382 -13.8373671291 1 + 1.0000000000e+08 5.4288352332e+09 0.6500000000 -4.2325503160 -13.3948530013 1 + 1.0000000000e+08 6.0106607628e+09 0.6500000000 -4.7651740951 -12.9000045811 1 + 1.0000000000e+08 6.6548423840e+09 0.6500000000 -5.3633523384 -12.3521350793 1 + 1.0000000000e+08 7.3680629973e+09 0.6500000000 -6.0349472150 -11.7846818528 1 + 1.0000000000e+08 8.1577217310e+09 0.6500000000 -6.7881653086 -11.4262642633 1 + 1.0000000000e+08 9.0320107014e+09 0.6500000000 -7.6263287309 -11.0865976478 1 + 1.0000000000e+08 1.0000000000e+10 0.6500000000 -8.5198307182 -10.7737584967 1 + 1.0000000000e+09 4.0000000000e+09 0.6500000000 -2.1686642542 -15.4135601959 0 + 1.0000000000e+09 4.4286927156e+09 0.6500000000 -2.4590379093 -15.1436760063 1 + 1.0000000000e+09 4.9033297922e+09 0.6500000000 -2.7867935242 -14.8391271318 1 + 1.0000000000e+09 5.4288352332e+09 0.6500000000 -3.1566062187 -14.4957844807 1 + 1.0000000000e+09 6.0106607628e+09 0.6500000000 -3.5737239146 -14.1090934401 1 + 1.0000000000e+09 6.6548423840e+09 0.6500000000 -4.0440362374 -13.6746794679 1 + 1.0000000000e+09 7.3680629973e+09 0.6500000000 -4.5741490489 -13.1886122154 1 + 1.0000000000e+09 8.1577217310e+09 0.6500000000 -5.1714414325 -12.6505348485 1 + 1.0000000000e+09 9.0320107014e+09 0.6500000000 -5.8439871511 -12.0827672096 1 + 1.0000000000e+09 1.0000000000e+10 0.6500000000 -6.5997012945 -11.6763905107 1 + 1.0000000000e+10 4.0000000000e+09 0.6500000000 -1.3794466033 -16.3246961572 0 + 1.0000000000e+10 4.4286927156e+09 0.6500000000 -1.5845558746 -16.1361734166 1 + 1.0000000000e+10 4.9033297922e+09 0.6500000000 -1.8179385128 -15.9216322620 1 + 1.0000000000e+10 5.4288352332e+09 0.6500000000 -2.0832924664 -15.6777872875 1 + 1.0000000000e+10 6.0106607628e+09 0.6500000000 -2.3847837732 -15.4008846806 1 + 1.0000000000e+10 6.6548423840e+09 0.6500000000 -2.7271044222 -15.0871030464 1 + 1.0000000000e+10 7.3680629973e+09 0.6500000000 -3.1155372402 -14.7321000031 1 + 1.0000000000e+10 8.1577217310e+09 0.6500000000 -3.5560283588 -14.3315813824 1 + 1.0000000000e+10 9.0320107014e+09 0.6500000000 -4.0552644947 -13.8817455175 1 + 1.0000000000e+10 1.0000000000e+10 0.6500000000 -4.6207246098 -13.3805467847 1 + 1.0000000000e+06 4.0000000000e+09 0.7000000000 -4.4484985951 -13.1143475495 0 + 1.0000000000e+06 4.4286927156e+09 0.7000000000 -4.9835286255 -12.6117296860 1 + 1.0000000000e+06 4.9033297922e+09 0.7000000000 -5.5821458226 -12.0561621904 1 + 1.0000000000e+06 5.4288352332e+09 0.7000000000 -6.2518304171 -11.5720026737 1 + 1.0000000000e+06 6.0106607628e+09 0.7000000000 -7.0008152024 -11.2561319512 1 + 1.0000000000e+06 6.6548423840e+09 0.7000000000 -7.8361158291 -10.9067031322 1 + 1.0000000000e+06 7.3680629973e+09 0.7000000000 -8.7392063208 -10.5687519073 1 + 1.0000000000e+06 8.1577217310e+09 0.7000000000 -9.6021715722 -10.5229664649 1 + 1.0000000000e+06 9.0320107014e+09 0.7000000000 -10.5729780670 -11.2517541081 1 + 1.0000000000e+06 1.0000000000e+10 0.7000000000 -11.8235708727 -12.5558647155 1 + 1.0000000000e+07 4.0000000000e+09 0.7000000000 -3.6549664086 -13.8914709006 0 + 1.0000000000e+07 4.4286927156e+09 0.7000000000 -4.1049270968 -13.4695590574 1 + 1.0000000000e+07 4.9033297922e+09 0.7000000000 -4.6093594165 -12.9975417362 1 + 1.0000000000e+07 5.4288352332e+09 0.7000000000 -5.1747701322 -12.4726172619 1 + 1.0000000000e+07 6.0106607628e+09 0.7000000000 -5.8084293227 -11.9132714976 1 + 1.0000000000e+07 6.6548423840e+09 0.7000000000 -6.5184076951 -11.5627682932 1 + 1.0000000000e+07 7.3680629973e+09 0.7000000000 -7.3127256426 -11.2391950016 1 + 1.0000000000e+07 8.1577217310e+09 0.7000000000 -8.1900933044 -10.8968169016 1 + 1.0000000000e+07 9.0320107014e+09 0.7000000000 -9.0920365426 -10.6614850215 1 + 1.0000000000e+07 1.0000000000e+10 0.7000000000 -9.9582648682 -10.8798139159 1 + 1.0000000000e+08 4.0000000000e+09 0.7000000000 -2.8617629403 -14.6856624315 0 + 1.0000000000e+08 4.4286927156e+09 0.7000000000 -3.2266384172 -14.3451718155 1 + 1.0000000000e+08 4.9033297922e+09 0.7000000000 -3.6368702837 -13.9626100890 1 + 1.0000000000e+08 5.4288352332e+09 0.7000000000 -4.0979886993 -13.5334235631 1 + 1.0000000000e+08 6.0106607628e+09 0.7000000000 -4.6161879732 -13.0535392989 1 + 1.0000000000e+08 6.6548423840e+09 0.7000000000 -5.1984011582 -12.5233515125 1 + 1.0000000000e+08 7.3680629973e+09 0.7000000000 -5.8523647926 -11.9934421813 1 + 1.0000000000e+08 8.1577217310e+09 0.7000000000 -6.5864090547 -11.6821331659 1 + 1.0000000000e+08 9.0320107014e+09 0.7000000000 -7.4067994755 -11.3654235744 1 + 1.0000000000e+08 1.0000000000e+10 0.7000000000 -8.3013506853 -11.0538132814 1 + 1.0000000000e+09 4.0000000000e+09 0.7000000000 -2.0695815696 -15.5177250154 0 + 1.0000000000e+09 4.4286927156e+09 0.7000000000 -2.3493240008 -15.2586547889 1 + 1.0000000000e+09 4.9033297922e+09 0.7000000000 -2.6653095613 -14.9660806276 1 + 1.0000000000e+09 5.4288352332e+09 0.7000000000 -3.0220912983 -14.6360018221 1 + 1.0000000000e+09 6.0106607628e+09 0.7000000000 -3.4247819832 -14.2640168398 1 + 1.0000000000e+09 6.6548423840e+09 0.7000000000 -3.8791217034 -13.8459596123 1 + 1.0000000000e+09 7.3680629973e+09 0.7000000000 -4.3915522296 -13.3782729765 1 + 1.0000000000e+09 8.1577217310e+09 0.7000000000 -4.9692869712 -12.8617482228 1 + 1.0000000000e+09 9.0320107014e+09 0.7000000000 -5.6203221431 -12.3289938297 1 + 1.0000000000e+09 1.0000000000e+10 0.7000000000 -6.3529920884 -11.9887960558 1 + 1.0000000000e+10 4.0000000000e+09 0.7000000000 -1.2805272580 -16.4326417920 0 + 1.0000000000e+10 4.4286927156e+09 0.7000000000 -1.4749982900 -16.2549301551 1 + 1.0000000000e+10 4.9033297922e+09 0.7000000000 -1.6966040866 -16.0523580418 1 + 1.0000000000e+10 5.4288352332e+09 0.7000000000 -1.9489205221 -15.8217642481 1 + 1.0000000000e+10 6.0106607628e+09 0.7000000000 -2.2359784782 -15.5595330834 1 + 1.0000000000e+10 6.6548423840e+09 0.7000000000 -2.5623203021 -15.2619989382 1 + 1.0000000000e+10 7.3680629973e+09 0.7000000000 -2.9330631536 -14.9250007236 1 + 1.0000000000e+10 8.1577217310e+09 0.7000000000 -3.3539699474 -14.5444727373 1 + 1.0000000000e+10 9.0320107014e+09 0.7000000000 -3.8315269899 -14.1169483246 1 + 1.0000000000e+10 1.0000000000e+10 0.7000000000 -4.3730130767 -13.6410243937 1 +AMReX (23.07-445-gc09c45c31007) finalized diff --git a/util/approx_math/approx_math.H b/util/approx_math/approx_math.H index 9f4a6a4835..97fa6b43b3 100644 --- a/util/approx_math/approx_math.H +++ b/util/approx_math/approx_math.H @@ -3,6 +3,7 @@ #include #include +#include #include #include #include @@ -38,6 +39,19 @@ amrex::Real fast_atan_1(const amrex::Real x) { constexpr amrex::Real A = 0.0776509570923569_rt; constexpr amrex::Real B = -0.287434475393028_rt; constexpr amrex::Real C = GCEM_PI / 4.0_rt - A - B; + constexpr amrex::Real third = 1.0_rt / 3.0_rt; + + /// + /// If x < 0.378, then using arctan(x) ~ x - x^3/3 + x^5/5 + /// gives better answer than the approximation below. + /// And accuracy increase as x << 0.378. + /// + + if (std::abs(x) < 0.378_rt) { + // return -0.3333333333_rt*amrex::Math::powi<3>(x) + x; + amrex::Real x2 = x*x; + return ((0.2_rt*x2 - third)*x2 + 1.0_rt)*x; + } amrex::Real x2 = x*x; return ((A*x2 + B)*x2 + C)*x; @@ -52,17 +66,6 @@ amrex::Real fast_atan(const amrex::Real x) { constexpr amrex::Real PI_2 = 0.5_rt * GCEM_PI; - /// - /// If x < 0.113, then using arctan(x) ~ x - /// gives better answer than the approximation below. - /// And accuracy increase as x << 0.113. - /// Also significantly faster. - /// - - if (std::abs(x) < 0.113_rt) { - return x; - } - // Check for large number, close to infinity. // Error is ~1e-8 rad by not checking actual inf