936     int iterations, error = 0, jN, iN, kN, lN = 0;  in icvSingularValueDecomposition()  local
958     for( iN = 0; iN < MN; iN += N )  in icvSingularValueDecomposition()
961             U[iN + j] = A[iN + j];  in icvSingularValueDecomposition()
973     for( i = 0, iN = 0; i < N; i++, iN += N )  in icvSingularValueDecomposition()
977         lN = iN + N;  in icvSingularValueDecomposition()
986         for( kN = iN; kN < MN; kN += N )  in icvSingularValueDecomposition()
992             for( kN = iN; kN < MN; kN += N )  in icvSingularValueDecomposition()
999             f = U[iN + i];  in icvSingularValueDecomposition()
1002             U[iN + i] = f - g;  in icvSingularValueDecomposition()
1009                 for( kN = iN; kN < MN; kN += N )  in icvSingularValueDecomposition()
1017                 for( kN = iN; kN < MN; kN += N )  in icvSingularValueDecomposition()
1024             for( kN = iN; kN < MN; kN += N )  in icvSingularValueDecomposition()
1037             scale += fabs( U[iN + k] );  in icvSingularValueDecomposition()
1045                 U[iN + k] /= scale;  in icvSingularValueDecomposition()
1046                 s += (U[iN + k]) * (U[iN + k]);  in icvSingularValueDecomposition()
1049             f = U[iN + l];  in icvSingularValueDecomposition()
1055                 rv1[k] = U[iN + k] / h;  in icvSingularValueDecomposition()
1063                     s += U[jN + k] * U[iN + k];  in icvSingularValueDecomposition()
1071                 U[iN + k] *= scale;  in icvSingularValueDecomposition()
1086         for( i = N - 1, iN = NN - N; i >= 0; i--, iN -= N )  in icvSingularValueDecomposition()
1097                         V[jN + i] = U[iN + j] / U[iN + l] / g;  in icvSingularValueDecomposition()
1105                             s += U[iN + k] * V[kN + j];  in icvSingularValueDecomposition()
1114                     V[iN + j] = 0;  in icvSingularValueDecomposition()
1119             V[iN + i] = 1;  in icvSingularValueDecomposition()
1122             lN = iN;  in icvSingularValueDecomposition()
1131         for( i = N - 1, iN = NN - N; i >= 0; i--, iN -= N )  in icvSingularValueDecomposition()
1135             lN = iN + N;  in icvSingularValueDecomposition()
1139                 U[iN + j] = 0;  in icvSingularValueDecomposition()
1153                     f = s / U[iN + i] / g;  in icvSingularValueDecomposition()
1155                     for( kN = iN; kN < MN; kN += N )  in icvSingularValueDecomposition()
1159                 for( jN = iN; jN < MN; jN += N )  in icvSingularValueDecomposition()
1165                 for( jN = iN; jN < MN; jN += N )  in icvSingularValueDecomposition()
1169             U[iN + i] += 1;  in icvSingularValueDecomposition()