Lines Matching refs:wa1
55 VectorType wa1, wa2;
62 wa1 = qtb;
63 wa1.tail(n-rank).setZero();
65 wa1.head(rank) = s.topLeftCorner(rank,rank).template triangularView<Upper>().solve(qtb.head(rank));
67 x = qr.colsPermutation()*wa1;
86 wa1 = qr.colsPermutation().inverse() * diag.cwiseProduct(wa2)/dxnorm;
87 s.topLeftCorner(n,n).transpose().template triangularView<Lower>().solveInPlace(wa1);
88 temp = wa1.blueNorm();
94 wa1[j] = s.col(j).head(j+1).dot(qtb.head(j+1)) / diag[qr.colsPermutation().indices()(j)];
96 gnorm = wa1.stableNorm();
115 wa1 = sqrt(par)* diag;
118 lmqrsolv(s, qr.colsPermutation(), wa1, qtb, x, sdiag);
132 wa1 = qr.colsPermutation().inverse() * diag.cwiseProduct(wa2/dxnorm);
135 wa1[j] /= sdiag[j];
136 temp = wa1[j];
138 wa1[i] -= s.coeff(i,j) * temp;
140 temp = wa1.blueNorm();
