Add (broken) validation test for thermalfluid benchmark

Author: serenity4

test/benchmark.jl

@@ -0,0 +1,24 @@
+using Pkg; Pkg.activate(joinpath(dirname(@__DIR__), "benchmark"))
+using DAECompiler
+using DAECompiler: compute_residual_vectors
+using SciMLBase, Sundials
+using Test
+
+include("../benchmark/thermalfluid.jl")
+
+@testset "Validation" begin
+  Benchmark{3}()()
+
+  u = zeros(68)
+  du = zeros(68)
+  residuals, expanded_residuals = compute_residual_vectors(Benchmark{3}(), u, du)
+  @test length(residuals) == length(expanded_residuals)
+  @test_broken residuals ≈ expanded_residuals
+  # indices = findall(i -> residuals[i] ≉ expanded_residuals[i], eachindex(residuals))
+  # residuals[indices]
+  # expanded_residuals[indices]
+end
+
+let sol = solve(DAECProblem(Benchmark{3}(), [1:9;] .=> 0.), IDA())
+    @test_broken sol.retcode == ReturnCode.Success
+end

test/runtests.jl

@@ -11,5 +11,5 @@ include("validation.jl")
 
 using Pkg
 Pkg.activate(joinpath(dirname(@__DIR__), "benchmark")) do
-  include("../benchmark/main.jl")
+  include("benchmark.jl")
 end