Skip to content

Format everything#590

Closed
ChrisRackauckas wants to merge 4 commits into
masterfrom
format
Closed

Format everything#590
ChrisRackauckas wants to merge 4 commits into
masterfrom
format

Conversation

@ChrisRackauckas
Copy link
Copy Markdown
Member

No description provided.

Copy link
Copy Markdown
Contributor

@github-actions github-actions Bot left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Remaining comments which cannot be posted as a review comment to avoid GitHub Rate Limit

JuliaFormatter

[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

@testset "TC: $(nameof(typeof(termination_condition)))" for termination_condition in
TERMINATION_CONDITIONS


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

u0=[-10.0, -1.0, 1.0, 2.0, 3.0, 4.0, 10.0]
p=[0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]
sol=solve_oop(newton_fails, u0, p; solver = TrustRegion())


[JuliaFormatter] reported by reviewdog 🐶

max_trust_radius=[10.0, 100.0, 1000.0]
initial_trust_radius=[10.0, 1.0, 0.1]
step_threshold=[0.0, 0.01, 0.25]
shrink_threshold=[0.25, 0.3, 0.5]
expand_threshold=[0.5, 0.8, 0.9]
shrink_factor=[0.1, 0.3, 0.5]
expand_factor=[1.5, 2.0, 3.0]
max_shrink_times=[10, 20, 30]
list_of_options=zip(


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

sol=solve_oop(quadratic_f, [1.0, 1.0], 2.0; solver = alg)


[JuliaFormatter] reported by reviewdog 🐶

err=maximum(abs, quadratic_f(sol.u, 2.0))


[JuliaFormatter] reported by reviewdog 🐶

maxiterations=[2, 3, 4, 5]
u0=[1.0, 1.0]


[JuliaFormatter] reported by reviewdog 🐶

@testset "TC: $(nameof(typeof(termination_condition)))" for termination_condition in
TERMINATION_CONDITIONS


[JuliaFormatter] reported by reviewdog 🐶

u0=[-10.0, -1.0, 1.0, 2.0, 3.0, 4.0, 10.0]
p=[0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]
sol=solve_oop(newton_fails, u0, p; solver = LevenbergMarquardt())


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

@testset "TC: $(nameof(typeof(termination_condition)))" for termination_condition in
TERMINATION_CONDITIONS


[JuliaFormatter] reported by reviewdog 🐶

damping_initial=[0.5, 2.0, 5.0]
damping_increase_factor=[1.5, 3.0, 10.0]
damping_decrease_factor=Float64[2, 5, 10.0]
finite_diff_step_geodesic=[0.02, 0.2, 0.3]
α_geodesic=[0.6, 0.8, 0.9]
b_uphill=Float64[0, 1, 2]
min_damping_D=[1e-12, 1e-9, 1e-4]
list_of_options=zip(


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

sol=solve_oop(quadratic_f, [1.0, 1.0], 2.0; solver = alg, maxiters = 10000)


[JuliaFormatter] reported by reviewdog 🐶

err=maximum(abs, quadratic_f(sol.u, 2.0))


[JuliaFormatter] reported by reviewdog 🐶

Δ=Tridiagonal(-ones(99), 2*ones(100), -ones(99))
return u+0.1*u .* Δ*u-p


[JuliaFormatter] reported by reviewdog 🐶

Δ=Tridiagonal(-ones(99), 2*ones(100), -ones(99))
du.=u+0.1*u .* Δ*u-p


[JuliaFormatter] reported by reviewdog 🐶

Δ=Tridiagonal(-ones(99), 2*ones(100), -ones(99))
return v+0.1*(u .* Δ*v+v .* Δ*u)


[JuliaFormatter] reported by reviewdog 🐶

Δ=Tridiagonal(-ones(99), 2*ones(100), -ones(99))
du.=v+0.1*(u .* Δ*v+v .* Δ*u)


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

prob=NonlinearProblem(NonlinearFunction{false}(F; jvp = JVP), u0, u0)
sol=solve(prob, NewtonRaphson(; linsolve = KrylovJL_GMRES()); abstol = 1e-13)
err=maximum(abs, sol.resid)


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

prob=NonlinearProblem(NonlinearFunction{true}(F!; jvp = JVP!), u0, u0)
sol=solve(prob, NewtonRaphson(; linsolve = KrylovJL_GMRES()); abstol = 1e-13)
err=maximum(abs, sol.resid)


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

j,
1] = alpha *
(u[im1, j, 1] + u[ip1, j, 1] + u[i, jp1, 1] + u[i, jm1, 1] -
4u[i, j, 1]) +
B +
u[i, j, 1]^2 * u[i, j, 2] - (A + 1) * u[i, j, 1] +
brusselator_f(x, y)


[JuliaFormatter] reported by reviewdog 🐶

j,
2] = alpha *
(u[im1, j, 2] + u[ip1, j, 2] + u[i, jp1, 2] + u[i, jm1, 2] -
4u[i, j, 2]) + A * u[i, j, 1] - u[i, j, 1]^2 * u[i, j, 2]


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

u0s=([1.0, 1.0], @SVector[1.0, 1.0], 1.0)


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

@testset "TC: $(nameof(typeof(termination_condition)))" for termination_condition in
TERMINATION_CONDITIONS


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

@testset "TC: $(nameof(typeof(termination_condition)))" for termination_condition in
TERMINATION_CONDITIONS


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

@testset "TC: $(nameof(typeof(termination_condition)))" for termination_condition in
TERMINATION_CONDITIONS


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

u0s=([1.0, 1.0], @SVector[1.0, 1.0], 1.0)


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

u0=[-10.0, -1.0, 1.0, 2.0, 3.0, 4.0, 10.0]
p=[0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]
sol=solve_oop(newton_fails, u0, p; solver = DFSane())


[JuliaFormatter] reported by reviewdog 🐶

σ_min=[1e-10, 1e-5, 1e-4]
σ_max=[1e10, 1e5, 1e4]
σ_1=[1.0, 0.5, 2.0]
M=[10, 1, 100]
γ=[1e-4, 1e-3, 1e-5]
τ_min=[0.1, 0.2, 0.3]
τ_max=[0.5, 0.8, 0.9]
nexp=[2, 1, 2]
η_strategy=[
(f_1, k, x, F)->f_1/k^2, (f_1, k, x, F)->f_1/k^3,
(f_1, k, x, F)->f_1/k^4


[JuliaFormatter] reported by reviewdog 🐶

list_of_options=zip(σ_min, σ_max, σ_1, M, γ, τ_min, τ_max, nexp, η_strategy)


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

probN=NonlinearProblem{false}(quadratic_f, [1.0, 1.0], 2.0)
sol=solve(probN, alg, abstol = 1e-11)


[JuliaFormatter] reported by reviewdog 🐶

@testset "TC: $(nameof(typeof(termination_condition)))" for termination_condition in
TERMINATION_CONDITIONS


[JuliaFormatter] reported by reviewdog 🐶

du[1]=cos(u[2])-u[1]
du[2]=sin(u[1]+u[2])+u[2]
du[3]=2u[4]+u[3]+1.0
du[4]=u[5]^2+u[4]
du[5]=u[3]^2+u[5]
du[6]=u[1]+u[2]+u[3]+u[4]+u[5]+2.0u[6]+2.5u[7]+1.5u[8]
du[7]=u[1]+u[2]+u[3]+2.0u[4]+u[5]+4.0u[6]-1.5u[7]+1.5u[8]
du[8]=u[1]+2.0u[2]+3.0u[3]+5.0u[4]+6.0u[5]+u[6]-u[7]-u[8]


[JuliaFormatter] reported by reviewdog 🐶

prob=NonlinearProblem(f, zeros(8))
sol=solve(prob, NewtonRaphson())


[JuliaFormatter] reported by reviewdog 🐶

u0=zeros(2)
p=zeros(3)


[JuliaFormatter] reported by reviewdog 🐶

du[1]=cos(u[2])-u[1]
du[2]=sin(u[1]+u[2])+u[2]


[JuliaFormatter] reported by reviewdog 🐶

explicitfun1(p, sols)=nothing
prob1=NonlinearProblem(


[JuliaFormatter] reported by reviewdog 🐶

sol1=solve(prob1, NewtonRaphson())


[JuliaFormatter] reported by reviewdog 🐶

du[1]=2u[2]+u[1]+1.0
du[2]=u[3]^2+u[2]
du[3]=u[1]^2+u[3]


[JuliaFormatter] reported by reviewdog 🐶

explicitfun2(p, sols)=nothing
prob2=NonlinearProblem(


[JuliaFormatter] reported by reviewdog 🐶

sol2=solve(prob2, NewtonRaphson())


[JuliaFormatter] reported by reviewdog 🐶

du[1]=p[1]+2.0u[1]+2.5u[2]+1.5u[3]
du[2]=p[2]+4.0u[1]-1.5u[2]+1.5u[3]
du[3]=p[3]++u[1]-u[2]-u[3]


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

p[1]=sols[1][1]+sols[1][2]+sols[2][1]+sols[2][2]+sols[2][3]
p[2]=sols[1][1]+sols[1][2]+sols[2][1]+2.0sols[2][2]+sols[2][3]
p[3]=sols[1][1]+2.0sols[1][2]+3.0sols[2][1]+5.0sols[2][2]+
6.0sols[2][3]


[JuliaFormatter] reported by reviewdog 🐶

sol3=solve(prob3, NewtonRaphson())
manualscc=[sol1; sol2; sol3]


[JuliaFormatter] reported by reviewdog 🐶

sccprob=SciMLBase.SCCNonlinearProblem([prob1, prob2, prob3],


[JuliaFormatter] reported by reviewdog 🐶

scc_sol=solve(sccprob, NewtonRaphson())


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

@testset "Termination Condition: $(nameof(typeof(termination_condition))) u0: $(nameof(typeof(u0)))" for termination_condition in
TERMINATION_CONDITIONS,


[JuliaFormatter] reported by reviewdog 🐶

@testset "Termination Condition: $(nameof(typeof(termination_condition))) u0: $(nameof(typeof(u0)))" for termination_condition in
TERMINATION_CONDITIONS,


[JuliaFormatter] reported by reviewdog 🐶

@testset "Termination Condition: $(nameof(typeof(termination_condition))) u0: $(nameof(typeof(u0)))" for termination_condition in
TERMINATION_CONDITIONS,


[JuliaFormatter] reported by reviewdog 🐶

u0=[-10.0, -1.0, 1.0, 2.0, 3.0, 4.0, 10.0]
p=[0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]


[JuliaFormatter] reported by reviewdog 🐶

prob=NonlinearProblem(quadratic_f, ones(4), 2.0; maxiters = 2)
sol=solve(prob, SimpleNewtonRaphson())


[JuliaFormatter] reported by reviewdog 🐶

u, p) -> vcat(u .* u .- p, u .* u .- p)),


[JuliaFormatter] reported by reviewdog 🐶

@testset "$idx: $(dict["title"]) | alg #$(alg_id)" for (alg_id, alg) in
enumerate(alg_ops)


[JuliaFormatter] reported by reviewdog 🐶

alg_ops=(RobustMultiNewton(), FastShortcutNonlinearPolyalg())


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[1]]=[]
broken_tests[alg_ops[2]]=[]


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[1]]=[1]


[JuliaFormatter] reported by reviewdog 🐶

alg_ops=(SimpleHalley(; autodiff = AutoForwardDiff()),)


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[1]]=[1, 5, 15, 16, 18]


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[1]]=[11, 21]
broken_tests[alg_ops[2]]=[11, 21]
broken_tests[alg_ops[3]]=[11, 21]
broken_tests[alg_ops[4]]=[8, 11, 21]
broken_tests[alg_ops[5]]=[21]
broken_tests[alg_ops[6]]=[11, 21]
broken_tests[alg_ops[7]]=[3, 15, 16, 21]
broken_tests[alg_ops[8]]=[15, 16]


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[1]]=[11, 21]
broken_tests[alg_ops[2]]=[11, 21]
broken_tests[alg_ops[3]]=[11, 21]


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[1]]=[1, 2, 3, 5, 21]


