chainlink_complex.h

Go to the documentation of this file.

/****************************************************************************
 This file is part of ChainLink
 Copyright (C) 2007 Jeremy Magland (Jeremy.Magland@gmail.com)

 ChainLink is free software; you can redistribute it and/or modify
 it under the terms of the GNU General Public License as published by
 the Free Software Foundation; either version 2 of the License, or
 (at your option) any later version.

 ChainLink is distributed in the hope that it will be useful,
 but WITHOUT ANY WARRANTY; without even the implied warranty of
 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 GNU General Public License for more details.

 You should have received a copy of the GNU General Public License
 along with ChainLink; if not, write to the Free Software
 Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA

*****************************************************************************/

#ifndef CHAINLINK_COMPLEX_H
#define CHAINLINK_COMPLEX_H

#include <math.h>
#include <stdlib.h>
#include <stdio.h>

//comment out one of the following lines
//choose floating point or double precision
        typedef float real;
        //typedef double real;
        
//comment out one of the following lines
//on 64-bit systems choose int
        typedef long int32;
        //typedef int int32;

struct complex {
        real re,im;
};

complex cplex(real re,real im);

//operatations on complex numbers
complex operator*(complex X1, complex X2);
complex operator/(complex X1,complex X2);
complex operator+(complex X1, complex X2);
complex operator-(complex X1,complex X2);
complex operator*(complex X1, real X2);
complex operator/(complex X1,real X2);
complex operator+(complex X1, real X2);
complex operator-(complex X1,real X2);
complex operator/(real X1, complex X2);
bool operator==(complex X1, complex X2);
bool operator==(complex X1, real X2);
bool operator==(real X1,complex X2);

real abs(complex X);

complex conj(complex X);

complex exp(complex X);

complex log(complex X);

void get_complex_string(char *ret,complex X);

complex power(const complex &x, const complex &y);

#endif

---