complex_analysis.h

Go to the documentation of this file.

 
 
#ifndef THEORETICA_COMPLEX_FUNCTIONS
#define THEORETICA_COMPLEX_FUNCTIONS
 
#include "./complex.h"
#include "../core/real_analysis.h"
 
 
namespace theoretica {
 
 
    template<typename T>
    inline complex<T> identity(complex<T> z) {
        return z;
    }
 
 
    template<typename T>
    inline complex<T> conjugate(complex<T> z) {
        return complex<T>(z.a, -z.b);
    }
 
 
    template<typename T>
    inline complex<T> inverse(complex<T> z) {
        return conjugate(z) / z.sqr_norm();
    }
 
 
    template<typename T>
    inline complex<T> square(complex<T> z) {
        return complex<T>(
            square(z.Re()) - square(z.Im()),
            2 * z.Re() * z.Im()
        );
    }
 
 
    template<typename T>
    inline complex<T> cube(complex<T> z) {
        return complex<T>(
            cube(z.Re()) - 3 * z.Re() * square(z.Im()),
            3 * square(z.Re()) * z.Im() - cube(z.Im()));
    }
    
 
    template<typename T>
    inline complex<T> exp(complex<T> z) {
        return complex<T>(cos(z.Im()), sin(z.Im())) * exp(z.Re());
    }
 
 
    template<typename T>
    inline real abs(complex<T> z) {
        return z.norm();
    }
 
 
    template<typename T>
    inline complex<T> sin(complex<T> z) {
 
        const complex<T> t = z * complex<T>(0, 1);
        return (exp(t) - exp(-t)) / complex<T>(0, 2);
    }
 
 
    template<typename T>
    inline complex<T> cos(complex<T> z) {
 
        const complex<T> t = z * complex<T>(0, 1);
        return (exp(t) + exp(-t)) / 2.0;
    }
 
 
    template<typename T>
    inline complex<T> tan(complex<T> z) {
 
        const complex<T> t = z * complex<T>(0, 2);
        return (exp(t) + complex<T>(-1, 0)) / (exp(t) + complex<T>(1, 0)) * complex<T>(0, -1);
    }
 
 
    template<typename T>
    inline complex<T> sqrt(complex<T> z) {
 
        if(abs(z.a) < MACH_EPSILON && abs(z.b) < MACH_EPSILON)
            return complex<T>(0);
 
        return complex<T>(
            INVSQR2 * sqrt((z.norm() + z.Re())),
            INVSQR2 * sqrt((z.norm() - z.Re()))  * sgn(z.b));
    }
 
 
    template<typename T>
    inline complex<T> ln(complex<T> z) {
        return complex<T>(ln(z.norm()), z.arg());
    }
 
 
    template<typename T>
    inline complex<T> asin(complex<T> z) {
        return ln(complex<T>(0, 1) * z + sqrt(complex<T>(1, 0) - square(z))) * complex<T>(0, -1);
    }
 
 
    template<typename T>
    inline complex<T> acos(complex<T> z) {
        return ln(z + sqrt(square(z) - 1)) * complex<T>(0, -1);
    }
 
 
    template<typename T>
    inline complex<T> atan(complex<T> z) {
        return ln((complex<T>(0, 1) - z) / (complex<T>(0, 1) + z)) * complex<T>(0, -0.5);
    }
 
 
}
 
#endif