/external/eigen/unsupported/Eigen/src/NonLinearOptimization/ |
lmpar.h | 36 Matrix< Scalar, Dynamic, 1 > wa1, wa2; local 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) 193 Matrix< Scalar, Dynamic, 1 > wa1, wa2; local [all...] |
dogleg.h | 30 Matrix< Scalar, Dynamic, 1 > wa1(n), wa2(n); 56 wa1.fill(0.); 58 wa1.tail(n-j) += qrfac.row(j).tail(n-j) * qtb[j]; 59 wa1[j] /= diag[j]; 64 gnorm = wa1.stableNorm(); 72 wa1.array() /= (diag*gnorm).array(); 74 // wa2 = qrfac.template triangularView<Upper>() * wa1; 78 sum += qrfac(j,i) * wa1[i]; 102 x = temp * wa1 + alpha * x;
|
fdjac1.h | 31 Matrix< Scalar, Dynamic, 1 > wa1(n); 44 iflag = Functor(x, wa1); 48 fjac.col(j) = (wa1-fvec)/h; 60 iflag = Functor(x, wa1); 70 fjac.col(j).segment(start, length) = ( wa1.segment(start, length)-fvec.segment(start, length))/h;
|
HybridNonLinearSolver.h | 113 FVectorType wa1, wa2, wa3, wa4; member in class:Eigen::HybridNonLinearSolver 147 wa1.resize(n); wa2.resize(n); wa3.resize(n); wa4.resize(n); 237 internal::dogleg<Scalar>(R, diag, qtf, delta, wa1); 240 wa1 = -wa1; 241 wa2 = x + wa1; 242 pnorm = diag.cwiseProduct(wa1).stableNorm(); 260 wa3 = R.template triangularView<Upper>()*wa1 + qtf; 326 wa1 = diag.cwiseProduct( diag.cwiseProduct(wa1)/pnorm ) [all...] |
LevenbergMarquardt.h | 115 FVectorType wa1, wa2, wa3, wa4; member in class:Eigen::LevenbergMarquardt 169 wa1.resize(n); wa2.resize(n); wa3.resize(n); 268 internal::lmpar2<Scalar>(qrfac, diag, qtf, delta, par, wa1); 271 wa1 = -wa1; 272 wa2 = x + wa1; 273 pnorm = diag.cwiseProduct(wa1).stableNorm(); 292 wa3 = fjac.template triangularView<Upper>() * (qrfac.colsPermutation().inverse() *wa1); 383 wa1.resize(n); wa2.resize(n); wa3.resize(n); 466 wa1 = fjac.diagonal() [all...] |
/external/eigen/unsupported/Eigen/src/LevenbergMarquardt/ |
LMpar.h | 55 VectorType wa1, wa2; local 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() [all...] |
/external/libvorbis/lib/ |
smallft.c | 114 static void dradf2(int ido,int l1,float *cc,float *ch,float *wa1){ 144 tr2=wa1[i-2]*cc[t3-1]+wa1[i-1]*cc[t3]; 145 ti2=wa1[i-2]*cc[t3]-wa1[i-1]*cc[t3-1]; 169 static void dradf4(int ido,int l1,float *cc,float *ch,float *wa1, 211 cr2=wa1[i-2]*cc[t3-1]+wa1[i-1]*cc[t3]; 212 ci2=wa1[i-2]*cc[t3]-wa1[i-1]*cc[t3-1] [all...] |
/external/speex/libspeex/ |
smallft.c | 116 static void dradf2(int ido,int l1,float *cc,float *ch,float *wa1){ 146 tr2=wa1[i-2]*cc[t3-1]+wa1[i-1]*cc[t3]; 147 ti2=wa1[i-2]*cc[t3]-wa1[i-1]*cc[t3-1]; 171 static void dradf4(int ido,int l1,float *cc,float *ch,float *wa1, 213 cr2=wa1[i-2]*cc[t3-1]+wa1[i-1]*cc[t3]; 214 ci2=wa1[i-2]*cc[t3]-wa1[i-1]*cc[t3-1] [all...] |