currently we have 1e-5, but if this is good or not depends strongly on the field and the system geometries. We should have some auto value that relates to the dimension involved. This is one reason why we would like to have derivatives in Magpylib