Program Listing for File system.hpp

↰ Return to documentation for file (/home/docs/checkouts/readthedocs.org/user_builds/oscode/checkouts/stable/include/system.hpp)

#pragma once
#include "interpolator.hpp"

class de_system
    {
    private:
        int even_;

    public:
        template<typename X, typename Y, typename Z, typename X_it>de_system(X
        &ts, Y &ws, Z &gs, X_it x_it, int size, bool isglogw=false, bool
        islogg=false, int even=0, int check_grid=0);
        de_system(std::complex<double> (*)(double, void*), std::complex<double> (*)(double, void*), void*);
        de_system(std::complex<double> (*)(double), std::complex<double> (*)(double));
        de_system();
        std::function<std::complex<double>(double)> w;
        std::function<std::complex<double>(double)> g;
        LinearInterpolator<> Winterp;
        LinearInterpolator<> Ginterp;
        bool islogg_, islogw_;
        bool grid_fine_enough = 1;
        bool is_interpolated;
};

de_system::de_system(){
}

template<typename X, typename Y, typename Z, typename X_it>
de_system::de_system(X &ts, Y &ws, Z &gs, X_it x_it, int size, bool islogw, bool
islogg, int even, int check_grid){

    is_interpolated = 1;
    even_ = even;
    islogg_ = islogg;
    islogw_ = islogw;

    LinearInterpolator<X,Y,X_it> winterp(ts,ws,even_);
    LinearInterpolator<X,Z,X_it> ginterp(ts,gs,even_);

    Winterp = winterp;
    Ginterp = ginterp;

    if(even_==0){
        Winterp.set_interp_start(x_it);
        Ginterp.set_interp_start(x_it);
        Winterp.set_interp_bounds(ts,ts+size-1);
        Ginterp.set_interp_bounds(ts,ts+size-1);
    }

    if(check_grid == 1){
        int w_is_fine = Winterp.check_grid_fineness(size);
        int g_is_fine = Ginterp.check_grid_fineness(size);
        if(w_is_fine==1 && g_is_fine==1)
            grid_fine_enough = 1;
        else
            grid_fine_enough = 0;
    }

    if(islogw)
        w = std::bind(&LinearInterpolator<X,Y,X_it>::expit, Winterp,
        std::placeholders::_1);
    else
        w = std::bind(&LinearInterpolator<X,Y,X_it>::operator(), Winterp,
        std::placeholders::_1);
    if(islogg)
        g = std::bind(&LinearInterpolator<X,Z,X_it>::expit, Ginterp,
        std::placeholders::_1);
    else
        g = std::bind(&LinearInterpolator<X,Z,X_it>::operator(), Ginterp,
        std::placeholders::_1);


}

de_system::de_system(std::complex<double> (*W)(double, void*),
std::complex<double> (*G)(double, void*), void *p){

    is_interpolated = 0;
    w = [W, p](double x){ return W(x, p); };
    g = [G, p](double x){ return G(x, p); };
};

de_system::de_system(std::complex<double> (*W)(double), std::complex<double> (*G)(double)){

    is_interpolated = 0;
    w = W;
    g = G;
};