Massive.h

#ifndef gridripper_phys_scalar_Massive_h
#define gridripper_phys_scalar_Massive_h

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

namespace gridripper {

class Parameters;

namespace phys { namespace scalar {

using namespace std;

class Massive: public ScalarPotential
{
private:
    GReal_t massSq;

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

    Massive(GReal_t m) {
        massSq = m*m;
    }

    GReal_t eval_V(GReal_t psi, GReal_t T) {
        return 0.5*massSq*psi*psi;
    }

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

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

#endif /* gridripper_phys_scalar_Massive_h */