> > I guess he expects (-2)**(-2) = 1/4. > > But defined this way, (**) would be rather discontinuous indeed, that's what I was expecting. For Float and Double (-2)**(-2) = 1/4, but for the interval arithmetic library that uses the default implementation, it doesn't.