Double Complex */
| 15353 | } /* zdotc_ */ |
| 15354 | |
| 15355 | /* Double Complex */ VOID zdotu_(doublecomplex * ret_val, integer *n, |
| 15356 | doublecomplex *zx, integer *incx, doublecomplex *zy, integer *incy) |
| 15357 | { |
| 15358 | /* System generated locals */ |
| 15359 | integer i__1, i__2, i__3; |
| 15360 | doublecomplex z__1, z__2; |
| 15361 | |
| 15362 | /* Local variables */ |
| 15363 | static integer i__, ix, iy; |
| 15364 | static doublecomplex ztemp; |
| 15365 | |
| 15366 | |
| 15367 | /* |
| 15368 | Purpose |
| 15369 | ======= |
| 15370 | |
| 15371 | ZDOTU forms the dot product of two vectors. |
| 15372 | |
| 15373 | Further Details |
| 15374 | =============== |
| 15375 | |
| 15376 | jack dongarra, 3/11/78. |
| 15377 | modified 12/3/93, array(1) declarations changed to array(*) |
| 15378 | |
| 15379 | ===================================================================== |
| 15380 | */ |
| 15381 | |
| 15382 | /* Parameter adjustments */ |
| 15383 | --zy; |
| 15384 | --zx; |
| 15385 | |
| 15386 | /* Function Body */ |
| 15387 | ztemp.r = 0., ztemp.i = 0.; |
| 15388 | ret_val->r = 0., ret_val->i = 0.; |
| 15389 | if (*n <= 0) { |
| 15390 | return ; |
| 15391 | } |
| 15392 | if (*incx == 1 && *incy == 1) { |
| 15393 | goto L20; |
| 15394 | } |
| 15395 | |
| 15396 | /* |
| 15397 | code for unequal increments or equal increments |
| 15398 | not equal to 1 |
| 15399 | */ |
| 15400 | |
| 15401 | ix = 1; |
| 15402 | iy = 1; |
| 15403 | if (*incx < 0) { |
| 15404 | ix = (-(*n) + 1) * *incx + 1; |
| 15405 | } |
| 15406 | if (*incy < 0) { |
| 15407 | iy = (-(*n) + 1) * *incy + 1; |
| 15408 | } |
| 15409 | i__1 = *n; |
| 15410 | for (i__ = 1; i__ <= i__1; ++i__) { |
| 15411 | i__2 = ix; |
| 15412 | i__3 = iy; |