00001 /* zladiv.f -- translated by f2c (version 20061008). 00002 You must link the resulting object file with libf2c: 00003 on Microsoft Windows system, link with libf2c.lib; 00004 on Linux or Unix systems, link with .../path/to/libf2c.a -lm 00005 or, if you install libf2c.a in a standard place, with -lf2c -lm 00006 -- in that order, at the end of the command line, as in 00007 cc *.o -lf2c -lm 00008 Source for libf2c is in /netlib/f2c/libf2c.zip, e.g., 00009 00010 http://www.netlib.org/f2c/libf2c.zip 00011 */ 00012 00013 #include "f2c.h" 00014 #include "blaswrap.h" 00015 00016 /* Double Complex */ VOID zladiv_(doublecomplex * ret_val, doublecomplex *x, 00017 doublecomplex *y) 00018 { 00019 /* System generated locals */ 00020 doublereal d__1, d__2, d__3, d__4; 00021 doublecomplex z__1; 00022 00023 /* Builtin functions */ 00024 double d_imag(doublecomplex *); 00025 00026 /* Local variables */ 00027 doublereal zi, zr; 00028 extern /* Subroutine */ int dladiv_(doublereal *, doublereal *, 00029 doublereal *, doublereal *, doublereal *, doublereal *); 00030 00031 00032 /* -- LAPACK auxiliary routine (version 3.2) -- */ 00033 /* Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. */ 00034 /* November 2006 */ 00035 00036 /* .. Scalar Arguments .. */ 00037 /* .. */ 00038 00039 /* Purpose */ 00040 /* ======= */ 00041 00042 /* ZLADIV := X / Y, where X and Y are complex. The computation of X / Y */ 00043 /* will not overflow on an intermediary step unless the results */ 00044 /* overflows. */ 00045 00046 /* Arguments */ 00047 /* ========= */ 00048 00049 /* X (input) COMPLEX*16 */ 00050 /* Y (input) COMPLEX*16 */ 00051 /* The complex scalars X and Y. */ 00052 00053 /* ===================================================================== */ 00054 00055 /* .. Local Scalars .. */ 00056 /* .. */ 00057 /* .. External Subroutines .. */ 00058 /* .. */ 00059 /* .. Intrinsic Functions .. */ 00060 /* .. */ 00061 /* .. Executable Statements .. */ 00062 00063 d__1 = x->r; 00064 d__2 = d_imag(x); 00065 d__3 = y->r; 00066 d__4 = d_imag(y); 00067 dladiv_(&d__1, &d__2, &d__3, &d__4, &zr, &zi); 00068 z__1.r = zr, z__1.i = zi; 00069 ret_val->r = z__1.r, ret_val->i = z__1.i; 00070 00071 return ; 00072 00073 /* End of ZLADIV */ 00074 00075 } /* zladiv_ */