[JuliaFormatter] reported by reviewdog 🐶

if VERSIONv"1.11-"
broken_tests[alg_ops[2]]=[1, 2, 3, 5, 6, 11, 21]


[JuliaFormatter] reported by reviewdog 🐶

broken_tests[alg_ops[2]]=[1, 2, 3, 5, 6, 21]


[JuliaFormatter] reported by reviewdog 🐶

broken_tests[alg_ops[2]]=[1, 2, 3, 5, 6, 11, 21]


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[2]]=[1, 5, 8, 11, 18]
broken_tests[alg_ops[4]]=[5, 6, 8, 11]


[JuliaFormatter] reported by reviewdog 🐶

broken_tests[alg_ops[1]]=[1, 5, 11]
broken_tests[alg_ops[3]]=[1, 5, 6, 9, 11]
if VERSIONv"1.11-"
broken_tests[alg_ops[5]]=[1, 4, 5, 11]


[JuliaFormatter] reported by reviewdog 🐶

broken_tests[alg_ops[5]]=[1, 5, 11]


[JuliaFormatter] reported by reviewdog 🐶

broken_tests[alg_ops[1]]=[1, 5, 11, 15]
broken_tests[alg_ops[3]]=[1, 5, 6, 9, 11, 16]
broken_tests[alg_ops[5]]=[1, 5, 11]


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[1]]=[1, 2, 4, 5, 11, 18, 22]
broken_tests[alg_ops[2]]=[2, 4, 5, 7, 18, 22]
broken_tests[alg_ops[3]]=[1, 2, 4, 5, 11, 22]


[JuliaFormatter] reported by reviewdog 🐶

alg_ops=(PseudoTransient(),)


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[1]]=[1, 2, 3, 11, 15, 16]


[JuliaFormatter] reported by reviewdog 🐶

compatible(u0, p) || continue


[JuliaFormatter] reported by reviewdog 🐶

