Double Complex */
| 9372 | } /* zlacrm_ */ |
| 9373 | |
| 9374 | /* Double Complex */ VOID zladiv_(doublecomplex * ret_val, doublecomplex *x, |
| 9375 | doublecomplex *y) |
| 9376 | { |
| 9377 | /* System generated locals */ |
| 9378 | doublereal d__1, d__2, d__3, d__4; |
| 9379 | doublecomplex z__1; |
| 9380 | |
| 9381 | /* Local variables */ |
| 9382 | static doublereal zi, zr; |
| 9383 | extern /* Subroutine */ int dladiv_(doublereal *, doublereal *, |
| 9384 | doublereal *, doublereal *, doublereal *, doublereal *); |
| 9385 | |
| 9386 | |
| 9387 | /* |
| 9388 | -- LAPACK auxiliary routine (version 3.2) -- |
| 9389 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 9390 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 9391 | November 2006 |
| 9392 | |
| 9393 | |
| 9394 | Purpose |
| 9395 | ======= |
| 9396 | |
| 9397 | ZLADIV := X / Y, where X and Y are complex. The computation of X / Y |
| 9398 | will not overflow on an intermediary step unless the results |
| 9399 | overflows. |
| 9400 | |
| 9401 | Arguments |
| 9402 | ========= |
| 9403 | |
| 9404 | X (input) COMPLEX*16 |
| 9405 | Y (input) COMPLEX*16 |
| 9406 | The complex scalars X and Y. |
| 9407 | |
| 9408 | ===================================================================== |
| 9409 | */ |
| 9410 | |
| 9411 | |
| 9412 | d__1 = x->r; |
| 9413 | d__2 = d_imag(x); |
| 9414 | d__3 = y->r; |
| 9415 | d__4 = d_imag(y); |
| 9416 | dladiv_(&d__1, &d__2, &d__3, &d__4, &zr, &zi); |
| 9417 | z__1.r = zr, z__1.i = zi; |
| 9418 | ret_val->r = z__1.r, ret_val->i = z__1.i; |
| 9419 | |
| 9420 | return ; |
| 9421 | |
| 9422 | /* End of ZLADIV */ |
| 9423 | |
| 9424 | } /* zladiv_ */ |
| 9425 | |
| 9426 | /* Subroutine */ int zlaed0_(integer *qsiz, integer *n, doublereal *d__, |
| 9427 | doublereal *e, doublecomplex *q, integer *ldq, doublecomplex *qstore, |