ExactOverlap.h

Go to the documentation of this file.

// This file is part of the MercuryDPM project (https://www.mercurydpm.org).
// Copyright (c), The MercuryDPM Developers Team. All rights reserved.
// License: BSD 3-Clause License; see the LICENSE file in the root directory.
 
#ifndef ExactOverlap_H
#define ExactOverlap_H
 
#include <GeneralDefine.h>
#include <iostream>
#include "IntegralType.h"
#include "CG/Coordinates/R.h"
#include "CG/Coordinates/X.h"
#include "CG/Coordinates/Y.h"
#include "CG/Coordinates/Z.h"
#include "CG/Coordinates/RZ.h"
#include "CG/Coordinates/XY.h"
#include "CG/Coordinates/XZ.h"
#include "CG/Coordinates/YZ.h"
#include "CG/Coordinates/XYZ.h"
#include "CG/Coordinates/O.h"
 
class BaseParticle;
 
class BaseInteraction;
 
class DPMBase;
 
class Coordinates;
 
namespace CGFunctions
{
 
template<class Coordinates>
class ExactOverlap
{
public:
    
    typedef Coordinates CoordinatesType;
    
    ExactOverlap();
    
    ExactOverlap(const ExactOverlap& c) = default;
    
    ~ExactOverlap() = default;
    
    void write(std::ostream& os) const;
    
    void setWidth(Mdouble width) {};
    
    Mdouble getWidth() const { return constants::NaN; }
    
    void setStandardDeviation(Mdouble std) {};
    
    void setWidthAndCutoff(Mdouble width, Mdouble cutoff) {};
    
    Mdouble getCutoff() const { return constants::NaN; }
 
    Vec3D getCutoffVector(Mdouble radius, Vec3D cellSize) const { return cellSize + Vec3D(radius, radius, radius); }
    
    // Mdouble getEps() const;
    
    void setEps(Mdouble eps) {};
    
    Mdouble evaluateCGFunction(const Vec3D& position, const Coordinates r) { return constants::NaN; }
    
    Mdouble evaluateCGFunction(const Vec3D& position, const Mdouble rad, const Coordinates r);
    
    Vec3D evaluateCGFunctionDerivatives(const Vec3D& position, const Coordinates& r) 
    { return Vec3D(constants::NaN, constants::NaN, constants::NaN); }
    
    Mdouble evaluateCGIntegral(const BaseInteraction& i, const Coordinates r, IntegralType type = IntegralType::I_TO_P) 
    { return constants::NaN; }
    
};
 
typedef CGFunctions::ExactOverlap<CGCoordinates::O> ExactOverlapO;
typedef CGFunctions::ExactOverlap<CGCoordinates::X> ExactOverlapX;
typedef CGFunctions::ExactOverlap<CGCoordinates::Y> ExactOverlapY;
typedef CGFunctions::ExactOverlap<CGCoordinates::Z> ExactOverlapZ;
typedef CGFunctions::ExactOverlap<CGCoordinates::YZ> ExactOverlapYZ;
typedef CGFunctions::ExactOverlap<CGCoordinates::XZ> ExactOverlapXZ;
typedef CGFunctions::ExactOverlap<CGCoordinates::XY> ExactOverlapXY;
typedef CGFunctions::ExactOverlap<CGCoordinates::XYZ> ExactOverlapXYZ;
typedef CGFunctions::ExactOverlap<CGCoordinates::R> ExactOverlapR;
    
} //namespace CGFunctions
#include "ExactOverlap.hcc"
 
#endif