This program is intended to reproduce the results from the preprint
"Exponential integrators for quasilinear wave-type equations" by Benjamin Dörich and Marlis Hochbruck
The codes generates the lines in Figure 1
#####################################################################################
The program is tested with
- Ubuntu 16.04.7 LTS and Python 3.7.6 and the following version of its modules:
numpy - 1.15.4
scipy - 1.4.1
dolfin - 2018.1.0
fenics - 2018.1.0
matplotlib - 3.2.1
tikzplotlib - 0.9.6
- Ubuntu 18.04.5 LTS and Python 3.6.9 and the following version of its modules:
numpy - 1.19.2
scipy - 1.5.1
dolfin - 2019.2.0.dev0
fenics - 2019.2.0.dev0
matplotlib - 3.3.2
tikzplotlib - 0.9.4
#####################################################################################
In the folder "DoeH21_quasilinear" open a console and run the following commands after each other.
- Run "python3 1_generate_initial_matrices.py"
- Run "python3 2_paper_run_reference_solution.py"
- Run "python3 3_paper_run_Euler.py"
- Run "python3 4_paper_run_midpoint.py"
#####################################################################################
After running the calculations, the errors can be found in
the folder "preprint_DoeH21_quasilinear/extracted_data/".
(a) settings_* contains the information of the config file
(b) infos_* contains the computed errors
(c) plot_* contains a tikz-file which gives the plot for the
Euler and the midpoint rule separately.