@@ -527,27 +527,29 @@ end
527527@parameters k1 k2 ω
528528@variables X (t) Y (t)
529529eqs_1st_order = [D (Y) + Y - ω ~ 0 ,
530- X + k1 ~ Y + k2]
531- eqs_2nd_order = [D (D (Y)) + 2 ω* D (Y) + (ω^ 2 )* Y ~ 0 ,
532- X + k1 ~ Y + k2]
530+ X + k1 ~ Y + k2]
531+ eqs_2nd_order = [D (D (Y)) + 2 ω * D (Y) + (ω^ 2 ) * Y ~ 0 ,
532+ X + k1 ~ Y + k2]
533533@mtkbuild sys_1st_order = ODESystem (eqs_1st_order, t)
534534@mtkbuild sys_2nd_order = ODESystem (eqs_2nd_order, t)
535535
536536u0_1st_order_1 = [X => 1.0 , Y => 2.0 ]
537537u0_1st_order_2 = [Y => 2.0 ]
538538u0_2nd_order_1 = [X => 1.0 , Y => 2.0 , D (Y) => 0.5 ]
539539u0_2nd_order_2 = [Y => 2.0 , D (Y) => 0.5 ]
540- tspan = (0.0 , 10. )
540+ tspan = (0.0 , 10.0 )
541541ps = [ω => 0.5 , k1 => 2.0 , k2 => 3.0 ]
542542
543543oprob_1st_order_1 = ODEProblem (sys_1st_order, u0_1st_order_1, tspan, ps)
544544oprob_1st_order_2 = ODEProblem (sys_1st_order, u0_1st_order_2, tspan, ps)
545545oprob_2nd_order_1 = ODEProblem (sys_2nd_order, u0_2nd_order_1, tspan, ps) # gives sys_2nd_order
546546oprob_2nd_order_2 = ODEProblem (sys_2nd_order, u0_2nd_order_2, tspan, ps)
547547
548- @test solve (oprob_1st_order_1, Rosenbrock23 ()). retcode == SciMLBase. ReturnCode. InitialFailure
548+ @test solve (oprob_1st_order_1, Rosenbrock23 ()). retcode ==
549+ SciMLBase. ReturnCode. InitialFailure
549550@test solve (oprob_1st_order_2, Rosenbrock23 ())[Y][1 ] == 2.0
550- @test solve (oprob_2nd_order_1, Rosenbrock23 ()). retcode == SciMLBase. ReturnCode. InitialFailure
551+ @test solve (oprob_2nd_order_1, Rosenbrock23 ()). retcode ==
552+ SciMLBase. ReturnCode. InitialFailure
551553sol = solve (oprob_2nd_order_2, Rosenbrock23 ()) # retcode: Success
552554@test sol[Y][1 ] == 2.0
553555@test sol[1 ][2 ] == 0.5
0 commit comments