diff --git a/docs/src/plots/P3ShapeSolverPlots.jl b/docs/src/plots/P3ShapeSolverPlots.jl index 9b6ac8f92f..878b4b47b9 100644 --- a/docs/src/plots/P3ShapeSolverPlots.jl +++ b/docs/src/plots/P3ShapeSolverPlots.jl @@ -10,7 +10,7 @@ const PSP3 = CMP.ParametersP3 p3 = CMP.ParametersP3(FT) function guess_value(λ::FT, p1::FT, p2::FT, q1::FT, q2::FT) - return q1 * (λ/p1)^((log(q1) - log(q2)) / (log(p1) - log(p2))) + return q1 * (λ / p1)^((log(q1) - log(q2)) / (log(p1) - log(p2))) end function lambda_guess_plot(F_r::FT, ρ_r::FT) where {FT} diff --git a/src/P3Scheme.jl b/src/P3Scheme.jl index f8e2d51a3d..b682ebfbcf 100644 --- a/src/P3Scheme.jl +++ b/src/P3Scheme.jl @@ -347,7 +347,8 @@ function get_bounds( ql = q_gamma(p3, F_r, N / N̂, log(left), th) qr = q_gamma(p3, F_r, N / N̂, log(right), th) - guess = left * (q / (N̂ * ql))^((log(right) - log(left)) / (log(qr) - log(ql))) + guess = + left * (q / (N̂ * ql))^((log(right) - log(left)) / (log(qr) - log(ql))) max = log(guess * exp(radius)) min = log(guess)