diff --git a/_assets/images/AdjointSensitivities.png b/_assets/images/AdjointSensitivities.png new file mode 100644 index 0000000..738ce56 Binary files /dev/null and b/_assets/images/AdjointSensitivities.png differ diff --git a/posts/julia-hackathon-v6-2024-news_1.md b/posts/julia-hackathon-v6-2024-news_1.md index b084470..cd9e344 100644 --- a/posts/julia-hackathon-v6-2024-news_1.md +++ b/posts/julia-hackathon-v6-2024-news_1.md @@ -73,3 +73,19 @@ iteration 1 1.43s 100.0% 1.43s 56.2MiB 100.0% 56.2MiB forward 10 73.0μs 0.0% 7.30μs 0.00B 0.0% 0.00B ─────────────────────────────────────────────────────────────────────── ``` + +## Using Enzyme.jl in order to calculate adjoint sensitivies with JustRelax.jl +*Christian Schuler, Valentin Churavy, Albert de Montserrat, Pascal Aellig* + +[ParallelStencil.jl](https://github.com/omlins/ParallelStencil.jl) was made compatible with the latest [Enzyme.jl](https://github.com/EnzymeAD/Enzyme.jl) version (PR [#169](https://github.com/omlins/ParallelStencil.jl/pull/169) and PR [#170](https://github.com/omlins/ParallelStencil.jl/pull/170)). With the help of [Enzyme.jl](https://github.com/EnzymeAD/Enzyme.jl) and [ParallelStencil.jl](https://github.com/omlins/ParallelStencil.jl) the neccessary vector-Jacobian products (VJP) for the adjoint solve in [JustRelax.jl](https://github.com/PTsolvers/JustRelax.jl) are calculated. + +The figure shows a viscoelastic falling block example which was run with JustRelax.jl + +~~~ +
+ FallingBlock +
+~~~ + + +