Hi everyone, I would like to discuss the best practices for automatic generation of Jacobians when solving PDEs. I understand that the ‘way to go’ is to (1) define a residual function (in loop-style - pure Julia), (2) use SparsityDetection to get the pattern (assuming we don’t know a priori), (3) g…

Best practices for automatic generation of Jacobians

ChrisRackauckas September 12, 2021, 11:51am 4

Just use a closure. (du,u) -> f(du,u,T,k,dx)

1 Like

Topic		Replies	Views
How to optimize ForwardDiff.jacobian for large system of equations? General Usage	4	1281	August 15, 2019
Calculating Sparse Jacobians with Autodifferentiation Numerics package	8	2331	February 28, 2017
DifferentialEquationsjl with ForwardDiff.jl New to Julia	4	1086	January 24, 2022
How to deal with sparse Jacobians? Optimization (Mathematical)	13	3775	January 13, 2019
Using ForwardDiff.jacobian to provide jac to ODE solvers New to Julia question , diffeq	1	686	December 13, 2021