wavefunctioncrosssection2d_tests.cpp

Go to the documentation of this file.

/**
 * @file wavefunctioncrosssection2d_tests.cpp
 * @brief Implements regression and invariant tests for wavefunctioncrosssection2d tests.
 * @ingroup qpe_tests
 *
 * Source: @c tests/wavefunctioncrosssection2d_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 "wavefunctioncrosssection2d_tests.h"
 
#include "core/math/complexfield2d.h"
#include "core/quantum/wavefunctioncrosssection2d.h"
 
#include <cmath>
 
namespace qpe {
namespace tests {
 
/**
 * @brief Runs wave function cross section2 d tests.
 * @param error Output diagnostic text on failure.
 */
bool runWaveFunctionCrossSection2DTests(std::string &error)
{
    error.clear();
    ComplexField2D field(5u, 5u, 1.0f, 1.0f);
    for (std::size_t y = 0u; y < field.height(); ++y) {
        for (std::size_t x = 0u; x < field.width(); ++x) {
            const float re = static_cast<float>(10u * y + x);
            const float im = static_cast<float>(x) * 0.5f;
            field(x, y) = ComplexField2D::Value(re, im);
        }
    }
 
    const WaveFunctionCrossSection2DData horizontal =
        WaveFunctionCrossSection2D::sample(
            field, CrossSectionAxis::HorizontalY, 0.2);
    if (horizontal.coordinate.size() != 5u ||
        std::fabs(horizontal.actualPosition) > 1.0e-9) {
        error = "horizontal section size or snapped position is wrong";
        return false;
    }
    if (std::fabs(horizontal.realPart[3] - 23.0) > 1.0e-9) {
        error = "horizontal section sampled the wrong row/value";
        return false;
    }
    const double expectedP = 23.0 * 23.0 + 1.5 * 1.5;
    if (std::fabs(horizontal.probability[3] - expectedP) > 1.0e-9) {
        error = "probability value is inconsistent with complex sample";
        return false;
    }
 
    const WaveFunctionCrossSection2DData vertical =
        WaveFunctionCrossSection2D::sample(
            field, CrossSectionAxis::VerticalX, -1.1);
    if (vertical.coordinate.size() != 5u ||
        std::fabs(vertical.actualPosition + 1.0) > 1.0e-9) {
        error = "vertical section size or snapped position is wrong";
        return false;
    }
    if (std::fabs(vertical.realPart[4] - 41.0) > 1.0e-9) {
        error = "vertical section sampled the wrong column/value";
        return false;
    }
 
    return true;
}
 
} // namespace tests
} // namespace qpe