probabilitycontours2d_tests.cpp

Go to the documentation of this file.

/**
 * @file probabilitycontours2d_tests.cpp
 * @brief Implements regression and invariant tests for probabilitycontours2d tests.
 * @ingroup qpe_tests
 *
 * Source: @c tests/probabilitycontours2d_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 "tests/probabilitycontours2d_tests.h"
 
#include "core/math/complexfield2d.h"
#include "core/quantum/probabilitycontours2d.h"
 
#include <cmath>
#include <sstream>
 
namespace qpe {
namespace tests {
namespace {
 
/**
 * @brief Performs fail.
 * @param errorMessage The error message value used by the operation.
 * @param text The text value used by the operation.
 */
bool fail(std::string &errorMessage, const char *text)
{
    errorMessage = text;
    return false;
}
 
} // namespace
 
/**
 * @brief Runs probability contours2 d tests.
 * @param errorMessage The error message value used by the operation.
 */
bool runProbabilityContours2DTests(std::string &errorMessage)
{
    errorMessage.clear();
 
    ComplexField2D field(65u, 65u, 0.1f, 0.1f);
    const double centre = 32.0;
    for (std::size_t y = 0u; y < field.height(); ++y) {
        for (std::size_t x = 0u; x < field.width(); ++x) {
            const double dx = static_cast<double>(x) - centre;
            const double dy = static_cast<double>(y) - centre;
            const double p = std::exp(-(dx * dx + dy * dy) / 160.0);
            field(x, y) = ComplexField2D::Value(static_cast<float>(std::sqrt(p)), 0.0f);
        }
    }
 
    const ProbabilityContours2DResult linear = ProbabilityContours2D::compute(
        field, 6u, ProbabilityContourSpacing2D::Linear, 4.0f, 384u);
    if (linear.levels.size() != 6u || linear.segments.empty()) {
        return fail(errorMessage, "ProbabilityContours2D: linear Gaussian contours missing");
    }
    for (std::size_t i = 1u; i < linear.levels.size(); ++i) {
        if (!(linear.levels[i] > linear.levels[i - 1u])) {
            return fail(errorMessage, "ProbabilityContours2D: linear levels are not increasing");
        }
    }
 
    const ProbabilityContours2DResult logarithmic = ProbabilityContours2D::compute(
        field, 8u, ProbabilityContourSpacing2D::Logarithmic, 5.0f, 384u);
    if (logarithmic.levels.size() != 8u || logarithmic.segments.empty()) {
        return fail(errorMessage, "ProbabilityContours2D: logarithmic Gaussian contours missing");
    }
    if (!(logarithmic.levels.front() < linear.levels.front())) {
        return fail(errorMessage, "ProbabilityContours2D: logarithmic spacing did not reach faint structure");
    }
 
    ComplexField2D large(1024u, 256u, 0.1f, 0.1f);
    large.fill(ComplexField2D::Value(1.0f, 0.0f));
    const ProbabilityContours2DResult sampled = ProbabilityContours2D::compute(
        large, 4u, ProbabilityContourSpacing2D::Linear, 4.0f, 384u);
    if (sampled.samplingStride < 2u) {
        return fail(errorMessage, "ProbabilityContours2D: automatic live-view downsampling not applied");
    }
    if (!sampled.segments.empty()) {
        return fail(errorMessage, "ProbabilityContours2D: constant field should not produce contours");
    }
 
    return true;
}
 
} // namespace tests
} // namespace qpe