langevin_base.hpp

Go to the documentation of this file.

 
#ifndef BASE_HPP
#define BASE_HPP
 
#include "langevin_coefficients.hpp"
#include "langevin_parameters.hpp"

class BaseLangevin
{
protected:
    int n_cells;
    grid_t density_grid;
    grid_wiring_t grid_wiring;
   
    double dt;
    double dx;
    double mean_density;

    poisson_dist_t poisson_sampler;
    gamma_dist_t gamma_sampler;
    gaussian_dist_t gaussian_sampler;

    double linear_coefficient;
    double noise_coefficient;
    double lambda;
    double lambda_on_explcdt;

    grid_t k1_grid;
    grid_t k2_grid;
    grid_t k3_grid;
    grid_t aux_grid1;
    grid_t aux_grid2;
 
public:
    BaseLangevin() = default;
    bool construct_grid(const Parameters parameters);
    bool construct_1D_grid(const Parameters parameters);
    bool construct_2D_grid(const Parameters parameters);
    bool initialize_grid(const Parameters parameters, rng_t& rng);
    void prepare(const Coefficients& coefficients);
    bool check_boundary_conditions(const Parameters parameters);
    void apply_boundary_conditions(const Parameters parameters, int i_epoch);
    void integrate_rungekutta(rng_t& rng);
    void integrate_euler(rng_t& rng);
    double get_density_grid_value(const int) const;
    double get_mean_density() const;
    double get_poisson_mean() const;

    virtual void set_nonlinear_coefficients(const Coefficients& coefficients) {};
    virtual double nonlinear_rhs(const int i_cell, const grid_t& field) const 
        { return 0; };
};
 
#endif