Quartic.h

#ifndef gridripper_phys_scalar_Quartic_h
#define gridripper_phys_scalar_Quartic_h

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

namespace gridripper {

class Parameters;

namespace phys { namespace scalar {

using namespace std;

class Quartic: public ScalarPotential
{
private:
    GReal_t massSq;
    GReal_t c4;

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

    Quartic(GReal_t m, GReal_t c4): massSq(m*m), c4(c4) { }

    GReal_t eval_V(GReal_t psi, GReal_t T) {
        double psisq = psi*psi;
        return 0.5*(massSq + 0.5*c4*psisq)*psisq;
    }

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

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

#endif /* gridripper_phys_scalar_Quartic_h */