Subroutine */
| 15106 | } /* strsm_ */ |
| 15107 | |
| 15108 | /* Subroutine */ int zaxpy_(integer *n, doublecomplex *za, doublecomplex *zx, |
| 15109 | integer *incx, doublecomplex *zy, integer *incy) |
| 15110 | { |
| 15111 | /* System generated locals */ |
| 15112 | integer i__1, i__2, i__3, i__4; |
| 15113 | doublecomplex z__1, z__2; |
| 15114 | |
| 15115 | /* Local variables */ |
| 15116 | static integer i__, ix, iy; |
| 15117 | extern doublereal dcabs1_(doublecomplex *); |
| 15118 | |
| 15119 | |
| 15120 | /* |
| 15121 | Purpose |
| 15122 | ======= |
| 15123 | |
| 15124 | ZAXPY constant times a vector plus a vector. |
| 15125 | |
| 15126 | Further Details |
| 15127 | =============== |
| 15128 | |
| 15129 | jack dongarra, 3/11/78. |
| 15130 | modified 12/3/93, array(1) declarations changed to array(*) |
| 15131 | |
| 15132 | ===================================================================== |
| 15133 | */ |
| 15134 | |
| 15135 | /* Parameter adjustments */ |
| 15136 | --zy; |
| 15137 | --zx; |
| 15138 | |
| 15139 | /* Function Body */ |
| 15140 | if (*n <= 0) { |
| 15141 | return 0; |
| 15142 | } |
| 15143 | if (dcabs1_(za) == 0.) { |
| 15144 | return 0; |
| 15145 | } |
| 15146 | if (*incx == 1 && *incy == 1) { |
| 15147 | goto L20; |
| 15148 | } |
| 15149 | |
| 15150 | /* |
| 15151 | code for unequal increments or equal increments |
| 15152 | not equal to 1 |
| 15153 | */ |
| 15154 | |
| 15155 | ix = 1; |
| 15156 | iy = 1; |
| 15157 | if (*incx < 0) { |
| 15158 | ix = (-(*n) + 1) * *incx + 1; |
| 15159 | } |
| 15160 | if (*incy < 0) { |
| 15161 | iy = (-(*n) + 1) * *incy + 1; |
| 15162 | } |
| 15163 | i__1 = *n; |
| 15164 | for (i__ = 1; i__ <= i__1; ++i__) { |
| 15165 | i__2 = iy; |