3 /* Subroutine */ int drot_(integer *n, doublereal *dx, integer *incx,
4 doublereal *dy, integer *incy, doublereal *c__, doublereal *s)
6 /* System generated locals */
13 /* .. Scalar Arguments .. */
15 /* .. Array Arguments .. */
21 /* applies a plane rotation. */
22 /* jack dongarra, linpack, 3/11/78. */
23 /* modified 12/3/93, array(1) declarations changed to array(*) */
26 /* .. Local Scalars .. */
28 /* Parameter adjustments */
36 if (*incx == 1 && *incy == 1) {
40 /* code for unequal increments or equal increments not equal */
46 ix = (-(*n) + 1) * *incx + 1;
49 iy = (-(*n) + 1) * *incy + 1;
52 for (i__ = 1; i__ <= i__1; ++i__) {
53 dtemp = *c__ * dx[ix] + *s * dy[iy];
54 dy[iy] = *c__ * dy[iy] - *s * dx[ix];
62 /* code for both increments equal to 1 */
66 for (i__ = 1; i__ <= i__1; ++i__) {
67 dtemp = *c__ * dx[i__] + *s * dy[i__];
68 dy[i__] = *c__ * dy[i__] - *s * dx[i__];