plaquette.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 "Observables/plaquette.h"
#include "Observables/observable.h"
#include "Math/su3.h"
#include "ParallelTools/parallel.h"
#include "Math/lattice.h"
#include <cstdio>
#include <cmath>
#include <string>
#include <iostream>

Plaquette::Plaquette() : Observable(){
    m_name = "Plaquette";
}


void Plaquette::initObservable(Lattice* lattice){
    m_lat = lattice;
    m_size = m_lat->getSize();
    m_norm = 1.0 / 18.0 / m_size[0] / m_size[1] / m_size[2] / m_size[3] / Parallel::numProcs();
}

void Plaquette::compute(){
    // loop through the sublattice
    m_value = 0;
    for(int mu = 0; mu < 4; mu++){
        for(int nu = mu + 1; nu < 4; nu++){
            for(int x = 0; x < m_size[0]; x++){
            for(int y = 0; y < m_size[1]; y++){
            for(int z = 0; z < m_size[2]; z++){
            for(int t = 0; t < m_size[3]; t++){
                m_value += (   (*m_lat)(x,y,z,t)[mu]
                             *   m_lat->shift(x,y,z,t,nu, mu, 1)
                             * ~(m_lat->shift(x,y,z,t,mu, nu, 1))
                             * ~(*m_lat)(x,y,z,t)[nu]
                           ).realTrace();
            }}}}
        }
    }
    m_value *= m_norm;
    gatherResults();
}