Home | History | Annotate | Download | only in blas

Lines Matching full:vector

26 // computes the sum of magnitudes of all vector elements or, for a complex vector x, the sum
27 // res = |Rex1| + |Imx1| + |Rex2| + |Imx2| + ... + |Rexn| + |Imxn|, where x is a vector of order n
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();
39 // computes a dot product of a conjugated vector with another vector.
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()));
58 // computes a vector-vector dot product without complex conjugation.
69 if(*incx==1 && *incy==1) *res = (vector(x,*n).cwiseProduct(vector(y,*n))).sum();
70 else if(*incx>0 && *incy>0) *res = (vector(x,*n,*incx).cwiseProduct(vector(y,*n,*incy))).sum();
71 else if(*incx<0 && *incy>0) *res = (vector(x,*n,-*incx).reverse().cwiseProduct(vector(y,*n,*incy))).sum();
72 else if(*incx>0 && *incy<0) *res = (vector(x,*n,*incx).cwiseProduct(vector(y,*n,-*incy).reverse())).sum();
73 else if(*incx<0 && *incy<0) *res = (vector(x,*n,-*incx).reverse().cwiseProduct(vector(y,*n,-*incy).reverse())).sum();
85 return vector(x,*n).stableNorm();
87 return vector(x,*n,*incx).stableNorm();
99 StridedVectorType vx(vector(x,*n,std::abs(*incx)));
100 StridedVectorType vy(vector(y,*n,std::abs(*incy)));
122 if(*incx==1) vector(x,*n) *= alpha;
123 else vector(x,*n,std::abs(*incx)) *= alpha;