Version 2.1 will allow the user to define the shape (structure) of the Jacobian matrix. I.e., instead of providing analytic Jacobian (or not providing it at all, which is slow), the user can specify the positions of non-zero elements in the Jacobian. The solver will use automatic differentiation for the specified elements only. This should work as fast as analytic Jacobian without requiring the user to differentiate the vector function manually.

Tasks:

• Upgrade Eigen, autodiff, googletest
• Add Jacobian shape class and a helper function to populate an array containing positions of non-zero elements
• Compute Jacobian using the given matrix shape
• Add utility function that compares fully automatic Jacobian with the one computed from the given shape
• Add utility function that compares fully automatic Jacobian with the one provided by the user
• Update tests
• Update examples
• Update docs, README, and CHANGELOG
• Add performance tests (added locally for now)