diff --git a/.gitignore b/.gitignore index c36e7f4..763513e 100644 --- a/.gitignore +++ b/.gitignore @@ -1,9 +1 @@ .ipynb_checkpoints -data.out -raw.out -plot.png - -main - -*.o -*.mod diff --git a/README.md b/README.md new file mode 100644 index 0000000..7305851 --- /dev/null +++ b/README.md @@ -0,0 +1,17 @@ +The generic_rk module is an attempt to modularize solvers for ordinary differential equations (ODEs). + +Below is the result of an example problem known as the Brusselator using an adaptive timestep control. The first subplot shows the raw results of $y_1$ and $y_2$ within the range $t \in [0, 20]$. The second subplot shows the acceptable $dt$ used throughout the simulation. + +$$ +\begin{array}{r l} +y_1' &=& 1 + y_1^2 y_2 -4y_1 \\ +y_2' &=& 3y_1 - y_1^2 y_2 \\ +\end{array} +$$ + +$$ y_1(0) = 1.5, ~~~ y_2(0) = 3, ~~~ t \in [0, 20]$$ + +![alt text] [sample plot] + + +[sample plot]: ./plot.png diff --git a/obj/.gitignore b/obj/.gitignore index e69de29..a0c6846 100644 --- a/obj/.gitignore +++ b/obj/.gitignore @@ -0,0 +1,2 @@ +*.o +*.mod diff --git a/plot.png b/plot.png new file mode 100644 index 0000000..60ab816 Binary files /dev/null and b/plot.png differ diff --git a/src/test/.gitignore b/src/test/.gitignore index e69de29..fa7de1c 100644 --- a/src/test/.gitignore +++ b/src/test/.gitignore @@ -0,0 +1,3 @@ +plot.png +data.out +raw.out diff --git a/src/test/bin/.gitignore b/src/test/bin/.gitignore index e69de29..9e98c40 100644 --- a/src/test/bin/.gitignore +++ b/src/test/bin/.gitignore @@ -0,0 +1,2 @@ +main +main.exe diff --git a/src/test/obj/.gitignore b/src/test/obj/.gitignore index e69de29..a0c6846 100644 --- a/src/test/obj/.gitignore +++ b/src/test/obj/.gitignore @@ -0,0 +1,2 @@ +*.o +*.mod