Subroutine */
| 2223 | } /* chemv_ */ |
| 2224 | |
| 2225 | /* Subroutine */ int cher2_(char *uplo, integer *n, complex *alpha, complex * |
| 2226 | x, integer *incx, complex *y, integer *incy, complex *a, integer *lda) |
| 2227 | { |
| 2228 | /* System generated locals */ |
| 2229 | integer a_dim1, a_offset, i__1, i__2, i__3, i__4, i__5, i__6; |
| 2230 | real r__1; |
| 2231 | complex q__1, q__2, q__3, q__4; |
| 2232 | |
| 2233 | /* Local variables */ |
| 2234 | static integer i__, j, ix, iy, jx, jy, kx, ky, info; |
| 2235 | static complex temp1, temp2; |
| 2236 | extern logical lsame_(char *, char *); |
| 2237 | extern /* Subroutine */ int xerbla_(char *, integer *); |
| 2238 | |
| 2239 | |
| 2240 | /* |
| 2241 | Purpose |
| 2242 | ======= |
| 2243 | |
| 2244 | CHER2 performs the hermitian rank 2 operation |
| 2245 | |
| 2246 | A := alpha*x*conjg( y' ) + conjg( alpha )*y*conjg( x' ) + A, |
| 2247 | |
| 2248 | where alpha is a scalar, x and y are n element vectors and A is an n |
| 2249 | by n hermitian matrix. |
| 2250 | |
| 2251 | Arguments |
| 2252 | ========== |
| 2253 | |
| 2254 | UPLO - CHARACTER*1. |
| 2255 | On entry, UPLO specifies whether the upper or lower |
| 2256 | triangular part of the array A is to be referenced as |
| 2257 | follows: |
| 2258 | |
| 2259 | UPLO = 'U' or 'u' Only the upper triangular part of A |
| 2260 | is to be referenced. |
| 2261 | |
| 2262 | UPLO = 'L' or 'l' Only the lower triangular part of A |
| 2263 | is to be referenced. |
| 2264 | |
| 2265 | Unchanged on exit. |
| 2266 | |
| 2267 | N - INTEGER. |
| 2268 | On entry, N specifies the order of the matrix A. |
| 2269 | N must be at least zero. |
| 2270 | Unchanged on exit. |
| 2271 | |
| 2272 | ALPHA - COMPLEX . |
| 2273 | On entry, ALPHA specifies the scalar alpha. |
| 2274 | Unchanged on exit. |
| 2275 | |
| 2276 | X - COMPLEX array of dimension at least |
| 2277 | ( 1 + ( n - 1 )*abs( INCX ) ). |
| 2278 | Before entry, the incremented array X must contain the n |
| 2279 | element vector x. |
| 2280 | Unchanged on exit. |
| 2281 | |
| 2282 | INCX - INTEGER. |