nlls_problems=[prob_oop, prob_iip]


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(prob, solver; maxiters = 10000, abstol = 1e-8)


[JuliaFormatter] reported by reviewdog 🐶

resid=zeros(length(p))
ForwardDiff.jacobian!(J, (resid, θ)->loss_function(resid, θ, p), resid, θ)


[JuliaFormatter] reported by reviewdog 🐶

jac(θ, p)=ForwardDiff.jacobian->loss_function(θ, p), θ)


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

solvers=Any[FastLevenbergMarquardtJL(linsolve) for linsolve in (:cholesky, :qr)]
Sys.isapple()||push!(solvers, CMINPACK())


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(prob, solver; maxiters = 10000, abstol = 1e-8)


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(prob, solver; maxiters = 10000, abstol = 1e-8)


[JuliaFormatter] reported by reviewdog 🐶

x_sa=SA[-1.0, -0.5, 0.0, 0.5, 1.0]


[JuliaFormatter] reported by reviewdog 🐶

const y_target_sa=true_function(x_sa, θ_true)


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

θ_init_sa=SVector{4}(θ_init)
prob_sa=NonlinearLeastSquaresProblem{false}(loss_function_sa, θ_init_sa, x)


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(prob_sa, FastLevenbergMarquardtJL())


[JuliaFormatter] reported by reviewdog 🐶

du[1]=2-2u[1]
du[2]=u[1]-4u[2]


[JuliaFormatter] reported by reviewdog 🐶

u0=zeros(2)
prob_iip=SteadyStateProblem(f_iip, u0)


[JuliaFormatter] reported by reviewdog 🐶

f_oop(u, p, t)=[2-2u[1], u[1]-4u[2]]
u0=zeros(2)
prob_oop=SteadyStateProblem(f_oop, u0)


[JuliaFormatter] reported by reviewdog 🐶

du[1]=2-2u[1]
du[2]=u[1]-4u[2]


[JuliaFormatter] reported by reviewdog 🐶

u0=zeros(2)
prob_iip=NonlinearProblem{true}(f_iip, u0)


[JuliaFormatter] reported by reviewdog 🐶

f_oop(u, p)=[2-2u[1], u[1]-4u[2]]
u0=zeros(2)
prob_oop=NonlinearProblem{false}(f_oop, u0)


[JuliaFormatter] reported by reviewdog 🐶

f_tol(u, p)=u^2-2
prob_tol=NonlinearProblem(f_tol, 1.0)


[JuliaFormatter] reported by reviewdog 🐶

alg isa CMINPACK&&Sys.isapple()&&continue
alg isa PETScSNES&&Sys.iswindows()&&continue
sol=solve(prob_tol, alg, abstol = tol)


[JuliaFormatter] reported by reviewdog 🐶

f_jfnk(u, p)=u^2-2
prob_jfnk=NonlinearProblem(f_jfnk, 1.0)


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(prob_jfnk, SIAMFANLEquationsJL(linsolve = :gmres), abstol = tol)


[JuliaFormatter] reported by reviewdog 🐶

fvec[1]=(x[1]+3)*(x[2]^3-7)+18
fvec[2]=sin(x[2]*exp(x[1])-1)


[JuliaFormatter] reported by reviewdog 🐶

prob=NonlinearProblem{true}(f!, [0.1; 1.2])
sol=solve(prob, NLsolveJL(autodiff = :central))


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(prob, SIAMFANLEquationsJL())


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(prob, NLsolveJL(autodiff = :forward))


[JuliaFormatter] reported by reviewdog 🐶

f_custom_jac!(F, u, p)=(F[1:152]=u .^ 2 .- p)
j_custom_jac!(J, u, p)=(J[1:152, 1:152]=diagm(2 .* u))


[JuliaFormatter] reported by reviewdog 🐶

init=ones(152)
A=ones(152)
A[6]=0.8


[JuliaFormatter] reported by reviewdog 🐶

f=NonlinearFunction(f_custom_jac!; jac = j_custom_jac!)
p=A


[JuliaFormatter] reported by reviewdog 🐶

ProbN=NonlinearProblem(f, init, p)


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(ProbN, NLsolveJL(); abstol = 1e-8)


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(ProbN, SIAMFANLEquationsJL(; method = :newton); abstol = 1e-8)


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(ProbN, SIAMFANLEquationsJL(; method = :pseudotransient); abstol = 1e-8)


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(ProbN, PETScSNES(); abstol = 1e-8)


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

u0=[1.0, 1.0]
probN=NonlinearProblem{false}(f, u0)


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(probN, PETScSNES(); abstol = 1e-8)


[JuliaFormatter] reported by reviewdog 🐶

u0=[1.0f0, 1.0f0]
probN=NonlinearProblem{false}(f, u0)


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(probN, PETScSNES(); abstol = 1e-5)


[JuliaFormatter] reported by reviewdog 🐶

u0=Float16[1.0, 1.0]
probN=NonlinearProblem{false}(f, u0)


fj_scalar = @closure (Jx, x) -> begin
fj_scalar = @closure (Jx,
x) -> begin
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
x) -> begin
x) -> begin

Comment on lines +53 to +54
cx,
user_ctx) -> begin
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
cx,
user_ctx) -> begin
cx,
user_ctx) -> begin

Comment on lines +93 to +95
J,
_,
user_ctx) -> begin
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
J,
_,
user_ctx) -> begin
J,
_,
user_ctx) -> begin

Comment on lines +113 to +115
J,
_,
user_ctx) -> begin
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
J,
_,
user_ctx) -> begin
J,
_,
user_ctx) -> begin

Comment on lines +11 to +12
V,
P} = IntervalNonlinearProblem{
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
V,
P} = IntervalNonlinearProblem{
V,
P} = IntervalNonlinearProblem{

using Zygote, Enzyme, ForwardDiff, FiniteDiff

u0s = ([1.0, 1.0], @SVector[1.0, 1.0], 1.0)
u0s=([1.0, 1.0], @SVector[1.0, 1.0], 1.0)
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
u0s=([1.0, 1.0], @SVector[1.0, 1.0], 1.0)
u0s = ([1.0, 1.0], @SVector[1.0, 1.0], 1.0)

Comment on lines +49 to +50
Diagonal(randn!(similar(A, size(A, 1)))), LinearAlgebra.I
)
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
Diagonal(randn!(similar(A, size(A, 1)))), LinearAlgebra.I
)
Diagonal(randn!(similar(A, size(A, 1)))), LinearAlgebra.I
)


@testitem "NewtonRaphson: Iterator Interface" setup=[CoreRootfindTesting] tags=[:core] begin
p = range(0.01, 2, length = 200)
p=range(0.01, 2, length = 200)
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
p=range(0.01, 2, length = 200)
p = range(0.01, 2, length = 200)

Comment on lines +83 to +85
@testset "TC: $(nameof(typeof(termination_condition)))" for termination_condition in
TERMINATION_CONDITIONS

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
@testset "TC: $(nameof(typeof(termination_condition)))" for termination_condition in
TERMINATION_CONDITIONS
@testset "TC: $(nameof(typeof(termination_condition)))" for termination_condition in TERMINATION_CONDITIONS

