Lines Matching refs:bstep
52 LUImpl(_Tp* A, size_t astep, int m, _Tp* b, size_t bstep, int n, _Tp eps) in LUImpl() argument
56 bstep /= sizeof(b[0]); in LUImpl()
75 std::swap(b[i*bstep + j], b[k*bstep + j]); in LUImpl()
90 b[j*bstep + k] += alpha*b[i*bstep + k]; in LUImpl()
101 _Tp s = b[i*bstep + j]; in LUImpl()
103 s -= A[i*astep + k]*b[k*bstep + j]; in LUImpl()
104 b[i*bstep + j] = s*A[i*astep + i]; in LUImpl()
112 int LU(float* A, size_t astep, int m, float* b, size_t bstep, int n) in LU() argument
114 return LUImpl(A, astep, m, b, bstep, n, FLT_EPSILON*10); in LU()
118 int LU(double* A, size_t astep, int m, double* b, size_t bstep, int n) in LU() argument
120 return LUImpl(A, astep, m, b, bstep, n, DBL_EPSILON*100); in LU()
125 CholImpl(_Tp* A, size_t astep, int m, _Tp* b, size_t bstep, int n) in CholImpl() argument
131 bstep /= sizeof(b[0]); in CholImpl()
176 s = b[i*bstep + j]; in CholImpl()
178 s -= L[i*astep + k]*b[k*bstep + j]; in CholImpl()
179 b[i*bstep + j] = (_Tp)(s*L[i*astep + i]); in CholImpl()
187 s = b[i*bstep + j]; in CholImpl()
189 s -= L[k*astep + i]*b[k*bstep + j]; in CholImpl()
190 b[i*bstep + j] = (_Tp)(s*L[i*astep + i]); in CholImpl()
198 bool Cholesky(float* A, size_t astep, int m, float* b, size_t bstep, int n) in Cholesky() argument
200 return CholImpl(A, astep, m, b, bstep, n); in Cholesky()
203 bool Cholesky(double* A, size_t astep, int m, double* b, size_t bstep, int n) in Cholesky() argument
205 return CholImpl(A, astep, m, b, bstep, n); in Cholesky()