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


(cherry picked from commit 7a018dc3)