Comment on lines +101 to +103
preconditioners=[
(u0)->nothing,
u0->((args...)->(Diagonal(rand!(similar(u0))), nothing))
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
preconditioners=[
(u0)->nothing,
u0->((args...)->(Diagonal(rand!(similar(u0))), nothing))
preconditioners = [
(u0) -> nothing,
u0 -> ((args...) -> (Diagonal(rand!(similar(u0))), nothing))

ChrisRackauckas and others added 3 commits April 25, 2025 18:24
Co-authored-by: github-actions[bot] <41898282+github-actions[bot]@users.noreply.github.com>
Co-authored-by: github-actions[bot] <41898282+github-actions[bot]@users.noreply.github.com>
Co-authored-by: github-actions[bot] <41898282+github-actions[bot]@users.noreply.github.com>
Copy link
Copy Markdown
Contributor

@github-actions github-actions Bot left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Remaining comments which cannot be posted as a review comment to avoid GitHub Rate Limit

JuliaFormatter

[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

j,
1] = alpha *
(u[im1, j, 1] + u[ip1, j, 1] + u[i, jp1, 1] + u[i, jm1, 1] -
4u[i, j, 1]) +
B +
u[i, j, 1]^2 * u[i, j, 2] - (A + 1) * u[i, j, 1] +
brusselator_f(x, y)


[JuliaFormatter] reported by reviewdog 🐶

j,
2] = alpha *
(u[im1, j, 2] + u[ip1, j, 2] + u[i, jp1, 2] + u[i, jm1, 2] -
4u[i, j, 2]) + A * u[i, j, 1] - u[i, j, 1]^2 * u[i, j, 2]


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

u0s=([1.0, 1.0], @SVector[1.0, 1.0], 1.0)


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

@testset "TC: $(nameof(typeof(termination_condition)))" for termination_condition in
TERMINATION_CONDITIONS


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

@testset "TC: $(nameof(typeof(termination_condition)))" for termination_condition in
TERMINATION_CONDITIONS


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

@testset "TC: $(nameof(typeof(termination_condition)))" for termination_condition in
TERMINATION_CONDITIONS


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

u0s=([1.0, 1.0], @SVector[1.0, 1.0], 1.0)


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

u0=[-10.0, -1.0, 1.0, 2.0, 3.0, 4.0, 10.0]
p=[0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]
sol=solve_oop(newton_fails, u0, p; solver = DFSane())


[JuliaFormatter] reported by reviewdog 🐶

σ_min=[1e-10, 1e-5, 1e-4]
σ_max=[1e10, 1e5, 1e4]
σ_1=[1.0, 0.5, 2.0]
M=[10, 1, 100]
γ=[1e-4, 1e-3, 1e-5]
τ_min=[0.1, 0.2, 0.3]
τ_max=[0.5, 0.8, 0.9]
nexp=[2, 1, 2]
η_strategy=[
(f_1, k, x, F)->f_1/k^2, (f_1, k, x, F)->f_1/k^3,
(f_1, k, x, F)->f_1/k^4


[JuliaFormatter] reported by reviewdog 🐶

list_of_options=zip(σ_min, σ_max, σ_1, M, γ, τ_min, τ_max, nexp, η_strategy)


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

probN=NonlinearProblem{false}(quadratic_f, [1.0, 1.0], 2.0)
sol=solve(probN, alg, abstol = 1e-11)


[JuliaFormatter] reported by reviewdog 🐶

@testset "TC: $(nameof(typeof(termination_condition)))" for termination_condition in
TERMINATION_CONDITIONS


[JuliaFormatter] reported by reviewdog 🐶

du[1]=cos(u[2])-u[1]
du[2]=sin(u[1]+u[2])+u[2]
du[3]=2u[4]+u[3]+1.0
du[4]=u[5]^2+u[4]
du[5]=u[3]^2+u[5]
du[6]=u[1]+u[2]+u[3]+u[4]+u[5]+2.0u[6]+2.5u[7]+1.5u[8]
du[7]=u[1]+u[2]+u[3]+2.0u[4]+u[5]+4.0u[6]-1.5u[7]+1.5u[8]
du[8]=u[1]+2.0u[2]+3.0u[3]+5.0u[4]+6.0u[5]+u[6]-u[7]-u[8]


[JuliaFormatter] reported by reviewdog 🐶

prob=NonlinearProblem(f, zeros(8))
sol=solve(prob, NewtonRaphson())


[JuliaFormatter] reported by reviewdog 🐶

u0=zeros(2)
p=zeros(3)


[JuliaFormatter] reported by reviewdog 🐶

du[1]=cos(u[2])-u[1]
du[2]=sin(u[1]+u[2])+u[2]


[JuliaFormatter] reported by reviewdog 🐶

explicitfun1(p, sols)=nothing
prob1=NonlinearProblem(


[JuliaFormatter] reported by reviewdog 🐶

sol1=solve(prob1, NewtonRaphson())


[JuliaFormatter] reported by reviewdog 🐶

du[1]=2u[2]+u[1]+1.0
du[2]=u[3]^2+u[2]
du[3]=u[1]^2+u[3]


[JuliaFormatter] reported by reviewdog 🐶

explicitfun2(p, sols)=nothing
prob2=NonlinearProblem(


[JuliaFormatter] reported by reviewdog 🐶

sol2=solve(prob2, NewtonRaphson())


[JuliaFormatter] reported by reviewdog 🐶

du[1]=p[1]+2.0u[1]+2.5u[2]+1.5u[3]
du[2]=p[2]+4.0u[1]-1.5u[2]+1.5u[3]
du[3]=p[3]++u[1]-u[2]-u[3]


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

p[1]=sols[1][1]+sols[1][2]+sols[2][1]+sols[2][2]+sols[2][3]
p[2]=sols[1][1]+sols[1][2]+sols[2][1]+2.0sols[2][2]+sols[2][3]
p[3]=sols[1][1]+2.0sols[1][2]+3.0sols[2][1]+5.0sols[2][2]+
6.0sols[2][3]


[JuliaFormatter] reported by reviewdog 🐶

sol3=solve(prob3, NewtonRaphson())
manualscc=[sol1; sol2; sol3]


[JuliaFormatter] reported by reviewdog 🐶

sccprob=SciMLBase.SCCNonlinearProblem([prob1, prob2, prob3],


[JuliaFormatter] reported by reviewdog 🐶

scc_sol=solve(sccprob, NewtonRaphson())


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

@testset "Termination Condition: $(nameof(typeof(termination_condition))) u0: $(nameof(typeof(u0)))" for termination_condition in
TERMINATION_CONDITIONS,


[JuliaFormatter] reported by reviewdog 🐶

@testset "Termination Condition: $(nameof(typeof(termination_condition))) u0: $(nameof(typeof(u0)))" for termination_condition in
TERMINATION_CONDITIONS,


[JuliaFormatter] reported by reviewdog 🐶

@testset "Termination Condition: $(nameof(typeof(termination_condition))) u0: $(nameof(typeof(u0)))" for termination_condition in
TERMINATION_CONDITIONS,


[JuliaFormatter] reported by reviewdog 🐶

u0=[-10.0, -1.0, 1.0, 2.0, 3.0, 4.0, 10.0]
p=[0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]


[JuliaFormatter] reported by reviewdog 🐶

prob=NonlinearProblem(quadratic_f, ones(4), 2.0; maxiters = 2)
sol=solve(prob, SimpleNewtonRaphson())


[JuliaFormatter] reported by reviewdog 🐶

u, p) -> vcat(u .* u .- p, u .* u .- p)),


[JuliaFormatter] reported by reviewdog 🐶

@testset "$idx: $(dict["title"]) | alg #$(alg_id)" for (alg_id, alg) in
enumerate(alg_ops)


[JuliaFormatter] reported by reviewdog 🐶

alg_ops=(RobustMultiNewton(), FastShortcutNonlinearPolyalg())


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[1]]=[]
broken_tests[alg_ops[2]]=[]


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[1]]=[1]


[JuliaFormatter] reported by reviewdog 🐶

alg_ops=(SimpleHalley(; autodiff = AutoForwardDiff()),)


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[1]]=[1, 5, 15, 16, 18]


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[1]]=[11, 21]
broken_tests[alg_ops[2]]=[11, 21]
broken_tests[alg_ops[3]]=[11, 21]
broken_tests[alg_ops[4]]=[8, 11, 21]
broken_tests[alg_ops[5]]=[21]
broken_tests[alg_ops[6]]=[11, 21]
broken_tests[alg_ops[7]]=[3, 15, 16, 21]
broken_tests[alg_ops[8]]=[15, 16]


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[1]]=[11, 21]
broken_tests[alg_ops[2]]=[11, 21]
broken_tests[alg_ops[3]]=[11, 21]


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[1]]=[1, 2, 3, 5, 21]


