E-017 — Cor 1.6 thinning remainder R_tau = W/F - R_0 direct measurement
08-18-2026 · 320 rows · script
research/numerics/scripts/E-017-thinning-measurement.py
Metadata
Validation checks
| C1_W_rescaled_vs_mp | PASS |
| C2_R0_skeleton_dominant | FAIL |
Key values
| xi_flint_rel | 0.00000000000000008821833278620677 |
| C1_worst | 0.0000000000000038224510130899045 |
| report | {'4.0_pi': {'max_Rtau_over_R0': 1.0471206540091407, 'arg_M_deg': 90.0, 'C2_R0_vs_asymptotic': 0.03663281379835252, 'ratio_profile': {'0.5': 1.0471206540091407, '4.5': 1.0089222529454642, '8.5': 1.004385471663986, '12.5': 1.002691172170338, '16.5': 1.0018343882474554}}, '4.0_zero': {'max_Rtau_over_R0': 1.0679333353283942, 'arg_M_deg': -90.0, 'C2_R0_vs_asymptotic': 0.07326332357745524, 'ratio_profile': {'0.5': 1.0679333353283942, '4.5': 1.007167312755912, '8.5': 1.0034687273606262, '12.5': 1.0021290544442027, '16.5': 1.001455414041486}}, '6.0_pi': {'max_Rtau_over_R0': 1.0610941706671497, 'arg_M_deg': 90.0, 'C2_R0_vs_asymptotic': 0.004957508160826192, 'ratio_profile': {'0.5': 1.0610941706671497, '4.5': 1.0092001036041862, '8.5': 1.0045129740120455, '12.5': 1.0027792548957009, '16.5': 1.0019027560815716}}, '6.0_zero': {'max_Rtau_over_R0': 1.076103329487372, 'arg_M_deg': -90.0, 'C2_R0_vs_asymptotic': 0.009915010610412719, 'ratio_profile': {'0.5': 1.076103329487372, '4.5': 1.0080288607530927, '8.5': 1.0038984561227535, '12.5': 1.0023976261332865, '16.5': 1.0016410970456837}}, '8.0_pi': {'max_Rtau_over_R0': 1.068688355229517, 'arg_M_deg': 90.0, 'C2_R0_vs_asymptotic': 0.0006709252652428602, 'ratio_profile': {'0.5': 1.068688355229517, '4.5': 1.0093700235310568, '8.5': 1.0046033456882202, '12.5': 1.0028429914959156, '16.5': 1.0019507330410926}}, '8.0_zero': {'max_Rtau_over_R0': 1.0801095258799194, 'arg_M_deg': -90.0, 'C2_R0_vs_asymptotic': 0.0013418505163289656, 'ratio_profile': {'0.5': 1.0801095258799194, '4.5': 1.0084481599920176, '8.5': 1.0041028348136058, '12.5': 1.0025253381993795, '16.5': 1.0017307942059315}}, '10.0_pi': {'max_Rtau_over_R0': 1.0733197718087135, 'arg_M_deg': 90.0, 'C2_R0_vs_asymptotic': 9.079985954836377e-05, 'ratio_profile': {'0.5': 1.0733197718087135, '4.5': 1.0094431817705503, '8.5': 1.0046315068937992, '12.5': 1.0028604241200307, '16.5': 1.0019637968007054}}, '10.0_zero': {'max_Rtau_over_R0': 1.0824561061830875, 'arg_M_deg': -90.0, 'C2_R0_vs_asymptotic': 0.00018159971906163639, 'ratio_profile': {'0.5': 1.0824561061830875, '4.5': 1.008718091685042, '8.5': 1.0042436582330878, '12.5': 1.0026166536755394, '16.5': 1.001795570011337}}} |
Parameters
| a_sweep | [4.0, 6.0, 8.0, 10.0] |
| thetas | ['pi', 'zero'] |
| h | 0.01 |
| lambda | -1.0 |
| y_grid | [0.5, 20.0, 40] |
Chart
Data preview (320 rows total)
| a | theta | y | W_re | W_im | M_re | M_im | R0_re | R0_im | Rtau_re | Rtau_im | ratio_Rtau_R0 |
|---|---|---|---|---|---|---|---|---|---|---|---|
| 4.0 | pi | 0.500 | -3.0773101912e-14 | -2.0972636474e+02 | 2.8429064230e-15 | 1.9375116330e+01 | -1.4507441631e+01 | -4.1068003556e+02 | 1.4507441631e+01 | 4.3005515189e+02 | 1.047121e+00 |
| 4.0 | pi | 1.000 | -2.6894093546e-13 | -2.1960693944e+03 | 1.2407839651e-14 | 1.0131769960e+02 | -1.0915966879e+02 | -2.9809576516e+03 | 1.0915966879e+02 | 3.0822753512e+03 | 1.033944e+00 |
| 4.0 | pi | 1.500 | -2.0889110486e-12 | -1.7937827376e+04 | 6.4121428500e-14 | 5.5062139499e+02 | -8.0685262948e+02 | -2.2019212029e+04 | 8.0685262948e+02 | 2.2569833424e+04 | 1.024973e+00 |
| 4.0 | pi | 2.000 | -1.5775695189e-11 | -1.3851084596e+05 | 3.6218736875e-13 | 3.1800106583e+03 | -5.9619153032e+03 | -1.6270021158e+05 | 5.9619153032e+03 | 1.6588022224e+05 | 1.019519e+00 |
| 4.0 | pi | 2.500 | -1.1786859454e-10 | -1.0476652130e+06 | 2.1572912761e-12 | 1.9174904335e+04 | -4.4052931499e+04 | -1.2022008578e+06 | 4.4052931499e+04 | 1.2213757622e+06 | 1.015929e+00 |
| 4.0 | pi | 3.000 | -8.7608148532e-10 | -7.8467604737e+06 | 1.3305723895e-11 | 1.1917479148e+05 | -3.2550958283e+05 | -8.8831295629e+06 | 3.2550958283e+05 | 9.0023043543e+06 | 1.013398e+00 |
| 4.0 | pi | 3.500 | -6.4927064777e-09 | -5.8451710278e+07 | 8.4108795076e-11 | 7.5720393930e+05 | -2.4052085683e+06 | -6.5637942672e+07 | 2.4052085683e+06 | 6.6395146611e+07 | 1.011521e+00 |
| 4.0 | pi | 4.000 | -4.8032185920e-08 | -4.3397650129e+08 | 5.4143367684e-10 | 4.8919175394e+06 | -1.7772221041e+07 | -4.8500244062e+08 | 1.7772221041e+07 | 4.8989435816e+08 | 1.010073e+00 |
| 4.0 | pi | 4.500 | -3.5492264232e-07 | -3.2151328714e+09 | 3.5344543712e-09 | 3.2017513329e+07 | -1.3131993827e+08 | -3.5837102418e+09 | 1.3131993827e+08 | 3.6157277552e+09 | 1.008922e+00 |
| 4.0 | pi | 5.000 | -2.6205371947e-06 | -2.3784372292e+10 | 2.3329640985e-08 | 2.1174317531e+08 | -9.7033039082e+08 | -2.6480236019e+10 | 9.7033039082e+08 | 2.6691979195e+10 | 1.007986e+00 |
| 4.0 | pi | 5.500 | -1.9337656001e-05 | -1.7576444976e+11 | 1.5537939068e-07 | 1.4122793945e+09 | -7.1698256923e+09 | -1.9566394946e+11 | 7.1698256923e+09 | 1.9707622885e+11 | 1.007208e+00 |
| 4.0 | pi | 6.000 | -1.4264011322e-04 | -1.2978971582e+12 | 1.0425698095e-06 | 9.4864506370e+09 | -5.2978244260e+10 | -1.4457718991e+12 | 5.2978244260e+10 | 1.4552583497e+12 | 1.006553e+00 |
| 4.0 | pi | 6.500 | -1.0518416704e-03 | -9.5786479859e+12 | 7.0393384322e-06 | 6.4104082194e+10 | -3.9145921886e+11 | -1.0682889669e+13 | 3.9145921886e+11 | 1.0746993751e+13 | 1.005993e+00 |
| 4.0 | pi | 7.000 | -7.7546637428e-03 | -7.0661524572e+13 | 4.7783379584e-05 | 4.3540849257e+11 | -2.8925141286e+12 | -7.8936471061e+13 | 2.8925141286e+12 | 7.9371879553e+13 | 1.005509e+00 |
| 4.0 | pi | 7.500 | -5.7161489196e-02 | -5.2109926535e+14 | 3.2585693544e-04 | 2.9705980732e+12 | -2.1372949163e+13 | -5.8326601292e+14 | 2.1372949163e+13 | 5.8623661099e+14 | 1.005086e+00 |
| 4.0 | pi | 8.000 | -4.2129852204e-01 | -3.8419274772e+15 | 2.2311609714e-03 | 2.0346519614e+13 | -1.5792592037e+14 | -4.3097852901e+15 | 1.5792592037e+14 | 4.3301318097e+15 | 1.004715e+00 |
| 4.0 | pi | 8.500 | -3.1048145139e+00 | -2.8320033139e+16 | 1.5331487577e-02 | 1.3984353471e+14 | -1.1669234851e+15 | -3.1845245283e+16 | 1.1669234851e+15 | 3.1985088817e+16 | 1.004385e+00 |
| 4.0 | pi | 9.000 | -2.2879736434e+01 | -2.0872437280e+17 | 1.0568632743e-01 | 9.6414189340e+14 | -8.6224630942e+15 | -2.3530630388e+17 | 8.6224630942e+15 | 2.3627044577e+17 | 1.004092e+00 |
| 4.0 | pi | 9.500 | -1.6859472030e+02 | -1.5381620124e+18 | 7.3062443733e-01 | 6.6658004050e+15 | -6.3711863514e+16 | -1.7386914798e+18 | 6.3711863514e+16 | 1.7453572802e+18 | 1.003829e+00 |
| 4.0 | pi | 10.000 | -1.2422835057e+03 | -1.1334230408e+19 | 5.0639847028e+00 | 4.6202311419e+16 | -4.7077053367e+17 | -1.2847288883e+19 | 4.7077053367e+17 | 1.2893491194e+19 | 1.003591e+00 |
| 4.0 | pi | 10.500 | -9.1534775778e+03 | -8.3512628069e+19 | 3.5181390695e+01 | 3.2098078256e+17 | -3.4785498830e+18 | -9.4929338274e+19 | 3.4785498830e+18 | 9.5250319057e+19 | 1.003377e+00 |
| 4.0 | pi | 11.000 | -6.7444106271e+04 | -6.1530392630e+20 | 2.4494576839e+02 | 2.2346814474e+18 | -2.5703200229e+19 | -7.0143820594e+20 | 2.5703200229e+19 | 7.0367288739e+20 | 1.003182e+00 |
| 4.0 | pi | 11.500 | -4.9693248217e+05 | -4.5332582641e+21 | 1.7087900775e+03 | 1.5588408926e+19 | -1.8992238841e+20 | -5.1829662536e+21 | 1.8992238841e+20 | 5.1985546626e+21 | 1.003004e+00 |
| 4.0 | pi | 12.000 | -3.6614124684e+06 | -3.3397900013e+22 | 1.1942751066e+04 | 1.0893686778e+20 | -1.4033471824e+21 | -3.8297228407e+22 | 1.4033471824e+21 | 3.8406165275e+22 | 1.002841e+00 |
| 4.0 | pi | 12.500 | -2.6977376157e+07 | -2.4604784434e+23 | 8.3610324814e+04 | 7.6257009079e+20 | -1.0369411057e+22 | -2.8298036913e+23 | 1.0369411057e+22 | 2.8374293922e+23 | 1.002691e+00 |