Implement zero-RK in Julia #1
Description
Abstract
Just a casual thread here. The acceleration of zero-RK is pretty impressive. Recently, there are a couple of undergoing Julia packages for chemical kinetic modeling in Julia, motivated by the powerful auto-differentiation and adjoint sensitivity in Julia eco-system.
One of the pain we have been experiencing is the relatively slow adjoint computation due to the stiffness. I am wondering preconditioner in zero-RK might be useful for the adjoint method as well.
None of I and Matt knows preconditioner very well. If you are interested in Julia's implementations, we could start some discussions on how to accelerate the development.
Reference
ReactionMechanismGenerator/ReactionMechanismSimulator.jl#71
https://github.com/ReactionMechanismGenerator/ReactionMechanismSimulator.jl
https://github.com/DENG-MIT/Arrhenius.jl