00001 /*************************************************************************** 00002 * Copyright (C) 2010 by Matthias Ihrke * 00003 * ihrke@nld.ds.mpg.de 00004 * * 00005 * This program is free software; you can redistribute it and/or modify * 00006 * it under the terms of the GNU General Public License as published by * 00007 * the Free Software Foundation; either version 2 of the License, or * 00008 * (at your option) any later version. * 00009 * * 00010 * This program is distributed in the hope that it will be useful, * 00011 * but WITHOUT ANY WARRANTY; without even the implied warranty of * 00012 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * 00013 * GNU General Public License for more details. * 00014 * * 00015 * You should have received a copy of the GNU General Public License * 00016 * aint with this program; if not, write to the * 00017 * Free Software Foundation, Inc., * 00018 * 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. * 00019 ***************************************************************************/ 00020 00027 #ifndef COMPLEX_H 00028 # define COMPLEX_H 00029 #include "definitions.h" 00030 00031 #ifdef __cplusplus 00032 extern "C" { 00033 #endif 00034 00035 typedef struct { 00036 double re; 00037 double im; 00038 } Complex; 00039 00040 00041 /* ---------------------------------------------------------------------------- 00042 -- Complex Arithmetic 00043 ---------------------------------------------------------------------------- */ 00044 00045 Complex complex ( double re, double im ); 00046 Complex complex_add ( Complex a, Complex b ); 00047 Complex complex_add_dbl ( Complex a, double b ); 00048 Complex complex_sub ( Complex a, Complex b ); 00049 Complex complex_mul ( Complex a, Complex b ); 00050 Complex complex_mul_double( Complex a, double b ); 00051 double complex_abs ( Complex a ); 00052 Complex complex_exp ( Complex a ); 00053 Complex complex_conj( Complex a ); 00054 Complex complex_neg ( Complex a ); 00055 Complex complex_div ( Complex a, Complex b); 00056 Complex complex_sqrt( Complex x ); 00057 Complex complex_bilinear_transform(Complex pz); 00058 00059 00060 #ifdef __cplusplus 00061 } 00062 #endif 00063 00064 #endif /* COMPLEX_H */