egfrd/GreensFunction3D.hpp at develop · Jintram/egfrd · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
#if !defined( __FREEPAIRGREENSFUNCTION )
#define __FREEPAIRGREENSFUNCTION

#include "compat.h"

#include <gsl/gsl_integration.h>

#include "Logger.hpp"
#include "PairGreensFunction.hpp"

/**
   Pair Green's function for the case where the pair never interact.

   Therefore, drawTime() always returns +INFINITY.
   kf == sigma == 0.
*/

class GreensFunction3D: public PairGreensFunction
{

private:

    static const Real TOLERANCE = 1e-8;

    static const Real H = 7;

public:

    GreensFunction3D(Real D, Real r0)
        : PairGreensFunction(D, 0.0, r0, 0.0)
    {
        ; // do nothing
    }


    virtual ~GreensFunction3D();

    virtual Real drawTime(Real rnd) const;

    Real drawR(Real rnd, Real t) const;

    Real drawTheta(Real rnd, Real r, Real t) const;

    Real p_r(Real r, Real t) const;

    Real ip_r(Real r, Real t ) const;


    Real p_theta(Real theta, Real r, Real t) const;

    Real ip_theta(Real theta, Real r, Real t ) const;

    std::string dump() const;

    const char* getName() const
    {
        return "GreensFunction3D";
    }

private:
    static Logger& log_;
};


#endif // __PLAINPAIRGREENSFUNCTION