sb_acc_i(3PI) Sun Interval Performance Library sb_acc_i(3PI)NAMEsb_acc_i - Symmetric band matrix accumulation and scale.
DESCRIPTION
B = ALPHA * A + BETA * B
or
B = ALPHA * A_transposed + BETA * B
The routine SB_ACC_I scales an symmetric band interval matrix A (or its
transpose), scales an symmetric band 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 or P is less than zero.
SYNOPSIS
SUBROUTINE SB_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 (P+1,N) where P = band width.
ARGUMENTS
B (input/output)
On entry, B is array of DIMENSION (P+1,N)
where P = band width.
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 (P+1,N)
where P = band width.
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.
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 sb_acc_i(3PI)