区间算术中的mod（模，fmod）

问：

如何计算作为区间给出的数字的模？

我尝试过使用 mpmath 库，但这似乎不起作用：

import mpmath
def mod(a: float, b: float, eps = 1e-9) -> float:
    a = mpmath.iv.mpf([a-eps, a+eps])
    b = mpmath.iv.mpf([b-eps, b+eps])
    return mpmath.fmod(a, b)

x = mod(0.2, 0.333333333))  # ValueError: can only create mpf from zero-width interval

我不拘泥于该库，只是想确保我的计算不会受到数值问题的影响。