Complex.h

Go to the documentation of this file.

/* Copyright 2019 Kristofer Björnson
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
 * You may obtain a copy of the License at
 *
 *     http://www.apache.org/licenses/LICENSE-2.0
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 */
 
#ifndef COM_DAFER45_TBTK_COMPLEX
#define COM_DAFER45_TBTK_COMPLEX
 
#include "TBTK/Boolean.h"
#include "TBTK/PseudoSerializable.h"
#include "TBTK/Real.h"
 
#include <complex>
#include <sstream>
 
namespace TBTK{
 
#if TBTK_WRAP_PRIMITIVE_TYPES
 
class Complex : PseudoSerializable{
public:
    Complex(){};
 
    constexpr Complex(std::complex<double> value) : value(value) {}
 
    constexpr Complex(Real real, Real imag) : value(real, imag) {}
 
    constexpr Complex(double value) : value(value) {}
 
    Complex(const std::string &serialization, Serializable::Mode mode);
 
    constexpr operator std::complex<double>() const{
        return value;
    };
 
    Complex& operator=(const std::complex<double> &rhs){
        value = rhs;
 
        return *this;
    }
 
    Complex& operator+=(const Complex &rhs){
        value += rhs.value;
 
        return *this;
    }
 
    friend std::complex<double>& operator+=(
        std::complex<double> &lhs,
        const Complex &rhs
    ){
        lhs += rhs.value;
 
        return lhs;
    }
 
    friend Complex operator+(const Complex &lhs, const double &rhs){
        Complex result(lhs);
        result.value += rhs;
 
        return result;
    }
 
    friend Complex operator+(const double &lhs, const Complex &rhs){
        Complex result(rhs);
        result.value += lhs;
 
        return result;
    }
 
    Complex& operator-=(const Complex &rhs){
        value -= rhs.value;
 
        return *this;
    }
 
    Complex operator-(const Complex &rhs) const{
        Complex result(*this);
        result -= rhs;
 
        return result;
    }
 
    friend Complex operator-(const Complex &lhs, const double &rhs){
        Complex result(lhs);
        result.value -= rhs;
 
        return result;
    }
 
    friend Complex operator-(const double &lhs, const Complex &rhs){
        Complex result(lhs - rhs.value.real(), -rhs.value.imag());
 
        return result;
    }
 
    Complex operator-() const{
        return Complex(-value);
    }
 
    Complex& operator*=(const Complex &rhs){
        value *= rhs.value;
 
        return *this;
    }
 
    friend Complex operator*(const Complex &lhs, const double &rhs){
        Complex result(lhs);
        result.value *= rhs;
 
        return result;
    }
 
    friend Complex operator*(const double &lhs, const Complex &rhs){
        Complex result(rhs);
        result.value *= lhs;
 
        return result;
    }
 
    Complex& operator/=(const Complex &rhs){
        value /= rhs.value;
 
        return *this;
    }
 
    friend Complex operator/(const Complex &lhs, const double &rhs){
        Complex result(lhs);
        result.value /= rhs;
 
        return result;
    }
 
    Boolean operator==(const Complex &rhs) const{
        return value == rhs.value;
    }
 
    Boolean operator!=(const Complex &rhs) const{
        return value != rhs.value;
    }
 
    friend std::ostream& operator<<(std::ostream &os, const Complex &complex){
        os << complex.value;
 
        return os;
    }
 
    friend std::istream& operator>>(std::istream &is, Complex &complex){
        is >> complex.value;
 
        return is;
    }
 
    std::string serialize(Serializable::Mode mode) const;
 
    friend constexpr Real real(const Complex &complex){
        return real(complex.value);
    }
 
    friend constexpr Real imag(const Complex &complex){
        return imag(complex.value);
    }
 
    friend Complex conj(const Complex &complex){
        return conj(complex.value);
    }
 
    friend Real abs(const Complex &complex){
        return abs(complex.value);
    }
private:
    std::complex<double> value;
};
 
inline Complex::Complex(
    const std::string &serialization,
    Serializable::Mode mode
){
    switch(mode){
    case Serializable::Mode::JSON:
    {
        std::stringstream ss(serialization);
        ss >> value;
 
        break;
    }
    default:
        TBTKExit(
            "Complex::Complex()",
            "Only Serializable::Mode::JSON is supported yet.",
            ""
        );
    }
}
 
inline std::string Complex::serialize(Serializable::Mode mode) const{
    switch(mode){
    case Serializable::Mode::JSON:
    {
        std::stringstream ss;
        ss << value;
 
        return ss.str();
    }
    default:
        TBTKExit(
            "Complex::serialize()",
            "Only Serializable::Mode::JSON is supported yet.",
            ""
        );
    }
}
 
#else
    typedef std::complex<double> Complex;
#endif
 
};  //End of namespace TBTK
 
#endif