Actually, I realize now that updating a sparse Jacobian inplace is not so straightforward… See, e.g. these many discussions around this topic:
- Update matrix block inplace from a sparse array,
- Efficiently updating values in a sparse diagonal matrix,
- Update bands of a CSC Matrix,
- Optimal way to do in-place modification of sparse matrix?,
- A version of sparse! that updates a previously existing SparseMatrixCSC?
- Example use of `sparse!` - #2 by bennedich
@ctkelley, Are there any plans to allow for out-of-place (sparse) jacobian functions? I ask because I already have available a fast-but-out-of-place jacobian function 