* Use more suitable numeric diff epsilon, leading to slightly better Newton convergence
* This will also make the test pass for the new staggered implementation