# How to tackle 2D reaction-diffusion equation?

**URL:** https://discourse.julialang.org/t/how-to-tackle-2d-reaction-diffusion-equation/54110
**Category:** Numerics
**Created:** [January 28, 2021, 10:41am UTC](https://discourse.julialang.org/t/how-to-tackle-2d-reaction-diffusion-equation/54110 "2021-01-28T10:41:04Z")
**Posts on this page:** 1
**Showing post:** 2

<div class="post-metadata">

### Author: ![samo](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/samo/32/35398_2.png) [@samo](https://discourse.julialang.org/u/samo)
#### Post date: [January 28, 2021, 10:45am UTC](https://discourse.julialang.org/t/how-to-tackle-2d-reaction-diffusion-equation/54110/2 "2021-01-28T10:45:37Z")

</div>

See here a diffusion example:

> **[GitHub - omlins/ParallelStencil.jl: Package for writing high-level code for...](https://github.com/omlins/ParallelStencil.jl#50-lines-example-deployable-on-gpu-and-cpu)**
>
> Package for writing high-level code for parallel high-performance stencil computations that can be deployed on both GPUs and CPUs - GitHub - omlins/ParallelStencil.jl: Package for writing high-leve...

You can add reaction to it. It is using the package [ParallelStencil.jl](https://github.com/omlins/ParallelStencil.jl) which empowers domain scientists to write architecture-agnostic high-level code for parallel high-performance stencil computations on GPUs and CPUs.

In addition, you might want to use implicit time stepping and convergence accelerators. [Here](https://www.youtube.com/watch?v=vPsfZUqI4_0) you find a presentation of a methodology using the pseudo-transient method and [ParallelStencil.jl](https://github.com/omlins/ParallelStencil.jl) (and [ImplicitGlobalGrid.jl](https://github.com/eth-cscs/ImplicitGlobalGrid.jl)).

---

_[View the full topic](https://discourse.julialang.org/t/how-to-tackle-2d-reaction-diffusion-equation/54110)._