functions.h

Go to the documentation of this file.

 
#pragma once
#include "../backend/config.h" // DG_DEVICE
 
 
namespace dg{
 
 
template<class T, class ...Ts>
DG_DEVICE
T zero(T x, Ts ...xs){return T(0);}
 
template<class T, class ...Ts>
DG_DEVICE
T one(T x, Ts ...xs){return T(1);}
 
struct ZERO
{
    template<class T, class ...Ts>
    DG_DEVICE
    T operator()(T x, Ts ...xs)const{return T(0);}
};
 
struct ONE
{
    template<class T, class ...Ts>
    DG_DEVICE
    T operator()(T x, Ts ...xs)const{return T(1);}
};
 
struct CONSTANT
{
    CONSTANT( double cte): m_value(cte){}
 
    template<class T, class ...Ts>
    DG_DEVICE
    T operator()(T x, Ts ...xs)const{return T(m_value);}
    private:
    double m_value;
};
 
DG_DEVICE inline double cooX1d( double x) {return x;}
DG_DEVICE inline double cooX2d( double x, double y) {return x;}
DG_DEVICE inline double cooX3d( double x, double y, double z) {return x;}
 
DG_DEVICE inline double cooY2d( double x, double y) {return y;}
DG_DEVICE inline double cooY3d( double x, double y, double z) {return y;}
DG_DEVICE inline double cooZ3d( double x, double y, double z) {return z;}
 
 
DG_DEVICE inline double cooRZP2X( double R, double Z, double P){ return R*sin(P);}
DG_DEVICE inline double cooRZP2Y( double R, double Z, double P){ return R*cos(P);}
DG_DEVICE inline double cooRZP2Z( double R, double Z, double P){ return Z;}
 
 
} //namespace dg