[JuliaFormatter] reported by reviewdog 🐶

if VERSIONv"1.11-"
broken_tests[alg_ops[2]]=[1, 2, 3, 5, 6, 11, 21]


[JuliaFormatter] reported by reviewdog 🐶

broken_tests[alg_ops[2]]=[1, 2, 3, 5, 6, 21]


[JuliaFormatter] reported by reviewdog 🐶

broken_tests[alg_ops[2]]=[1, 2, 3, 5, 6, 11, 21]


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[2]]=[1, 5, 8, 11, 18]
broken_tests[alg_ops[4]]=[5, 6, 8, 11]


[JuliaFormatter] reported by reviewdog 🐶

broken_tests[alg_ops[1]]=[1, 5, 11]
broken_tests[alg_ops[3]]=[1, 5, 6, 9, 11]
if VERSIONv"1.11-"
broken_tests[alg_ops[5]]=[1, 4, 5, 11]


[JuliaFormatter] reported by reviewdog 🐶

broken_tests[alg_ops[5]]=[1, 5, 11]


[JuliaFormatter] reported by reviewdog 🐶

broken_tests[alg_ops[1]]=[1, 5, 11, 15]
broken_tests[alg_ops[3]]=[1, 5, 6, 9, 11, 16]
broken_tests[alg_ops[5]]=[1, 5, 11]


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[1]]=[1, 2, 4, 5, 11, 18, 22]
broken_tests[alg_ops[2]]=[2, 4, 5, 7, 18, 22]
broken_tests[alg_ops[3]]=[1, 2, 4, 5, 11, 22]


[JuliaFormatter] reported by reviewdog 🐶

alg_ops=(PseudoTransient(),)


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[1]]=[1, 2, 3, 11, 15, 16]


[JuliaFormatter] reported by reviewdog 🐶

compatible(u0, p) || continue


[JuliaFormatter] reported by reviewdog 🐶

nlls_problems=[prob_oop, prob_iip]


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(prob, solver; maxiters = 10000, abstol = 1e-8)


[JuliaFormatter] reported by reviewdog 🐶

resid=zeros(length(p))
ForwardDiff.jacobian!(J, (resid, θ)->loss_function(resid, θ, p), resid, θ)


[JuliaFormatter] reported by reviewdog 🐶

jac(θ, p)=ForwardDiff.jacobian->loss_function(θ, p), θ)


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

solvers=Any[FastLevenbergMarquardtJL(linsolve) for linsolve in (:cholesky, :qr)]
Sys.isapple()||push!(solvers, CMINPACK())


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(prob, solver; maxiters = 10000, abstol = 1e-8)


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(prob, solver; maxiters = 10000, abstol = 1e-8)


[JuliaFormatter] reported by reviewdog 🐶

x_sa=SA[-1.0, -0.5, 0.0, 0.5, 1.0]


[JuliaFormatter] reported by reviewdog 🐶

const y_target_sa=true_function(x_sa, θ_true)


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

θ_init_sa=SVector{4}(θ_init)
prob_sa=NonlinearLeastSquaresProblem{false}(loss_function_sa, θ_init_sa, x)


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(prob_sa, FastLevenbergMarquardtJL())


[JuliaFormatter] reported by reviewdog 🐶

du[1]=2-2u[1]
du[2]=u[1]-4u[2]


[JuliaFormatter] reported by reviewdog 🐶

u0=zeros(2)
prob_iip=SteadyStateProblem(f_iip, u0)


[JuliaFormatter] reported by reviewdog 🐶

f_oop(u, p, t)=[2-2u[1], u[1]-4u[2]]
u0=zeros(2)
prob_oop=SteadyStateProblem(f_oop, u0)


[JuliaFormatter] reported by reviewdog 🐶

du[1]=2-2u[1]
du[2]=u[1]-4u[2]


[JuliaFormatter] reported by reviewdog 🐶

u0=zeros(2)
prob_iip=NonlinearProblem{true}(f_iip, u0)


[JuliaFormatter] reported by reviewdog 🐶

f_oop(u, p)=[2-2u[1], u[1]-4u[2]]
u0=zeros(2)
prob_oop=NonlinearProblem{false}(f_oop, u0)


[JuliaFormatter] reported by reviewdog 🐶

f_tol(u, p)=u^2-2
prob_tol=NonlinearProblem(f_tol, 1.0)


[JuliaFormatter] reported by reviewdog 🐶

alg isa CMINPACK&&Sys.isapple()&&continue
alg isa PETScSNES&&Sys.iswindows()&&continue
sol=solve(prob_tol, alg, abstol = tol)


[JuliaFormatter] reported by reviewdog 🐶

f_jfnk(u, p)=u^2-2
prob_jfnk=NonlinearProblem(f_jfnk, 1.0)


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(prob_jfnk, SIAMFANLEquationsJL(linsolve = :gmres), abstol = tol)


[JuliaFormatter] reported by reviewdog 🐶

fvec[1]=(x[1]+3)*(x[2]^3-7)+18
fvec[2]=sin(x[2]*exp(x[1])-1)


[JuliaFormatter] reported by reviewdog 🐶

