Remove the local copy of Jacobian matrix when doing LU decomposition

[sjsprecious]

This PR:

• removes the generation of local copy of Jacobian matrix
• reconstructs the Jacobian matrix when a substepping is needed for the LU decomposition

fix #644

[codecov-commenter]

Codecov Report

All modified and coverable lines are covered by tests ✅

Project coverage is 93.51%. Comparing base (7c5f4bf) to head (821beb5).

Additional details and impacted files

@@           Coverage Diff           @@
##             main     #646   +/-   ##
=======================================
  Coverage   93.51%   93.51%           
=======================================
  Files          49       49           
  Lines        3500     3502    +2     
=======================================
+ Hits         3273     3275    +2     
  Misses        227      227           

☔ View full report in Codecov by Sentry.
📢 Have feedback on the report? Share it here.