Lines Matching refs:wa1
34 Matrix< Scalar, Dynamic, 1 > wa1, wa2;
39 wa1 = qtb;
44 wa1[j] = 0.;
47 wa1[j] /= r(j,j);
48 temp = wa1[j];
50 wa1[i] -= r(i,j) * temp;
54 x[ipvt[j]] = wa1[j];
75 wa1[j] = diag[l] * (wa2[l] / dxnorm);
82 sum += r(i,j) * wa1[i];
83 wa1[j] = (wa1[j] - sum) / r(j,j);
85 temp = wa1.blueNorm();
91 wa1[j] = r.col(j).head(j+1).dot(qtb.head(j+1)) / diag[ipvt[j]];
93 gnorm = wa1.stableNorm();
112 wa1 = sqrt(par)* diag;
115 qrsolv<Scalar>(r, ipvt, wa1, qtb, x, sdiag);
131 wa1[j] = diag[l] * (wa2[l] / dxnorm);
134 wa1[j] /= sdiag[j];
135 temp = wa1[j];
137 wa1[i] -= r(i,j) * temp;
139 temp = wa1.blueNorm();
189 Matrix< Scalar, Dynamic, 1 > wa1, wa2;
196 wa1 = qtb;
197 wa1.tail(n-rank).setZero();
198 qr.matrixQR().topLeftCorner(rank, rank).template triangularView<Upper>().solveInPlace(wa1.head(rank));
200 x = qr.colsPermutation()*wa1;
219 wa1 = qr.colsPermutation().inverse() * diag.cwiseProduct(wa2)/dxnorm;
220 qr.matrixQR().topLeftCorner(n, n).transpose().template triangularView<Lower>().solveInPlace(wa1);
221 temp = wa1.blueNorm();
227 wa1[j] = qr.matrixQR().col(j).head(j+1).dot(qtb.head(j+1)) / diag[qr.colsPermutation().indices()(j)];
229 gnorm = wa1.stableNorm();
249 wa1 = sqrt(par)* diag;
252 qrsolv<Scalar>(s, qr.colsPermutation().indices(), wa1, qtb, x, sdiag);
266 wa1 = qr.colsPermutation().inverse() * diag.cwiseProduct(wa2/dxnorm);
268 // qr.matrixQR().topLeftCorner(n, n).transpose().template triangularView<Lower>().solveInPlace(wa1);
270 wa1[j] /= sdiag[j];
271 temp = wa1[j];
273 wa1[i] -= s(i,j) * temp;
275 temp = wa1.blueNorm();