prob=NonlinearProblem{true}(f!, [0.1; 1.2])
sol=solve(prob, NLsolveJL(autodiff = :central))


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(prob, SIAMFANLEquationsJL())


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(prob, NLsolveJL(autodiff = :forward))


[JuliaFormatter] reported by reviewdog 🐶

f_custom_jac!(F, u, p)=(F[1:152]=u .^ 2 .- p)
j_custom_jac!(J, u, p)=(J[1:152, 1:152]=diagm(2 .* u))


[JuliaFormatter] reported by reviewdog 🐶

init=ones(152)
A=ones(152)
A[6]=0.8


[JuliaFormatter] reported by reviewdog 🐶

f=NonlinearFunction(f_custom_jac!; jac = j_custom_jac!)
p=A


[JuliaFormatter] reported by reviewdog 🐶

ProbN=NonlinearProblem(f, init, p)


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(ProbN, NLsolveJL(); abstol = 1e-8)


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(ProbN, SIAMFANLEquationsJL(; method = :newton); abstol = 1e-8)


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(ProbN, SIAMFANLEquationsJL(; method = :pseudotransient); abstol = 1e-8)


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(ProbN, PETScSNES(); abstol = 1e-8)


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

u0=[1.0, 1.0]
probN=NonlinearProblem{false}(f, u0)


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(probN, PETScSNES(); abstol = 1e-8)


[JuliaFormatter] reported by reviewdog 🐶

u0=[1.0f0, 1.0f0]
probN=NonlinearProblem{false}(f, u0)


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(probN, PETScSNES(); abstol = 1e-5)


[JuliaFormatter] reported by reviewdog 🐶

u0=Float16[1.0, 1.0]
probN=NonlinearProblem{false}(f, u0)


@testitem "PseudoTransient: Iterator Interface" setup=[CoreRootfindTesting] tags=[:core] begin
p = range(0.01, 2, length = 200)
p=range(0.01, 2, length = 200)
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
p=range(0.01, 2, length = 200)
p = range(0.01, 2, length = 200)

Comment on lines +169 to +171
@testset "TC: $(nameof(typeof(termination_condition)))" for termination_condition in
TERMINATION_CONDITIONS

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
@testset "TC: $(nameof(typeof(termination_condition)))" for termination_condition in
TERMINATION_CONDITIONS
@testset "TC: $(nameof(typeof(termination_condition)))" for termination_condition in TERMINATION_CONDITIONS

using Zygote, Enzyme, ForwardDiff, FiniteDiff

