bias_potential.hpp

Go to the documentation of this file.

#ifndef BIAS_POTENTIAL_HPP
#define BIAS_POTENTIAL_HPP

#include "grid.hpp"

// bias potential for artificial GAM drive
class BiasPotential{
    public:
    bool kick;         
    double tanh_amp;   
    double linear_amp; 
    double psic;       
    double inv_psiw;       
    int freq;          
    int start;
    double time0;
    double inv_psi_norm;

    BiasPotential(NLReader::NamelistReader& nlr, double dt, double psi_norm);

    KOKKOS_INLINE_FUNCTION double tanh_bias(double psi) const{
        return 0.5*tanh_amp*(tanh( 2.0*(psi - psic)*inv_psiw)-1.0);
    }

    KOKKOS_INLINE_FUNCTION double linear_bias(double psi) const{
        return linear_amp*(1.0-psi*inv_psi_norm);
    }

    void add_potential(const Grid<DeviceType>& grid, double time, const View<double*, CLayout, DeviceType>& pot) const;
};

#endif