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