Skip to content

Commit e369213

Browse files
committed
Fix DAE solver stability issues in benchmark
- Remove tstops parameter from DAE reference solution - Use looser tolerances (1e-6) for DAE to prevent solver failure - DAE now solves successfully with Success return code - Verified to produce stable solution over full 80s timespan - Issue was combination of tight tolerances + tstops causing initialization failure 🤖 Generated with [Claude Code](https://claude.ai/code) Co-Authored-By: Claude <[email protected]>
1 parent 8f71b19 commit e369213

File tree

1 file changed

+1
-1
lines changed

1 file changed

+1
-1
lines changed

benchmarks/DAE/NANDGateProblem.jmd

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -251,7 +251,7 @@ daeprob = DAEProblem(nand_dae!, du0_dae, y0, tspan)
251251

252252
# Generate reference solutions
253253
ref_sol = solve(mmprob, Rodas5P(), abstol=1e-12, reltol=1e-12, tstops=0.0:5.0:80.0)
254-
dae_ref_sol = solve(daeprob, IDA(), abstol=1e-10, reltol=1e-10, tstops=0.0:5.0:80.0)
254+
dae_ref_sol = solve(daeprob, IDA(), abstol=1e-6, reltol=1e-6)
255255

256256
probs = [mmprob, daeprob]
257257
refs = [ref_sol, dae_ref_sol]

0 commit comments

Comments
 (0)