/external/eigen/blas/ |
D | level1_cplx_impl.h | 28 …lar EIGEN_CAT(EIGEN_CAT(REAL_SCALAR_SUFFIX,SCALAR_SUFFIX),asum_)(int *n, RealScalar *px, int *incx) 35 if(*incx==1) return vector(x,*n).unaryExpr<scalar_norm1_op>().sum(); 36 else return vector(x,*n,std::abs(*incx)).unaryExpr<scalar_norm1_op>().sum(); 40 int EIGEN_BLAS_FUNC(dotcw)(int *n, RealScalar *px, int *incx, RealScalar *py, int *incy, RealScalar… 50 if(*incx==1 && *incy==1) *res = (vector(x,*n).dot(vector(y,*n))); 51 else if(*incx>0 && *incy>0) *res = (vector(x,*n,*incx).dot(vector(y,*n,*incy))); 52 else if(*incx<0 && *incy>0) *res = (vector(x,*n,-*incx).reverse().dot(vector(y,*n,*incy))); 53 else if(*incx>0 && *incy<0) *res = (vector(x,*n,*incx).dot(vector(y,*n,-*incy).reverse())); 54 …else if(*incx<0 && *incy<0) *res = (vector(x,*n,-*incx).reverse().dot(vector(y,*n,-*incy).reverse(… 59 int EIGEN_BLAS_FUNC(dotuw)(int *n, RealScalar *px, int *incx, RealScalar *py, int *incy, RealScalar… [all …]
|
D | level1_impl.h | 12 int EIGEN_BLAS_FUNC(axpy)(int *n, RealScalar *palpha, RealScalar *px, int *incx, RealScalar *py, in… in EIGEN_BLAS_FUNC() 20 if(*incx==1 && *incy==1) vector(y,*n) += alpha * vector(x,*n); in EIGEN_BLAS_FUNC() 21 else if(*incx>0 && *incy>0) vector(y,*n,*incy) += alpha * vector(x,*n,*incx); in EIGEN_BLAS_FUNC() 22 else if(*incx>0 && *incy<0) vector(y,*n,-*incy).reverse() += alpha * vector(x,*n,*incx); in EIGEN_BLAS_FUNC() 23 else if(*incx<0 && *incy>0) vector(y,*n,*incy) += alpha * vector(x,*n,-*incx).reverse(); in EIGEN_BLAS_FUNC() 24 …else if(*incx<0 && *incy<0) vector(y,*n,-*incy).reverse() += alpha * vector(x,*n,-*incx).reverse(); in EIGEN_BLAS_FUNC() 29 int EIGEN_BLAS_FUNC(copy)(int *n, RealScalar *px, int *incx, RealScalar *py, int *incy) in EIGEN_BLAS_FUNC() 37 if(*incx==1 && *incy==1) in EIGEN_BLAS_FUNC() 41 if(*incx<0) x = x - (*n-1)*(*incx); in EIGEN_BLAS_FUNC() 46 x += *incx; in EIGEN_BLAS_FUNC() [all …]
|
D | level1_real_impl.h | 14 RealScalar EIGEN_BLAS_FUNC(asum)(int *n, RealScalar *px, int *incx) in EIGEN_BLAS_FUNC() 22 if(*incx==1) return vector(x,*n).cwiseAbs().sum(); in EIGEN_BLAS_FUNC() 23 else return vector(x,*n,std::abs(*incx)).cwiseAbs().sum(); in EIGEN_BLAS_FUNC() 27 Scalar EIGEN_BLAS_FUNC(dot)(int *n, RealScalar *px, int *incx, RealScalar *py, int *incy) in EIGEN_BLAS_FUNC() 36 if(*incx==1 && *incy==1) return (vector(x,*n).cwiseProduct(vector(y,*n))).sum(); in EIGEN_BLAS_FUNC() 37 else if(*incx>0 && *incy>0) return (vector(x,*n,*incx).cwiseProduct(vector(y,*n,*incy))).sum(); in EIGEN_BLAS_FUNC() 38 …else if(*incx<0 && *incy>0) return (vector(x,*n,-*incx).reverse().cwiseProduct(vector(y,*n,*incy))… in EIGEN_BLAS_FUNC() 39 …else if(*incx>0 && *incy<0) return (vector(x,*n,*incx).cwiseProduct(vector(y,*n,-*incy).reverse())… in EIGEN_BLAS_FUNC() 40 …else if(*incx<0 && *incy<0) return (vector(x,*n,-*incx).reverse().cwiseProduct(vector(y,*n,-*incy)… in EIGEN_BLAS_FUNC() 46 Scalar EIGEN_BLAS_FUNC(nrm2)(int *n, RealScalar *px, int *incx) in EIGEN_BLAS_FUNC() [all …]
|
D | double.cpp | 22 double BLASFUNC(dsdot)(int* n, float* x, int* incx, float* y, int* incy) in BLASFUNC() 26 …if(*incx==1 && *incy==1) return (vector(x,*n).cast<double>().cwiseProduct(vector(y,*n).cast<dou… in BLASFUNC() 27 …else if(*incx>0 && *incy>0) return (vector(x,*n,*incx).cast<double>().cwiseProduct(vector(y,*n,*in… in BLASFUNC() 28 …else if(*incx<0 && *incy>0) return (vector(x,*n,-*incx).reverse().cast<double>().cwiseProduct(vect… in BLASFUNC() 29 …else if(*incx>0 && *incy<0) return (vector(x,*n,*incx).cast<double>().cwiseProduct(vector(y,*n,-*i… in BLASFUNC() 30 …else if(*incx<0 && *incy<0) return (vector(x,*n,-*incx).reverse().cast<double>().cwiseProduct(vect… in BLASFUNC()
|
D | level2_cplx_impl.h | 19 …nt *n, RealScalar *palpha, RealScalar *pa, int *lda, RealScalar *px, int *incx, RealScalar *pbeta,… in EIGEN_BLAS_FUNC() 47 else if(*incx==0) info = 7; in EIGEN_BLAS_FUNC() 55 Scalar* actual_x = get_compact_vector(x,*n,*incx); in EIGEN_BLAS_FUNC() 111 int EIGEN_BLAS_FUNC(hpr)(char *uplo, int *n, RealScalar *palpha, RealScalar *px, int *incx, RealSca… in EIGEN_BLAS_FUNC() 135 else if(*incx==0) info = 5; in EIGEN_BLAS_FUNC() 142 Scalar* x_cpy = get_compact_vector(x, *n, *incx); in EIGEN_BLAS_FUNC() 162 int EIGEN_BLAS_FUNC(hpr2)(char *uplo, int *n, RealScalar *palpha, RealScalar *px, int *incx, RealSc… in EIGEN_BLAS_FUNC() 187 else if(*incx==0) info = 5; in EIGEN_BLAS_FUNC() 195 Scalar* x_cpy = get_compact_vector(x, *n, *incx); in EIGEN_BLAS_FUNC() 217 int EIGEN_BLAS_FUNC(her)(char *uplo, int *n, RealScalar *palpha, RealScalar *px, int *incx, RealSca… in EIGEN_BLAS_FUNC() [all …]
|
D | level2_real_impl.h | 13 …nt *n, RealScalar *palpha, RealScalar *pa, int *lda, RealScalar *px, int *incx, RealScalar *pbeta,… in EIGEN_BLAS_FUNC() 41 else if(*incx==0) info = 7; in EIGEN_BLAS_FUNC() 49 Scalar* actual_x = get_compact_vector(x,*n,*incx); in EIGEN_BLAS_FUNC() 71 int EIGEN_BLAS_FUNC(syr)(char *uplo, int *n, RealScalar *palpha, RealScalar *px, int *incx, RealSca… in EIGEN_BLAS_FUNC() 110 else if(*incx==0) info = 5; in EIGEN_BLAS_FUNC() 118 Scalar* x_cpy = get_compact_vector(x,*n,*incx); in EIGEN_BLAS_FUNC() 132 int EIGEN_BLAS_FUNC(syr2)(char *uplo, int *n, RealScalar *palpha, RealScalar *px, int *incx, RealSc… in EIGEN_BLAS_FUNC() 171 else if(*incx==0) info = 5; in EIGEN_BLAS_FUNC() 180 Scalar* x_cpy = get_compact_vector(x,*n,*incx); in EIGEN_BLAS_FUNC() 234 int EIGEN_BLAS_FUNC(spr)(char *uplo, int *n, Scalar *palpha, Scalar *px, int *incx, Scalar *pap) in EIGEN_BLAS_FUNC() [all …]
|
D | common.h | 121 T* get_compact_vector(T* x, int n, int incx) in get_compact_vector() argument 123 if(incx==1) in get_compact_vector() 127 if(incx<0) vector(ret,n) = vector(x,n,-incx).reverse(); in get_compact_vector() 128 else vector(ret,n) = vector(x,n, incx); in get_compact_vector() 133 T* copy_back(T* x_cpy, T* x, int n, int incx) in copy_back() argument 138 if(incx<0) vector(x,n,-incx).reverse() = vector(x_cpy,n); in copy_back() 139 else vector(x,n, incx) = vector(x_cpy,n); in copy_back()
|
D | level2_impl.h | 201 RealScalar *px, int *incx, RealScalar *pbeta, RealScalar *py, int *incy) in EIGEN_BLAS_FUNC() 217 else if(*incx==0) info = 10; in EIGEN_BLAS_FUNC() 230 Scalar* actual_x = get_compact_vector(x,actual_n,*incx); in EIGEN_BLAS_FUNC() 270 … *uplo, char *opa, char *diag, int *n, int *k, RealScalar *pa, int *lda, RealScalar *px, int *incx) 283 else if(*incx==0) info = 9; 292 Scalar* actual_x = get_compact_vector(x,actual_n,*incx); 332 …r *uplo, char *op, char *diag, int *n, int *k, RealScalar *pa, int *lda, RealScalar *px, int *incx) in EIGEN_BLAS_FUNC() 373 else if(*incx==0) info = 9; in EIGEN_BLAS_FUNC() 382 Scalar* actual_x = get_compact_vector(x,actual_n,*incx); in EIGEN_BLAS_FUNC() 390 if(actual_x!=x) delete[] copy_back(actual_x,x,actual_n,*incx); in EIGEN_BLAS_FUNC() [all …]
|
/external/eigen/blas/f2c/ |
D | stbmv.c | 16 integer *k, real *a, integer *lda, real *x, integer *incx, ftnlen in stbmv_() argument 203 } else if (*incx == 0) { in stbmv_() 222 if (*incx <= 0) { in stbmv_() 223 kx = 1 - (*n - 1) * *incx; in stbmv_() 224 } else if (*incx != 1) { in stbmv_() 237 if (*incx == 1) { in stbmv_() 269 ix += *incx; in stbmv_() 276 jx += *incx; in stbmv_() 278 kx += *incx; in stbmv_() 284 if (*incx == 1) { in stbmv_() [all …]
|
D | dtbmv.c | 16 integer *k, doublereal *a, integer *lda, doublereal *x, integer *incx, in dtbmv_() argument 203 } else if (*incx == 0) { in dtbmv_() 222 if (*incx <= 0) { in dtbmv_() 223 kx = 1 - (*n - 1) * *incx; in dtbmv_() 224 } else if (*incx != 1) { in dtbmv_() 237 if (*incx == 1) { in dtbmv_() 269 ix += *incx; in dtbmv_() 276 jx += *incx; in dtbmv_() 278 kx += *incx; in dtbmv_() 284 if (*incx == 1) { in dtbmv_() [all …]
|
D | drotm.c | 15 /* Subroutine */ int drotm_(integer *n, doublereal *dx, integer *incx, in drotm_() argument 99 if (! (*incx == *incy && *incx > 0)) { in drotm_() 103 nsteps = *n * *incx; in drotm_() 115 i__2 = *incx; in drotm_() 128 i__1 = *incx; in drotm_() 143 i__2 = *incx; in drotm_() 155 if (*incx < 0) { in drotm_() 156 kx = (1 - *n) * *incx + 1; in drotm_() 178 kx += *incx; in drotm_() 192 kx += *incx; in drotm_() [all …]
|
D | srotm.c | 15 /* Subroutine */ int srotm_(integer *n, real *sx, integer *incx, real *sy, in srotm_() argument 100 if (! (*incx == *incy && *incx > 0)) { in srotm_() 104 nsteps = *n * *incx; in srotm_() 116 i__2 = *incx; in srotm_() 129 i__1 = *incx; in srotm_() 144 i__2 = *incx; in srotm_() 156 if (*incx < 0) { in srotm_() 157 kx = (1 - *n) * *incx + 1; in srotm_() 179 kx += *incx; in srotm_() 193 kx += *incx; in srotm_() [all …]
|
D | complexdots.c | 22 *incx, complex *cy, integer *incy) in cdotc_() 33 cdotcw_(n, &cx[1], incx, &cy[1], incy, &res); in cdotc_() 38 *incx, complex *cy, integer *incy) in cdotu_() 49 cdotuw_(n, &cx[1], incx, &cy[1], incy, &res); in cdotu_() 53 doublecomplex zdotc_(integer *n, doublecomplex *cx, integer *incx, in zdotc_() argument 65 zdotcw_(n, &cx[1], incx, &cy[1], incy, &res); in zdotc_() 69 doublecomplex zdotu_(integer *n, doublecomplex *cx, integer *incx, in zdotu_() argument 81 zdotuw_(n, &cx[1], incx, &cy[1], incy, &res); in zdotu_()
|
D | dspmv.c | 16 doublereal *ap, doublereal *x, integer *incx, doublereal *beta, in dspmv_() argument 145 } else if (*incx == 0) { in dspmv_() 163 if (*incx > 0) { in dspmv_() 166 kx = 1 - (*n - 1) * *incx; in dspmv_() 221 if (*incx == 1 && *incy == 1) { in dspmv_() 251 ix += *incx; in dspmv_() 256 jx += *incx; in dspmv_() 266 if (*incx == 1 && *incy == 1) { in dspmv_() 296 ix += *incx; in dspmv_() 303 jx += *incx; in dspmv_()
|
D | sspmv.c | 16 real *x, integer *incx, real *beta, real *y, integer *incy, ftnlen in sspmv_() argument 145 } else if (*incx == 0) { in sspmv_() 163 if (*incx > 0) { in sspmv_() 166 kx = 1 - (*n - 1) * *incx; in sspmv_() 221 if (*incx == 1 && *incy == 1) { in sspmv_() 251 ix += *incx; in sspmv_() 256 jx += *incx; in sspmv_() 266 if (*incx == 1 && *incy == 1) { in sspmv_() 296 ix += *incx; in sspmv_() 303 jx += *incx; in sspmv_()
|
D | dsbmv.c | 16 alpha, doublereal *a, integer *lda, doublereal *x, integer *incx, in dsbmv_() argument 186 } else if (*incx == 0) { in dsbmv_() 204 if (*incx > 0) { in dsbmv_() 207 kx = 1 - (*n - 1) * *incx; in dsbmv_() 262 if (*incx == 1 && *incy == 1) { in dsbmv_() 295 ix += *incx; in dsbmv_() 301 jx += *incx; in dsbmv_() 304 kx += *incx; in dsbmv_() 314 if (*incx == 1 && *incy == 1) { in dsbmv_() 347 ix += *incx; in dsbmv_() [all …]
|
D | ssbmv.c | 16 real *a, integer *lda, real *x, integer *incx, real *beta, real *y, in ssbmv_() argument 188 } else if (*incx == 0) { in ssbmv_() 206 if (*incx > 0) { in ssbmv_() 209 kx = 1 - (*n - 1) * *incx; in ssbmv_() 264 if (*incx == 1 && *incy == 1) { in ssbmv_() 297 ix += *incx; in ssbmv_() 303 jx += *incx; in ssbmv_() 306 kx += *incx; in ssbmv_() 316 if (*incx == 1 && *incy == 1) { in ssbmv_() 349 ix += *incx; in ssbmv_() [all …]
|
D | ztbmv.c | 17 *incx, ftnlen uplo_len, ftnlen trans_len, ftnlen diag_len) in ztbmv_() 207 } else if (*incx == 0) { in ztbmv_() 227 if (*incx <= 0) { in ztbmv_() 228 kx = 1 - (*n - 1) * *incx; in ztbmv_() 229 } else if (*incx != 1) { in ztbmv_() 242 if (*incx == 1) { in ztbmv_() 300 ix += *incx; in ztbmv_() 313 jx += *incx; in ztbmv_() 315 kx += *incx; in ztbmv_() 321 if (*incx == 1) { in ztbmv_() [all …]
|
D | ctbmv.c | 16 integer *k, complex *a, integer *lda, complex *x, integer *incx, in ctbmv_() argument 207 } else if (*incx == 0) { in ctbmv_() 227 if (*incx <= 0) { in ctbmv_() 228 kx = 1 - (*n - 1) * *incx; in ctbmv_() 229 } else if (*incx != 1) { in ctbmv_() 242 if (*incx == 1) { in ctbmv_() 300 ix += *incx; in ctbmv_() 313 jx += *incx; in ctbmv_() 315 kx += *incx; in ctbmv_() 321 if (*incx == 1) { in ctbmv_() [all …]
|
/external/cblas/testing/ |
D | c_zblas2.c | 13 CBLAS_TEST_ZOMPLEX *a, int *lda, const void *x, int *incx, in F77_zgemv() argument 29 cblas_zgemv( CblasRowMajor, trans, *m, *n, alpha, A, LDA, x, *incx, in F77_zgemv() 35 *m, *n, alpha, a, *lda, x, *incx, beta, y, *incy ); in F77_zgemv() 38 *m, *n, alpha, a, *lda, x, *incx, beta, y, *incy ); in F77_zgemv() 43 CBLAS_TEST_ZOMPLEX *x, int *incx, in F77_zgbmv() argument 77 *incx, beta, y, *incy ); in F77_zgbmv() 82 *incx, beta, y, *incy ); in F77_zgbmv() 85 *incx, beta, y, *incy ); in F77_zgbmv() 89 CBLAS_TEST_ZOMPLEX *x, int *incx, CBLAS_TEST_ZOMPLEX *y, int *incy, in F77_zgeru() argument 103 cblas_zgeru( CblasRowMajor, *m, *n, alpha, x, *incx, y, *incy, A, LDA ); in F77_zgeru() [all …]
|
D | c_cblas2.c | 13 CBLAS_TEST_COMPLEX *a, int *lda, const void *x, int *incx, in F77_cgemv() argument 29 cblas_cgemv( CblasRowMajor, trans, *m, *n, alpha, A, LDA, x, *incx, in F77_cgemv() 35 *m, *n, alpha, a, *lda, x, *incx, beta, y, *incy ); in F77_cgemv() 38 *m, *n, alpha, a, *lda, x, *incx, beta, y, *incy ); in F77_cgemv() 43 CBLAS_TEST_COMPLEX *x, int *incx, in F77_cgbmv() argument 77 *incx, beta, y, *incy ); in F77_cgbmv() 82 *incx, beta, y, *incy ); in F77_cgbmv() 85 *incx, beta, y, *incy ); in F77_cgbmv() 89 CBLAS_TEST_COMPLEX *x, int *incx, CBLAS_TEST_COMPLEX *y, int *incy, in F77_cgeru() argument 103 cblas_cgeru( CblasRowMajor, *m, *n, alpha, x, *incx, y, *incy, A, LDA ); in F77_cgeru() [all …]
|
D | c_dblas2.c | 12 double *a, int *lda, double *x, int *incx, double *beta, in F77_dgemv() argument 27 *m, *n, *alpha, A, LDA, x, *incx, *beta, y, *incy ); in F77_dgemv() 32 *m, *n, *alpha, a, *lda, x, *incx, *beta, y, *incy ); in F77_dgemv() 35 *m, *n, *alpha, a, *lda, x, *incx, *beta, y, *incy ); in F77_dgemv() 38 void F77_dger(int *order, int *m, int *n, double *alpha, double *x, int *incx, in F77_dger() argument 53 cblas_dger(CblasRowMajor, *m, *n, *alpha, x, *incx, y, *incy, A, LDA ); in F77_dger() 60 cblas_dger( CblasColMajor, *m, *n, *alpha, x, *incx, y, *incy, a, *lda ); in F77_dger() 64 int *n, double *a, int *lda, double *x, int *incx) { in F77_dtrmv() argument 81 cblas_dtrmv(CblasRowMajor, uplo, trans, diag, *n, A, LDA, x, *incx); in F77_dtrmv() 85 cblas_dtrmv(CblasColMajor, uplo, trans, diag, *n, a, *lda, x, *incx); in F77_dtrmv() [all …]
|
D | c_sblas2.c | 12 float *a, int *lda, float *x, int *incx, float *beta, in F77_sgemv() argument 27 *m, *n, *alpha, A, LDA, x, *incx, *beta, y, *incy ); in F77_sgemv() 32 *m, *n, *alpha, a, *lda, x, *incx, *beta, y, *incy ); in F77_sgemv() 35 *m, *n, *alpha, a, *lda, x, *incx, *beta, y, *incy ); in F77_sgemv() 38 void F77_sger(int *order, int *m, int *n, float *alpha, float *x, int *incx, in F77_sger() argument 53 cblas_sger(CblasRowMajor, *m, *n, *alpha, x, *incx, y, *incy, A, LDA ); in F77_sger() 60 cblas_sger( CblasColMajor, *m, *n, *alpha, x, *incx, y, *incy, a, *lda ); in F77_sger() 64 int *n, float *a, int *lda, float *x, int *incx) { in F77_strmv() argument 81 cblas_strmv(CblasRowMajor, uplo, trans, diag, *n, A, LDA, x, *incx); in F77_strmv() 85 cblas_strmv(CblasColMajor, uplo, trans, diag, *n, a, *lda, x, *incx); in F77_strmv() [all …]
|
/external/eigen/Eigen/src/Core/products/ |
D | TriangularMatrixVector_MKL.h | 108 MKL_INT m, n, lda, incx, incy; \ 117 incx = 1; \ 126 MKLPREFIX##trmv(&uplo, &trans, &diag, &n, (const MKLTYPE*)_lhs, &lda, (MKLTYPE*)x, &incx); \ 129 MKLPREFIX##axpy(&n, &alpha_,(const MKLTYPE*)x, &incx, (MKLTYPE*)_res, &incy); \ 148 …MKLPREFIX##gemv(&trans, &m, &n, &alpha_, (const MKLTYPE*)a, &lda, (const MKLTYPE*)x, &incx, &beta_… 193 MKL_INT m, n, lda, incx, incy; \ 202 incx = 1; \ 211 MKLPREFIX##trmv(&uplo, &trans, &diag, &n, (const MKLTYPE*)_lhs, &lda, (MKLTYPE*)x, &incx); \ 214 MKLPREFIX##axpy(&n, &alpha_,(const MKLTYPE*)x, &incx, (MKLTYPE*)_res, &incy); \ 233 …MKLPREFIX##gemv(&trans, &n, &m, &alpha_, (const MKLTYPE*)a, &lda, (const MKLTYPE*)x, &incx, &beta_…
|
/external/cblas/src/ |
D | cblas_zhpr2.c | 28 #define F77_incX incx in cblas_zhpr2() 31 int n, i, j, incx=incX, incy=incY; in cblas_zhpr2() local 119 incx = 1; in cblas_zhpr2() 121 incx = -1; in cblas_zhpr2()
|