c/hvapprox.c (compute_hua_wang_direction): Pre-compute sinl().

Author: MLopez-Ibanez

c/hvapprox.c

@@ -573,15 +573,21 @@ compute_hua_wang_direction(double * restrict direction, dimension_t dim,
 {
     ASSUME(2 <= dim && dim <= 32);
     dimension_t k, j;
-    direction[0] = STATIC_CAST(double, sinl(theta[0]));
+    for (j = 1; j < dim; j++)
+        direction[j] = STATIC_CAST(double, cosl(theta[dim - j - 1]));
+
+    double * sinl_theta = malloc((dim - 1) * sizeof(*theta));
+    for (k = 0; k < dim - 1; k++)
+        sinl_theta[k] = STATIC_CAST(double, sinl(theta[k]));
+    direction[0] = sinl_theta[0];
     for (k = 1; k < dim - 1; k++)
-        direction[0] *= STATIC_CAST(double, sinl(theta[k]));
+        direction[0] *= sinl_theta[k];
     for (j = 1; j < dim; j++) {
-        direction[j] = STATIC_CAST(double, cosl(theta[dim - j - 1]));
         for (k = 0; k < dim - j - 1; k++) {
-            direction[j] *= STATIC_CAST(double, sinl(theta[k]));
+            direction[j] *= sinl_theta[k];
         }
     }
+    free(sinl_theta);
     for (k = 0; k < dim; k++) {
         // FIXME: Can direction[k] be negative? If not, then we don't need fabs().
         direction[k] = (fabs(direction[k]) <= ALMOST_ZERO_WEIGHT)