inline void her2( const Order, const UpLo, const int n, const double alpha,
const double* x, const int incx, const double* y, const int incy,
double* a, const int lda ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasDsyr2( blas_option< UpLo >::value, n, alpha, x, incx, y, incy, a,
lda );
}
void magma_dsyr2(
magma_uplo_t uplo,
magma_int_t n,
double alpha,
const double *dx, magma_int_t incx,
const double *dy, magma_int_t incy,
double *dA, magma_int_t ldda )
{
cublasDsyr2(
cublas_uplo_const( uplo ),
n,
alpha, dx, incx,
dy, incy,
dA, ldda );
}
void d_syr2(SEXP ruplo, SEXP ralpha, SEXP rx, SEXP rincx, SEXP ry, SEXP rincy,
SEXP ra, SEXP rlda)
{
char
uplo = getSymLoc(ruplo);
double
alpha = asReal(ralpha),
* a, * x, * y;
int
rowsa, colsa, lda = asInteger(rlda),
nx, incx = asInteger(rincx),
ny, incy = asInteger(rincy);
unpackVector(rx, &nx, &x);
unpackVector(ry, &ny, &y);
unpackMatrix(ra, &rowsa, &colsa, &a);
cublasDsyr2(uplo, rowsa, alpha, x, incx, y, incy, a, lda);
checkCublasError("d_syr2");
}