StatisticsPoint.h

Go to the documentation of this file.

enum CG {HeavisideSphere, Gaussian, polynomial};

template <StatType T>
class StatisticsVector;

template <StatType T>
struct StatisticsPoint;

template <StatType T>
std::ostream& operator<<(std::ostream& os, const StatisticsPoint<T> &stat);

template <StatType T>
struct StatisticsPoint {
public:
        
        //Constructors and Destructors  
        
        StatisticsPoint() {
                this->Position.set_zero(); 
                mirrorParticle=-1;
                //~ this->set_zero(); 
        }
        
        StatisticsPoint(const StatisticsPoint &other) {*this = other; mirrorParticle=-1;}
        
        //set and get access functions
        static void set_gb(StatisticsVector<T>* new_gb) {gb = new_gb;}
        void set_CG_type(const char* CG_type) {this->gb->set_CG_type(CG_type);}
        CG get_CG_type() const {return this->gb->get_CG_type();}

        void set_w2(Mdouble new_) {this->gb->set_w2(new_);}
        Mdouble get_w2() const {return this->gb->get_w2();}
        Mdouble get_w() const {return this->gb->get_w();}
        Mdouble get_cutoff() {return this->gb->get_cutoff();}
        Mdouble get_cutoff2() {return this->gb->get_cutoff2();}
        Mdouble get_xmaxStat() {return this->gb->get_xmaxStat();}
        Mdouble get_ymaxStat() {return this->gb->get_ymaxStat();}
        Mdouble get_zmaxStat() {return this->gb->get_zmaxStat();}
        Mdouble get_xminStat() {return this->gb->get_xminStat();}
        Mdouble get_yminStat() {return this->gb->get_yminStat();}
        Mdouble get_zminStat() {return this->gb->get_zminStat();}
        void get_n(int& nx_, int& ny_, int& nz_) {this->gb->get_n(nx_,ny_,nz_);}
        Mdouble evaluatePolynomial(Mdouble r) {return gb->evaluatePolynomial(r);}
        Mdouble evaluateIntegral(Mdouble n1, Mdouble n2, Mdouble t) {return gb->evaluateIntegral(n1,n2,t);}
                        
        void set_CG_invvolume();
        Mdouble get_CG_invvolume() {return this->CG_invvolume;}
        void set_Gaussian_invvolume(int dim);
        void set_Heaviside_invvolume();
        
        void set_Polynomial_invvolume(int dim);
                
        void set_Position(Vec3D new_) {this->Position = new_;}
        Vec3D get_Position() const {return this->Position;}

        friend std::ostream& operator<< <T>(std::ostream& os, const StatisticsPoint<T> &stat);

        void set_zero();

        StatisticsPoint<T> getSquared();

        inline StatisticsPoint<T>& operator=(const StatisticsPoint<T> &P);

        inline StatisticsPoint<T>& operator+=(const StatisticsPoint<T> &P);

        inline StatisticsPoint<T>& operator-=(const StatisticsPoint<T> &P);

        inline StatisticsPoint<T>& operator/=(const Mdouble a);

        inline void timeAverage(const int n);

        Mdouble get_distance2(Vec3D &P);

        Mdouble dot(Vec3D &P, Vec3D &Q);

        Mdouble CG_function(Vec3D &PI);

        Mdouble CG_function_2D(Vec3D &PI);

        Mdouble CG_function_1D(Vec3D &PI);

        Vec3D CG_gradient(Vec3D &P, Mdouble phi);

        Vec3D CG_integral_gradient(Vec3D &P1, Vec3D &P2, Vec3D &P1_P2_normal);

        Mdouble CG_integral(Vec3D &P1, Vec3D &P2, Vec3D &P1_P2_normal, Mdouble P1_P2_distance);

        Mdouble CG_integral_2D(Vec3D &P1, Vec3D &P2, Vec3D &P1_P2_normal, Mdouble P1_P2_distance);
        Mdouble CG_integral_3D(Vec3D &P1, Vec3D &P2, Vec3D &P1_P2_normal, Mdouble P1_P2_distance);


        Mdouble CG_integral_1D(Vec3D &P1, Vec3D &P2, Vec3D &P1_P2_normal, Mdouble P1_P2_distance);


        string print() const;
        string print_sqrt() const;
        string write_variable_names();
        string write() const;

        
public:
        //Values of macroscopical fields at Position
        Mdouble Nu; 
        Mdouble Density; 
        Vec3D Momentum; 
        Vec3D DisplacementMomentum; 
        MatrixSymmetric3D Displacement; 
        MatrixSymmetric3D MomentumFlux; 
        MatrixSymmetric3D DisplacementMomentumFlux; 
        Vec3D EnergyFlux; 
        Matrix3D NormalStress; 
        Matrix3D TangentialStress; 
        Vec3D NormalTraction; 
        Vec3D TangentialTraction; 
        MatrixSymmetric3D Fabric; 
        Vec3D CollisionalHeatFlux; 
        Mdouble Dissipation; 
        Mdouble Potential; 
        
        Mdouble CG_invvolume;

        int mirrorParticle;

private:
        static StatisticsVector<T>* gb; 
        Vec3D Position; 
        
};