X-Git-Url: https://pintos-os.org/cgi-bin/gitweb.cgi?a=blobdiff_plain;f=src%2Fmath%2Fmoments.c;h=d129f6ab038b3eee06104e716b15fbcf432c61c5;hb=d75247c28e0dce9c21070e4ee14fdc6a2338fb77;hp=a5be12e9cd0d53a108cd91ea98b3d85bd6daf4f3;hpb=9f087e7aa4cdff1d5d46d5e188c0017a9d2d0029;p=pspp-builds.git diff --git a/src/math/moments.c b/src/math/moments.c index a5be12e9..d129f6ab 100644 --- a/src/math/moments.c +++ b/src/math/moments.c @@ -20,6 +20,7 @@ #include #include #include +#include #include #include "xalloc.h" @@ -54,7 +55,7 @@ calc_moments (enum moment max_moment, { double s3 = s2 * sqrt (s2); double g1 = (w * d3) / ((w - 1.0) * (w - 2.0) * s3); - if (finite (g1)) + if (isfinite (g1)) *skewness = g1; } if (max_moment >= MOMENT_KURTOSIS && kurtosis != NULL && w > 3.) @@ -62,7 +63,7 @@ calc_moments (enum moment max_moment, double den = (w - 2.) * (w - 3.) * pow2 (s2); double g2 = (w * (w + 1) * d4 / (w - 1.) / den - 3. * pow2 (d2) / den); - if (finite (g2)) + if (isfinite (g2)) *kurtosis = g2; } }