From d44da6a6bd3db38b4d1e13f6ca968beeb51e2d41 Mon Sep 17 00:00:00 2001 From: Caleb Munger Date: Sat, 11 Jul 2026 17:03:16 -0500 Subject: [PATCH 1/2] make PSReadIC test for more IC coverage --- tests/PSReadIC/2d.xyzr | 72 ++++++++++++ tests/PSReadIC/input | 136 +++++++++++++++++++++++ tests/PSReadIC/reference/2d-parallel.dat | 109 ++++++++++++++++++ tests/PSReadIC/reference/2d-serial.dat | 109 ++++++++++++++++++ tests/PSReadIC/test | 27 +++++ 5 files changed, 453 insertions(+) create mode 100644 tests/PSReadIC/2d.xyzr create mode 100644 tests/PSReadIC/input create mode 100644 tests/PSReadIC/reference/2d-parallel.dat create mode 100644 tests/PSReadIC/reference/2d-serial.dat create mode 100755 tests/PSReadIC/test diff --git a/tests/PSReadIC/2d.xyzr b/tests/PSReadIC/2d.xyzr new file mode 100644 index 0000000000..8afb1072a7 --- /dev/null +++ b/tests/PSReadIC/2d.xyzr @@ -0,0 +1,72 @@ +0.000764027 0.000458356 0.000547547 9.157919480757395e-05 +0.000706756 0.000963831 0.000531243 5.7207217962404245e-05 +0.000144156 0.000105138 0.00052914 4.557273362100805e-05 +0.00048385 0.000200307 0.000542433 3.565844848912189e-05 +0.000365064 0.000315486 0.000582259 5.521603158904549e-05 +0.000176361 0.000989813 0.000515152 3.358780886213228e-05 +0.000702453 0.000610648 0.000540604 7.092845092280766e-05 +2.27987e-05 0.000990404 0.000577217 4.955587920186915e-05 +0.00040073 0.000907542 0.000571807 5.9589985477084446e-05 +0.000849761 2.32098e-05 0.000507006 6.878530906354948e-05 +0.000622208 0.00049455 0.000552457 5.44304758332052e-05 +0.000882236 0.000325443 0.000635867 3.80243553329561e-05 +0.00042026 0.00013373 0.000530357 3.68397104603657e-05 +0.00097258 0.000197985 0.000542534 7.88339056201301e-05 +0.000877619 0.000916387 0.000565639 3.784097512083696e-05 +0.000228351 0.000958798 0.000567998 2.0658690712068698e-05 +0.000695329 0.000716354 0.000533267 2.5292938146583816e-05 +7.71642e-05 0.000414946 0.000573195 5.3025923019796264e-05 +0.000583362 0.000596005 0.000606495 3.9406608163877306e-05 +0.000782716 0.000247894 0.000441496 8.270165993176396e-05 +0.00023922 1.95622e-05 0.00054558 1.934974117134453e-05 +0.000610581 0.000733289 0.00059792 6.093866457515767e-05 +0.00015809 0.000811974 0.000542924 5.703805926087346e-05 +0.000509939 0.000111422 0.000608132 5.428359461927924e-05 +0.000283844 0.000127578 0.000595859 4.17727048767922e-05 +0.000983321 0.00039492 0.00052771 5.982674719761284e-06 +0.000632833 0.0008414 0.000472655 4.059336990125662e-05 +0.000251781 0.000396611 0.000546737 6.218654093615097e-05 +0.000911487 0.000436719 0.000546635 5.603933765302946e-05 +0.000379316 3.60579e-05 0.000516485 4.8087646660969464e-05 +4.8336e-06 0.000475747 0.000538566 4.1403269970796874e-05 +0.000503531 0.000728128 0.000542115 4.463550534493841e-05 +0.000557648 0.000229484 0.000538757 3.788750982061122e-05 +9.26847e-05 0.000177594 0.000569166 4.2680715609215044e-05 +0.000339821 0.00047318 0.000570772 5.422516414506402e-05 +0.000105376 1.81598e-05 0.000572021 1.939862445141182e-05 +0.000586606 0.000920505 0.000674707 5.085406362133194e-05 +0.000348653 0.000241817 0.000544105 1.967585089669914e-05 +0.000362486 0.000397412 0.000562888 1.9792653948097408e-05 +0.000444208 0.00060632 0.000549128 7.518244545417664e-05 +0.000682256 8.2097e-05 0.000594876 6.30707145225211e-05 +0.000997479 0.000831419 0.00057031 0.00010459649034807909 +0.000504422 0.000472898 0.000489653 4.6592666391793035e-05 +0.000512838 0.00039833 0.000563108 1.902505782917835e-05 +0.000336609 0.00076567 0.000591 7.836551174001665e-05 +0.000437898 0.000413714 0.000601825 4.202559991258272e-05 +0.000293909 0.000985566 0.000560448 4.404782806277768e-05 +0.000444061 0.000807689 0.0004977 2.282329863358987e-05 +0.000986337 0.00040912 0.000587197 8.45418212656519e-06 +0.000785199 0.000805494 0.000599348 9.990220243107881e-05 +4.82455e-05 0.000109893 0.000516956 2.016649635801031e-06 +0.000839015 0.000139968 0.000651987 3.4051663855941284e-05 +0.000619205 0.000342904 0.00059458 9.115766017160855e-05 +0.00040966 0.00022367 0.000550016 3.9770253983183376e-05 +0.000198606 0.000600709 0.000448806 0.00013190870802449455 +0.000185071 0.000905208 0.000554471 3.9309641367060494e-05 +7.95143e-05 4.02304e-05 0.000529048 1.2603503884203527e-05 +0.000531178 0.000811227 0.000554845 4.184926820020237e-05 +0.000970769 0.000323016 0.000534723 3.967879509340839e-05 +0.000341543 0.000175547 0.00052579 3.323859800531813e-05 +0.000481826 1.32864e-06 0.000550954 5.142493452236119e-05 +0.000911282 0.00061958 0.000525757 0.00012376734021058412 +0.00059811 0.00014595 0.000520425 3.990573420397158e-05 +0.000424591 0.000712924 0.000526166 1.9630960360026034e-05 +0.000393845 0.000978356 0.000558189 1.0342702531158138e-05 +0.000962533 7.93134e-05 0.00058867 3.449657693707887e-05 +0.000422701 0.000504412 0.000569123 2.666703679164433e-05 +0.000474654 0.000309178 0.000484258 4.769952210077653e-05 +0.000223923 0.000246501 0.000529212 8.69647349017483e-05 +0.000113959 0.000947068 0.000583195 3.805190998218708e-05 +0.000490793 0.000868199 0.000559346 2.6147893612174206e-05 +7.9948e-05 0.000292751 0.000572211 6.422356688917259e-05 diff --git a/tests/PSReadIC/input b/tests/PSReadIC/input new file mode 100644 index 0000000000..912fe354a5 --- /dev/null +++ b/tests/PSReadIC/input @@ -0,0 +1,136 @@ +#@ +#@ [2d-serial-short] +#@ dim=2 +#@ nprocs=1 +#@ args = stop_time=1e-8 +#@ check-file = reference/2d-serial.dat +#@ +#@ [2d-parallel-long] +#@ dim=2 +#@ nprocs=2 +#@ args = stop_time=1e-8 +#@ check-file = reference/2d-parallel.dat +#@ + +alamo.program = flame +plot_file = tests/PSReadIC/output + +# AMR Parameters +amr.plot_int = 1 +amr.blocking_factor = 4 +amr.grid_eff = 1 +amr.refinement_criterion = 0.1 +amr.phi_refinement_criterion = 0.5 +amr.cell.all = 1 +amr.node.all = 0 + + +# Geometry and AMR level specification +amr.n_cell = 8 8 8 +geometry.prob_lo = 0.0000_m 0.0000_m 0.0000_m +geometry.prob_hi = 0.0005_m 0.0005_m 0.0005_m +amr.max_level = 4 +geometry.is_periodic = 0 0 0 + +# Timestep and duration +timestep = 1e-8 +stop_time = 1e-10 + +# Phase field parameters +pf.eps = 0.000001_m +pf.lambda = 0.001_J/m^2 +pf.eta.ic.type = psread +pf.eta.ic.psread.file.unit = m +pf.eta.ic.psread.file.name = tests/PSReadIC/2d.xyzr +pf.eta.ic.psread.eps = 2.0e-5_m +pf.eta.ic.psread.invert = 1 +pf.eta.ic.psread.x0 = 1e-4 1e-4 0.0 +pf.eta.ic.psread.mult = 1.1 + +# Eta boundary condition +pf.eta.bc.type = constant +pf.eta.bc.constant.type.xlo = neumann +pf.eta.bc.constant.type.xhi = neumann +pf.eta.bc.constant.type.ylo = neumann +pf.eta.bc.constant.type.yhi = neumann +pf.eta.bc.constant.type.zlo = neumann +pf.eta.bc.constant.type.zhi = neumann +pf.eta.bc.constant.val.xlo = 0.0 +pf.eta.bc.constant.val.xhi = 1.0 +pf.eta.bc.constant.val.ylo = 0.0 +pf.eta.bc.constant.val.yhi = 0.0 +pf.eta.bc.constant.val.zlo = 0.0 +pf.eta.bc.constant.val.zhi = 0.0 + +# Phi initial condition (spheres) +phi.ic.type = psread +phi.ic.psread.file.unit = m +phi.ic.psread.file.name = tests/PSReadIC/2d.xyzr +phi.ic.psread.eps = 1.13e-5_m + +# Use the full-feedback propellant model +propellant.type = fullfeedback +# Arrhenius rate law params +propellant.fullfeedback.m_ap = 1.45e5_m/s +propellant.fullfeedback.m_htpb = 1.04e2_m/s +propellant.fullfeedback.E_ap = 11000.0_K +propellant.fullfeedback.E_htpb = 7500.0_K +propellant.fullfeedback.mob_ap = true +propellant.fullfeedback.bound = 300.0_K +# Thermal transport params +propellant.fullfeedback.rho_ap = 1950.0_kg/m^3 +propellant.fullfeedback.rho_htpb = 920.0_kg/m^3 +propellant.fullfeedback.k_ap = 0.4186e0_W/m/K +propellant.fullfeedback.k_htpb = 0.13_W/m/K +propellant.fullfeedback.cp_ap = 1297.9_J/kg/K +propellant.fullfeedback.cp_htpb = 2418.29_J/kg/K +# Mass-to-heat flux estimation params +propellant.fullfeedback.a1 = 1.11400_1 +propellant.fullfeedback.a2 = 0.46000_1 +propellant.fullfeedback.a3 = 2.79700_1 +propellant.fullfeedback.b1 = 0.32300_1 +propellant.fullfeedback.b2 = 0.42000_1 +propellant.fullfeedback.b3 = 0.32250_1 +propellant.fullfeedback.c1 = -0.09906_1 +propellant.fullfeedback.mlocal_ap = 1000.0_kg/m^3/s +propellant.fullfeedback.mlocal_comb = 0.0_kg/m^3/s +propellant.fullfeedback.phi.zeta_0 = 1.13e-5_m +propellant.fullfeedback.phi.zeta = 1.13e-5 _m +propellant.fullfeedback.Pref = 1.0_MPa + + +# Activate thermal transport +thermal.on = 1 + +# Temperature boundary conditions +thermal.temp.bc.type = constant +thermal.temp.bc.constant.type.xlo = neumann +thermal.temp.bc.constant.type.xhi = neumann +thermal.temp.bc.constant.type.ylo = neumann +thermal.temp.bc.constant.type.yhi = neumann +thermal.temp.bc.constant.type.zlo = neumann +thermal.temp.bc.constant.type.zhi = neumann +thermal.temp.bc.constant.val.xlo = 0.0 +thermal.temp.bc.constant.val.xhi = 0.0 +thermal.temp.bc.constant.val.ylo = 0.0 +thermal.temp.bc.constant.val.yhi = 0.0 +thermal.temp.bc.constant.val.zlo = 0.0 +thermal.temp.bc.constant.val.zhi = 0.0 + +# Temperature initial conditions +temp.ic.type="constant" +temp.ic.constant.value=500_K + +# Laser source term +laser.ic.type = constant +laser.ic.constant.value = 0.0_W/m^2 +thermal.hc = 1.0e7_W/m^2 + +# Pressure +chamber.pressure = 4.0_MPa # [MPa] + +elastic.type = disable + + + + diff --git a/tests/PSReadIC/reference/2d-parallel.dat b/tests/PSReadIC/reference/2d-parallel.dat new file mode 100644 index 0000000000..87a1591685 --- /dev/null +++ b/tests/PSReadIC/reference/2d-parallel.dat @@ -0,0 +1,109 @@ +,x,y,z,eta,phi +0,3.90625e-06,0.00043359375,0.5,1.0,0.4582811715952242 +1,1.171875e-05,0.00043359375,0.5,1.0,0.4533531771300272 +2,1.953125e-05,0.00043359375,0.5,1.0,0.48465814735371426 +3,2.734375e-05,0.00043359375,0.5,1.0,0.6187702232104019 +4,3.5156250000000004e-05,0.00043359375,0.5,1.0,0.8099244424445855 +5,4.296875e-05,0.00043359375,0.5,1.0,0.9445592890653867 +6,5.0781250000000004e-05,0.00043359375,0.5,1.0,0.9899106690231667 +7,5.859375e-05,0.00043359375,0.5,1.0,0.9987226708136097 +8,6.640625e-05,0.00043359375000000006,0.5,1.0,0.9998294728446361 +9,7.421875e-05,0.00043359375000000006,0.5,1.0,0.9999535186180926 +10,8.0078125e-05,0.00043164062500000006,0.5,1.0,0.9999943161163413 +11,8.398437500000001e-05,0.00043164062500000006,0.5,1.0,0.9999895642889608 +12,8.7890625e-05,0.00043164062500000006,0.5,1.0,0.9999725150198857 +13,9.1796875e-05,0.00043164062500000006,0.5,1.0,0.9999099977284864 +14,9.5703125e-05,0.00043164062500000006,0.5,1.0,0.9996756221836025 +15,9.960937500000001e-05,0.00043164062500000006,0.5,1.0,0.9988238561477296 +16,0.00010351562500000001,0.00043164062500000006,0.5,1.0,0.9959673688757547 +17,0.000107421875,0.00043164062500000006,0.5,0.8835343699143661,0.9874509381755153 +18,0.000111328125,0.00043164062500000006,0.5,0.8208671646157901,0.9654712883579316 +19,0.000115234375,0.00043164062500000006,0.5,0.740447256770379,0.9172788186274969 +20,0.000119140625,0.00043164062500000006,0.5,0.6446918295319187,0.8286858464419236 +21,0.000123046875,0.00043164062500000006,0.5,0.5388871141854548,0.6934180757154493 +22,0.000126953125,0.00043164062500000006,0.5,0.4303774999022171,0.5230522787196434 +23,0.000130859375,0.00043164062500000006,0.5,0.327065287188968,0.3469737918465656 +24,0.000134765625,0.00043164062500000006,0.5,0.23572152800605672,0.19825339418877008 +25,0.00013867187500000001,0.00043164062500000006,0.5,0.16069643584910842,0.06722993939132167 +26,0.00014257812500000002,0.00043164062500000006,0.5,0.10342525048783302,0.0 +27,0.00014648437500000002,0.00043164062500000006,0.5,0.06276864574204888,0.0 +28,0.000150390625,0.00043164062500000006,0.5,0.035906463899466834,0.0 +29,0.000154296875,0.00043164062500000006,0.5,0.019370086331836456,0.0 +30,0.00016015625,0.00043359375000000006,0.5,0.007910337399139489,0.0 +31,0.00016796875,0.00043359375000000006,0.5,0.0018120266525033513,0.0 +32,0.00017578125,0.00043359375000000006,0.5,0.00039759964545038606,0.0 +33,0.00018359375,0.00043359375000000006,0.5,0.00012249730266478132,0.025747918070621933 +34,0.00019140625,0.00043359375000000006,0.5,0.00011325211219714859,0.15443901372791802 +35,0.00019921875000000001,0.00043359375000000006,0.5,0.0003381402778630527,0.40233391751008196 +36,0.00020703125,0.00043359375000000006,0.5,0.0015231130960052797,0.672315709249885 +37,0.00021484375,0.00043359375000000006,0.5,0.0067424768577353245,0.8634516875769216 +38,0.000220703125,0.00043164062500000006,0.5,0.016706464664568266,0.9687428716690855 +39,0.000224609375,0.00043164062500000006,0.5,0.031377608252043276,0.9848584579809165 +40,0.000228515625,0.00043164062500000006,0.5,0.05559037927896038,0.9927392119335315 +41,0.00023242187500000002,0.00043164062500000006,0.5,0.09283330995526129,0.9964418101299328 +42,0.000236328125,0.00043164062500000006,0.5,0.14616106890974956,0.998153390257225 +43,0.000240234375,0.00043164062500000006,0.5,0.21718402864992847,0.9989463945822822 +44,0.000244140625,0.00043164062500000006,0.5,0.305106466416148,0.9993147640037052 +45,0.000248046875,0.00043164062500000006,0.5,0.4062277011976442,0.9994761631709994 +46,0.000251953125,0.00043164062500000006,0.5,0.514236973096488,0.9995194778817386 +47,0.000255859375,0.00043164062500000006,0.5,0.6213458666453286,0.9994674167066254 +48,0.000259765625,0.00043164062500000006,0.5,0.6523439128133405,0.9992926164031527 +49,0.000263671875,0.00043164062500000006,0.5,0.6937417449035357,0.998898237327045 +50,0.000267578125,0.00043164062500000006,0.5,0.7039234649993117,0.9980498737388569 +51,0.000271484375,0.00043164062500000006,0.5,0.6841301611019979,0.9962176589430582 +52,0.00027539062500000003,0.00043164062500000006,0.5,0.6707762212007404,0.9922562938388506 +53,0.00027929687500000003,0.00043164062500000006,0.5,0.5926683111782173,0.9838461985307727 +54,0.00028320312500000003,0.00043164062500000006,0.5,0.512820071885169,0.9667287443991215 +55,0.00028710937500000004,0.00043164062500000006,0.5,0.43477946748927565,0.9370021576890688 +56,0.00029101562500000004,0.00043164062500000006,0.5,0.3616458842026158,0.8926582739231124 +57,0.00029492187500000005,0.00043164062500000006,0.5,0.2957269026442766,0.8347080272160838 +58,0.00029882812500000005,0.00043164062500000006,0.5,0.23838775313694593,0.7745539452947771 +59,0.00030273437500000005,0.00043164062500000006,0.5,0.19008567696192225,0.7315880621219365 +60,0.00030664062500000006,0.00043164062500000006,0.5,0.15053738366430325,0.720084201568417 +61,0.00031054687500000006,0.00043164062500000006,0.5,0.11895049545853575,0.7408545427050972 +62,0.00031445312500000006,0.00043164062500000006,0.5,0.09425678677264238,0.7797203886904559 +63,0.000318359375,0.00043164062500000006,0.5,0.075305904320766,0.8232279314546171 +64,0.000322265625,0.00043164062500000006,0.5,0.06100181346413858,0.8652624889770042 +65,0.000326171875,0.00043164062500000006,0.5,0.05038274537594123,0.8981342938629702 +66,0.000330078125,0.00043164062500000006,0.5,0.042656032003524746,0.9191696364636631 +67,0.000333984375,0.00043164062500000006,0.5,0.03720270220531097,0.9314811888879476 +68,0.00033789062500000003,0.00043164062500000006,0.5,0.033565555548943915,0.937122601176575 +69,0.00034179687500000003,0.00043164062500000006,0.5,0.03143118716752369,0.9370905002307901 +70,0.00034570312500000004,0.00043164062500000006,0.5,0.03061291751220796,0.9313791198983925 +71,0.00034960937500000004,0.00043164062500000006,0.5,0.03103869931202019,0.9189798823410011 +72,0.00035351562500000004,0.00043164062500000006,0.5,0.032746052998986364,0.897826535687336 +73,0.00035742187500000005,0.00043164062500000006,0.5,0.035884732215552884,0.8647948813507121 +74,0.00036132812500000005,0.00043164062500000006,0.5,0.04072676073536963,0.8159879414985258 +75,0.00036523437500000005,0.00043164062500000006,0.5,0.04768226911073992,0.7476192072705634 +76,0.000369140625,0.00043164062500000006,0.5,0.05731778166048307,0.6576883711775607 +77,0.000373046875,0.00043164062500000006,0.5,0.07037099334506203,0.5481426878831556 +78,0.000376953125,0.00043164062500000006,0.5,0.0877526963275912,0.42642257545575435 +79,0.000380859375,0.00043164062500000006,0.5,0.1105230353942922,0.3048592220905742 +80,0.000384765625,0.00043164062500000006,0.5,0.1398272117768098,0.19722116560831943 +81,0.000388671875,0.00043164062500000006,0.5,0.17677749010723032,0.1678610307161649 +82,0.000392578125,0.00043164062500000006,0.5,0.22227653795347058,0.2282379134165803 +83,0.000396484375,0.00043164062500000006,0.5,0.2767932751785763,0.3512146614648787 +84,0.00040039062500000003,0.00043164062500000006,0.5,0.3401249927568327,0.5197524683186292 +85,0.00040429687500000003,0.00043164062500000006,0.5,0.4112023248556993,0.681807525502341 +86,0.00040820312500000004,0.00043164062500000006,0.5,0.4880062268070724,0.811976863203694 +87,0.00041210937500000004,0.00043164062500000006,0.5,0.5676571744429646,0.9003208354367438 +88,0.00041601562500000004,0.00043164062500000006,0.5,0.6467006915761834,0.9517714808776861 +89,0.00041992187500000005,0.00043164062500000006,0.5,0.7215561004276289,0.978031920424588 +90,0.00042382812500000005,0.00043164062500000006,0.5,0.7890362461056692,0.9901089469743323 +91,0.000427734375,0.00043164062500000006,0.5,0.8468105962133419,0.9952790268496151 +92,0.000431640625,0.00043164062500000006,0.5,0.7835297915920545,0.997386420875755 +93,0.000435546875,0.00043164062500000006,0.5,0.813112582945398,0.9981528368536158 +94,0.000439453125,0.00043164062500000006,0.5,0.7326544224958522,0.9982138373458963 +95,0.000443359375,0.00043164062500000006,0.5,0.6380592020564385,0.997619584311696 +96,0.000447265625,0.00043164062500000006,0.5,0.5346112883878035,0.9958781596259787 +97,0.000451171875,0.00043164062500000006,0.5,0.42931656875456514,0.9915555696939881 +98,0.000455078125,0.00043164062500000006,0.5,0.32948353810703557,0.9813476003025368 +99,0.00045898437500000003,0.00043164062500000006,0.5,0.24122230431057734,0.9587643157545075 +100,0.00046289062500000003,0.00043164062500000006,0.5,0.16836626219655204,0.917917532292608 +101,0.00046679687500000004,0.00043164062500000006,0.5,0.11211884629916469,0.856387492723463 +102,0.00047070312500000004,0.00043164062500000006,0.5,0.07141359412910775,0.7763375016584368 +103,0.00047460937500000004,0.00043164062500000006,0.5,0.04371994443225746,0.6943510972969217 +104,0.000478515625,0.00043164062500000006,0.5,0.0259394580603316,0.6354585405288947 +105,0.000482421875,0.00043164062500000006,0.5,0.015114926193632972,0.6139670453370898 +106,0.00048828125,0.00043359375000000006,0.5,0.005116696568355672,0.7004280271117741 +107,0.00049609375,0.00043359375000000006,0.5,0.0020004992962890933,0.7562051131588987 diff --git a/tests/PSReadIC/reference/2d-serial.dat b/tests/PSReadIC/reference/2d-serial.dat new file mode 100644 index 0000000000..87a1591685 --- /dev/null +++ b/tests/PSReadIC/reference/2d-serial.dat @@ -0,0 +1,109 @@ +,x,y,z,eta,phi +0,3.90625e-06,0.00043359375,0.5,1.0,0.4582811715952242 +1,1.171875e-05,0.00043359375,0.5,1.0,0.4533531771300272 +2,1.953125e-05,0.00043359375,0.5,1.0,0.48465814735371426 +3,2.734375e-05,0.00043359375,0.5,1.0,0.6187702232104019 +4,3.5156250000000004e-05,0.00043359375,0.5,1.0,0.8099244424445855 +5,4.296875e-05,0.00043359375,0.5,1.0,0.9445592890653867 +6,5.0781250000000004e-05,0.00043359375,0.5,1.0,0.9899106690231667 +7,5.859375e-05,0.00043359375,0.5,1.0,0.9987226708136097 +8,6.640625e-05,0.00043359375000000006,0.5,1.0,0.9998294728446361 +9,7.421875e-05,0.00043359375000000006,0.5,1.0,0.9999535186180926 +10,8.0078125e-05,0.00043164062500000006,0.5,1.0,0.9999943161163413 +11,8.398437500000001e-05,0.00043164062500000006,0.5,1.0,0.9999895642889608 +12,8.7890625e-05,0.00043164062500000006,0.5,1.0,0.9999725150198857 +13,9.1796875e-05,0.00043164062500000006,0.5,1.0,0.9999099977284864 +14,9.5703125e-05,0.00043164062500000006,0.5,1.0,0.9996756221836025 +15,9.960937500000001e-05,0.00043164062500000006,0.5,1.0,0.9988238561477296 +16,0.00010351562500000001,0.00043164062500000006,0.5,1.0,0.9959673688757547 +17,0.000107421875,0.00043164062500000006,0.5,0.8835343699143661,0.9874509381755153 +18,0.000111328125,0.00043164062500000006,0.5,0.8208671646157901,0.9654712883579316 +19,0.000115234375,0.00043164062500000006,0.5,0.740447256770379,0.9172788186274969 +20,0.000119140625,0.00043164062500000006,0.5,0.6446918295319187,0.8286858464419236 +21,0.000123046875,0.00043164062500000006,0.5,0.5388871141854548,0.6934180757154493 +22,0.000126953125,0.00043164062500000006,0.5,0.4303774999022171,0.5230522787196434 +23,0.000130859375,0.00043164062500000006,0.5,0.327065287188968,0.3469737918465656 +24,0.000134765625,0.00043164062500000006,0.5,0.23572152800605672,0.19825339418877008 +25,0.00013867187500000001,0.00043164062500000006,0.5,0.16069643584910842,0.06722993939132167 +26,0.00014257812500000002,0.00043164062500000006,0.5,0.10342525048783302,0.0 +27,0.00014648437500000002,0.00043164062500000006,0.5,0.06276864574204888,0.0 +28,0.000150390625,0.00043164062500000006,0.5,0.035906463899466834,0.0 +29,0.000154296875,0.00043164062500000006,0.5,0.019370086331836456,0.0 +30,0.00016015625,0.00043359375000000006,0.5,0.007910337399139489,0.0 +31,0.00016796875,0.00043359375000000006,0.5,0.0018120266525033513,0.0 +32,0.00017578125,0.00043359375000000006,0.5,0.00039759964545038606,0.0 +33,0.00018359375,0.00043359375000000006,0.5,0.00012249730266478132,0.025747918070621933 +34,0.00019140625,0.00043359375000000006,0.5,0.00011325211219714859,0.15443901372791802 +35,0.00019921875000000001,0.00043359375000000006,0.5,0.0003381402778630527,0.40233391751008196 +36,0.00020703125,0.00043359375000000006,0.5,0.0015231130960052797,0.672315709249885 +37,0.00021484375,0.00043359375000000006,0.5,0.0067424768577353245,0.8634516875769216 +38,0.000220703125,0.00043164062500000006,0.5,0.016706464664568266,0.9687428716690855 +39,0.000224609375,0.00043164062500000006,0.5,0.031377608252043276,0.9848584579809165 +40,0.000228515625,0.00043164062500000006,0.5,0.05559037927896038,0.9927392119335315 +41,0.00023242187500000002,0.00043164062500000006,0.5,0.09283330995526129,0.9964418101299328 +42,0.000236328125,0.00043164062500000006,0.5,0.14616106890974956,0.998153390257225 +43,0.000240234375,0.00043164062500000006,0.5,0.21718402864992847,0.9989463945822822 +44,0.000244140625,0.00043164062500000006,0.5,0.305106466416148,0.9993147640037052 +45,0.000248046875,0.00043164062500000006,0.5,0.4062277011976442,0.9994761631709994 +46,0.000251953125,0.00043164062500000006,0.5,0.514236973096488,0.9995194778817386 +47,0.000255859375,0.00043164062500000006,0.5,0.6213458666453286,0.9994674167066254 +48,0.000259765625,0.00043164062500000006,0.5,0.6523439128133405,0.9992926164031527 +49,0.000263671875,0.00043164062500000006,0.5,0.6937417449035357,0.998898237327045 +50,0.000267578125,0.00043164062500000006,0.5,0.7039234649993117,0.9980498737388569 +51,0.000271484375,0.00043164062500000006,0.5,0.6841301611019979,0.9962176589430582 +52,0.00027539062500000003,0.00043164062500000006,0.5,0.6707762212007404,0.9922562938388506 +53,0.00027929687500000003,0.00043164062500000006,0.5,0.5926683111782173,0.9838461985307727 +54,0.00028320312500000003,0.00043164062500000006,0.5,0.512820071885169,0.9667287443991215 +55,0.00028710937500000004,0.00043164062500000006,0.5,0.43477946748927565,0.9370021576890688 +56,0.00029101562500000004,0.00043164062500000006,0.5,0.3616458842026158,0.8926582739231124 +57,0.00029492187500000005,0.00043164062500000006,0.5,0.2957269026442766,0.8347080272160838 +58,0.00029882812500000005,0.00043164062500000006,0.5,0.23838775313694593,0.7745539452947771 +59,0.00030273437500000005,0.00043164062500000006,0.5,0.19008567696192225,0.7315880621219365 +60,0.00030664062500000006,0.00043164062500000006,0.5,0.15053738366430325,0.720084201568417 +61,0.00031054687500000006,0.00043164062500000006,0.5,0.11895049545853575,0.7408545427050972 +62,0.00031445312500000006,0.00043164062500000006,0.5,0.09425678677264238,0.7797203886904559 +63,0.000318359375,0.00043164062500000006,0.5,0.075305904320766,0.8232279314546171 +64,0.000322265625,0.00043164062500000006,0.5,0.06100181346413858,0.8652624889770042 +65,0.000326171875,0.00043164062500000006,0.5,0.05038274537594123,0.8981342938629702 +66,0.000330078125,0.00043164062500000006,0.5,0.042656032003524746,0.9191696364636631 +67,0.000333984375,0.00043164062500000006,0.5,0.03720270220531097,0.9314811888879476 +68,0.00033789062500000003,0.00043164062500000006,0.5,0.033565555548943915,0.937122601176575 +69,0.00034179687500000003,0.00043164062500000006,0.5,0.03143118716752369,0.9370905002307901 +70,0.00034570312500000004,0.00043164062500000006,0.5,0.03061291751220796,0.9313791198983925 +71,0.00034960937500000004,0.00043164062500000006,0.5,0.03103869931202019,0.9189798823410011 +72,0.00035351562500000004,0.00043164062500000006,0.5,0.032746052998986364,0.897826535687336 +73,0.00035742187500000005,0.00043164062500000006,0.5,0.035884732215552884,0.8647948813507121 +74,0.00036132812500000005,0.00043164062500000006,0.5,0.04072676073536963,0.8159879414985258 +75,0.00036523437500000005,0.00043164062500000006,0.5,0.04768226911073992,0.7476192072705634 +76,0.000369140625,0.00043164062500000006,0.5,0.05731778166048307,0.6576883711775607 +77,0.000373046875,0.00043164062500000006,0.5,0.07037099334506203,0.5481426878831556 +78,0.000376953125,0.00043164062500000006,0.5,0.0877526963275912,0.42642257545575435 +79,0.000380859375,0.00043164062500000006,0.5,0.1105230353942922,0.3048592220905742 +80,0.000384765625,0.00043164062500000006,0.5,0.1398272117768098,0.19722116560831943 +81,0.000388671875,0.00043164062500000006,0.5,0.17677749010723032,0.1678610307161649 +82,0.000392578125,0.00043164062500000006,0.5,0.22227653795347058,0.2282379134165803 +83,0.000396484375,0.00043164062500000006,0.5,0.2767932751785763,0.3512146614648787 +84,0.00040039062500000003,0.00043164062500000006,0.5,0.3401249927568327,0.5197524683186292 +85,0.00040429687500000003,0.00043164062500000006,0.5,0.4112023248556993,0.681807525502341 +86,0.00040820312500000004,0.00043164062500000006,0.5,0.4880062268070724,0.811976863203694 +87,0.00041210937500000004,0.00043164062500000006,0.5,0.5676571744429646,0.9003208354367438 +88,0.00041601562500000004,0.00043164062500000006,0.5,0.6467006915761834,0.9517714808776861 +89,0.00041992187500000005,0.00043164062500000006,0.5,0.7215561004276289,0.978031920424588 +90,0.00042382812500000005,0.00043164062500000006,0.5,0.7890362461056692,0.9901089469743323 +91,0.000427734375,0.00043164062500000006,0.5,0.8468105962133419,0.9952790268496151 +92,0.000431640625,0.00043164062500000006,0.5,0.7835297915920545,0.997386420875755 +93,0.000435546875,0.00043164062500000006,0.5,0.813112582945398,0.9981528368536158 +94,0.000439453125,0.00043164062500000006,0.5,0.7326544224958522,0.9982138373458963 +95,0.000443359375,0.00043164062500000006,0.5,0.6380592020564385,0.997619584311696 +96,0.000447265625,0.00043164062500000006,0.5,0.5346112883878035,0.9958781596259787 +97,0.000451171875,0.00043164062500000006,0.5,0.42931656875456514,0.9915555696939881 +98,0.000455078125,0.00043164062500000006,0.5,0.32948353810703557,0.9813476003025368 +99,0.00045898437500000003,0.00043164062500000006,0.5,0.24122230431057734,0.9587643157545075 +100,0.00046289062500000003,0.00043164062500000006,0.5,0.16836626219655204,0.917917532292608 +101,0.00046679687500000004,0.00043164062500000006,0.5,0.11211884629916469,0.856387492723463 +102,0.00047070312500000004,0.00043164062500000006,0.5,0.07141359412910775,0.7763375016584368 +103,0.00047460937500000004,0.00043164062500000006,0.5,0.04371994443225746,0.6943510972969217 +104,0.000478515625,0.00043164062500000006,0.5,0.0259394580603316,0.6354585405288947 +105,0.000482421875,0.00043164062500000006,0.5,0.015114926193632972,0.6139670453370898 +106,0.00048828125,0.00043359375000000006,0.5,0.005116696568355672,0.7004280271117741 +107,0.00049609375,0.00043359375000000006,0.5,0.0020004992962890933,0.7562051131588987 diff --git a/tests/PSReadIC/test b/tests/PSReadIC/test new file mode 100755 index 0000000000..80ca0e2591 --- /dev/null +++ b/tests/PSReadIC/test @@ -0,0 +1,27 @@ +#!/usr/bin/env python3 +import sys +sys.path.insert(0,"../../scripts") +import testlib +import glob + +outdir = sys.argv[1] +reference = sys.argv[2] + +path = sorted(glob.glob("{}/*cell/".format(outdir)))[-1] + +testlib.validate(path=path, + outdir=outdir, + start=[0.0000,0.00043,0], + end=[0.0005,0.00043,0], + vars=["eta","phi"], + tolerance=[0.01,0.01], + abs_tolerance=[0.01,0.01], + reference=reference, + generate_ref_data = False, +) +exit(0) + + + + + From 8485b14b61b58e65864e41d99d2687f3e3605f3e Mon Sep 17 00:00:00 2001 From: Caleb Munger Date: Mon, 13 Jul 2026 10:50:59 -0500 Subject: [PATCH 2/2] move psread test to SCPSpheresElastic --- tests/PSReadIC/2d.xyzr | 72 ---------- tests/PSReadIC/input | 136 ------------------ tests/PSReadIC/reference/2d-parallel.dat | 109 -------------- tests/PSReadIC/reference/2d-serial.dat | 109 -------------- tests/PSReadIC/test | 27 ---- tests/SCPSpheresElastic/input | 17 ++- .../reference/2d-psread-coverage.dat | 109 ++++++++++++++ 7 files changed, 125 insertions(+), 454 deletions(-) delete mode 100644 tests/PSReadIC/2d.xyzr delete mode 100644 tests/PSReadIC/input delete mode 100644 tests/PSReadIC/reference/2d-parallel.dat delete mode 100644 tests/PSReadIC/reference/2d-serial.dat delete mode 100755 tests/PSReadIC/test create mode 100644 tests/SCPSpheresElastic/reference/2d-psread-coverage.dat diff --git a/tests/PSReadIC/2d.xyzr b/tests/PSReadIC/2d.xyzr deleted file mode 100644 index 8afb1072a7..0000000000 --- a/tests/PSReadIC/2d.xyzr +++ /dev/null @@ -1,72 +0,0 @@ -0.000764027 0.000458356 0.000547547 9.157919480757395e-05 -0.000706756 0.000963831 0.000531243 5.7207217962404245e-05 -0.000144156 0.000105138 0.00052914 4.557273362100805e-05 -0.00048385 0.000200307 0.000542433 3.565844848912189e-05 -0.000365064 0.000315486 0.000582259 5.521603158904549e-05 -0.000176361 0.000989813 0.000515152 3.358780886213228e-05 -0.000702453 0.000610648 0.000540604 7.092845092280766e-05 -2.27987e-05 0.000990404 0.000577217 4.955587920186915e-05 -0.00040073 0.000907542 0.000571807 5.9589985477084446e-05 -0.000849761 2.32098e-05 0.000507006 6.878530906354948e-05 -0.000622208 0.00049455 0.000552457 5.44304758332052e-05 -0.000882236 0.000325443 0.000635867 3.80243553329561e-05 -0.00042026 0.00013373 0.000530357 3.68397104603657e-05 -0.00097258 0.000197985 0.000542534 7.88339056201301e-05 -0.000877619 0.000916387 0.000565639 3.784097512083696e-05 -0.000228351 0.000958798 0.000567998 2.0658690712068698e-05 -0.000695329 0.000716354 0.000533267 2.5292938146583816e-05 -7.71642e-05 0.000414946 0.000573195 5.3025923019796264e-05 -0.000583362 0.000596005 0.000606495 3.9406608163877306e-05 -0.000782716 0.000247894 0.000441496 8.270165993176396e-05 -0.00023922 1.95622e-05 0.00054558 1.934974117134453e-05 -0.000610581 0.000733289 0.00059792 6.093866457515767e-05 -0.00015809 0.000811974 0.000542924 5.703805926087346e-05 -0.000509939 0.000111422 0.000608132 5.428359461927924e-05 -0.000283844 0.000127578 0.000595859 4.17727048767922e-05 -0.000983321 0.00039492 0.00052771 5.982674719761284e-06 -0.000632833 0.0008414 0.000472655 4.059336990125662e-05 -0.000251781 0.000396611 0.000546737 6.218654093615097e-05 -0.000911487 0.000436719 0.000546635 5.603933765302946e-05 -0.000379316 3.60579e-05 0.000516485 4.8087646660969464e-05 -4.8336e-06 0.000475747 0.000538566 4.1403269970796874e-05 -0.000503531 0.000728128 0.000542115 4.463550534493841e-05 -0.000557648 0.000229484 0.000538757 3.788750982061122e-05 -9.26847e-05 0.000177594 0.000569166 4.2680715609215044e-05 -0.000339821 0.00047318 0.000570772 5.422516414506402e-05 -0.000105376 1.81598e-05 0.000572021 1.939862445141182e-05 -0.000586606 0.000920505 0.000674707 5.085406362133194e-05 -0.000348653 0.000241817 0.000544105 1.967585089669914e-05 -0.000362486 0.000397412 0.000562888 1.9792653948097408e-05 -0.000444208 0.00060632 0.000549128 7.518244545417664e-05 -0.000682256 8.2097e-05 0.000594876 6.30707145225211e-05 -0.000997479 0.000831419 0.00057031 0.00010459649034807909 -0.000504422 0.000472898 0.000489653 4.6592666391793035e-05 -0.000512838 0.00039833 0.000563108 1.902505782917835e-05 -0.000336609 0.00076567 0.000591 7.836551174001665e-05 -0.000437898 0.000413714 0.000601825 4.202559991258272e-05 -0.000293909 0.000985566 0.000560448 4.404782806277768e-05 -0.000444061 0.000807689 0.0004977 2.282329863358987e-05 -0.000986337 0.00040912 0.000587197 8.45418212656519e-06 -0.000785199 0.000805494 0.000599348 9.990220243107881e-05 -4.82455e-05 0.000109893 0.000516956 2.016649635801031e-06 -0.000839015 0.000139968 0.000651987 3.4051663855941284e-05 -0.000619205 0.000342904 0.00059458 9.115766017160855e-05 -0.00040966 0.00022367 0.000550016 3.9770253983183376e-05 -0.000198606 0.000600709 0.000448806 0.00013190870802449455 -0.000185071 0.000905208 0.000554471 3.9309641367060494e-05 -7.95143e-05 4.02304e-05 0.000529048 1.2603503884203527e-05 -0.000531178 0.000811227 0.000554845 4.184926820020237e-05 -0.000970769 0.000323016 0.000534723 3.967879509340839e-05 -0.000341543 0.000175547 0.00052579 3.323859800531813e-05 -0.000481826 1.32864e-06 0.000550954 5.142493452236119e-05 -0.000911282 0.00061958 0.000525757 0.00012376734021058412 -0.00059811 0.00014595 0.000520425 3.990573420397158e-05 -0.000424591 0.000712924 0.000526166 1.9630960360026034e-05 -0.000393845 0.000978356 0.000558189 1.0342702531158138e-05 -0.000962533 7.93134e-05 0.00058867 3.449657693707887e-05 -0.000422701 0.000504412 0.000569123 2.666703679164433e-05 -0.000474654 0.000309178 0.000484258 4.769952210077653e-05 -0.000223923 0.000246501 0.000529212 8.69647349017483e-05 -0.000113959 0.000947068 0.000583195 3.805190998218708e-05 -0.000490793 0.000868199 0.000559346 2.6147893612174206e-05 -7.9948e-05 0.000292751 0.000572211 6.422356688917259e-05 diff --git a/tests/PSReadIC/input b/tests/PSReadIC/input deleted file mode 100644 index 912fe354a5..0000000000 --- a/tests/PSReadIC/input +++ /dev/null @@ -1,136 +0,0 @@ -#@ -#@ [2d-serial-short] -#@ dim=2 -#@ nprocs=1 -#@ args = stop_time=1e-8 -#@ check-file = reference/2d-serial.dat -#@ -#@ [2d-parallel-long] -#@ dim=2 -#@ nprocs=2 -#@ args = stop_time=1e-8 -#@ check-file = reference/2d-parallel.dat -#@ - -alamo.program = flame -plot_file = tests/PSReadIC/output - -# AMR Parameters -amr.plot_int = 1 -amr.blocking_factor = 4 -amr.grid_eff = 1 -amr.refinement_criterion = 0.1 -amr.phi_refinement_criterion = 0.5 -amr.cell.all = 1 -amr.node.all = 0 - - -# Geometry and AMR level specification -amr.n_cell = 8 8 8 -geometry.prob_lo = 0.0000_m 0.0000_m 0.0000_m -geometry.prob_hi = 0.0005_m 0.0005_m 0.0005_m -amr.max_level = 4 -geometry.is_periodic = 0 0 0 - -# Timestep and duration -timestep = 1e-8 -stop_time = 1e-10 - -# Phase field parameters -pf.eps = 0.000001_m -pf.lambda = 0.001_J/m^2 -pf.eta.ic.type = psread -pf.eta.ic.psread.file.unit = m -pf.eta.ic.psread.file.name = tests/PSReadIC/2d.xyzr -pf.eta.ic.psread.eps = 2.0e-5_m -pf.eta.ic.psread.invert = 1 -pf.eta.ic.psread.x0 = 1e-4 1e-4 0.0 -pf.eta.ic.psread.mult = 1.1 - -# Eta boundary condition -pf.eta.bc.type = constant -pf.eta.bc.constant.type.xlo = neumann -pf.eta.bc.constant.type.xhi = neumann -pf.eta.bc.constant.type.ylo = neumann -pf.eta.bc.constant.type.yhi = neumann -pf.eta.bc.constant.type.zlo = neumann -pf.eta.bc.constant.type.zhi = neumann -pf.eta.bc.constant.val.xlo = 0.0 -pf.eta.bc.constant.val.xhi = 1.0 -pf.eta.bc.constant.val.ylo = 0.0 -pf.eta.bc.constant.val.yhi = 0.0 -pf.eta.bc.constant.val.zlo = 0.0 -pf.eta.bc.constant.val.zhi = 0.0 - -# Phi initial condition (spheres) -phi.ic.type = psread -phi.ic.psread.file.unit = m -phi.ic.psread.file.name = tests/PSReadIC/2d.xyzr -phi.ic.psread.eps = 1.13e-5_m - -# Use the full-feedback propellant model -propellant.type = fullfeedback -# Arrhenius rate law params -propellant.fullfeedback.m_ap = 1.45e5_m/s -propellant.fullfeedback.m_htpb = 1.04e2_m/s -propellant.fullfeedback.E_ap = 11000.0_K -propellant.fullfeedback.E_htpb = 7500.0_K -propellant.fullfeedback.mob_ap = true -propellant.fullfeedback.bound = 300.0_K -# Thermal transport params -propellant.fullfeedback.rho_ap = 1950.0_kg/m^3 -propellant.fullfeedback.rho_htpb = 920.0_kg/m^3 -propellant.fullfeedback.k_ap = 0.4186e0_W/m/K -propellant.fullfeedback.k_htpb = 0.13_W/m/K -propellant.fullfeedback.cp_ap = 1297.9_J/kg/K -propellant.fullfeedback.cp_htpb = 2418.29_J/kg/K -# Mass-to-heat flux estimation params -propellant.fullfeedback.a1 = 1.11400_1 -propellant.fullfeedback.a2 = 0.46000_1 -propellant.fullfeedback.a3 = 2.79700_1 -propellant.fullfeedback.b1 = 0.32300_1 -propellant.fullfeedback.b2 = 0.42000_1 -propellant.fullfeedback.b3 = 0.32250_1 -propellant.fullfeedback.c1 = -0.09906_1 -propellant.fullfeedback.mlocal_ap = 1000.0_kg/m^3/s -propellant.fullfeedback.mlocal_comb = 0.0_kg/m^3/s -propellant.fullfeedback.phi.zeta_0 = 1.13e-5_m -propellant.fullfeedback.phi.zeta = 1.13e-5 _m -propellant.fullfeedback.Pref = 1.0_MPa - - -# Activate thermal transport -thermal.on = 1 - -# Temperature boundary conditions -thermal.temp.bc.type = constant -thermal.temp.bc.constant.type.xlo = neumann -thermal.temp.bc.constant.type.xhi = neumann -thermal.temp.bc.constant.type.ylo = neumann -thermal.temp.bc.constant.type.yhi = neumann -thermal.temp.bc.constant.type.zlo = neumann -thermal.temp.bc.constant.type.zhi = neumann -thermal.temp.bc.constant.val.xlo = 0.0 -thermal.temp.bc.constant.val.xhi = 0.0 -thermal.temp.bc.constant.val.ylo = 0.0 -thermal.temp.bc.constant.val.yhi = 0.0 -thermal.temp.bc.constant.val.zlo = 0.0 -thermal.temp.bc.constant.val.zhi = 0.0 - -# Temperature initial conditions -temp.ic.type="constant" -temp.ic.constant.value=500_K - -# Laser source term -laser.ic.type = constant -laser.ic.constant.value = 0.0_W/m^2 -thermal.hc = 1.0e7_W/m^2 - -# Pressure -chamber.pressure = 4.0_MPa # [MPa] - -elastic.type = disable - - - - diff --git a/tests/PSReadIC/reference/2d-parallel.dat b/tests/PSReadIC/reference/2d-parallel.dat deleted file mode 100644 index 87a1591685..0000000000 --- a/tests/PSReadIC/reference/2d-parallel.dat +++ /dev/null @@ -1,109 +0,0 @@ -,x,y,z,eta,phi -0,3.90625e-06,0.00043359375,0.5,1.0,0.4582811715952242 -1,1.171875e-05,0.00043359375,0.5,1.0,0.4533531771300272 -2,1.953125e-05,0.00043359375,0.5,1.0,0.48465814735371426 -3,2.734375e-05,0.00043359375,0.5,1.0,0.6187702232104019 -4,3.5156250000000004e-05,0.00043359375,0.5,1.0,0.8099244424445855 -5,4.296875e-05,0.00043359375,0.5,1.0,0.9445592890653867 -6,5.0781250000000004e-05,0.00043359375,0.5,1.0,0.9899106690231667 -7,5.859375e-05,0.00043359375,0.5,1.0,0.9987226708136097 -8,6.640625e-05,0.00043359375000000006,0.5,1.0,0.9998294728446361 -9,7.421875e-05,0.00043359375000000006,0.5,1.0,0.9999535186180926 -10,8.0078125e-05,0.00043164062500000006,0.5,1.0,0.9999943161163413 -11,8.398437500000001e-05,0.00043164062500000006,0.5,1.0,0.9999895642889608 -12,8.7890625e-05,0.00043164062500000006,0.5,1.0,0.9999725150198857 -13,9.1796875e-05,0.00043164062500000006,0.5,1.0,0.9999099977284864 -14,9.5703125e-05,0.00043164062500000006,0.5,1.0,0.9996756221836025 -15,9.960937500000001e-05,0.00043164062500000006,0.5,1.0,0.9988238561477296 -16,0.00010351562500000001,0.00043164062500000006,0.5,1.0,0.9959673688757547 -17,0.000107421875,0.00043164062500000006,0.5,0.8835343699143661,0.9874509381755153 -18,0.000111328125,0.00043164062500000006,0.5,0.8208671646157901,0.9654712883579316 -19,0.000115234375,0.00043164062500000006,0.5,0.740447256770379,0.9172788186274969 -20,0.000119140625,0.00043164062500000006,0.5,0.6446918295319187,0.8286858464419236 -21,0.000123046875,0.00043164062500000006,0.5,0.5388871141854548,0.6934180757154493 -22,0.000126953125,0.00043164062500000006,0.5,0.4303774999022171,0.5230522787196434 -23,0.000130859375,0.00043164062500000006,0.5,0.327065287188968,0.3469737918465656 -24,0.000134765625,0.00043164062500000006,0.5,0.23572152800605672,0.19825339418877008 -25,0.00013867187500000001,0.00043164062500000006,0.5,0.16069643584910842,0.06722993939132167 -26,0.00014257812500000002,0.00043164062500000006,0.5,0.10342525048783302,0.0 -27,0.00014648437500000002,0.00043164062500000006,0.5,0.06276864574204888,0.0 -28,0.000150390625,0.00043164062500000006,0.5,0.035906463899466834,0.0 -29,0.000154296875,0.00043164062500000006,0.5,0.019370086331836456,0.0 -30,0.00016015625,0.00043359375000000006,0.5,0.007910337399139489,0.0 -31,0.00016796875,0.00043359375000000006,0.5,0.0018120266525033513,0.0 -32,0.00017578125,0.00043359375000000006,0.5,0.00039759964545038606,0.0 -33,0.00018359375,0.00043359375000000006,0.5,0.00012249730266478132,0.025747918070621933 -34,0.00019140625,0.00043359375000000006,0.5,0.00011325211219714859,0.15443901372791802 -35,0.00019921875000000001,0.00043359375000000006,0.5,0.0003381402778630527,0.40233391751008196 -36,0.00020703125,0.00043359375000000006,0.5,0.0015231130960052797,0.672315709249885 -37,0.00021484375,0.00043359375000000006,0.5,0.0067424768577353245,0.8634516875769216 -38,0.000220703125,0.00043164062500000006,0.5,0.016706464664568266,0.9687428716690855 -39,0.000224609375,0.00043164062500000006,0.5,0.031377608252043276,0.9848584579809165 -40,0.000228515625,0.00043164062500000006,0.5,0.05559037927896038,0.9927392119335315 -41,0.00023242187500000002,0.00043164062500000006,0.5,0.09283330995526129,0.9964418101299328 -42,0.000236328125,0.00043164062500000006,0.5,0.14616106890974956,0.998153390257225 -43,0.000240234375,0.00043164062500000006,0.5,0.21718402864992847,0.9989463945822822 -44,0.000244140625,0.00043164062500000006,0.5,0.305106466416148,0.9993147640037052 -45,0.000248046875,0.00043164062500000006,0.5,0.4062277011976442,0.9994761631709994 -46,0.000251953125,0.00043164062500000006,0.5,0.514236973096488,0.9995194778817386 -47,0.000255859375,0.00043164062500000006,0.5,0.6213458666453286,0.9994674167066254 -48,0.000259765625,0.00043164062500000006,0.5,0.6523439128133405,0.9992926164031527 -49,0.000263671875,0.00043164062500000006,0.5,0.6937417449035357,0.998898237327045 -50,0.000267578125,0.00043164062500000006,0.5,0.7039234649993117,0.9980498737388569 -51,0.000271484375,0.00043164062500000006,0.5,0.6841301611019979,0.9962176589430582 -52,0.00027539062500000003,0.00043164062500000006,0.5,0.6707762212007404,0.9922562938388506 -53,0.00027929687500000003,0.00043164062500000006,0.5,0.5926683111782173,0.9838461985307727 -54,0.00028320312500000003,0.00043164062500000006,0.5,0.512820071885169,0.9667287443991215 -55,0.00028710937500000004,0.00043164062500000006,0.5,0.43477946748927565,0.9370021576890688 -56,0.00029101562500000004,0.00043164062500000006,0.5,0.3616458842026158,0.8926582739231124 -57,0.00029492187500000005,0.00043164062500000006,0.5,0.2957269026442766,0.8347080272160838 -58,0.00029882812500000005,0.00043164062500000006,0.5,0.23838775313694593,0.7745539452947771 -59,0.00030273437500000005,0.00043164062500000006,0.5,0.19008567696192225,0.7315880621219365 -60,0.00030664062500000006,0.00043164062500000006,0.5,0.15053738366430325,0.720084201568417 -61,0.00031054687500000006,0.00043164062500000006,0.5,0.11895049545853575,0.7408545427050972 -62,0.00031445312500000006,0.00043164062500000006,0.5,0.09425678677264238,0.7797203886904559 -63,0.000318359375,0.00043164062500000006,0.5,0.075305904320766,0.8232279314546171 -64,0.000322265625,0.00043164062500000006,0.5,0.06100181346413858,0.8652624889770042 -65,0.000326171875,0.00043164062500000006,0.5,0.05038274537594123,0.8981342938629702 -66,0.000330078125,0.00043164062500000006,0.5,0.042656032003524746,0.9191696364636631 -67,0.000333984375,0.00043164062500000006,0.5,0.03720270220531097,0.9314811888879476 -68,0.00033789062500000003,0.00043164062500000006,0.5,0.033565555548943915,0.937122601176575 -69,0.00034179687500000003,0.00043164062500000006,0.5,0.03143118716752369,0.9370905002307901 -70,0.00034570312500000004,0.00043164062500000006,0.5,0.03061291751220796,0.9313791198983925 -71,0.00034960937500000004,0.00043164062500000006,0.5,0.03103869931202019,0.9189798823410011 -72,0.00035351562500000004,0.00043164062500000006,0.5,0.032746052998986364,0.897826535687336 -73,0.00035742187500000005,0.00043164062500000006,0.5,0.035884732215552884,0.8647948813507121 -74,0.00036132812500000005,0.00043164062500000006,0.5,0.04072676073536963,0.8159879414985258 -75,0.00036523437500000005,0.00043164062500000006,0.5,0.04768226911073992,0.7476192072705634 -76,0.000369140625,0.00043164062500000006,0.5,0.05731778166048307,0.6576883711775607 -77,0.000373046875,0.00043164062500000006,0.5,0.07037099334506203,0.5481426878831556 -78,0.000376953125,0.00043164062500000006,0.5,0.0877526963275912,0.42642257545575435 -79,0.000380859375,0.00043164062500000006,0.5,0.1105230353942922,0.3048592220905742 -80,0.000384765625,0.00043164062500000006,0.5,0.1398272117768098,0.19722116560831943 -81,0.000388671875,0.00043164062500000006,0.5,0.17677749010723032,0.1678610307161649 -82,0.000392578125,0.00043164062500000006,0.5,0.22227653795347058,0.2282379134165803 -83,0.000396484375,0.00043164062500000006,0.5,0.2767932751785763,0.3512146614648787 -84,0.00040039062500000003,0.00043164062500000006,0.5,0.3401249927568327,0.5197524683186292 -85,0.00040429687500000003,0.00043164062500000006,0.5,0.4112023248556993,0.681807525502341 -86,0.00040820312500000004,0.00043164062500000006,0.5,0.4880062268070724,0.811976863203694 -87,0.00041210937500000004,0.00043164062500000006,0.5,0.5676571744429646,0.9003208354367438 -88,0.00041601562500000004,0.00043164062500000006,0.5,0.6467006915761834,0.9517714808776861 -89,0.00041992187500000005,0.00043164062500000006,0.5,0.7215561004276289,0.978031920424588 -90,0.00042382812500000005,0.00043164062500000006,0.5,0.7890362461056692,0.9901089469743323 -91,0.000427734375,0.00043164062500000006,0.5,0.8468105962133419,0.9952790268496151 -92,0.000431640625,0.00043164062500000006,0.5,0.7835297915920545,0.997386420875755 -93,0.000435546875,0.00043164062500000006,0.5,0.813112582945398,0.9981528368536158 -94,0.000439453125,0.00043164062500000006,0.5,0.7326544224958522,0.9982138373458963 -95,0.000443359375,0.00043164062500000006,0.5,0.6380592020564385,0.997619584311696 -96,0.000447265625,0.00043164062500000006,0.5,0.5346112883878035,0.9958781596259787 -97,0.000451171875,0.00043164062500000006,0.5,0.42931656875456514,0.9915555696939881 -98,0.000455078125,0.00043164062500000006,0.5,0.32948353810703557,0.9813476003025368 -99,0.00045898437500000003,0.00043164062500000006,0.5,0.24122230431057734,0.9587643157545075 -100,0.00046289062500000003,0.00043164062500000006,0.5,0.16836626219655204,0.917917532292608 -101,0.00046679687500000004,0.00043164062500000006,0.5,0.11211884629916469,0.856387492723463 -102,0.00047070312500000004,0.00043164062500000006,0.5,0.07141359412910775,0.7763375016584368 -103,0.00047460937500000004,0.00043164062500000006,0.5,0.04371994443225746,0.6943510972969217 -104,0.000478515625,0.00043164062500000006,0.5,0.0259394580603316,0.6354585405288947 -105,0.000482421875,0.00043164062500000006,0.5,0.015114926193632972,0.6139670453370898 -106,0.00048828125,0.00043359375000000006,0.5,0.005116696568355672,0.7004280271117741 -107,0.00049609375,0.00043359375000000006,0.5,0.0020004992962890933,0.7562051131588987 diff --git a/tests/PSReadIC/reference/2d-serial.dat b/tests/PSReadIC/reference/2d-serial.dat deleted file mode 100644 index 87a1591685..0000000000 --- a/tests/PSReadIC/reference/2d-serial.dat +++ /dev/null @@ -1,109 +0,0 @@ -,x,y,z,eta,phi -0,3.90625e-06,0.00043359375,0.5,1.0,0.4582811715952242 -1,1.171875e-05,0.00043359375,0.5,1.0,0.4533531771300272 -2,1.953125e-05,0.00043359375,0.5,1.0,0.48465814735371426 -3,2.734375e-05,0.00043359375,0.5,1.0,0.6187702232104019 -4,3.5156250000000004e-05,0.00043359375,0.5,1.0,0.8099244424445855 -5,4.296875e-05,0.00043359375,0.5,1.0,0.9445592890653867 -6,5.0781250000000004e-05,0.00043359375,0.5,1.0,0.9899106690231667 -7,5.859375e-05,0.00043359375,0.5,1.0,0.9987226708136097 -8,6.640625e-05,0.00043359375000000006,0.5,1.0,0.9998294728446361 -9,7.421875e-05,0.00043359375000000006,0.5,1.0,0.9999535186180926 -10,8.0078125e-05,0.00043164062500000006,0.5,1.0,0.9999943161163413 -11,8.398437500000001e-05,0.00043164062500000006,0.5,1.0,0.9999895642889608 -12,8.7890625e-05,0.00043164062500000006,0.5,1.0,0.9999725150198857 -13,9.1796875e-05,0.00043164062500000006,0.5,1.0,0.9999099977284864 -14,9.5703125e-05,0.00043164062500000006,0.5,1.0,0.9996756221836025 -15,9.960937500000001e-05,0.00043164062500000006,0.5,1.0,0.9988238561477296 -16,0.00010351562500000001,0.00043164062500000006,0.5,1.0,0.9959673688757547 -17,0.000107421875,0.00043164062500000006,0.5,0.8835343699143661,0.9874509381755153 -18,0.000111328125,0.00043164062500000006,0.5,0.8208671646157901,0.9654712883579316 -19,0.000115234375,0.00043164062500000006,0.5,0.740447256770379,0.9172788186274969 -20,0.000119140625,0.00043164062500000006,0.5,0.6446918295319187,0.8286858464419236 -21,0.000123046875,0.00043164062500000006,0.5,0.5388871141854548,0.6934180757154493 -22,0.000126953125,0.00043164062500000006,0.5,0.4303774999022171,0.5230522787196434 -23,0.000130859375,0.00043164062500000006,0.5,0.327065287188968,0.3469737918465656 -24,0.000134765625,0.00043164062500000006,0.5,0.23572152800605672,0.19825339418877008 -25,0.00013867187500000001,0.00043164062500000006,0.5,0.16069643584910842,0.06722993939132167 -26,0.00014257812500000002,0.00043164062500000006,0.5,0.10342525048783302,0.0 -27,0.00014648437500000002,0.00043164062500000006,0.5,0.06276864574204888,0.0 -28,0.000150390625,0.00043164062500000006,0.5,0.035906463899466834,0.0 -29,0.000154296875,0.00043164062500000006,0.5,0.019370086331836456,0.0 -30,0.00016015625,0.00043359375000000006,0.5,0.007910337399139489,0.0 -31,0.00016796875,0.00043359375000000006,0.5,0.0018120266525033513,0.0 -32,0.00017578125,0.00043359375000000006,0.5,0.00039759964545038606,0.0 -33,0.00018359375,0.00043359375000000006,0.5,0.00012249730266478132,0.025747918070621933 -34,0.00019140625,0.00043359375000000006,0.5,0.00011325211219714859,0.15443901372791802 -35,0.00019921875000000001,0.00043359375000000006,0.5,0.0003381402778630527,0.40233391751008196 -36,0.00020703125,0.00043359375000000006,0.5,0.0015231130960052797,0.672315709249885 -37,0.00021484375,0.00043359375000000006,0.5,0.0067424768577353245,0.8634516875769216 -38,0.000220703125,0.00043164062500000006,0.5,0.016706464664568266,0.9687428716690855 -39,0.000224609375,0.00043164062500000006,0.5,0.031377608252043276,0.9848584579809165 -40,0.000228515625,0.00043164062500000006,0.5,0.05559037927896038,0.9927392119335315 -41,0.00023242187500000002,0.00043164062500000006,0.5,0.09283330995526129,0.9964418101299328 -42,0.000236328125,0.00043164062500000006,0.5,0.14616106890974956,0.998153390257225 -43,0.000240234375,0.00043164062500000006,0.5,0.21718402864992847,0.9989463945822822 -44,0.000244140625,0.00043164062500000006,0.5,0.305106466416148,0.9993147640037052 -45,0.000248046875,0.00043164062500000006,0.5,0.4062277011976442,0.9994761631709994 -46,0.000251953125,0.00043164062500000006,0.5,0.514236973096488,0.9995194778817386 -47,0.000255859375,0.00043164062500000006,0.5,0.6213458666453286,0.9994674167066254 -48,0.000259765625,0.00043164062500000006,0.5,0.6523439128133405,0.9992926164031527 -49,0.000263671875,0.00043164062500000006,0.5,0.6937417449035357,0.998898237327045 -50,0.000267578125,0.00043164062500000006,0.5,0.7039234649993117,0.9980498737388569 -51,0.000271484375,0.00043164062500000006,0.5,0.6841301611019979,0.9962176589430582 -52,0.00027539062500000003,0.00043164062500000006,0.5,0.6707762212007404,0.9922562938388506 -53,0.00027929687500000003,0.00043164062500000006,0.5,0.5926683111782173,0.9838461985307727 -54,0.00028320312500000003,0.00043164062500000006,0.5,0.512820071885169,0.9667287443991215 -55,0.00028710937500000004,0.00043164062500000006,0.5,0.43477946748927565,0.9370021576890688 -56,0.00029101562500000004,0.00043164062500000006,0.5,0.3616458842026158,0.8926582739231124 -57,0.00029492187500000005,0.00043164062500000006,0.5,0.2957269026442766,0.8347080272160838 -58,0.00029882812500000005,0.00043164062500000006,0.5,0.23838775313694593,0.7745539452947771 -59,0.00030273437500000005,0.00043164062500000006,0.5,0.19008567696192225,0.7315880621219365 -60,0.00030664062500000006,0.00043164062500000006,0.5,0.15053738366430325,0.720084201568417 -61,0.00031054687500000006,0.00043164062500000006,0.5,0.11895049545853575,0.7408545427050972 -62,0.00031445312500000006,0.00043164062500000006,0.5,0.09425678677264238,0.7797203886904559 -63,0.000318359375,0.00043164062500000006,0.5,0.075305904320766,0.8232279314546171 -64,0.000322265625,0.00043164062500000006,0.5,0.06100181346413858,0.8652624889770042 -65,0.000326171875,0.00043164062500000006,0.5,0.05038274537594123,0.8981342938629702 -66,0.000330078125,0.00043164062500000006,0.5,0.042656032003524746,0.9191696364636631 -67,0.000333984375,0.00043164062500000006,0.5,0.03720270220531097,0.9314811888879476 -68,0.00033789062500000003,0.00043164062500000006,0.5,0.033565555548943915,0.937122601176575 -69,0.00034179687500000003,0.00043164062500000006,0.5,0.03143118716752369,0.9370905002307901 -70,0.00034570312500000004,0.00043164062500000006,0.5,0.03061291751220796,0.9313791198983925 -71,0.00034960937500000004,0.00043164062500000006,0.5,0.03103869931202019,0.9189798823410011 -72,0.00035351562500000004,0.00043164062500000006,0.5,0.032746052998986364,0.897826535687336 -73,0.00035742187500000005,0.00043164062500000006,0.5,0.035884732215552884,0.8647948813507121 -74,0.00036132812500000005,0.00043164062500000006,0.5,0.04072676073536963,0.8159879414985258 -75,0.00036523437500000005,0.00043164062500000006,0.5,0.04768226911073992,0.7476192072705634 -76,0.000369140625,0.00043164062500000006,0.5,0.05731778166048307,0.6576883711775607 -77,0.000373046875,0.00043164062500000006,0.5,0.07037099334506203,0.5481426878831556 -78,0.000376953125,0.00043164062500000006,0.5,0.0877526963275912,0.42642257545575435 -79,0.000380859375,0.00043164062500000006,0.5,0.1105230353942922,0.3048592220905742 -80,0.000384765625,0.00043164062500000006,0.5,0.1398272117768098,0.19722116560831943 -81,0.000388671875,0.00043164062500000006,0.5,0.17677749010723032,0.1678610307161649 -82,0.000392578125,0.00043164062500000006,0.5,0.22227653795347058,0.2282379134165803 -83,0.000396484375,0.00043164062500000006,0.5,0.2767932751785763,0.3512146614648787 -84,0.00040039062500000003,0.00043164062500000006,0.5,0.3401249927568327,0.5197524683186292 -85,0.00040429687500000003,0.00043164062500000006,0.5,0.4112023248556993,0.681807525502341 -86,0.00040820312500000004,0.00043164062500000006,0.5,0.4880062268070724,0.811976863203694 -87,0.00041210937500000004,0.00043164062500000006,0.5,0.5676571744429646,0.9003208354367438 -88,0.00041601562500000004,0.00043164062500000006,0.5,0.6467006915761834,0.9517714808776861 -89,0.00041992187500000005,0.00043164062500000006,0.5,0.7215561004276289,0.978031920424588 -90,0.00042382812500000005,0.00043164062500000006,0.5,0.7890362461056692,0.9901089469743323 -91,0.000427734375,0.00043164062500000006,0.5,0.8468105962133419,0.9952790268496151 -92,0.000431640625,0.00043164062500000006,0.5,0.7835297915920545,0.997386420875755 -93,0.000435546875,0.00043164062500000006,0.5,0.813112582945398,0.9981528368536158 -94,0.000439453125,0.00043164062500000006,0.5,0.7326544224958522,0.9982138373458963 -95,0.000443359375,0.00043164062500000006,0.5,0.6380592020564385,0.997619584311696 -96,0.000447265625,0.00043164062500000006,0.5,0.5346112883878035,0.9958781596259787 -97,0.000451171875,0.00043164062500000006,0.5,0.42931656875456514,0.9915555696939881 -98,0.000455078125,0.00043164062500000006,0.5,0.32948353810703557,0.9813476003025368 -99,0.00045898437500000003,0.00043164062500000006,0.5,0.24122230431057734,0.9587643157545075 -100,0.00046289062500000003,0.00043164062500000006,0.5,0.16836626219655204,0.917917532292608 -101,0.00046679687500000004,0.00043164062500000006,0.5,0.11211884629916469,0.856387492723463 -102,0.00047070312500000004,0.00043164062500000006,0.5,0.07141359412910775,0.7763375016584368 -103,0.00047460937500000004,0.00043164062500000006,0.5,0.04371994443225746,0.6943510972969217 -104,0.000478515625,0.00043164062500000006,0.5,0.0259394580603316,0.6354585405288947 -105,0.000482421875,0.00043164062500000006,0.5,0.015114926193632972,0.6139670453370898 -106,0.00048828125,0.00043359375000000006,0.5,0.005116696568355672,0.7004280271117741 -107,0.00049609375,0.00043359375000000006,0.5,0.0020004992962890933,0.7562051131588987 diff --git a/tests/PSReadIC/test b/tests/PSReadIC/test deleted file mode 100755 index 80ca0e2591..0000000000 --- a/tests/PSReadIC/test +++ /dev/null @@ -1,27 +0,0 @@ -#!/usr/bin/env python3 -import sys -sys.path.insert(0,"../../scripts") -import testlib -import glob - -outdir = sys.argv[1] -reference = sys.argv[2] - -path = sorted(glob.glob("{}/*cell/".format(outdir)))[-1] - -testlib.validate(path=path, - outdir=outdir, - start=[0.0000,0.00043,0], - end=[0.0005,0.00043,0], - vars=["eta","phi"], - tolerance=[0.01,0.01], - abs_tolerance=[0.01,0.01], - reference=reference, - generate_ref_data = False, -) -exit(0) - - - - - diff --git a/tests/SCPSpheresElastic/input b/tests/SCPSpheresElastic/input index 36fabbbe2f..092be8f86b 100644 --- a/tests/SCPSpheresElastic/input +++ b/tests/SCPSpheresElastic/input @@ -1,4 +1,18 @@ -#@ +#@ [2d-psread-coverage] +#@ dim=2 +#@ nprocs=1 +#@ check-file = reference/2d-psread-coverage.dat +#@ args = timestep = 1e-8 +#@ args = stop_time = 1e-10 +#@ args = pf.eta.ic.type = psread +#@ args = pf.eta.ic.psread.file.unit = m +#@ args = pf.eta.ic.psread.file.name = tests/SCPSpheresElastic/2d.xyzr +#@ args = pf.eta.ic.psread.eps = 2.0e-5_m +#@ args = pf.eta.ic.psread.invert = 1 +#@ args = pf.eta.ic.psread.x0 = 1e-4 1e-4 0.0 +#@ args = pf.eta.ic.psread.mult = 1.1 +#@ args = allow_unused=True +#@ #@ [2d-serial-short] #@ dim=2 #@ nprocs=1 @@ -12,6 +26,7 @@ #@ check-file = reference/2d-parallel-long.dat #@ + alamo.program = flame plot_file = tests/SCPSpheresElastic/output diff --git a/tests/SCPSpheresElastic/reference/2d-psread-coverage.dat b/tests/SCPSpheresElastic/reference/2d-psread-coverage.dat new file mode 100644 index 0000000000..3e61544cac --- /dev/null +++ b/tests/SCPSpheresElastic/reference/2d-psread-coverage.dat @@ -0,0 +1,109 @@ +,x,y,z,eta,disp_x,temp +0,1.953125e-06,0.00043164062500000006,0.5,0.999999947991569,-2.7254790089824363e-06,500.06623196097485 +1,5.859375e-06,0.00043164062500000006,0.5,0.9999999999999989,-2.7137430174841725e-06,500.0000018494134 +2,9.765625e-06,0.00043164062500000006,0.5,1.0,-2.702362719322129e-06,500.00000000003183 +3,1.3671875e-05,0.00043164062500000006,0.5,1.0,-2.6908031920153993e-06,500.0 +4,1.953125e-05,0.00043359374999999995,0.5,1.0,-2.68079219753425e-06,500.0 +5,2.734375e-05,0.00043359374999999995,0.5,1.0,-2.6520240304241237e-06,500.0 +6,3.5156250000000004e-05,0.00043359374999999995,0.5,1.0,-2.619656100986421e-06,500.0 +7,4.296875e-05,0.00043359374999999995,0.5,1.0,-2.5844691059492305e-06,500.0 +8,5.0781250000000004e-05,0.00043359374999999995,0.5,1.0,-2.5482408835635892e-06,500.0 +9,5.859375e-05,0.00043359374999999995,0.5,1.0,-2.51208427739249e-06,500.0 +10,6.640625e-05,0.00043359374999999995,0.5,1.0,-2.4762262520732006e-06,500.0 +11,7.421875e-05,0.00043359374999999995,0.5,1.0,-2.440597389693026e-06,500.0 +12,8.203125e-05,0.00043359374999999995,0.5,1.0,-2.405127452284172e-06,500.00000000000006 +13,8.984375e-05,0.00043359374999999995,0.5,1.0,-2.3697874885448235e-06,500.0000000000002 +14,9.5703125e-05,0.00043164062500000006,0.5,1.0,-2.337107791726974e-06,500.0000000000014 +15,9.960937500000001e-05,0.00043164062500000006,0.5,0.9999999999999994,-2.319544127242415e-06,500.0000000385357 +16,0.000103515625,0.00043164062500000006,0.5,0.9999999876909482,-2.301998698182939e-06,500.0006079567512 +17,0.000107421875,0.00043164062500000006,0.5,0.8835343699143661,-2.2845436333155715e-06,476.7079963470066 +18,0.000111328125,0.00043164062500000006,0.5,0.8208671646157901,-2.2672926775797562e-06,464.17341564195004 +19,0.000115234375,0.00043164062500000006,0.5,0.7404472567703791,-2.2504311855058375e-06,448.08956929773217 +20,0.000119140625,0.00043164062500000006,0.5,0.6446918295319187,-2.234329214182718e-06,428.9385820418053 +21,0.000123046875,0.00043164062500000006,0.5,0.5388871141571729,-2.2194960674844604e-06,407.77770734317716 +22,0.000126953125,0.00043164062500000006,0.5,0.43037749975064415,-2.2064018186103624e-06,386.07587147217345 +23,0.000130859375,0.00043164062500000006,0.5,0.3270652870057237,-2.1950300593480044e-06,365.41326441958654 +24,0.000134765625,0.00043164062500000006,0.5,0.23572152785504413,-2.184092184756887e-06,347.14438488417534 +25,0.00013867187500000001,0.00043164062500000006,0.5,0.1606964357557567,-2.1705329294958795e-06,332.139314831103 +26,0.00014257812500000002,0.00043164062500000006,0.5,0.10342525043425982,-2.1517496687478092e-06,320.68506026888065 +27,0.000146484375,0.00043164062500000006,0.5,0.06276864570870505,-2.128753347183505e-06,312.5537328191314 +28,0.000150390625,0.00043164062500000006,0.5,0.03590646388169183,-2.1040930144312018e-06,307.1812939787235 +29,0.000154296875,0.00043164062500000006,0.5,0.01937008632350626,-2.0792282153420267e-06,303.8740175047069 +30,0.00016015625,0.00043359375000000006,0.5,0.007910337397425115,-2.0401965913539734e-06,301.58206755577817 +31,0.00016796875,0.00043359375000000006,0.5,0.0018120266525033513,-1.9987204751415984e-06,300.3624053351687 +32,0.00017578125,0.00043359375000000006,0.5,0.00039759964545038606,-1.9528298691386644e-06,300.07951992935733 +33,0.00018359375,0.00043359375000000006,0.5,0.00012249730266478132,-1.8697733482083474e-06,300.0244994605613 +34,0.00019140625,0.00043359375000000006,0.5,0.00011325211219714859,-1.750791032820562e-06,300.0226504224622 +35,0.00019921875000000001,0.00043359375000000006,0.5,0.0003381402778630527,-1.6456100552077926e-06,300.06762805574095 +36,0.00020703125000000002,0.00043359375000000006,0.5,0.0015231130960052797,-1.5730661886157942e-06,300.304622622082 +37,0.00021484375,0.00043359375000000006,0.5,0.0067424768539680584,-1.5188040654383442e-06,301.3484954306152 +38,0.000220703125,0.000431640625,0.5,0.01670646463079328,-1.4893224028980504e-06,303.34129321663227 +39,0.000224609375,0.000431640625,0.5,0.0313776081764083,-1.467844253560395e-06,306.2755241577657 +40,0.000228515625,0.000431640625,0.5,0.055590379131016596,-1.4470231060420163e-06,311.1180848429554 +41,0.00023242187500000002,0.000431640625,0.5,0.09283330970864305,-1.4267278236619203e-06,318.56668696045995 +42,0.000236328125,0.000431640625,0.5,0.14616106855979386,-1.406904278097399e-06,329.2322731812362 +43,0.000240234375,0.000431640625,0.5,0.21718402823418073,-1.387521001564421e-06,343.4369277772623 +44,0.000244140625,0.000431640625,0.5,0.3051064660180982,-1.3685607470797267e-06,361.02150755087393 +45,0.000248046875,0.000431640625,0.5,0.406227700918689,-1.3500164752910623e-06,381.2458592687379 +46,0.000251953125,0.000431640625,0.5,0.5142369730103853,-1.3318798435519198e-06,402.8477458080302 +47,0.000255859375,0.000431640625,0.5,0.6213458657886387,-1.3141283052463806e-06,424.26887067786583 +48,0.000259765625,0.000431640625,0.5,0.6523439128133406,-1.296727305909328e-06,430.46943061732736 +49,0.000263671875,0.000431640625,0.5,0.6937417449035358,-1.279656796634185e-06,438.74810173141924 +50,0.000267578125,0.000431640625,0.5,0.7039234649993117,-1.262912145168623e-06,440.78443701166975 +51,0.000271484375,0.000431640625,0.5,0.684130161101998,-1.2464414704268876e-06,436.8267943357596 +52,0.00027539062500000003,0.000431640625,0.5,0.6707762208729219,-1.2301629706559404e-06,434.1551516520934 +53,0.00027929687500000003,0.000431640625,0.5,0.5926683111782173,-1.2140761127430877e-06,418.5338425000262 +54,0.00028320312500000003,0.000431640625,0.5,0.512820071785883,-1.1981709483237515e-06,402.5643634299852 +55,0.000287109375,0.000431640625,0.5,0.43477946725572936,-1.182377576997528e-06,386.9562438248246 +56,0.000291015625,0.000431640625,0.5,0.361645883878257,-1.1666450078447632e-06,372.3294723672921 +57,0.000294921875,0.000431640625,0.5,0.29572690228220383,-1.1509086246093955e-06,359.14561976872443 +58,0.00029882812500000005,0.00043164062500000006,0.5,0.23838775278223467,-1.1350145481706077e-06,347.6777355712519 +59,0.00030273437500000005,0.00043164062500000006,0.5,0.1900856766390146,-1.1186280227620514e-06,338.01726847356736 +60,0.00030664062500000006,0.00043164062500000006,0.5,0.1505373833796441,-1.1012465921783682e-06,330.10756454465127 +61,0.00031054687500000006,0.00043164062500000006,0.5,0.11895049521054253,-1.082349386756508e-06,323.79015360387865 +62,0.00031445312500000006,0.00043164062500000006,0.5,0.09425678655923347,-1.0615845236268234e-06,318.8513900761477 +63,0.00031835937500000007,0.00043164062500000006,0.5,0.07530590413954072,-1.0388188935454794e-06,315.06120042007853 +64,0.00032226562500000007,0.00043164062500000006,0.5,0.061001813311297975,-1.0140627022249251e-06,312.20037463862695 +65,0.00032617187500000007,0.00043164062500000006,0.5,0.05038274524728651,-9.874341112518242e-07,310.07655671426795 +66,0.0003300781250000001,0.00043164062500000006,0.5,0.042656031894459955,-9.591246170647837e-07,308.5312115851417 +67,0.000333984375,0.00043164062500000006,0.5,0.037202702111004494,-9.293674055090761e-07,307.44054420315916 +68,0.00033789062500000003,0.00043164062500000006,0.5,0.03356555546491881,-8.984427849693326e-07,306.7131140562027 +69,0.00034179687500000003,0.00043164062500000006,0.5,0.03143118708981076,-8.666829291362102e-07,306.2862399469343 +70,0.00034570312500000004,0.00043164062500000006,0.5,0.030612917437278005,-8.344730782282548e-07,306.1225858538605 +71,0.00034960937500000004,0.00043164062500000006,0.5,0.031038699236640443,-8.022373948600779e-07,306.20774228284756 +72,0.00035351562500000004,0.00043164062500000006,0.5,0.03274605292009599,-7.704066010558463e-07,306.54921334331067 +73,0.00035742187500000005,0.00043164062500000006,0.5,0.03588473213023869,-7.393896636679949e-07,307.1769498627597 +74,0.00036132812500000005,0.00043164062500000006,0.5,0.04072676064102961,-7.095736821324763e-07,308.1453568102413 +75,0.00036523437500000005,0.00043164062500000006,0.5,0.0476822690055171,-6.813289903003362e-07,309.53646067053575 +76,0.00036914062500000006,0.00043164062500000006,0.5,0.0573177815439676,-6.550093208608057e-07,311.4635667126494 +77,0.00037304687500000006,0.00043164062500000006,0.5,0.07037099321903507,-6.309478269550206e-07,314.07421372247956 +78,0.00037695312500000006,0.00043164062500000006,0.5,0.08775269619627445,-6.094441661520415e-07,317.5505587681298 +79,0.00038085937500000007,0.00043164062500000006,0.5,0.11052303526349849,-5.907884548925204e-07,322.10462994397153 +80,0.000384765625,0.00043164062500000006,0.5,0.13982721165198983,-5.75190213470398e-07,327.9654693597169 +81,0.000388671875,0.00043164062500000006,0.5,0.17677748997427092,-5.62487947411793e-07,335.35553753645075 +82,0.000392578125,0.00043164062500000006,0.5,0.22227653778924153,-5.520049099785188e-07,344.45538287965326 +83,0.000396484375,0.00043164062500000006,0.5,0.2767932749726171,-5.42483491973293e-07,355.3588155591067 +84,0.00040039062500000003,0.00043164062500000006,0.5,0.34012499252167594,-5.323310336621909e-07,368.0253130345521 +85,0.00040429687500000003,0.00043164062500000006,0.5,0.41120232464198775,-5.205093655852377e-07,382.2408736264229 +86,0.00040820312500000004,0.00043164062500000006,0.5,0.488006226678582,-5.06863297134483e-07,397.6016285383439 +87,0.00041210937500000004,0.00043164062500000006,0.5,0.5676571744429646,-4.917112614170922e-07,413.5316252604507 +88,0.00041601562500000004,0.00043164062500000006,0.5,0.6467006915761834,-4.754920531251931e-07,429.3402899958386 +89,0.00041992187500000005,0.00043164062500000006,0.5,0.721556100427629,-4.5856286356331133e-07,444.3112958119607 +90,0.00042382812500000005,0.00043164062500000006,0.5,0.7890362461056692,-4.4115745722168133e-07,457.80723201289914 +91,0.000427734375,0.00043164062500000006,0.5,0.8468105918084019,-4.234169406496953e-07,469.35935251585624 +92,0.000431640625,0.00043164062500000006,0.5,0.7835297915920545,-4.054407938267584e-07,456.7071032708084 +93,0.000435546875,0.00043164062500000006,0.5,0.8131125799370595,-3.872773151693242e-07,462.62134594783686 +94,0.000439453125,0.00043164062500000006,0.5,0.7326544224958522,-3.6893087394371027e-07,446.5309661533894 +95,0.000443359375,0.00043164062500000006,0.5,0.6380592020564385,-3.503728724469086e-07,427.61201729967087 +96,0.000447265625,0.00043164062500000006,0.5,0.5346112883500763,-3.31539468343993e-07,406.9225530674465 +97,0.000451171875,0.00043164062500000006,0.5,0.4293165685219856,-3.1235178425446937e-07,385.8636623633551 +98,0.000455078125,0.00043164062500000006,0.5,0.32948353774317657,-2.92697074855792e-07,365.8969576664204 +99,0.00045898437500000003,0.00043164062500000006,0.5,0.24122230391126098,-2.724273829162645e-07,348.24461613272973 +100,0.00046289062500000003,0.00043164062500000006,0.5,0.16836626184926198,-2.5135632230054375e-07,333.6733382316478 +101,0.00046679687500000004,0.00043164062500000006,0.5,0.11211884604878315,-2.2925029329840716e-07,322.4238123327418 +102,0.00047070312500000004,0.00043164062500000006,0.5,0.07141359397570844,-2.0586423519156267e-07,314.282736470658 +103,0.00047460937500000004,0.00043164062500000006,0.5,0.04371994435012459,-1.810780077578208e-07,308.74399402219865 +104,0.000478515625,0.00043164062500000006,0.5,0.02593945802059488,-1.5506845496131043e-07,305.1878928509446 +105,0.000482421875,0.00043164062500000006,0.5,0.01511492617335381,-1.284012381860753e-07,303.0229854127716 +106,0.00048828125,0.00043359375000000006,0.5,0.005116696568355561,-9.047745040917093e-08,301.0233393430549 +107,0.00049609375,0.00043359375000000006,0.5,0.0020004992962892043,-3.04525130035454e-08,300.40010033098173