Add a flag `insert(..., verify=True)` that checks for linearity of `compute_current`

[michaeldeistler]

Internally, the implicit Euler solver requires that compute_current is affine in v. We should add an automated check for this, as it has caused confusion and NaN e.g. here.