SawComplex.cpp

/***************************************************************************
*   Copyright (C) 2006 by Nestor Aguirre                                  *
*   nfaguirrec@unal.edu.co                                                *
*                                                                         *
*   This program 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.                                   *
*                                                                         *
*   This program 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 this program; if not, write to the                         *
*   Free Software Foundation, Inc.,                                       *
*   59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.             *
***************************************************************************/
#include "SawComplex.h"

SawComplex::SawComplex()
:SawObject()
{
        complex = gsl_complex_rect( 0.0, 0.0 ) ;
}

SawComplex::SawComplex( double a, double b,  int representation )
{
        if( representation == CARTESIAN_REPRESENTATION )
                complex = gsl_complex_rect( a, b ) ;
        if( representation == POLAR_REPRESENTATION ) 
                complex = gsl_complex_polar( a, b ) ;
}

SawComplex::SawComplex( const SawComplex& c )
{
        this->complex = c.getGSLComplex() ;
}

SawComplex::SawComplex( double d )
{
        GSL_SET_REAL( &complex, d ) ;
        GSL_SET_IMAG( &complex, 0.0 ) ;
}

SawComplex& SawComplex::operator = ( double d )
{
//      GSL_SET_REAL( &complex, d ) ;
//      GSL_SET_IMAG( &complex, 0.0 ) ;
        this->setReal( d ) ;
        this->setImg( 0.0 ) ;
}

SawComplex& SawComplex::operator = ( const SawComplex& c )
{
//      GSL_SET_REAL( &complex, c.getReal() ) ;
//      GSL_SET_IMAG( &complex, c.getImg() ) ;
        this->setReal( c.getReal() ) ;
        this->setImg( c.getImg() ) ;
}

SawComplex::~SawComplex()
{
}

double SawComplex::getReal() const
{
        return GSL_REAL( complex ) ;
}

double SawComplex::getImg() const
{
        return GSL_IMAG( complex ) ;
}

void SawComplex::setReal( double real )
{
        GSL_SET_REAL( &complex, real ) ;
}

void SawComplex::setImg( double img )
{
        GSL_SET_IMAG( &complex, img ) ;
}

gsl_complex SawComplex::getGSLComplex() const
{
        return complex ;
}

void SawComplex::setGSLComplex( gsl_complex complex )
{
        this->complex = complex ;
}

double SawComplex::getArg() const 
{
        return gsl_complex_arg( complex ) ;
}

SawComplex SawComplex::getConjugate() const 
{
}

SawComplex SawComplex::getInverse() const 
{
}

SawComplex SawComplex::operator + ( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_add( this->getGSLComplex(), complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::operator - ( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_sub( this->getGSLComplex(), complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::operator * ( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_mul( this->getGSLComplex(), complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::operator / ( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_div( this->getGSLComplex(), complex.getGSLComplex() ) ) ;
        return output ;
}

ostream& operator << ( ostream& os, const SawComplex& e )
{
        os.precision( e.precisionOutputStream ) ;
        
        if( e.getImg() != 0.0 )
                os << e.getReal() << " + " << e.getImg() << "i" ;
        else
                os << e.getReal() ;
                                
        return os ;
}

istream& operator >> ( istream& is, SawComplex& e )
{
        double r, i ;
        is >> r ;
        is >> i ;
        e.setReal( r ) ;
        e.setImg( i ) ;
        
        return is ;
}

SawComplex SawComplex::sqrt( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_sqrt( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::pow( SawComplex complex, SawComplex a )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_pow( complex.getGSLComplex(), a.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::pow( SawComplex complex, double a )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_pow_real( complex.getGSLComplex(), a ) ) ;
        return output ;
}

SawComplex SawComplex::exp( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_exp( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::log( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_log( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::log10( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_log10( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::sin( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_sin( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::cos( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_cos( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::tan( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_tan( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::sec( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_sec( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::csc( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_csc( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::cot( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_cot( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::asin( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_arcsin( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::acos( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_arccos( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::atan( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_arctan( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::asec( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_arcsec( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::acsc( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_arccsc( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::acot( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_arccot( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::sinh( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_sinh( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::cosh( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_cosh( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::tanh( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_tanh( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::sech( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_sech( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::csch( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_csch( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::coth( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_coth( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::asinh( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_arcsinh( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::acosh( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_arccosh( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::atanh( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_arctanh( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::asech( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_arcsech( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::acsch( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_arccsch( complex.getGSLComplex() ) ) ;
        return output ;
}

SawComplex SawComplex::acoth( SawComplex complex )
{
        SawComplex output ;
        output.setGSLComplex( gsl_complex_arccoth( complex.getGSLComplex() ) ) ;
        return output ;
}

---