Skip to content

Equil Developement - #13

Merged
jhalpern30 merged 109 commits into
mainfrom
equil_dev
Oct 9, 2025
Merged

Equil Developement#13
jhalpern30 merged 109 commits into
mainfrom
equil_dev

Conversation

@logan-nc

Copy link
Copy Markdown
Collaborator

This is a placeholder for the equil development. More documentation will be added later

@logan-nclogan-nc added the feature New capability label Jul 28, 2025
@logan-nc

Copy link
Copy Markdown
CollaboratorAuthor

@jmlmir369 the equil_test notebook gives me the following error after the second cell:

Details
Starting equilibrium reconstruction...
--- Julia Equilibrium Setup ---
Equilibrium file: beta_1.00
Type = efit, Jac_type = boozer
----------------------------------------
--> Processing EFIT g-file: beta_1.00
--> Parsed from header: nw=129, nh=128
--> All main data blocks parsed successfully.
--> Creating 1D profile splines...
--> 1D Spline fitting complete.
--> Creating 2D psi spline...
--> 2D Spline fitting complete.
--- Starting Direct Equilibrium Processing ---
Finding magnetic axis...
Iter 1: R = 3.115238, Z = 0.000008, |ΔR|=1.15e-01, |ΔZ|=1.11e-05
Iter 2: R = 3.117635, Z = -0.000003, |ΔR|=2.40e-03, |ΔZ|=1.10e-05
Iter 3: R = 3.117636, Z = -0.000003, |ΔR|=1.07e-06, |ΔZ|=1.04e-09
Iter 4: R = 3.117636, Z = -0.000003, |ΔR|=2.12e-13, |ΔZ|=5.62e-16
Magnetic axis found at R=3.1176357129137466, Z=-2.9184696943092166e-6.
Finding inboard separatrix crossing...
Restart attempt 1/6 with initial R = 2.675802
inboard separatrix found at R=1.5000175239830411.
Finding outboard separatrix crossing...
Restart attempt 1/6 with initial R = 3.500684
outboard separatrix found at R=4.499995001114418.
Starting loop over flux surfaces...
--> Processing surface ipsi = 100 / 100 (ψ_norm = 1.0000)
--> Processing surface ipsi = 99 / 100 (ψ_norm = 0.9998)
--> Processing surface ipsi = 98 / 100 (ψ_norm = 0.9990)
--> Processing surface ipsi = 97 / 100 (ψ_norm = 0.9978)
--> Processing surface ipsi = 96 / 100 (ψ_norm = 0.9961)
--> Processing surface ipsi = 95 / 100 (ψ_norm = 0.9938)
--> Processing surface ipsi = 94 / 100 (ψ_norm = 0.9911)
--> Processing surface ipsi = 93 / 100 (ψ_norm = 0.9880)
--> Processing surface ipsi = 92 / 100 (ψ_norm = 0.9843)
--> Processing surface ipsi = 91 / 100 (ψ_norm = 0.9801)
--> Processing surface ipsi = 90 / 100 (ψ_norm = 0.9755)
--> Processing surface ipsi = 89 / 100 (ψ_norm = 0.9704)
--> Processing surface ipsi = 88 / 100 (ψ_norm = 0.9649)
--> Processing surface ipsi = 87 / 100 (ψ_norm = 0.9589)
--> Processing surface ipsi = 86 / 100 (ψ_norm = 0.9524)
--> Processing surface ipsi = 85 / 100 (ψ_norm = 0.9455)
--> Processing surface ipsi = 84 / 100 (ψ_norm = 0.9382)
--> Processing surface ipsi = 83 / 100 (ψ_norm = 0.9304)
--> Processing surface ipsi = 82 / 100 (ψ_norm = 0.9222)
--> Processing surface ipsi = 81 / 100 (ψ_norm = 0.9135)
--> Processing surface ipsi = 80 / 100 (ψ_norm = 0.9045)
--> Processing surface ipsi = 79 / 100 (ψ_norm = 0.8951)
--> Processing surface ipsi = 78 / 100 (ψ_norm = 0.8853)
--> Processing surface ipsi = 77 / 100 (ψ_norm = 0.8751)
--> Processing surface ipsi = 76 / 100 (ψ_norm = 0.8645)
--> Processing surface ipsi = 75 / 100 (ψ_norm = 0.8536)
--> Processing surface ipsi = 74 / 100 (ψ_norm = 0.8423)
--> Processing surface ipsi = 73 / 100 (ψ_norm = 0.8307)
--> Processing surface ipsi = 72 / 100 (ψ_norm = 0.8187)
--> Processing surface ipsi = 71 / 100 (ψ_norm = 0.8065)
--> Processing surface ipsi = 70 / 100 (ψ_norm = 0.7939)
--> Processing surface ipsi = 69 / 100 (ψ_norm = 0.7810)
--> Processing surface ipsi = 68 / 100 (ψ_norm = 0.7679)
--> Processing surface ipsi = 67 / 100 (ψ_norm = 0.7545)
--> Processing surface ipsi = 66 / 100 (ψ_norm = 0.7409)
--> Processing surface ipsi = 65 / 100 (ψ_norm = 0.7270)
--> Processing surface ipsi = 64 / 100 (ψ_norm = 0.7129)
--> Processing surface ipsi = 63 / 100 (ψ_norm = 0.6986)
--> Processing surface ipsi = 62 / 100 (ψ_norm = 0.6841)
--> Processing surface ipsi = 61 / 100 (ψ_norm = 0.6694)
--> Processing surface ipsi = 60 / 100 (ψ_norm = 0.6545)
--> Processing surface ipsi = 59 / 100 (ψ_norm = 0.6395)
--> Processing surface ipsi = 58 / 100 (ψ_norm = 0.6243)
--> Processing surface ipsi = 57 / 100 (ψ_norm = 0.6091)
--> Processing surface ipsi = 56 / 100 (ψ_norm = 0.5937)
--> Processing surface ipsi = 55 / 100 (ψ_norm = 0.5782)
--> Processing surface ipsi = 54 / 100 (ψ_norm = 0.5627)
--> Processing surface ipsi = 53 / 100 (ψ_norm = 0.5471)
--> Processing surface ipsi = 52 / 100 (ψ_norm = 0.5314)
--> Processing surface ipsi = 51 / 100 (ψ_norm = 0.5157)
--> Processing surface ipsi = 50 / 100 (ψ_norm = 0.5000)
--> Processing surface ipsi = 49 / 100 (ψ_norm = 0.4843)
--> Processing surface ipsi = 48 / 100 (ψ_norm = 0.4686)
--> Processing surface ipsi = 47 / 100 (ψ_norm = 0.4529)
--> Processing surface ipsi = 46 / 100 (ψ_norm = 0.4373)
--> Processing surface ipsi = 45 / 100 (ψ_norm = 0.4218)
--> Processing surface ipsi = 44 / 100 (ψ_norm = 0.4063)
--> Processing surface ipsi = 43 / 100 (ψ_norm = 0.3909)
--> Processing surface ipsi = 42 / 100 (ψ_norm = 0.3757)
--> Processing surface ipsi = 41 / 100 (ψ_norm = 0.3605)
--> Processing surface ipsi = 40 / 100 (ψ_norm = 0.3455)
--> Processing surface ipsi = 39 / 100 (ψ_norm = 0.3306)
--> Processing surface ipsi = 38 / 100 (ψ_norm = 0.3159)
--> Processing surface ipsi = 37 / 100 (ψ_norm = 0.3014)
--> Processing surface ipsi = 36 / 100 (ψ_norm = 0.2871)
--> Processing surface ipsi = 35 / 100 (ψ_norm = 0.2730)
--> Processing surface ipsi = 34 / 100 (ψ_norm = 0.2591)
--> Processing surface ipsi = 33 / 100 (ψ_norm = 0.2455)
--> Processing surface ipsi = 32 / 100 (ψ_norm = 0.2321)
--> Processing surface ipsi = 31 / 100 (ψ_norm = 0.2190)
--> Processing surface ipsi = 30 / 100 (ψ_norm = 0.2061)
--> Processing surface ipsi = 29 / 100 (ψ_norm = 0.1935)
--> Processing surface ipsi = 28 / 100 (ψ_norm = 0.1813)
--> Processing surface ipsi = 27 / 100 (ψ_norm = 0.1693)
--> Processing surface ipsi = 26 / 100 (ψ_norm = 0.1577)
--> Processing surface ipsi = 25 / 100 (ψ_norm = 0.1464)
--> Processing surface ipsi = 24 / 100 (ψ_norm = 0.1355)
--> Processing surface ipsi = 23 / 100 (ψ_norm = 0.1249)
--> Processing surface ipsi = 22 / 100 (ψ_norm = 0.1147)
--> Processing surface ipsi = 21 / 100 (ψ_norm = 0.1049)
--> Processing surface ipsi = 20 / 100 (ψ_norm = 0.0955)
--> Processing surface ipsi = 19 / 100 (ψ_norm = 0.0865)
--> Processing surface ipsi = 18 / 100 (ψ_norm = 0.0778)
--> Processing surface ipsi = 17 / 100 (ψ_norm = 0.0696)
--> Processing surface ipsi = 16 / 100 (ψ_norm = 0.0618)
--> Processing surface ipsi = 15 / 100 (ψ_norm = 0.0545)
--> Processing surface ipsi = 14 / 100 (ψ_norm = 0.0476)
--> Processing surface ipsi = 13 / 100 (ψ_norm = 0.0411)
--> Processing surface ipsi = 12 / 100 (ψ_norm = 0.0351)
--> Processing surface ipsi = 11 / 100 (ψ_norm = 0.0296)
--> Processing surface ipsi = 10 / 100 (ψ_norm = 0.0245)
--> Processing surface ipsi = 9 / 100 (ψ_norm = 0.0199)
--> Processing surface ipsi = 8 / 100 (ψ_norm = 0.0157)
--> Processing surface ipsi = 7 / 100 (ψ_norm = 0.0120)
--> Processing surface ipsi = 6 / 100 (ψ_norm = 0.0089)
--> Processing surface ipsi = 5 / 100 (ψ_norm = 0.0062)
--> Processing surface ipsi = 4 / 100 (ψ_norm = 0.0039)
--> Processing surface ipsi = 3 / 100 (ψ_norm = 0.0022)
--> Processing surface ipsi = 2 / 100 (ψ_norm = 0.0010)
--> Processing surface ipsi = 1 / 100 (ψ_norm = 0.0002)
--> Processing surface ipsi = 0 / 100 (ψ_norm = 0.0000)
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1934374679157955e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1934374948360171e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1934375193476013e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1934375295502649e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1934375531994367e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1934376554413518e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.193437828039486e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1934380310274918e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1934384499927036e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1934386938553648e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.193439056469854e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.193439292570479e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1934410283937524e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.193441394475206e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1934416090629412e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1934427306779514e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1934428292428542e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1934430110875732e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1934437724816438e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.193445154973529e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1934463939723531e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1934469865584334e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1934472743031836e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1934501517506865e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1934543255631748e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1934617577783578e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1934702707256078e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1934892441181676e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1935017521879349e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.193626832885607e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1948776398623283e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1949124806082059e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1949351145212304e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1950266010298486e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1950798644100508e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1951320173385062e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1954209035931711e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1955979224805278e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1958216255467463e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1978303755935059e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1987165213617355e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.2026989035462863e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.2054460246601643e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.2135980083870523e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.217350825254459e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.2548789939285259e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.3240605433593921e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.3326140727851548e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.3364748500066515e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.3750826222216193e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.4657756332177526e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5289098064000347e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.699876777783192e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.8443803151307222e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=2.0609951463430802e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=2.2866658016625694e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=2.5558367030486857e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=2.893913495299692e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=3.2687398819289046e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=3.537398180795848e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=4.136372026042266e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=4.575859473776159e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=5.2390366396711705e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=5.990106673398998e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=6.326976640354607e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=6.574324394945173e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=7.21793671917911e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=7.92459286422923e-8.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.0552781230352525e-7.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1735371594770229e-7.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.3759064020781498e-7.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5777044773460637e-7.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.8277484316741093e-7.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=2.1036477471458033e-7.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=2.424767003536876e-7.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=2.7768150424954267e-7.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=3.1460411362485093e-7.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=3.603762293230009e-7.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=3.6888324370964935e-7.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=3.7457433983124644e-7.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=4.1741194164380976e-7.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=4.5686811910959415e-7.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=5.126833411821898e-7.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=5.745926707680156e-7.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=6.505858605404238e-7.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=7.363917434448673e-7.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=8.303391391637222e-7.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=9.425258686684986e-7.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.0668247288059904e-6.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.1752833580567294e-6.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.2001084552173785e-6.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.2165527557898307e-6.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.3265168113674344e-6.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.398800072399472e-6.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.577707972511546e-6.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.7296453915946846e-6.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.9308650926585927e-6.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=2.2394346818227485e-6.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=2.504298265669411e-6.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=2.834840418653443e-6.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=3.2418055890493277e-6.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=3.7025580871698877e-6.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=4.242483868100705e-6.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=4.850755782515066e-6.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=5.5066856735766475e-6.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=6.327476352309914e-6.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=7.160498773141204e-6.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=8.271743729612105e-6.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=9.452906822096365e-6.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.0799644079177515e-5.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.2392219498170774e-5.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.3464877891978464e-5.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.4387719623518411e-5.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.651583631622297e-5.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.8187587462191317e-5.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=2.0667578345989457e-5.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=2.339255391097283e-5.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=2.664168169607294e-5.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=3.0136588016150942e-5.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=3.462774578525419e-5.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=3.9579746576569746e-5.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=4.544813721353825e-5.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=5.2047913706281105e-5.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=5.913385425328114e-5.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=6.824712555036012e-5.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=7.81755360878256e-5.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=8.9701040207718e-5.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.00010299097696124239.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.00011851664553937796.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.0001363015354844364.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.0001568335477784417.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.0001802242953210489.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.00020590660708798213.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.00023763007081471096.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.0002729467727211139.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.00031433594127417666.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.00036128294633940537.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.0004109886728530024.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.0004740649000404196.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.0005442333139948664.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.0006266634517413995.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.0007207035639563165.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.0008290585237755077.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.0009500180484805502.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.0010953587423253254.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.0012593362396505579.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.0014500973913959972.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.0016671638277117897.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.0019181833692584259.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.0022032977948149984.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.00253820512566912.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.002918584696195748.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.003361526808547803.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.003868306795682205.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.004449329708968075.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.005117167821384341.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.005862454081521323.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.006758955089663796.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.0077640051164053845.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.00894471881790833.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.010291829291077893.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.011851502616641803.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.013643469834152537.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.0157100555541902.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.018088041104358115.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.018953291176736286.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.019540708128638697.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.022207997714070806.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.024394475588162875.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.02761197262562182.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.03106083614153266.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.035340214167752154.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.04021807420047132.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.04597931562265935.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.05262471884798083.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.06035128525724905.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.06926999972643211.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.07956863638115155.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.09146723815462528.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.10517794784055935.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.12097978481632937.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.13916304705769048.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.16011180572836708.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.18419841809402737.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.21189562599100414.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.24370928040337322.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.2802254187742643.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.32208673807608956.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.36999365183955724.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.42479392444668757.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.4873383111900738.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.5585992803628508.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.6396420219484465.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.7321506797698337.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.8391331421135546.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=0.9766446327001435.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.0947079676827345.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.2072021873341978.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.30587418828741.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.3736304596797844.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.4243530873010282.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.4619512671414754.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.4899254234544093.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.510712319003834.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5261577338448067.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.537621366646649.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5460888913625994.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.552401642746174.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5571170033359831.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5606167993630424.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5632315115118685.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5651768482989477.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5666217717834252.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5676944401433965.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5684871166257095.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5690794079057602.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5695172105883726.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5698427918929447.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5700874979432522.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5702697386244195.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5704050526611875.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5705051897389326.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5705775916496751.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5706333842893387.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.570675274596314.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5707063830328256.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5707294917851171.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5707465660834423.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5707588066500737.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5707683503155072.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5707755397819043.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5707808702707897.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5707847880402162.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.570787504146472.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5707896720774257.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5707912019393502.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5707924268982218.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.570793278202827.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5707940371387843.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5707946257476952.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.570795061706333.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5707953847039517.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: direct_refine did not converge after 10 iterations at eta=1.5707956161105714.
└ @ JPEC.Equilibrium.Direct /Users/nlogan/Code/JPEC_equil/src/Equilibrium/direct.jl:407
┌ Warning: At t=1.5707963267948117, dt was forced below floating point epsilon 2.220446049250313e-16, and step error estimate = 9.29469967593341. Aborting. There is either an error in your model specification or the true solution is unstable (or the true solution can not be represented in the precision of Float64).
└ @ SciMLBase /Users/nlogan/.julia/packages/SciMLBase/DbVzk/src/integrator_interface.jl:623
</details>

