From efe19783b2b7c6657cc8036387f37ecf81b5fa94 Mon Sep 17 00:00:00 2001 From: Nuwan Yapa Date: Fri, 25 Apr 2025 18:38:43 -0400 Subject: [PATCH] Sqrt fix (still not working) --- calculations/ECCC.jl | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/calculations/ECCC.jl b/calculations/ECCC.jl index 1a9035d..8b4b00a 100644 --- a/calculations/ECCC.jl +++ b/calculations/ECCC.jl @@ -44,7 +44,7 @@ target_c = range(0.78, 0.45, 7) # original: range(0.75, 0.40, 8) target_E = [quick_pole_E(V_system(c)) for c in target_c] # Solve coefficients as a linear system -M_left_element(c, i) = complex(c - c0)^(i/2) +M_left_element(c, i) = alt_sqrt(c - c0)^i M_left = M_left_element.(training_c, (0:order)') a = zeros(ComplexF64, dim, order+1) b = zeros(ComplexF64, dim, order+1) @@ -58,7 +58,7 @@ for i in 1:dim end # Pade approximant -polynomial(a, c) = sum(i -> a[:, i+1] .* complex(c - c0)^(i/2), 0:order) +polynomial(a, c) = sum(i -> a[:, i+1] .* alt_sqrt(c - c0)^i, 0:order) pade_approx(c) = polynomial(a, c) ./ polynomial(b, c) # Extrapolate