style cleanup: comments
[blender-staging.git] / source / blender / blenlib / intern / math_matrix.c
index 1f61b37a1afae116e054a1bb87df5ea6135549fb..3e19ea4f999dff55c872795ae54a1a957ed0209a 100644 (file)
@@ -1408,8 +1408,8 @@ void svd_m4(float U[4][4], float s[4], float V[4][4], float A_[4][4])
 
        int i = 0, j = 0, k = 0, p, pp, iter;
 
-       // Reduce A to bidiagonal form, storing the diagonal elements
-       // in s and the super-diagonal elements in e.
+       /* Reduce A to bidiagonal form, storing the diagonal elements
+        * in s and the super-diagonal elements in e. */
 
        int nct = minf(m - 1, n);
        int nrt = maxf(0, minf(n - 2, m));
@@ -1421,9 +1421,9 @@ void svd_m4(float U[4][4], float s[4], float V[4][4], float A_[4][4])
        for (k = 0; k < maxf(nct, nrt); k++) {
                if (k < nct) {
 
-                       // Compute the transformation for the k-th column and
-                       // place the k-th diagonal in s[k].
-                       // Compute 2-norm of k-th column without under/overflow.
+                       /* Compute the transformation for the k-th column and
+                        * place the k-th diagonal in s[k].
+                        * Compute 2-norm of k-th column without under/overflow. */
                        s[k] = 0;
                        for (i = k; i < m; i++) {
                                s[k] = hypotf(s[k], A[i][k]);
@@ -1444,7 +1444,7 @@ void svd_m4(float U[4][4], float s[4], float V[4][4], float A_[4][4])
                for (j = k + 1; j < n; j++) {
                        if ((k < nct) && (s[k] != 0.0f)) {
 
-                               // Apply the transformation.
+                               /* Apply the transformation. */
 
                                float t = 0;
                                for (i = k; i < m; i++) {
@@ -1456,24 +1456,24 @@ void svd_m4(float U[4][4], float s[4], float V[4][4], float A_[4][4])
                                }
                        }
 
-                       // Place the k-th row of A into e for the
-                       // subsequent calculation of the row transformation.
+                       /* Place the k-th row of A into e for the */
+                       /* subsequent calculation of the row transformation. */
 
                        e[j] = A[k][j];
                }
                if (k < nct) {
 
-                       // Place the transformation in U for subsequent back
-                       // multiplication.
+                       /* Place the transformation in U for subsequent back
+                        * multiplication. */
 
                        for (i = k; i < m; i++)
                                U[i][k] = A[i][k];
                }
                if (k < nrt) {
 
-                       // Compute the k-th row transformation and place the
-                       // k-th super-diagonal in e[k].
-                       // Compute 2-norm without under/overflow.
+                       /* Compute the k-th row transformation and place the
+                        * k-th super-diagonal in e[k].
+                        * Compute 2-norm without under/overflow. */
                        e[k] = 0;
                        for (i = k + 1; i < n; i++) {
                                e[k] = hypotf(e[k], e[i]);
@@ -1493,7 +1493,7 @@ void svd_m4(float U[4][4], float s[4], float V[4][4], float A_[4][4])
                        if ((k + 1 < m) & (e[k] != 0.0f)) {
                                float invek1;
 
-                               // Apply the transformation.
+                               /* Apply the transformation. */
 
                                for (i = k + 1; i < m; i++) {
                                        work[i] = 0.0f;
@@ -1512,15 +1512,15 @@ void svd_m4(float U[4][4], float s[4], float V[4][4], float A_[4][4])
                                }
                        }
 
-                       // Place the transformation in V for subsequent
-                       // back multiplication.
+                       /* Place the transformation in V for subsequent
+                        * back multiplication. */
 
                        for (i = k + 1; i < n; i++)
                                V[i][k] = e[i];
                }
        }
 
-       // Set up the final bidiagonal matrix or order p.
+       /* Set up the final bidiagonal matrix or order p. */
 
        p = minf(n, m + 1);
        if (nct < n) {
@@ -1534,7 +1534,7 @@ void svd_m4(float U[4][4], float s[4], float V[4][4], float A_[4][4])
        }
        e[p - 1] = 0.0f;
 
-       // If required, generate U.
+       /* If required, generate U. */
 
        for (j = nct; j < nu; j++) {
                for (i = 0; i < m; i++) {
@@ -1570,7 +1570,7 @@ void svd_m4(float U[4][4], float s[4], float V[4][4], float A_[4][4])
                }
        }
 
-       // If required, generate V.
+       /* If required, generate V. */
 
        for (k = n - 1; k >= 0; k--) {
                if ((k < nrt) & (e[k] != 0.0f)) {
@@ -1591,7 +1591,7 @@ void svd_m4(float U[4][4], float s[4], float V[4][4], float A_[4][4])
                V[k][k] = 1.0f;
        }
 
-       // Main iteration loop for the singular values.
+       /* Main iteration loop for the singular values. */
 
        pp = p - 1;
        iter = 0;
@@ -1599,20 +1599,20 @@ void svd_m4(float U[4][4], float s[4], float V[4][4], float A_[4][4])
        while (p > 0) {
                int kase = 0;
 
-               // Test for maximum iterations to avoid infinite loop
+               /* Test for maximum iterations to avoid infinite loop */
                if (maxiter == 0)
                        break;
                maxiter--;
 
-               // This section of the program inspects for
-               // negligible elements in the s and e arrays.  On
-               // completion the variables kase and k are set as follows.
-
-               // kase = 1       if s(p) and e[k - 1] are negligible and k<p
-               // kase = 2       if s(k) is negligible and k<p
-               // kase = 3       if e[k - 1] is negligible, k<p, and
-               //               s(k), ..., s(p) are not negligible (qr step).
-               // kase = 4       if e(p - 1) is negligible (convergence).
+               /* This section of the program inspects for
+                * negligible elements in the s and e arrays.  On
+                * completion the variables kase and k are set as follows.
+                *
+                * kase = 1       if s(p) and e[k - 1] are negligible and k<p
+                * kase = 2       if s(k) is negligible and k<p
+                * kase = 3       if e[k - 1] is negligible, k<p, and
+                *               s(k), ..., s(p) are not negligible (qr step).
+                * kase = 4       if e(p - 1) is negligible (convergence). */
 
                for (k = p - 2; k >= -1; k--) {
                        if (k == -1) {
@@ -1653,11 +1653,11 @@ void svd_m4(float U[4][4], float s[4], float V[4][4], float A_[4][4])
                }
                k++;
 
-               // Perform the task indicated by kase.
+               /* Perform the task indicated by kase. */
 
                switch (kase) {
 
-                       // Deflate negligible s(p).
+                       /* Deflate negligible s(p). */
 
                        case 1:
                        {
@@ -1683,7 +1683,7 @@ void svd_m4(float U[4][4], float s[4], float V[4][4], float A_[4][4])
                                break;
                        }
 
-                       // Split at negligible s(k).
+                       /* Split at negligible s(k). */
 
                        case 2:
                        {
@@ -1707,12 +1707,12 @@ void svd_m4(float U[4][4], float s[4], float V[4][4], float A_[4][4])
                                break;
                        }
 
-                       // Perform one qr step.
+                       /* Perform one qr step. */
 
                        case 3:
                        {
 
-                               // Calculate the shift.
+                               /* Calculate the shift. */
 
                                float scale = maxf(maxf(maxf(maxf(
                                                   fabsf(s[p - 1]), fabsf(s[p - 2])), fabsf(e[p - 2])),
@@ -1737,7 +1737,7 @@ void svd_m4(float U[4][4], float s[4], float V[4][4], float A_[4][4])
                                f = (sk + sp) * (sk - sp) + shift;
                                g = sk * ek;
 
-                               // Chase zeros.
+                               /* Chase zeros. */
 
                                for (j = k; j < p - 1; j++) {
                                        float t = hypotf(f, g);
@@ -1779,12 +1779,12 @@ void svd_m4(float U[4][4], float s[4], float V[4][4], float A_[4][4])
                                iter = iter + 1;
                                break;
                        }
-                       // Convergence.
+                       /* Convergence. */
 
                        case 4:
                        {
 
-                               // Make the singular values positive.
+                               /* Make the singular values positive. */
 
                                if (s[k] <= 0.0f) {
                                        s[k] = (s[k] < 0.0f ? -s[k] : 0.0f);
@@ -1793,7 +1793,7 @@ void svd_m4(float U[4][4], float s[4], float V[4][4], float A_[4][4])
                                                V[i][k] = -V[i][k];
                                }
 
-                               // Order the singular values.
+                               /* Order the singular values. */
 
                                while (k < pp) {
                                        float t;