## Codes

Contact me for the following programs. Alternatively, you may look at public repositories on bitbucket, github, gitlab accounts.

• ### BifurcationKit.jl

Bifurcation Analysis code in Julia. Allows to perform continuation and (automatic) bifurcation analysis in large dimensions on CPUs / GPUs. ODEs are also well within the scope of the package.

• ### PiecewiseDeterministicMarkovProcesses.jl

Simulation of PDMPs in Julia based on this paper with speed on par with C. PDMP means piecewise deterministic Markov process.

• ### LSODA

A Julia wrapper to call the LSODA algorithm by Linda Petzold and Alan Hindmarsh. It solves ODE by switching automatically between stiff and non-stiff methods.

• ### Neuron Code

Simulation code (asic_model_container.pdf) in Neuron associated with our paper posted on BioRxiv. Download the file (right click, Save Link As...), change the name into asic_model_container.zip and unzip it.

• ### Pytrilinos

I have been working with Bill Spotz from the Sandia Labs to develop PyTrilinos, e.g. the nonlinear / numerical continuation part, in Trilinos > v.12.

• ### Neural fields equations in PETSc

Simulation of Neural Fields equations using petsc4py. This allows parallel and fast simulations of models $$\frac{d}{dt}V(x,t) = -V(x,t) + \int_\Omega J(x,y)S(V(y,t))dy,\ y\in\mathbb R^2$$

together with the numerical bifurcation analysis. The number of unknowns can be >$10^6$.

• ### Hopf curves for DDE

Computing Hopf bifurcation curves for delay differential equations (paper), basically finding $\lambda\in i\mathbb R$ such that there is a solution in $U$ to

$$\left(\lambda+l\right) U_i(x)=\sum\limits_{j=1}^p\int\limits_\Omega J_{ij}(x,y)e^{-\lambda\tau(x,y)}U_j(y)dy,\ 1\leq i\leq p,\ x,y\in\mathbb R^d$$

or