Skip to content

Commit 88be6f7

Browse files
fixup! fixup! test: test initialization on NonlinearProblem and NonlinearLeastSquaresProblem
1 parent c7bc9ac commit 88be6f7

File tree

1 file changed

+14
-7
lines changed

1 file changed

+14
-7
lines changed

test/initializationsystem.jl

+14-7
Original file line numberDiff line numberDiff line change
@@ -957,33 +957,38 @@ end
957957
@mtkbuild sys = System(
958958
[D(x) ~ x * p + y^2 * q + rhss; alge_eqs],
959959
t; guesses = [x => 0.0, y => 0.0, z => 0.0, p => 0.0, q => 0.0])
960-
prob = Problem(sys, [x => 1.0], (0.0, 1.0), [p => 1.0, q => missing]; initialization_eqs = is_nlsolve ? alge_eqs : [])
960+
prob = Problem(sys, [x => 1.0], (0.0, 1.0), [p => 1.0, q => missing];
961+
initialization_eqs = is_nlsolve ? alge_eqs : [])
961962
@test is_variable(prob.f.initialization_data.initializeprob, q)
962963
ps = prob.p
963964
newps = SciMLStructures.replace(Tunable(), ps, ForwardDiff.Dual.(ps.tunable))
964965
prob2 = remake(prob; p = newps)
965-
@test eltype(state_values(prob2.f.initialization_data.initializeprob)) <: ForwardDiff.Dual
966+
@test eltype(state_values(prob2.f.initialization_data.initializeprob)) <:
967+
ForwardDiff.Dual
966968
@test eltype(prob2.f.initialization_data.initializeprob.p.tunable) <:
967969
ForwardDiff.Dual
968970
@test state_values(prob2.f.initialization_data.initializeprob)
969971
state_values(prob.f.initialization_data.initializeprob)
970972

971973
prob2 = remake(prob; u0 = ForwardDiff.Dual.(prob.u0))
972-
@test eltype(state_values(prob2.f.initialization_data.initializeprob)) <: ForwardDiff.Dual
974+
@test eltype(state_values(prob2.f.initialization_data.initializeprob)) <:
975+
ForwardDiff.Dual
973976
@test eltype(prob2.f.initialization_data.initializeprob.p.tunable) <: Float64
974977
@test state_values(prob2.f.initialization_data.initializeprob)
975978
state_values(prob.f.initialization_data.initializeprob)
976979

977980
prob2 = remake(prob; u0 = ForwardDiff.Dual.(prob.u0), p = newps)
978-
@test eltype(state_values(prob2.f.initialization_data.initializeprob)) <: ForwardDiff.Dual
981+
@test eltype(state_values(prob2.f.initialization_data.initializeprob)) <:
982+
ForwardDiff.Dual
979983
@test eltype(prob2.f.initialization_data.initializeprob.p.tunable) <:
980984
ForwardDiff.Dual
981985
@test state_values(prob2.f.initialization_data.initializeprob)
982986
state_values(prob.f.initialization_data.initializeprob)
983987

984988
prob2 = remake(prob; u0 = [x => ForwardDiff.Dual(1.0)],
985989
p = [p => ForwardDiff.Dual(1.0), q => missing])
986-
@test eltype(state_values(prob2.f.initialization_data.initializeprob)) <: ForwardDiff.Dual
990+
@test eltype(state_values(prob2.f.initialization_data.initializeprob)) <:
991+
ForwardDiff.Dual
987992
@test eltype(prob2.f.initialization_data.initializeprob.p.tunable) <:
988993
ForwardDiff.Dual
989994
@test state_values(prob2.f.initialization_data.initializeprob)
@@ -1027,12 +1032,14 @@ end
10271032
@mtkbuild sys = System(
10281033
[D(x) ~ x + p * y^2 + rhss; alge_eqs], t; guesses = [
10291034
y => 1.0, p => 1.0])
1030-
prob = Problem(sys, [x => 1.0], (0.0, 1.0), [p => 1.0]; initialization_eqs = is_nlsolve ? alge_eqs : [])
1035+
prob = Problem(sys, [x => 1.0], (0.0, 1.0), [p => 1.0];
1036+
initialization_eqs = is_nlsolve ? alge_eqs : [])
10311037
@test is_variable(prob.f.initialization_data.initializeprob, y)
10321038
prob2 = @test_nowarn remake(prob; p = [p => 3.0]) # ensure no over/under-determined warning
10331039
@test is_variable(prob.f.initialization_data.initializeprob, y)
10341040

1035-
prob = Problem(sys, [y => 1.0, x => 2.0], (0.0, 1.0), [p => missing]; initialization_eqs = is_nlsolve ? alge_eqs : [])
1041+
prob = Problem(sys, [y => 1.0, x => 2.0], (0.0, 1.0), [p => missing];
1042+
initialization_eqs = is_nlsolve ? alge_eqs : [])
10361043
@test is_variable(prob.f.initialization_data.initializeprob, p)
10371044
prob2 = @test_nowarn remake(prob; u0 = [y => 0.5])
10381045
@test is_variable(prob.f.initialization_data.initializeprob, p)

0 commit comments

Comments
 (0)