double[] y = new double[S != null ? S.m2 : 1]; /* get workspace */ Dcs_utsolve.cs_utsolve(N.U, y); /* y = R'\y */ for (int k = m - 1; k >= 0; k--) /* apply Householder refl. to x */
double[] y = new double[S != null ? S.m2 : 1]; /* get workspace */ Dcs_utsolve.cs_utsolve(N.U, y); /* y = R'\y */ for (int k = m - 1; k >= 0; k--) /* apply Householder refl. to x */
double[] y = new double[S != null ? S.m2 : 1]; /* get workspace */ Dcs_utsolve.cs_utsolve(N.U, y); /* y = R'\y */ for (int k = m - 1; k >= 0; k--) {
if (ok) { Dcs_utsolve.cs_utsolve(N.U, x); /* x = R'\x */ for (k = m - 1; k >= 0; k--) /* apply Householder refl. to x */