SineGordon.h

#ifndef gridripper_phys_scalar_SineGordon_h
#define gridripper_phys_scalar_SineGordon_h

#include "ScalarPotential.h"
#include <gridripper/lang/IllegalArgumentException.h>
#include <cmath>

namespace gridripper {

class Parameters;

namespace phys { namespace scalar {

using namespace std;

class SineGordon: public ScalarPotential
{
private:
    GReal_t mass;
    GReal_t massSq;

public:
    SineGordon(const string& str, const Parameters& p)
        throw(IllegalArgumentException&);

    GReal_t eval_V(GReal_t psi, GReal_t T) {
        return massSq*(1 - cos(psi));
    }

    GReal_t eval_dVdpsi(GReal_t psi, GReal_t T) {
        return massSq*sin(psi);
    }
};

} } } // namespace gridripper::phys::scalar

#endif /* gridripper_phys_scalar_SineGordon_h */