Update gpu_ode_modelingtoolkit_dae.jl

Author: ChrisRackauckas

test/gpu_kernel_de/stiff_ode/gpu_ode_modelingtoolkit_dae.jl

@@ -1,6 +1,8 @@
-using ModelingToolkit, DiffEqGPU, OrdinaryDiffEq, LinearAlgebra, Test
+using ModelingToolkit, DiffEqGPU, OrdinaryDiffEq, LinearAlgebra, StaticArrays, Test
 using ModelingToolkit: t_nounits as t, D_nounits as D
 
+include("../../utils.jl")
+
 # Define the cartesian pendulum DAE system
 @parameters g = 9.81 L = 1.0
 @variables x(t) y(t) [state_priority = 10] λ(t)
@@ -16,16 +18,16 @@ eqs = [D(D(x)) ~ λ * x / L
 @named pendulum = ODESystem(eqs, t, [x, y, λ], [g, L])
 
 # Perform structural simplification with index reduction
-pendulum_sys = structural_simplify(dae_index_lowering(pendulum))
+pendulum_sys = mtkcompile(pendulum)
 
 # Initial conditions: pendulum starts at bottom right position
-u0 = [x => 1.0, y => 0.0, λ => -g]  # λ initial guess for tension
+u0 = SA[x => 1.0, y => 0.0, λ => -g]  # λ initial guess for tension
 
 # Time span  
 tspan = (0.0f0, 1.0f0)
 
 # Create the ODE problem
-prob = ODEProblem(pendulum_sys, u0, tspan, Float32[])
+prob = ODEProblem(pendulum_sys, u0, tspan)
 
 # Verify DAE properties
 @test SciMLBase.has_initialization_data(prob.f) == true
@@ -54,4 +56,4 @@ sol2 = solve(monteprob, GPURodas4(), EnsembleGPUKernel(backend),
     dt = 0.01f0,
     adaptive = false)
 
-@test length(sol2.u) == 2
+@test length(sol2.u) == 2