htswanalysis/MACS/lib/gsl/gsl-1.11/complex/gsl_complex_math.h

   1 /* complex/gsl_complex_math.h
   2  *
   3  * Copyright (C) 1996, 1997, 1998, 1999, 2000, 2004, 2007 Jorma Olavi Tähtinen, Brian Gough
   4  *
   5  * This program is free software; you can redistribute it and/or modify
   6  * it under the terms of the GNU General Public License as published by
   7  * the Free Software Foundation; either version 3 of the License, or (at
   8  * your option) any later version.
   9  *
  10  * This program is distributed in the hope that it will be useful, but
  11  * WITHOUT ANY WARRANTY; without even the implied warranty of
  12  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  13  * General Public License for more details.
  14  *
  15  * You should have received a copy of the GNU General Public License
  16  * along with this program; if not, write to the Free Software
  17  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
  18  */
  19
  20 #ifndef __GSL_COMPLEX_MATH_H__
  21 #define __GSL_COMPLEX_MATH_H__
  22 #include <gsl/gsl_complex.h>
  23
  24 #undef __BEGIN_DECLS
  25 #undef __END_DECLS
  26 #ifdef __cplusplus
  27 #define __BEGIN_DECLS extern "C" {
  28 #define __END_DECLS }
  29 #else
  30 #define __BEGIN_DECLS           /* empty */
  31 #define __END_DECLS             /* empty */
  32 #endif
  33
  34 __BEGIN_DECLS
  35
  36 /* Complex numbers */
  37
  38 gsl_complex gsl_complex_rect (double x, double y);  /* r= real+i*imag */
  39 gsl_complex gsl_complex_polar (double r, double theta); /* r= r e^(i theta) */
  40
  41 #ifdef HAVE_INLINE
  42 extern inline gsl_complex
  43 gsl_complex_rect (double x, double y)
  44 {                               /* return z = x + i y */
  45   gsl_complex z;
  46   GSL_SET_COMPLEX (&z, x, y);
  47   return z;
  48 }
  49 #endif
  50
  51 #define GSL_COMPLEX_ONE (gsl_complex_rect(1.0,0.0))
  52 #define GSL_COMPLEX_ZERO (gsl_complex_rect(0.0,0.0))
  53 #define GSL_COMPLEX_NEGONE (gsl_complex_rect(-1.0,0.0))
  54
  55 /* Properties of complex numbers */
  56
  57 double gsl_complex_arg (gsl_complex z); /* return arg(z), -pi< arg(z) <=+pi */
  58 double gsl_complex_abs (gsl_complex z);   /* return |z|   */
  59 double gsl_complex_abs2 (gsl_complex z);  /* return |z|^2 */
  60 double gsl_complex_logabs (gsl_complex z); /* return log|z| */
  61
  62 /* Complex arithmetic operators */
  63
  64 gsl_complex gsl_complex_add (gsl_complex a, gsl_complex b);  /* r=a+b */
  65 gsl_complex gsl_complex_sub (gsl_complex a, gsl_complex b);  /* r=a-b */
  66 gsl_complex gsl_complex_mul (gsl_complex a, gsl_complex b);  /* r=a*b */
  67 gsl_complex gsl_complex_div (gsl_complex a, gsl_complex b);  /* r=a/b */
  68
  69 gsl_complex gsl_complex_add_real (gsl_complex a, double x);  /* r=a+x */
  70 gsl_complex gsl_complex_sub_real (gsl_complex a, double x);  /* r=a-x */
  71 gsl_complex gsl_complex_mul_real (gsl_complex a, double x);  /* r=a*x */
  72 gsl_complex gsl_complex_div_real (gsl_complex a, double x);  /* r=a/x */
  73
  74 gsl_complex gsl_complex_add_imag (gsl_complex a, double y);  /* r=a+iy */
  75 gsl_complex gsl_complex_sub_imag (gsl_complex a, double y);  /* r=a-iy */
  76 gsl_complex gsl_complex_mul_imag (gsl_complex a, double y);  /* r=a*iy */
  77 gsl_complex gsl_complex_div_imag (gsl_complex a, double y);  /* r=a/iy */
  78
  79 gsl_complex gsl_complex_conjugate (gsl_complex z);  /* r=conj(z) */
  80 gsl_complex gsl_complex_inverse (gsl_complex a);    /* r=1/a */
  81 gsl_complex gsl_complex_negative (gsl_complex a);    /* r=-a */
  82
  83 /* Elementary Complex Functions */
  84
  85 gsl_complex gsl_complex_sqrt (gsl_complex z);  /* r=sqrt(z) */
  86 gsl_complex gsl_complex_sqrt_real (double x);  /* r=sqrt(x) (x<0 ok) */
  87
  88 gsl_complex gsl_complex_pow (gsl_complex a, gsl_complex b);  /* r=a^b */
  89 gsl_complex gsl_complex_pow_real (gsl_complex a, double b);  /* r=a^b */
  90
  91 gsl_complex gsl_complex_exp (gsl_complex a);    /* r=exp(a) */
  92 gsl_complex gsl_complex_log (gsl_complex a);    /* r=log(a) (base e) */
  93 gsl_complex gsl_complex_log10 (gsl_complex a);  /* r=log10(a) (base 10) */
  94 gsl_complex gsl_complex_log_b (gsl_complex a, gsl_complex b);   /* r=log_b(a) (base=b) */
  95
  96 /* Complex Trigonometric Functions */
  97
  98 gsl_complex gsl_complex_sin (gsl_complex a);  /* r=sin(a) */
  99 gsl_complex gsl_complex_cos (gsl_complex a);  /* r=cos(a) */
 100 gsl_complex gsl_complex_sec (gsl_complex a);  /* r=sec(a) */
 101 gsl_complex gsl_complex_csc (gsl_complex a);  /* r=csc(a) */
 102 gsl_complex gsl_complex_tan (gsl_complex a);  /* r=tan(a) */
 103 gsl_complex gsl_complex_cot (gsl_complex a);  /* r=cot(a) */
 104
 105 /* Inverse Complex Trigonometric Functions */
 106
 107 gsl_complex gsl_complex_arcsin (gsl_complex a);  /* r=arcsin(a) */
 108 gsl_complex gsl_complex_arcsin_real (double a);  /* r=arcsin(a) */
 109 gsl_complex gsl_complex_arccos (gsl_complex a);  /* r=arccos(a) */
 110 gsl_complex gsl_complex_arccos_real (double a);  /* r=arccos(a) */
 111 gsl_complex gsl_complex_arcsec (gsl_complex a);  /* r=arcsec(a) */
 112 gsl_complex gsl_complex_arcsec_real (double a);  /* r=arcsec(a) */
 113 gsl_complex gsl_complex_arccsc (gsl_complex a);  /* r=arccsc(a) */
 114 gsl_complex gsl_complex_arccsc_real (double a);  /* r=arccsc(a) */
 115 gsl_complex gsl_complex_arctan (gsl_complex a);  /* r=arctan(a) */
 116 gsl_complex gsl_complex_arccot (gsl_complex a);  /* r=arccot(a) */
 117
 118 /* Complex Hyperbolic Functions */
 119
 120 gsl_complex gsl_complex_sinh (gsl_complex a);  /* r=sinh(a) */
 121 gsl_complex gsl_complex_cosh (gsl_complex a);  /* r=coshh(a) */
 122 gsl_complex gsl_complex_sech (gsl_complex a);  /* r=sech(a) */
 123 gsl_complex gsl_complex_csch (gsl_complex a);  /* r=csch(a) */
 124 gsl_complex gsl_complex_tanh (gsl_complex a);  /* r=tanh(a) */
 125 gsl_complex gsl_complex_coth (gsl_complex a);  /* r=coth(a) */
 126
 127 /* Inverse Complex Hyperbolic Functions */
 128
 129 gsl_complex gsl_complex_arcsinh (gsl_complex a);  /* r=arcsinh(a) */
 130 gsl_complex gsl_complex_arccosh (gsl_complex a);  /* r=arccosh(a) */
 131 gsl_complex gsl_complex_arccosh_real (double a);  /* r=arccosh(a) */
 132 gsl_complex gsl_complex_arcsech (gsl_complex a);  /* r=arcsech(a) */
 133 gsl_complex gsl_complex_arccsch (gsl_complex a);  /* r=arccsch(a) */
 134 gsl_complex gsl_complex_arctanh (gsl_complex a);  /* r=arctanh(a) */
 135 gsl_complex gsl_complex_arctanh_real (double a);  /* r=arctanh(a) */
 136 gsl_complex gsl_complex_arccoth (gsl_complex a);  /* r=arccoth(a) */
 137
 138 __END_DECLS
 139
 140 #endif /* __GSL_COMPLEX_MATH_H__ */