Lines Matching refs:wa1
36 Matrix< Scalar, Dynamic, 1 > wa1, wa2;
41 wa1 = qtb;
46 wa1[j] = 0.;
49 wa1[j] /= r(j,j);
50 temp = wa1[j];
52 wa1[i] -= r(i,j) * temp;
56 x[ipvt[j]] = wa1[j];
77 wa1[j] = diag[l] * (wa2[l] / dxnorm);
84 sum += r(i,j) * wa1[i];
85 wa1[j] = (wa1[j] - sum) / r(j,j);
87 temp = wa1.blueNorm();
93 wa1[j] = r.col(j).head(j+1).dot(qtb.head(j+1)) / diag[ipvt[j]];
95 gnorm = wa1.stableNorm();
114 wa1 = sqrt(par)* diag;
117 qrsolv<Scalar>(r, ipvt, wa1, qtb, x, sdiag);
133 wa1[j] = diag[l] * (wa2[l] / dxnorm);
136 wa1[j] /= sdiag[j];
137 temp = wa1[j];
139 wa1[i] -= r(i,j) * temp;
141 temp = wa1.blueNorm();
193 Matrix< Scalar, Dynamic, 1 > wa1, wa2;
200 wa1 = qtb;
201 wa1.tail(n-rank).setZero();
202 qr.matrixQR().topLeftCorner(rank, rank).template triangularView<Upper>().solveInPlace(wa1.head(rank));
204 x = qr.colsPermutation()*wa1;
223 wa1 = qr.colsPermutation().inverse() * diag.cwiseProduct(wa2)/dxnorm;
224 qr.matrixQR().topLeftCorner(n, n).transpose().template triangularView<Lower>().solveInPlace(wa1);
225 temp = wa1.blueNorm();
231 wa1[j] = qr.matrixQR().col(j).head(j+1).dot(qtb.head(j+1)) / diag[qr.colsPermutation().indices()(j)];
233 gnorm = wa1.stableNorm();
253 wa1 = sqrt(par)* diag;
256 qrsolv<Scalar>(s, qr.colsPermutation().indices(), wa1, qtb, x, sdiag);
270 wa1 = qr.colsPermutation().inverse() * diag.cwiseProduct(wa2/dxnorm);
272 // qr.matrixQR().topLeftCorner(n, n).transpose().template triangularView<Lower>().solveInPlace(wa1);
274 wa1[j] /= sdiag[j];
275 temp = wa1[j];
277 wa1[i] -= s(i,j) * temp;
279 temp = wa1.blueNorm();
