spr_i(3PI) Sun Interval Performance Library spr_i(3PI)NAMEspr_i - Symmetric rank one update.
DESCRIPTION
A = ALPHA * X * X_transposed + BETA * A_transposed
The routine SPR_I multiplies an interval vectors X by a transposed
interval vector X, scales the resulting interval matrix by interval
ALPHA and adds it to the scaled by BETA transposed symmetric packed
interval matrix A. The resulting interval is copied in the interval
matrix A. The routine returns immediately if N is equal to zero.
SYNOPSIS
SUBROUTINE SPR_I( A, X [, UPLO] [, ALPHA] [, BETA] )
INTERVAL, INTENT(INOUT) :: A(:)
INTERVAL, INTENT(IN) :: X(:)
TYPE(blas_uplo_type), INTENT(IN), OPTIONAL :: UPLO
INTERVAL, INTENT(IN), OPTIONAL :: ALPHA, BETA
where
X has have shape (N),
A has shape (N*(N+1)/2).
ARGUMENTS
A (input/output)
On entry, array of DIMENSION (N*(N+1)/2).
On exit,
A = ALPHA * X * X_transposed + BETA * A_transposed
X (input)
Array of DIMENSION N.
Unchanged on exit.
UPLO(input)
Scalar element.
Specifies if the matrix operand is an upper or
lower tringular part of the matrix.
Possible values:
UPLO = blas_upper (default)
UPLO = blas_lower
Unchanged on exit.
ALPHA (input)
On entry, ALPHA specifies the scalar ALPHA.
Default value is ALPHA = [1,1].
Unchanged on exit.
BETA (input)
On entry, BETA specifies the scalar BETA.
Default value is BETA = [0,0].
Unchanged on exit.
ERRORS
blas_error_alloc
A or X is unallocated array
blas_error_shape
A or X has unappropriate shape
6 Mar 2009 spr_i(3PI)