puregauge.cpp

Go to the documentation of this file.

/******************************************************************************
* 
* MIT License
* 
* Copyright (c) 2018 Giovanni Pederiva
* 
* Permission is hereby granted, free of charge, to any person obtaining a copy
* of this software and associated documentation files (the "Software"), to deal
* in the Software without restriction, including without limitation the rights
* to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
* copies of the Software, and to permit persons to whom the Software is
* furnished to do so, subject to the following conditions:
* 
* The above copyright notice and this permission notice shall be included in 
* all copies or substantial portions of the Software.
* 
* THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
* IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
* FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
* AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
* LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
* OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
* SOFTWARE.
******************************************************************************/

#include <cstdio>
#include <cmath>
#include "Actions/puregauge.h"
#include "Math/latticemath.h"
#include "Utils/latticeunits.h"
#include "ParallelTools/parallel.h"
#include "Utils/clusterspecifier.h"
#ifndef LACONIA
#include <mpi/mpi.h>
#else
#include <mpi.h>
#endif


PureGauge::PureGauge(double beta){
    m_beta = beta;
    m_staplesLattice = new GluonField(Parallel::latticeSubSize());
    m_constPart = new GluonField(Parallel::latticeSubSize());
    m_omega.allocate(Parallel::latticeSubSize());
    trace.allocate(Parallel::latticeSubSize());
    LatticeUnits::initialize(m_beta);
}

double PureGauge::compute(int x, int y, int z, int t, int mu, SU3& newLink){
    //return -m_beta / 3.0 * getMultSumTrace(newLink - (*m_lat)(x,y,z,t)[mu], m_constPart);
}



void PureGauge::computeStaples(int mu){
    setToZero(m_omega);
    Lattice<SU3> shiftmu, tempStaple;
    for(int nu = 0; nu < 4; nu++){
        if(nu != mu){
            shiftmu = shift((*m_lat)[nu],mu,1);
            tempStaple  = shiftmu;
            tempStaple *= adj( shift((*m_lat)[mu],nu,1));
            tempStaple *= adj( (*m_lat)[nu]);
            m_omega += tempStaple;

            tempStaple  = adj(shiftmu);
            tempStaple *= adj( (*m_lat)[mu]);
            tempStaple *= (*m_lat)[nu];

            m_omega += shift(tempStaple, nu, -1);
        }
    }
}



Lattice<SU3> PureGauge::computeDerivative(int mu){
    computeStaples(mu);
    m_omega = adj(m_omega) * (adj((*m_lat)[mu]));
    m_omega = (adj(m_omega)-m_omega);

    for(int site = 0; site < m_omega.sites; site++){
        double tr = m_omega.at(site).imagTrace()/3.0;
        for(int i = 1; i < 18; i+=8)
            m_omega.at(site).mat[i] -= tr;

        for(int i = 0; i < 18; i+=2){
            double temp = m_omega.at(site).mat[i];
            m_omega.at(site).mat[i] = -m_omega.at(site).mat[i+1]*0.5;
            m_omega.at(site).mat[i+1] = temp*0.5;
        }
    }
    return m_omega;
}



void PureGauge::computeStaplez(GluonField* lattice){
}