diff options
| author | shmel1k <[email protected]> | 2022-09-02 12:44:59 +0300 | 
|---|---|---|
| committer | shmel1k <[email protected]> | 2022-09-02 12:44:59 +0300 | 
| commit | 90d450f74722da7859d6f510a869f6c6908fd12f (patch) | |
| tree | 538c718dedc76cdfe37ad6d01ff250dd930d9278 /contrib/libs/clapack/zladiv.c | |
| parent | 01f64c1ecd0d4ffa9e3a74478335f1745f26cc75 (diff) | |
[] add metering mode to CLI
Diffstat (limited to 'contrib/libs/clapack/zladiv.c')
| -rw-r--r-- | contrib/libs/clapack/zladiv.c | 75 | 
1 files changed, 75 insertions, 0 deletions
diff --git a/contrib/libs/clapack/zladiv.c b/contrib/libs/clapack/zladiv.c new file mode 100644 index 00000000000..d92be5a1cbf --- /dev/null +++ b/contrib/libs/clapack/zladiv.c @@ -0,0 +1,75 @@ +/* zladiv.f -- translated by f2c (version 20061008). +   You must link the resulting object file with libf2c: +	on Microsoft Windows system, link with libf2c.lib; +	on Linux or Unix systems, link with .../path/to/libf2c.a -lm +	or, if you install libf2c.a in a standard place, with -lf2c -lm +	-- in that order, at the end of the command line, as in +		cc *.o -lf2c -lm +	Source for libf2c is in /netlib/f2c/libf2c.zip, e.g., + +		http://www.netlib.org/f2c/libf2c.zip +*/ + +#include "f2c.h" +#include "blaswrap.h" + +/* Double Complex */ VOID zladiv_(doublecomplex * ret_val, doublecomplex *x,  +	doublecomplex *y) +{ +    /* System generated locals */ +    doublereal d__1, d__2, d__3, d__4; +    doublecomplex z__1; + +    /* Builtin functions */ +    double d_imag(doublecomplex *); + +    /* Local variables */ +    doublereal zi, zr; +    extern /* Subroutine */ int dladiv_(doublereal *, doublereal *,  +	    doublereal *, doublereal *, doublereal *, doublereal *); + + +/*  -- LAPACK auxiliary routine (version 3.2) -- */ +/*     Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. */ +/*     November 2006 */ + +/*     .. Scalar Arguments .. */ +/*     .. */ + +/*  Purpose */ +/*  ======= */ + +/*  ZLADIV := X / Y, where X and Y are complex.  The computation of X / Y */ +/*  will not overflow on an intermediary step unless the results */ +/*  overflows. */ + +/*  Arguments */ +/*  ========= */ + +/*  X       (input) COMPLEX*16 */ +/*  Y       (input) COMPLEX*16 */ +/*          The complex scalars X and Y. */ + +/*  ===================================================================== */ + +/*     .. Local Scalars .. */ +/*     .. */ +/*     .. External Subroutines .. */ +/*     .. */ +/*     .. Intrinsic Functions .. */ +/*     .. */ +/*     .. Executable Statements .. */ + +    d__1 = x->r; +    d__2 = d_imag(x); +    d__3 = y->r; +    d__4 = d_imag(y); +    dladiv_(&d__1, &d__2, &d__3, &d__4, &zr, &zi); +    z__1.r = zr, z__1.i = zi; +     ret_val->r = z__1.r,  ret_val->i = z__1.i; + +    return ; + +/*     End of ZLADIV */ + +} /* zladiv_ */  | 
