In this repository we provide a small code written in Julia to test the high-order, positivity preserving, conservative DeC Patankar scheme. Details can be found here https://arxiv.org/abs/1905.09237
\ No newline at end of file
In this repository we provide a small code written in Julia to test the high-order, positivity preserving, conservative DeC Patankar scheme. Details can be found here https://arxiv.org/abs/1905.09237
## Installation
This code was tested on Julia-1.0.4 and requires the packages LinearAlgebra, Plots and SymPy. To download Julia just visit [Julia webpage](https://julialang.org/)
##Usage
In notebook folder one can find the library files. The files "*_model.jl" contain the production destruction function for the specific model. These are used by the "DeC_procedure.jl" function, that is the core of the algorithm and can be used by
```julia
dts,U=dec_correction(test,M_sub,K_iter,u0,ts)#model, subtimesteps, corrections, initial conditions, times
```
In "Time_integration.jl" one can run some prepared testcases, just changing the test name in
```julia
test="robertson"#"linear" #"non_linear" #choose a model
```
A figure will be produced in the folder Figures
## Authors
*mPDeC* has been developed at University of Zurich by Philipp \"Offner and Davide Torlo
## License
*mPDeC* code is freely available, but the authors do not
