scatteringobservables2d_tests.cpp

Go to the documentation of this file.

/**
 * @file scatteringobservables2d_tests.cpp
 * @brief Implements regression and invariant tests for scatteringobservables2d tests.
 * @ingroup qpe_tests
 *
 * Source: @c tests/scatteringobservables2d_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 "scatteringobservables2d_tests.h"
 
#include "core/math/complexfield2d.h"
#include "core/quantum/scatteringobservables2d.h"
 
#include <cmath>
#include <sstream>
 
namespace qpe { namespace tests {
namespace {
 
/**
 * @brief Performs near.
 * @param a The a value used by the operation.
 * @param b The b value used by the operation.
 * @param tolerance The tolerance value used by the operation.
 */
bool near(double a, double b, double tolerance)
{
    return std::abs(a - b) <= tolerance;
}
 
} // namespace
 
/**
 * @brief Runs scattering observables2 d tests.
 * @param error Output diagnostic text on failure.
 */
bool runScatteringObservables2DTests(std::string &error)
{
    error.clear();
 
    // 6 x 2 grid, dx=dy=1. x coordinates are
    // -2.5, -1.5, -0.5, +0.5, +1.5, +2.5.
    ComplexField2D field(6, 2, 1.0f, 1.0f);
    field.fill(ComplexField2D::Value(1.0f, 0.0f));
 
    const ScatteringRegions2D m = ScatteringObservables2D::measureX(field, -1.0, 1.0);
 
    // Two x columns per region, two y rows: 4 samples in every region.
    if (!near(m.total, 12.0, 1.0e-12) ||
        !near(m.left, 4.0, 1.0e-12) ||
        !near(m.interaction, 4.0, 1.0e-12) ||
        !near(m.right, 4.0, 1.0e-12) ||
        !near(m.leftFraction, 1.0 / 3.0, 1.0e-12) ||
        !near(m.interactionFraction, 1.0 / 3.0, 1.0e-12) ||
        !near(m.rightFraction, 1.0 / 3.0, 1.0e-12)) {
        std::ostringstream ss;
        ss << "region integration mismatch: total=" << m.total
           << " left=" << m.left
           << " interaction=" << m.interaction
           << " right=" << m.right;
        error = ss.str();
        return false;
    }
 
    // Boundary order is intentionally reversible for convenience.
    const ScatteringRegions2D reversed = ScatteringObservables2D::measureX(field, 1.0, -1.0);
    if (!near(reversed.left, m.left, 1.0e-12) ||
        !near(reversed.interaction, m.interaction, 1.0e-12) ||
        !near(reversed.right, m.right, 1.0e-12)) {
        error = "reversed boundaries changed the result";
        return false;
    }
 
    return true;
}
 
} } // namespace qpe::tests