🚀 go-pugleaf

         David Hinds said:
Really, it seems to me that if support for interval arithmetic is built
into a compiler, it should be able to optimize away most of the rounding
mode changes without too much trouble.  This should be trivial for the
cases where the mode changes would otherwise do the most damage, such as
in tight inner loops with no dependencies, i.e., vector operations.  Just
split and do two loops, one for the lower and one for the upper vectors.

         I believe if you try to do this with
      DO 10 J=1,N
   10 X(J)=Y(J)*Z(J)
         you will discover things are not so simple.
                       James B. Shearer