b0_unscaled `nearly_ge` 0 &&
b1_unscaled `nearly_ge` 0 &&
b2_unscaled `nearly_ge` 0 &&
b0_unscaled `nearly_ge` 0 &&
b1_unscaled `nearly_ge` 0 &&
b2_unscaled `nearly_ge` 0 &&
-- would do if we used the regular b0,..b3 functions.
b0_unscaled :: Double
b0_unscaled = volume inner_tetrahedron
-- would do if we used the regular b0,..b3 functions.
b0_unscaled :: Double
b0_unscaled = volume inner_tetrahedron