@matt-pharr

Copy link
Copy Markdown
Collaborator

What is the status of this PR? It looks like only the last two commits did not make it into #21. Should we be merging this into main or just wait for #21?

@jmlmir369

Copy link
Copy Markdown
Collaborator

Collaborator

What is the status of this PR? It looks like only the last two commits did not make it into #21. Should we be merging this into main or just wait for #21?

Most of the development has been completed except for inverse.f. It's also already been merged with the ode_team branch. So unless there's a need to rewrite the vacuum code in Julia (since Brendan Lyons is currently modernizing the spline package, which is based on the spline used in this code, it would be better to merge the spline from this branch before proceeding with the vacuum code), there's probably no need to do anything at this point. It might be better to merge everything after inverse.jl and the ODE team’s work are completed.

@matt-pharr

matt-pharr commented Aug 5, 2025 via email

Copy link
Copy Markdown
Collaborator

@jaesun57

Copy link
Copy Markdown
Collaborator

I’ve pushed all the work I’ve done so far on converting inverse.f.
Would you prefer that I open a new pull request specifically for this, or would it be better to merge these changes into the current one?

The Jupyter notebook (.ipynb) is not ready yet — @ebursch might be able to work on that.
On the other hand, lar.f is working well as expected.
Let me know what you think!

