egfrd/GreensFunction3DSym.hpp at develop · TheOtherOne/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
66
67
#if !defined( __FREEGREENSFUNCTION )
#define __FREEGREENSFUNCTION

#include "compat.h"

#include <gsl/gsl_integration.h>

#include "Logger.hpp"
#include "GreensFunction.hpp"

/**
  Green's Function for a free diffusion particle.
*/

class GreensFunction3DSym
    :
    public GreensFunction
{

private:

    static const Real TOLERANCE = 1e-8;

    static const Real H = 6;

public:

    GreensFunction3DSym( const Real D )
        :
        GreensFunction( D )
    {
        ; // do nothing
    }


    ~GreensFunction3DSym()
    {
        ; // do nothing
    }

    Real drawTime( const Real ) const
    {
        return INFINITY;
    }

    Real drawR( const Real rnd, const Real t ) const;

    Real p_r( const Real r, const Real t ) const;

    Real ip_r( const Real r, const Real t ) const;


    std::string dump() const;

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

private:

    static Logger& log_;
};


#endif // __FREEGREENSFUNCTION