Barry Smith <[email protected]> writes: >So make all the “extra” equations be simply f.v[i][j] = x.v[i][j] (or >x.p or x.u depending on where) and put a 1 on the diagonal of that >row/column of the Jacobian).
This would typically be f[j][i].v = x[j][i].v, etc.
pgpvsQhsyJUgd.pgp
Description: PGP signature
