Subroutine */
| 8002 | } /* dswap_ */ |
| 8003 | |
| 8004 | /* Subroutine */ int dsymv_(char *uplo, integer *n, doublereal *alpha, |
| 8005 | doublereal *a, integer *lda, doublereal *x, integer *incx, doublereal |
| 8006 | *beta, doublereal *y, integer *incy) |
| 8007 | { |
| 8008 | /* System generated locals */ |
| 8009 | integer a_dim1, a_offset, i__1, i__2; |
| 8010 | |
| 8011 | /* Local variables */ |
| 8012 | static integer i__, j, ix, iy, jx, jy, kx, ky, info; |
| 8013 | static doublereal temp1, temp2; |
| 8014 | extern logical lsame_(char *, char *); |
| 8015 | extern /* Subroutine */ int xerbla_(char *, integer *); |
| 8016 | |
| 8017 | |
| 8018 | /* |
| 8019 | Purpose |
| 8020 | ======= |
| 8021 | |
| 8022 | DSYMV performs the matrix-vector operation |
| 8023 | |
| 8024 | y := alpha*A*x + beta*y, |
| 8025 | |
| 8026 | where alpha and beta are scalars, x and y are n element vectors and |
| 8027 | A is an n by n symmetric matrix. |
| 8028 | |
| 8029 | Arguments |
| 8030 | ========== |
| 8031 | |
| 8032 | UPLO - CHARACTER*1. |
| 8033 | On entry, UPLO specifies whether the upper or lower |
| 8034 | triangular part of the array A is to be referenced as |
| 8035 | follows: |
| 8036 | |
| 8037 | UPLO = 'U' or 'u' Only the upper triangular part of A |
| 8038 | is to be referenced. |
| 8039 | |
| 8040 | UPLO = 'L' or 'l' Only the lower triangular part of A |
| 8041 | is to be referenced. |
| 8042 | |
| 8043 | Unchanged on exit. |
| 8044 | |
| 8045 | N - INTEGER. |
| 8046 | On entry, N specifies the order of the matrix A. |
| 8047 | N must be at least zero. |
| 8048 | Unchanged on exit. |
| 8049 | |
| 8050 | ALPHA - DOUBLE PRECISION. |
| 8051 | On entry, ALPHA specifies the scalar alpha. |
| 8052 | Unchanged on exit. |
| 8053 | |
| 8054 | A - DOUBLE PRECISION array of DIMENSION ( LDA, n ). |
| 8055 | Before entry with UPLO = 'U' or 'u', the leading n by n |
| 8056 | upper triangular part of the array A must contain the upper |
| 8057 | triangular part of the symmetric matrix and the strictly |
| 8058 | lower triangular part of A is not referenced. |
| 8059 | Before entry with UPLO = 'L' or 'l', the leading n by n |
| 8060 | lower triangular part of the array A must contain the lower |
| 8061 | triangular part of the symmetric matrix and the strictly |