luscheraction.h

Go to the documentation of this file.

#ifndef LUSCHERACTION_H
#define LUSCHERACTION_H

#include "action.h"

class LuscherAction : public Action
{
private:
    int m_muIndex[4];

    int m_nuIndex[4];

    double m_beta;

    double m_multiplicationFactor;

    SU3 m_staple, m_staple1, m_staple2;

    Lattice<SU3> m_latticeStaple,m_tempStaple1,m_tempStaple2;

public:
    LuscherAction();
    ~LuscherAction();
    double getDeltaAction(SU3 U, SU3 UPrime);
    void computeStaple(Lattice<SU3> *lattice, unsigned int i, unsigned int j, unsigned int k, unsigned int l, int mu);
    Lattice<SU3> getActionDerivative(Lattice<SU3> *lattice, int mu);

    inline void updateMuIndex(int mu) {
        for (int i = 0; i < 4; i++)
        {
            m_muIndex[i] = 0;
        }
        m_muIndex[mu] = 1;
    }

    inline void updateNuIndex(int nu) {
        for (int i = 0; i < 4; i++)
        {
            m_nuIndex[i] = 0;
        }
        m_nuIndex[nu] = 1;
    }
};

#endif // LUSCHERACTION_H