sp_acc_i(3PI) Sun Interval Performance Library sp_acc_i(3PI)NAMEsp_acc_i - Symmetric packed matrix accumulation and scale.
DESCRIPTION
B = ALPHA * A + BETA * B
or
B = ALPHA * A_transposed + BETA * B
The routine SP_ACC_I scales an symmetric packed interval matrix A (or
its transpose) and scales an symmetric packed interval matrix B and
accumulates the result in the interval matrix B. Matrices A (or
A_transposed) and B have the same storage format. This routine returns
immediately if N is equal to zero.
SYNOPSIS
SUBROUTINE SP_ACC_I( A, B [, UPLO] [, TRANSA] [, ALPHA] [, BETA] )
INTERVAL, INTENT(INOUT) :: B(:)
INTERVAL, INTENT(IN) :: A(:)
INTERVAL, INTENT(IN), OPTIONAL :: ALPHA, BETA
TYPE(blas_trans_type), INTENT(IN), OPTIONAL :: TRANSA
TYPE(blas_uplo_type), INTENT(IN), OPTIONAL :: UPLO
where
A and B have shape (N*(N+1)/2).
ARGUMENTS
B (input/output)
On entry, B is array of DIMENSION (N*(N+1)/2).
On exit, if TRANSA = blas_no_trans:
B = ALPHA * A + BETA * B
If TRANSA /= blas_no_trans:
B = ALPHA * A_transposed + BETA * B
A (input)
On entry, A is array of DIMENSION (N*(N+1)/2).
Unchanged on exit.
ALPHA (input)
On entry, ALPHA specifies the scalar ALPHA.
Unchanged on exit.
BETA (input)
On entry, BETA specifies the scalar BETA.
Unchanged on exit.
TRANSA (input)
Scalar element.
Specifies the matrix A transpose.
Possible values:
TRANSA = blas_no_trans (default)
TRANSA = blas_trans
Unchanged on exit.
UPLO(input)
Scalar element.
Specifies if the matrix operand is an upper or
lower tringular part of the symmetric matrix.
Possible values:
UPLO = blas_upper (default)
UPLO = blas_lower
Unchanged on exit.
ERRORS
blas_error_alloc
A or B is unallocated array
blas_error_shape
shape(A) /= shape(B)
6 Mar 2009 sp_acc_i(3PI)