diracspinor2d_tests.cpp

Go to the documentation of this file.

/**
 * @file diracspinor2d_tests.cpp
 * @brief Implements regression and invariant tests for diracspinor2d tests.
 * @ingroup qpe_tests
 *
 * Source: @c tests/diracspinor2d_tests.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 "diracspinor2d_tests.h"
#include "core/quantum/diracspinor2d.h"
#include <cmath>
namespace qpe { namespace tests {
/**
 * @brief Runs dirac spinor2 d tests.
 * @param error Output diagnostic text on failure.
 */
bool runDiracSpinor2DTests(std::string &error)
{
    error.clear();
    DiracSpinor2D s(16u, 8u, 0.25f, 0.5f);
    s.fill(ComplexField2D::Value(1.0f,0.0f), ComplexField2D::Value(0.5f,0.0f));
    if (!s.normalize() || std::abs(s.norm2()-1.0)>2.0e-6) { error="Dirac spinor normalization failed"; return false; }
    ComplexField2D proxy;
    s.writeDensityProxy(proxy);
    if (!s.matchesGrid(proxy) || std::abs(proxy.norm2()-1.0)>2.0e-6) { error="Dirac density proxy does not preserve density"; return false; }
    if (std::abs(s.upperPopulation()+s.lowerPopulation()-1.0)>2.0e-6) { error="Dirac component populations do not sum to one"; return false; }
    return true;
}
} }