Differential algebraic equations solvers HBT10DAE and DP87DAE tested on four problems:

- car axis
- 8-node transistor amplifier 
- pendulum
- 5-node transistor amplifier 

********************************

The test files to be run are:

HBT10DAE_caraxis.cpp //tests HBT10DAE on the car axis problem
HBT10DAE_transamCWI.cpp //tests HBT10DAE on the 8-node transistor amplifier problem
HBT10DAE_pend.cpp //tests HBT10DAE on the pendulum problem
HBT10DAE_transam.cpp //tests HBT10DAE on the 5-node transistor amplifier problem

DP87DAE_caraxis.cpp //tests DP87DAE on the car axis problem
DP87DAE_transamCWI.cpp //tests DP87DAE on the 8-node transistor amplifier problem
DP87DAE_pend.cpp //tests DP87DAE on the pendulum problem
DP87DAE_transam.cpp //tests DP87DAE on the 5-node transistor amplifier problem

********************************

Results are written in (in the same folder):

DAE tests HBT10 DP87.txt

********************************

Files containing the two solvers:

solver_HBT10DAE.cpp //contains the 1-step 9-stage Hermite-Birkhoff-Taylor DAE solver of order 10
solver_DP87DAE.cpp //contains the Dormand-Prince DP(8,7)13M DAE solver

********************************

Files containing the function y'=f to be evaluated for the four test problems:

prob_caraxis.cpp //car axis
prob_transamCWI.cpp //8-node transistor amplifier
prob_pendulum.cpp //pendulum
prob_transam.cpp //5-node transistor amplifier

********************************

Other files:

autoDiff.cpp //contains routines needed for automatic differentiation
linSolver.cpp //contains a routine to solver a linear system Ax=b with A square

********************************

Output:

Example: "Solver: HBT10DAE Prb=pendulum TOL=0.0001 NFE=3454 SCD=1.66824 n=243 nreje=58 CPU time=0.00351648********"

//Solver: either HBT10DAE or DP87DAE, 
//Prb: problem name, 
//TOL: absolute tolerance, 
//NFE: number of function evaluations
//SCD: number of significant correct digits
//n: number of success steps
//nreje: number of rejected steps
//CPU time: CPU time in seconds