dspr - perform the symmetric rank 1 operation A := alpha*x*x' + A
SUBROUTINE DSPR( UPLO, N, ALPHA, X, INCX, A) CHARACTER * 1 UPLO INTEGER N, INCX DOUBLE PRECISION ALPHA DOUBLE PRECISION X(*), A(*)
SUBROUTINE DSPR_64( UPLO, N, ALPHA, X, INCX, A) CHARACTER * 1 UPLO INTEGER*8 N, INCX DOUBLE PRECISION ALPHA DOUBLE PRECISION X(*), A(*)
SUBROUTINE SPR( UPLO, N, ALPHA, X, [INCX], A) CHARACTER(LEN=1) :: UPLO INTEGER :: N, INCX REAL(8) :: ALPHA REAL(8), DIMENSION(:) :: X, A
SUBROUTINE SPR_64( UPLO, N, ALPHA, X, [INCX], A) CHARACTER(LEN=1) :: UPLO INTEGER(8) :: N, INCX REAL(8) :: ALPHA REAL(8), DIMENSION(:) :: X, A
#include <sunperf.h>
void dspr(char uplo, int n, double alpha, double *x, int incx, double *a);
void dspr_64(char uplo, long n, double alpha, double *x, long incx, double *a);
dspr performs the symmetric rank 1 operation A := alpha*x*x' + A, where alpha is a real scalar, x is an n element vector and A is an n by n symmetric matrix, supplied in packed form.
UPLO = 'U' or 'u' The upper triangular part of A is supplied in A.
UPLO = 'L' or 'l' The lower triangular part of A is supplied in A.
Unchanged on exit.