@matt-pharr

Copy link
Copy Markdown
Collaborator

@satelite2517 Is there more work to be done on EQUIL beyond the jupyter notebook example? It should be pretty quick to make an example right?

@jaesun57

jaesun57 commented Aug 6, 2025

Copy link
Copy Markdown
Collaborator

I think there aren't many works left. Adding some methods to grid_type andipynbexample would be the only thing left.
Or ... error correction? Oh and have to work on some optimization comments that Brendan left.

@logan-nc

Copy link
Copy Markdown
CollaboratorAuthor

I think there aren't many works left. Adding some methods to grid_type andipynbexample would be the only thing left. Or ... error correction? Oh and have to work on some optimization comments that Brendan left.

@satelite2517 I believe this is the only remaining thing this branch is waiting on. I think we should wrap it up for now, since we have the basic functionality needed to run. Future branches can add features with more equilibrium interfaces.

I'd like to merge this branch to main tomorrow, Aug 29 NYC time. Please let me know if you are finished or if I should wait longer.

@matt-pharr

Copy link
Copy Markdown
Collaborator

@satelite2517 @logan-nc What is the status of this branch?

@logan-nc

Copy link
Copy Markdown
CollaboratorAuthor

As per the comment above, I think it was just waiting on @satelite2517 to wrap up a few things in progress.
If not being actively developed, I think we should merge and close the branch.
Remaining work can be done on smaller feature branches

@jhalpern30
jhalpern30 merged commit 1222e34 into mainOct 9, 2025
@matt-pharr
matt-pharr deleted the equil_dev branch October 9, 2025 21:50
Sign up for freeto join this conversation on GitHub. Already have an account? Sign in to comment

Labels

featureNew capability

Projects

None yet

Development

Successfully merging this pull request may close these issues.

Refactor EQUIL

8 participants

@logan-nc@matt-pharr@jmlmir369@jaesun57@bclyons12@ebursch@jhalpern30@d-burg