I want to modify the Jacobian inside an IJacobian call. Right now I just
have
TSComputeIJacobian(ts, ..., A, B, ...)
but I want to do
TSComputeIJacobian(ts, ..., A0, B, ...)
MatCreateLRC(A0, U, NULL, V, &A)
However, the MatCreateLRC() needs a Mat *, so what I really want is to do
something like
TSComputeIJacobian(ts, ..., A, B, ...)
MatCreateLRC(A, U, NULL, V, &Anew)
<copy implementation of Anew into A>
Isn't there a function for this?
Thanks,
Matt
--
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>
