In each case NAME.F will denote the drivers and NAME.RES will denote the numerical results obtained. The runs were done on an IBM RS6000. We used the Fortran 77 compiler with optimization: f77 -O3 MEBDFDAE.F NAME.F.
- Van der Pol equation: A stiff problem of dimension 2
- Robertson problem: of dimension 3
- HIRES problem: of dimension 8
- Cusp catastrophe model: of dimension 96
- Elastic beam Problem: of dimension 80
- Movement of a rectangular plate: of dimension 80
- Brusselator problem with 1D diffusion: of dimension 1000
- Kuramoto-Sivashinsky equation: of dimension 1024.
With this equation a fast fourier transform is needed and this is provided in FFT.F
- Oregonator, Belusov-Zhabotinskii reaction: of dimension 3
- E5: of dimension 4
- B5: The B5 problem
of DETEST with ALPHA=50,200,500,800,1000.
This is a severe test for non A-stable methods.
The drivers and the results are respectively
The Ring Modulator problem is described in [5]. With CS=1e-9 it is a moderately easy problem but with CS=2e-12 it represents a very challenging problem.
- Ring Modulator problem with CS=1e-9
- Ring Modulator problem with CS=2e-12
This is the end of the results we present for stiff ordinary differential
equations using MEBDFDAE.
In what follows we present the
Results for MEBDFDAE on some linearly implicit DAEs of the form MY´ =F(T,Y)
Note here that the mass matrix M is constant.- The pendulum problem [5] with INDEX=0, 1, 2, 3.
- Stiff pendulum problem [5] : index 3 DAE of dimension 5
The problem PENDUL0.F runs as an ODE since it is an index 0 DAE.
We also give results for some of the challenging problems in the
Amsterdam test set. In particular we consider:
Top of the page
Last updated April 07, 2000. Jeff Cash (j.cash@ma.ic.ac.uk)
people have visited this page since April 1997.