More robust rejection calculation
currently a ZeroDivisionError
is raised in the rejection calculation here.
In !370 (merged) a new save_divide
function is introduced which would allow to return infinity for these cases.
However, doing this will cause the plotting to fail, which needs to be adapted accordingly then.