math: fix long double ulperror calculation

diff --git a/src/math/util.c b/src/math/util.c
index f02f5a0..a44b5f6 100644 (file)
--- a/src/math/util.c
+++ b/src/math/util.c
@@ -87,7 +87,7 @@ float ulperrl(long double got, long double want, float dwant)
                got = copysignl(0x1p16383L, got);
                want *= 0.5;
        }
-       return scalbn(got - want, -eulpl(want)) + dwant;
+       return scalbnl(got - want, -eulpl(want)) + dwant;
 #else
        // TODO
        return inf;