radius_update_schemes = [
radius_update_schemes=[
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
radius_update_schemes=[
radius_update_schemes = [


@testitem "TrustRegion: Iterator Interface" setup=[CoreRootfindTesting] tags=[:core] begin
p = range(0.01, 2, length = 200)
p=range(0.01, 2, length = 200)
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
p=range(0.01, 2, length = 200)
p = range(0.01, 2, length = 200)

Comment on lines +234 to +236
u0=[-10.0, -1.0, 1.0, 2.0, 3.0, 4.0, 10.0]
p=[0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]
sol=solve_oop(newton_fails, u0, p; solver = TrustRegion())
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
u0=[-10.0, -1.0, 1.0, 2.0, 3.0, 4.0, 10.0]
p=[0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]
sol=solve_oop(newton_fails, u0, p; solver = TrustRegion())
u0 = [-10.0, -1.0, 1.0, 2.0, 3.0, 4.0, 10.0]
p = [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]
sol = solve_oop(newton_fails, u0, p; solver = TrustRegion())

Comment on lines +457 to +459
prob=NonlinearProblem(NonlinearFunction{false}(F; jvp = JVP), u0, u0)
sol=solve(prob, NewtonRaphson(; linsolve = KrylovJL_GMRES()); abstol = 1e-13)
err=maximum(abs, sol.resid)
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
prob=NonlinearProblem(NonlinearFunction{false}(F; jvp = JVP), u0, u0)
sol=solve(prob, NewtonRaphson(; linsolve = KrylovJL_GMRES()); abstol = 1e-13)
err=maximum(abs, sol.resid)
prob = NonlinearProblem(NonlinearFunction{false}(F; jvp = JVP), u0, u0)
sol = solve(prob, NewtonRaphson(; linsolve = KrylovJL_GMRES()); abstol = 1e-13)
err = maximum(abs, sol.resid)

@test err < 1e-6

sol = solve(
sol=solve(
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
sol=solve(
sol = solve(

abstol = 1e-13
)
err = maximum(abs, sol.resid)
err=maximum(abs, sol.resid)
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
err=maximum(abs, sol.resid)
err = maximum(abs, sol.resid)

Comment on lines +469 to +471
prob=NonlinearProblem(NonlinearFunction{true}(F!; jvp = JVP!), u0, u0)
sol=solve(prob, NewtonRaphson(; linsolve = KrylovJL_GMRES()); abstol = 1e-13)
err=maximum(abs, sol.resid)
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
prob=NonlinearProblem(NonlinearFunction{true}(F!; jvp = JVP!), u0, u0)
sol=solve(prob, NewtonRaphson(; linsolve = KrylovJL_GMRES()); abstol = 1e-13)
err=maximum(abs, sol.resid)
prob = NonlinearProblem(NonlinearFunction{true}(F!; jvp = JVP!), u0, u0)
sol = solve(prob, NewtonRaphson(; linsolve = KrylovJL_GMRES()); abstol = 1e-13)
err = maximum(abs, sol.resid)

@test err < 1e-6

sol = solve(
sol=solve(
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
sol=solve(
sol = solve(

Copy link
Copy Markdown
Contributor

@github-actions github-actions Bot left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Remaining comments which cannot be posted as a review comment to avoid GitHub Rate Limit

JuliaFormatter

[JuliaFormatter] reported by reviewdog 🐶

sol2=solve(prob2, NewtonRaphson())


[JuliaFormatter] reported by reviewdog 🐶

du[1]=p[1]+2.0u[1]+2.5u[2]+1.5u[3]
du[2]=p[2]+4.0u[1]-1.5u[2]+1.5u[3]
du[3]=p[3]++u[1]-u[2]-u[3]


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

p[1]=sols[1][1]+sols[1][2]+sols[2][1]+sols[2][2]+sols[2][3]
p[2]=sols[1][1]+sols[1][2]+sols[2][1]+2.0sols[2][2]+sols[2][3]
p[3]=sols[1][1]+2.0sols[1][2]+3.0sols[2][1]+5.0sols[2][2]+
6.0sols[2][3]


[JuliaFormatter] reported by reviewdog 🐶

sol3=solve(prob3, NewtonRaphson())
manualscc=[sol1; sol2; sol3]


[JuliaFormatter] reported by reviewdog 🐶

sccprob=SciMLBase.SCCNonlinearProblem([prob1, prob2, prob3],


[JuliaFormatter] reported by reviewdog 🐶

scc_sol=solve(sccprob, NewtonRaphson())


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

@testset "Termination Condition: $(nameof(typeof(termination_condition))) u0: $(nameof(typeof(u0)))" for termination_condition in
TERMINATION_CONDITIONS,


[JuliaFormatter] reported by reviewdog 🐶

@testset "Termination Condition: $(nameof(typeof(termination_condition))) u0: $(nameof(typeof(u0)))" for termination_condition in
TERMINATION_CONDITIONS,


[JuliaFormatter] reported by reviewdog 🐶

@testset "Termination Condition: $(nameof(typeof(termination_condition))) u0: $(nameof(typeof(u0)))" for termination_condition in
TERMINATION_CONDITIONS,


[JuliaFormatter] reported by reviewdog 🐶

u0=[-10.0, -1.0, 1.0, 2.0, 3.0, 4.0, 10.0]
p=[0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]


[JuliaFormatter] reported by reviewdog 🐶

prob=NonlinearProblem(quadratic_f, ones(4), 2.0; maxiters = 2)
sol=solve(prob, SimpleNewtonRaphson())


[JuliaFormatter] reported by reviewdog 🐶

u, p) -> vcat(u .* u .- p, u .* u .- p)),


[JuliaFormatter] reported by reviewdog 🐶

@testset "$idx: $(dict["title"]) | alg #$(alg_id)" for (alg_id, alg) in
enumerate(alg_ops)


[JuliaFormatter] reported by reviewdog 🐶

alg_ops=(RobustMultiNewton(), FastShortcutNonlinearPolyalg())


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[1]]=[]
broken_tests[alg_ops[2]]=[]


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[1]]=[1]


[JuliaFormatter] reported by reviewdog 🐶

alg_ops=(SimpleHalley(; autodiff = AutoForwardDiff()),)


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[1]]=[1, 5, 15, 16, 18]


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[1]]=[11, 21]
broken_tests[alg_ops[2]]=[11, 21]
broken_tests[alg_ops[3]]=[11, 21]
broken_tests[alg_ops[4]]=[8, 11, 21]
broken_tests[alg_ops[5]]=[21]
broken_tests[alg_ops[6]]=[11, 21]
broken_tests[alg_ops[7]]=[3, 15, 16, 21]
broken_tests[alg_ops[8]]=[15, 16]


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[1]]=[11, 21]
broken_tests[alg_ops[2]]=[11, 21]
broken_tests[alg_ops[3]]=[11, 21]


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[1]]=[1, 2, 3, 5, 21]


[JuliaFormatter] reported by reviewdog 🐶

if VERSIONv"1.11-"
broken_tests[alg_ops[2]]=[1, 2, 3, 5, 6, 11, 21]


[JuliaFormatter] reported by reviewdog 🐶

broken_tests[alg_ops[2]]=[1, 2, 3, 5, 6, 21]


[JuliaFormatter] reported by reviewdog 🐶

broken_tests[alg_ops[2]]=[1, 2, 3, 5, 6, 11, 21]


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[2]]=[1, 5, 8, 11, 18]
broken_tests[alg_ops[4]]=[5, 6, 8, 11]


[JuliaFormatter] reported by reviewdog 🐶

broken_tests[alg_ops[1]]=[1, 5, 11]
broken_tests[alg_ops[3]]=[1, 5, 6, 9, 11]
if VERSIONv"1.11-"
broken_tests[alg_ops[5]]=[1, 4, 5, 11]


[JuliaFormatter] reported by reviewdog 🐶

broken_tests[alg_ops[5]]=[1, 5, 11]


[JuliaFormatter] reported by reviewdog 🐶

broken_tests[alg_ops[1]]=[1, 5, 11, 15]
broken_tests[alg_ops[3]]=[1, 5, 6, 9, 11, 16]
broken_tests[alg_ops[5]]=[1, 5, 11]


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[1]]=[1, 2, 4, 5, 11, 18, 22]
broken_tests[alg_ops[2]]=[2, 4, 5, 7, 18, 22]
broken_tests[alg_ops[3]]=[1, 2, 4, 5, 11, 22]


[JuliaFormatter] reported by reviewdog 🐶

alg_ops=(PseudoTransient(),)


[JuliaFormatter] reported by reviewdog 🐶

broken_tests=Dict(alg=>Int[] for alg in alg_ops)
broken_tests[alg_ops[1]]=[1, 2, 3, 11, 15, 16]


[JuliaFormatter] reported by reviewdog 🐶

compatible(u0, p) || continue


[JuliaFormatter] reported by reviewdog 🐶

nlls_problems=[prob_oop, prob_iip]


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(prob, solver; maxiters = 10000, abstol = 1e-8)


[JuliaFormatter] reported by reviewdog 🐶

resid=zeros(length(p))
ForwardDiff.jacobian!(J, (resid, θ)->loss_function(resid, θ, p), resid, θ)


[JuliaFormatter] reported by reviewdog 🐶

jac(θ, p)=ForwardDiff.jacobian->loss_function(θ, p), θ)


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

solvers=Any[FastLevenbergMarquardtJL(linsolve) for linsolve in (:cholesky, :qr)]
Sys.isapple()||push!(solvers, CMINPACK())


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(prob, solver; maxiters = 10000, abstol = 1e-8)


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(prob, solver; maxiters = 10000, abstol = 1e-8)


[JuliaFormatter] reported by reviewdog 🐶

x_sa=SA[-1.0, -0.5, 0.0, 0.5, 1.0]


[JuliaFormatter] reported by reviewdog 🐶

const y_target_sa=true_function(x_sa, θ_true)


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

θ_init_sa=SVector{4}(θ_init)
prob_sa=NonlinearLeastSquaresProblem{false}(loss_function_sa, θ_init_sa, x)


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(prob_sa, FastLevenbergMarquardtJL())


[JuliaFormatter] reported by reviewdog 🐶

du[1]=2-2u[1]
du[2]=u[1]-4u[2]


[JuliaFormatter] reported by reviewdog 🐶

u0=zeros(2)
prob_iip=SteadyStateProblem(f_iip, u0)


[JuliaFormatter] reported by reviewdog 🐶

f_oop(u, p, t)=[2-2u[1], u[1]-4u[2]]
u0=zeros(2)
prob_oop=SteadyStateProblem(f_oop, u0)


[JuliaFormatter] reported by reviewdog 🐶

du[1]=2-2u[1]
du[2]=u[1]-4u[2]


[JuliaFormatter] reported by reviewdog 🐶

u0=zeros(2)
prob_iip=NonlinearProblem{true}(f_iip, u0)


[JuliaFormatter] reported by reviewdog 🐶

f_oop(u, p)=[2-2u[1], u[1]-4u[2]]
u0=zeros(2)
prob_oop=NonlinearProblem{false}(f_oop, u0)


[JuliaFormatter] reported by reviewdog 🐶

f_tol(u, p)=u^2-2
prob_tol=NonlinearProblem(f_tol, 1.0)


[JuliaFormatter] reported by reviewdog 🐶

alg isa CMINPACK&&Sys.isapple()&&continue
alg isa PETScSNES&&Sys.iswindows()&&continue
sol=solve(prob_tol, alg, abstol = tol)


[JuliaFormatter] reported by reviewdog 🐶

f_jfnk(u, p)=u^2-2
prob_jfnk=NonlinearProblem(f_jfnk, 1.0)


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(prob_jfnk, SIAMFANLEquationsJL(linsolve = :gmres), abstol = tol)


[JuliaFormatter] reported by reviewdog 🐶

fvec[1]=(x[1]+3)*(x[2]^3-7)+18
fvec[2]=sin(x[2]*exp(x[1])-1)


[JuliaFormatter] reported by reviewdog 🐶

prob=NonlinearProblem{true}(f!, [0.1; 1.2])
sol=solve(prob, NLsolveJL(autodiff = :central))


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(prob, SIAMFANLEquationsJL())


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(prob, NLsolveJL(autodiff = :forward))


[JuliaFormatter] reported by reviewdog 🐶

f_custom_jac!(F, u, p)=(F[1:152]=u .^ 2 .- p)
j_custom_jac!(J, u, p)=(J[1:152, 1:152]=diagm(2 .* u))


[JuliaFormatter] reported by reviewdog 🐶

init=ones(152)
A=ones(152)
A[6]=0.8


[JuliaFormatter] reported by reviewdog 🐶

f=NonlinearFunction(f_custom_jac!; jac = j_custom_jac!)
p=A


[JuliaFormatter] reported by reviewdog 🐶

ProbN=NonlinearProblem(f, init, p)


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(ProbN, NLsolveJL(); abstol = 1e-8)


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(ProbN, SIAMFANLEquationsJL(; method = :newton); abstol = 1e-8)


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(ProbN, SIAMFANLEquationsJL(; method = :pseudotransient); abstol = 1e-8)


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(ProbN, PETScSNES(); abstol = 1e-8)


[JuliaFormatter] reported by reviewdog 🐶


[JuliaFormatter] reported by reviewdog 🐶

u0=[1.0, 1.0]
probN=NonlinearProblem{false}(f, u0)


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(probN, PETScSNES(); abstol = 1e-8)


[JuliaFormatter] reported by reviewdog 🐶

u0=[1.0f0, 1.0f0]
probN=NonlinearProblem{false}(f, u0)


[JuliaFormatter] reported by reviewdog 🐶

sol=solve(probN, PETScSNES(); abstol = 1e-5)


[JuliaFormatter] reported by reviewdog 🐶

u0=Float16[1.0, 1.0]
probN=NonlinearProblem{false}(f, u0)

Comment thread docs/src/tutorials/large_systems.md Outdated
Comment on lines +87 to +89
j,
2] = alpha * (u[im1, j, 2] + u[ip1, j, 2] + u[i, jp1, 2] + u[i, jm1, 2] -
4u[i, j, 2]) + A * u[i, j, 1] - u[i, j, 1]^2 * u[i, j, 2]
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
j,
2] = alpha * (u[im1, j, 2] + u[ip1, j, 2] + u[i, jp1, 2] + u[i, jm1, 2] -
4u[i, j, 2]) + A * u[i, j, 1] - u[i, j, 1]^2 * u[i, j, 2]
j,
2] = alpha * (u[im1, j, 2] + u[ip1, j, 2] + u[i, jp1, 2] + u[i, jm1, 2] -
4u[i, j, 2]) + A * u[i, j, 1] - u[i, j, 1]^2 * u[i, j, 2]

