langevin_prepare.cpp

Go to the documentation of this file.

 
#include "langevin_types.hpp"
#include "langevin_coefficients.hpp"
#include "langevin_base.hpp"

void BaseLangevin::prepare(const Coefficients& coefficients)
{
    // Set "linear" coefficients
    linear_coefficient = coefficients.linear;
    noise_coefficient = coefficients.noise;
    const auto explcdt = exp(-linear_coefficient*dt);
    lambda = (
        (2*linear_coefficient*explcdt)
            / 
        ((1-explcdt)*(noise_coefficient*noise_coefficient))
    );
    lambda_on_explcdt = lambda / explcdt;
    // Set "nonlinear" coefficients in function supplied by child class
    set_nonlinear_coefficients(coefficients);
}