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