potentialcontours2d_tests.cpp

Go to the documentation of this file.

/**
 * @file potentialcontours2d_tests.cpp
 * @brief Implements regression and invariant tests for potentialcontours2d tests.
 * @ingroup qpe_tests
 *
 * Source: @c tests/potentialcontours2d_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/potentialcontours2d_tests.h"
#include "core/quantum/potentialcontours2d.h"
#include "core/quantum/potentialfield2d.h"
#include <cmath>
 
namespace qpe { namespace tests {
 
/**
 * @brief Runs potential contours2 d tests.
 * @param error Output diagnostic text on failure.
 */
bool runPotentialContours2DTests(std::string &error)
{
    error.clear();
    PotentialField2D binary(32, 32, 0.1f, 0.1f);
    binary.clear();
    binary.addRectangle(0.0f, 0.0f, 1.0f, 0.8f, 8.0f);
    const PotentialContours2DResult c1 = PotentialContours2D::compute(binary, 6u);
    if (!c1.binaryField || c1.levels.size() != 1u || c1.segments.empty()) { error = "Binary contour extraction failed"; return false; }
    if (std::fabs(c1.levels[0] - 4.0f) > 1.0e-5f) { error = "Binary midpoint level is wrong"; return false; }
 
    PotentialField2D ramp(8, 8, 0.25f, 0.25f);
    for (std::size_t y = 0; y < ramp.height(); ++y) {
        for (std::size_t x = 0; x < ramp.width(); ++x) {
            ramp.set(x, y, static_cast<float>(x));
        }
    }
    const PotentialContours2DResult c2 = PotentialContours2D::compute(ramp, 3u);
    if (c2.binaryField || c2.levels.size() != 3u || c2.segments.empty()) { error = "Continuous contour extraction failed"; return false; }
    return true;
}
 
} }