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