Improve testing on jacobians

Description

The current test for jacobians only does checks against a reference solution. This means that test should have been run before for each case.

Proposal

Instead of running the whole problem to check the final result, generate a jacobian for numeric and analytic local jacobian operators and compare their values. This would allow us to test basically any case to see if both jacobians are the same up to some error threshold.

How to test the implementation?

Errors between the two matrices are low enough.

Related issues

See #