u1gaugefield2d.cpp

Go to the documentation of this file.

/**
 * @file u1gaugefield2d.cpp
 * @brief Implements the U(1) lattice gauge field and gauge-covariant link construction.
 * @ingroup qpe_quantum
 *
 * Source: @c core/quantum/u1gaugefield2d.cpp. This file is part of Quantum Phenomena Explorer (QPE),
 * maintained by Officina Turini. See LICENSE_IT.txt / LICENSE.txt for the
 * QPE freeware source-available terms and THIRD_PARTY_NOTICES.txt for
 * third-party components.
 */
 
#include "u1gaugefield2d.h"
#include "core/math/complexfield2d.h"
 
namespace qpe {
 
U1GaugeField2D::U1GaugeField2D(std::size_t width, std::size_t height,
                               float dx, float dy)
{
    resize(width, height, dx, dy);
}
 
void U1GaugeField2D::resize(std::size_t width, std::size_t height,
                            float dx, float dy)
{
    width_ = width;
    height_ = height;
    dx_ = dx;
    dy_ = dy;
    xLinks_.assign(width_ * height_, Link(1.0f, 0.0f));
    yLinks_.assign(width_ * height_, Link(1.0f, 0.0f));
}
 
void U1GaugeField2D::clear() noexcept
{
    width_ = height_ = 0;
    xLinks_.clear();
    yLinks_.clear();
}
 
void U1GaugeField2D::fillIdentity() noexcept
{
    const Link one(1.0f, 0.0f);
    for (std::size_t i = 0; i < xLinks_.size(); ++i) xLinks_[i] = one;
    for (std::size_t i = 0; i < yLinks_.size(); ++i) yLinks_[i] = one;
}
 
bool U1GaugeField2D::empty() const noexcept
{
    return width_ == 0 || height_ == 0 || xLinks_.empty() || yLinks_.empty();
}
 
bool U1GaugeField2D::matchesGrid(const ComplexField2D &field) const noexcept
{
    return !empty() && width_ == field.width() && height_ == field.height() &&
           dx_ == field.dx() && dy_ == field.dy();
}
 
U1GaugeField2D::Link U1GaugeField2D::linkX(std::size_t x, std::size_t y) const noexcept
{
    return xLinks_[index(x, y)];
}
U1GaugeField2D::Link U1GaugeField2D::linkY(std::size_t x, std::size_t y) const noexcept
{
    return yLinks_[index(x, y)];
}
void U1GaugeField2D::setLinkX(std::size_t x, std::size_t y, Link value) noexcept
{
    xLinks_[index(x, y)] = value;
}
void U1GaugeField2D::setLinkY(std::size_t x, std::size_t y, Link value) noexcept
{
    yLinks_[index(x, y)] = value;
}
 
} // namespace qpe