@Elrod I like your idea best. I suppose we loose a tiny bit of accuracy.
It also does not break any symmetry and is fast:
1-((2-eps(Float64)) - oneopenconst(Float64)) = 1 - oneopenconst(Float64)
@Elrod I like your idea best. I suppose we loose a tiny bit of accuracy.
It also does not break any symmetry and is fast:
1-((2-eps(Float64)) - oneopenconst(Float64)) = 1 - oneopenconst(Float64)