Problem with sparsity pattern about cppad

[bradbell]

The sparsity pattern 'pattern_in' should be an n2 by n2 identity matrix:
sparse_rc<s_vector> pattern_in;
pattern_in.resize(n2, n2, n2);
for(size_t k = 0; k < n2; k++)
pattern_in.set(k, k, k);

In addition, the pattern passed into sparse_jac_for should be the pattern for the Jacobian:
auto n_sweep = g.sparse_jac_for(
group_max, x2, subset, pattern_jac, coloring, work
);

At the end, I would also print the values as follows:
std::cout << "n_sweep: " << n_sweep << std::endl;
const auto& c = subset.col();
const auto& r = subset.row();
const auto& v = subset.val();
for (int l = 0; l < subset.nnz(); ++l) {
std::cout << "J( " << r[l] << "," << c[l] << ")=" << v[l] << std::endl;
}

It this solves the problem for you, please close this issue.

[bradbell]

I have not head back, so I am closing this issue. Please reopen it with further details if this solution did not work for you.