JCublasJava/src/main/java/jcuda/jcublas/JCublas.java
public static void cublasDspr(char uplo, int n, double alpha, Pointer x, int incx, Pointer AP)
    {
        cublasDsprNative(uplo, n, alpha, x, incx, AP);
        checkResultBLAS();
    }
public static void cublasDspr(char uplo, int n, double alpha, Pointer x, int incx, Pointer AP)
    {
        cublasDsprNative(uplo, n, alpha, x, incx, AP);
        checkResultBLAS();
    }
public static void cublasDspr ( char uplo , int n , double alpha , Pointer x , int incx , Pointer AP ) { cublasDsprNative ( uplo , n , alpha , x , incx , AP ) ; checkResultBLAS ( ) ; }
<pre>
void
cublasDspr (char uplo, int n, double alpha, const double *x, int incx,
double *AP)

performs the symmetric rank 1 operation

A = alpha * x * transpose(x) + A,

where alpha is a double precision scalar and x is an n element double
precision vector. A is a symmetric n x n matrix consisting of double
precision elements that is supplied in packed form.

Input
-----
uplo   specifies whether the matrix data is stored in the upper or the lower
triangular part of array AP. If uplo == 'U' or 'u', then the upper
triangular part of A is supplied in AP. If uplo == 'L' or 'l', then
the lower triangular part of A is supplied in AP.
n      specifies the number of rows and columns of the matrix A. It must be
at least zero.
alpha  double precision scalar multiplier applied to x * transpose(x).
x      double precision array of length at least (1 + (n - 1) * abs(incx)).
incx   storage spacing between elements of x. incx must not be zero.
AP     double precision array with at least ((n * (n + 1)) / 2) elements. If
uplo == 'U' or 'u', the array AP contains the upper triangular part
of the symmetric matrix A, packed sequentially, column by column;
that is, if i <= j, then A[i,j] is stored is AP[i+(j*(j+1)/2)]. If
uplo == 'L' or 'L', the array AP contains the lower triangular part
of the symmetric matrix A, packed sequentially, column by column;
that is, if i >= j, then A[i,j] is stored in AP[i+((2*n-j+1)*j)/2].

Output
------
A      updated according to A = alpha * x * transpose(x) + A

Reference: http://www.netlib.org/blas/dspr.f

Error status for this function can be retrieved via cublasGetError().

Error Status
------------
CUBLAS_STATUS_NOT_INITIALIZED  if CUBLAS library has not been initialized
CUBLAS_STATUS_INVALID_VALUE    if n < 0, or incx == 0
CUBLAS_STATUS_ARCH_MISMATCH    if invoked on device without DP support
CUBLAS_STATUS_EXECUTION_FAILED if function failed to launch on GPU
</pre>
<pre > void cublasDspr ( char uplo int n double alpha const double * x int incx double * AP )