Complex */
| 194 | } /* ccopy_ */ |
| 195 | |
| 196 | /* Complex */ VOID cdotc_(complex * ret_val, integer *n, complex *cx, integer |
| 197 | *incx, complex *cy, integer *incy) |
| 198 | { |
| 199 | /* System generated locals */ |
| 200 | integer i__1, i__2; |
| 201 | complex q__1, q__2, q__3; |
| 202 | |
| 203 | /* Local variables */ |
| 204 | static integer i__, ix, iy; |
| 205 | static complex ctemp; |
| 206 | |
| 207 | |
| 208 | /* |
| 209 | Purpose |
| 210 | ======= |
| 211 | |
| 212 | forms the dot product of two vectors, conjugating the first |
| 213 | vector. |
| 214 | |
| 215 | Further Details |
| 216 | =============== |
| 217 | |
| 218 | jack dongarra, linpack, 3/11/78. |
| 219 | modified 12/3/93, array(1) declarations changed to array(*) |
| 220 | |
| 221 | ===================================================================== |
| 222 | */ |
| 223 | |
| 224 | /* Parameter adjustments */ |
| 225 | --cy; |
| 226 | --cx; |
| 227 | |
| 228 | /* Function Body */ |
| 229 | ctemp.r = 0.f, ctemp.i = 0.f; |
| 230 | ret_val->r = 0.f, ret_val->i = 0.f; |
| 231 | if (*n <= 0) { |
| 232 | return ; |
| 233 | } |
| 234 | if (*incx == 1 && *incy == 1) { |
| 235 | goto L20; |
| 236 | } |
| 237 | |
| 238 | /* |
| 239 | code for unequal increments or equal increments |
| 240 | not equal to 1 |
| 241 | */ |
| 242 | |
| 243 | ix = 1; |
| 244 | iy = 1; |
| 245 | if (*incx < 0) { |
| 246 | ix = (-(*n) + 1) * *incx + 1; |
| 247 | } |
| 248 | if (*incy < 0) { |
| 249 | iy = (-(*n) + 1) * *incy + 1; |
| 250 | } |
| 251 | i__1 = *n; |
| 252 | for (i__ = 1; i__ <= i__1; ++i__) { |
| 253 | r_cnjg(&q__3, &cx[ix]); |