Subroutine */
| 7198 | } /* dgemm_ */ |
| 7199 | |
| 7200 | /* Subroutine */ int dgemv_(char *trans, integer *m, integer *n, doublereal * |
| 7201 | alpha, doublereal *a, integer *lda, doublereal *x, integer *incx, |
| 7202 | doublereal *beta, doublereal *y, integer *incy) |
| 7203 | { |
| 7204 | /* System generated locals */ |
| 7205 | integer a_dim1, a_offset, i__1, i__2; |
| 7206 | |
| 7207 | /* Local variables */ |
| 7208 | static integer i__, j, ix, iy, jx, jy, kx, ky, info; |
| 7209 | static doublereal temp; |
| 7210 | static integer lenx, leny; |
| 7211 | extern logical lsame_(char *, char *); |
| 7212 | extern /* Subroutine */ int xerbla_(char *, integer *); |
| 7213 | |
| 7214 | |
| 7215 | /* |
| 7216 | Purpose |
| 7217 | ======= |
| 7218 | |
| 7219 | DGEMV performs one of the matrix-vector operations |
| 7220 | |
| 7221 | y := alpha*A*x + beta*y, or y := alpha*A'*x + beta*y, |
| 7222 | |
| 7223 | where alpha and beta are scalars, x and y are vectors and A is an |
| 7224 | m by n matrix. |
| 7225 | |
| 7226 | Arguments |
| 7227 | ========== |
| 7228 | |
| 7229 | TRANS - CHARACTER*1. |
| 7230 | On entry, TRANS specifies the operation to be performed as |
| 7231 | follows: |
| 7232 | |
| 7233 | TRANS = 'N' or 'n' y := alpha*A*x + beta*y. |
| 7234 | |
| 7235 | TRANS = 'T' or 't' y := alpha*A'*x + beta*y. |
| 7236 | |
| 7237 | TRANS = 'C' or 'c' y := alpha*A'*x + beta*y. |
| 7238 | |
| 7239 | Unchanged on exit. |
| 7240 | |
| 7241 | M - INTEGER. |
| 7242 | On entry, M specifies the number of rows of the matrix A. |
| 7243 | M must be at least zero. |
| 7244 | Unchanged on exit. |
| 7245 | |
| 7246 | N - INTEGER. |
| 7247 | On entry, N specifies the number of columns of the matrix A. |
| 7248 | N must be at least zero. |
| 7249 | Unchanged on exit. |
| 7250 | |
| 7251 | ALPHA - DOUBLE PRECISION. |
| 7252 | On entry, ALPHA specifies the scalar alpha. |
| 7253 | Unchanged on exit. |
| 7254 | |
| 7255 | A - DOUBLE PRECISION array of DIMENSION ( LDA, n ). |
| 7256 | Before entry, the leading m by n part of the array A must |
| 7257 | contain the matrix of coefficients. |