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

performs the symmetric rank 2 operation

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

where alpha is a single precision scalar, and x and y are n element single
precision vectors. A is a symmetric n x n matrix consisting of single
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 A. If uplo == 'U' or 'u', then only the
upper triangular part of A may be referenced and the lower triangular
part of A is inferred. If uplo == 'L' or 'l', then only the lower
triangular part of A may be referenced and the upper triangular part
of A is inferred.
n      specifies the number of rows and columns of the matrix A. It must be
at least zero.
alpha  single precision scalar multiplier applied to x * transpose(y) +
y * transpose(x).
x      single precision array of length at least (1 + (n - 1) * abs (incx)).
incx   storage spacing between elements of x. incx must not be zero.
y      single precision array of length at least (1 + (n - 1) * abs (incy)).
incy   storage spacing between elements of y. incy must not be zero.
AP     single 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(y)+alpha*y*transpose(x)+A

Reference: http://www.netlib.org/blas/sspr2.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, incx == 0, incy == 0
CUBLAS_STATUS_EXECUTION_FAILED if function failed to launch on GPU
</pre>