Used pow2(x) instead of x * x where appropriate.
[pspp-builds.git] / src / math / ts / innovations.c
index 553e20e8f7953a879595c3ddb0dce2e0e715b1dc..b9a7cf223442a42ce597748320647a23a0739939 100644 (file)
@@ -160,7 +160,7 @@ innovations_update_scale (struct innovations_estimate *est, double *theta,
       for (j = 0; j < i; j++)
        {
          k = i - j - 1;
-         result -= theta[k] * theta[k] * est->scale[j];
+         result -= pow2 (theta[k]) * est->scale[j];
        }
       est->scale[i] = result;
     }