Possible to use Jacobian pattern (Jacobian sparcity matrix)?

headmyshoulder / odeint-v2

odeint - solving ordinary differential equations in c++ v2

http://headmyshoulder.github.com/odeint-v2/

Other

337 stars 102 forks source link

Possible to use Jacobian pattern (Jacobian sparcity matrix)? #263

Open wimverbruggen opened 8 months ago

wimverbruggen commented 8 months ago

I don't find much about this in the documentation, but I'm trying to solve a stiff system for which I would need to specify the Jacobian. However, it would be much more convenient for me to specify the Jacobian sparsity matrix instead (which is in this case a tridiagonal matrix of ones). Is this possible with this library?