abstol = 1e-13
)
err = maximum(abs, sol.resid)
err=maximum(abs, sol.resid)
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
err=maximum(abs, sol.resid)
err = maximum(abs, sol.resid)

Comment on lines +21 to +27
j,
1] = alpha *
(u[im1, j, 1] + u[ip1, j, 1] + u[i, jp1, 1] + u[i, jm1, 1] -
4u[i, j, 1]) +
B +
u[i, j, 1]^2 * u[i, j, 2] - (A + 1) * u[i, j, 1] +
brusselator_f(x, y)
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
j,
1] = alpha *
(u[im1, j, 1] + u[ip1, j, 1] + u[i, jp1, 1] + u[i, jm1, 1] -
4u[i, j, 1]) +
B +
u[i, j, 1]^2 * u[i, j, 2] - (A + 1) * u[i, j, 1] +
brusselator_f(x, y)
j,
1] = alpha *
(u[im1, j, 1] + u[ip1, j, 1] + u[i, jp1, 1] + u[i, jm1, 1] -
4u[i, j, 1]) +
B +
u[i, j, 1]^2 * u[i, j, 2] - (A + 1) * u[i, j, 1] +
brusselator_f(x, y)

Comment on lines +29 to +32
j,
2] = alpha *
(u[im1, j, 2] + u[ip1, j, 2] + u[i, jp1, 2] + u[i, jm1, 2] -
4u[i, j, 2]) + A * u[i, j, 1] - u[i, j, 1]^2 * u[i, j, 2]
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
j,
2] = alpha *
(u[im1, j, 2] + u[ip1, j, 2] + u[i, jp1, 2] + u[i, jm1, 2] -
4u[i, j, 2]) + A * u[i, j, 1] - u[i, j, 1]^2 * u[i, j, 2]
j,
2] = alpha *
(u[im1, j, 2] + u[ip1, j, 2] + u[i, jp1, 2] + u[i, jm1, 2] -
4u[i, j, 2]) + A * u[i, j, 1] - u[i, j, 1]^2 * u[i, j, 2]


@compile_workload begin
@sync for prob in nonlinear_problems, alg in algs

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change

Comment on lines +26 to +27
du[1]=cos(u[2])-u[1]
du[2]=sin(u[1]+u[2])+u[2]
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
du[1]=cos(u[2])-u[1]
du[2]=sin(u[1]+u[2])+u[2]
du[1] = cos(u[2]) - u[1]
du[2] = sin(u[1] + u[2]) + u[2]

Comment on lines +29 to +30
explicitfun1(p, sols)=nothing
prob1=NonlinearProblem(
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
explicitfun1(p, sols)=nothing
prob1=NonlinearProblem(
explicitfun1(p, sols) = nothing
prob1 = NonlinearProblem(

prob1=NonlinearProblem(
NonlinearFunction{true, SciMLBase.NoSpecialize}(f1), zeros(2), p)
sol1 = solve(prob1, NewtonRaphson())
sol1=solve(prob1, NewtonRaphson())
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
sol1=solve(prob1, NewtonRaphson())
sol1 = solve(prob1, NewtonRaphson())

Comment on lines +35 to +37
du[1]=2u[2]+u[1]+1.0
du[2]=u[3]^2+u[2]
du[3]=u[1]^2+u[3]
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
du[1]=2u[2]+u[1]+1.0
du[2]=u[3]^2+u[2]
du[3]=u[1]^2+u[3]
du[1] = 2u[2] + u[1] + 1.0
du[2] = u[3]^2 + u[2]
du[3] = u[1]^2 + u[3]

Comment on lines +39 to +40
explicitfun2(p, sols)=nothing
prob2=NonlinearProblem(
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[JuliaFormatter] reported by reviewdog 🐶

Suggested change
explicitfun2(p, sols)=nothing
prob2=NonlinearProblem(
explicitfun2(p, sols) = nothing
prob2 = NonlinearProblem(

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant