quantumkapitzadirac2d_tests.cpp

Go to the documentation of this file.

/**
 * @file quantumkapitzadirac2d_tests.cpp
 * @brief Implements regression and invariant tests for quantumkapitzadirac2d tests.
 * @ingroup qpe_tests
 *
 * Source: @c tests/quantumkapitzadirac2d_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 "quantumkapitzadirac2d_tests.h"
 
#include "core/math/complexfield2d.h"
#include "core/quantum/momentumspace2d.h"
#include "core/quantum/potentialfield2d.h"
#include "core/quantum/potentialscene2d.h"
#include "core/quantum/quantumkapitzadirac2d.h"
 
#include <cmath>
 
namespace qpe {
namespace tests {
 
/**
 * @brief Runs quantum kapitza dirac2 d tests.
 * @param error Output diagnostic text on failure.
 */
bool runQuantumKapitzaDirac2DTests(std::string &error)
{
    ComplexField2D field(256u, 256u, 0.16f, 0.16f);
    QuantumKapitzaDiracInitialState2DConfig initial;
    if (!QuantumKapitzaDirac2D::initialize(field, initial, 1.0, nullptr) ||
        std::fabs(field.norm2() - 1.0) > 2.0e-5) {
        error = "Kapitza-Dirac initial beam failed normalization";
        return false;
    }
 
    QuantumKapitzaDirac2DConfig config;
    config.enabled = true;
    config.latticePeriod = 4.096;
    config.latticeDepth = 16.0;
    config.pulseStartTime = 0.0;
    config.pulseDuration = 0.30;
    config.analysisEndTime = 4.0;
    config.maximumOrder = 4;
    if (!QuantumKapitzaDirac2D::validateConfig(config, 1.0, nullptr)) {
        error = "Kapitza-Dirac default config is invalid";
        return false;
    }
    const double expectedG = 2.0 * 3.14159265358979323846 / 4.096;
    if (std::fabs(QuantumKapitzaDirac2D::reciprocalMomentum(config, 1.0) - expectedG) > 1.0e-12 ||
        std::fabs(QuantumKapitzaDirac2D::pulseArea(config, 1.0) - 2.4) > 1.0e-12) {
        error = "Kapitza-Dirac reciprocal momentum or pulse area is inconsistent";
        return false;
    }
 
    PotentialScene2D scene;
    PotentialField2D potential(field.width(), field.height(), field.dx(), field.dy());
    const QuantumKapitzaDiracPotentialState2D during =
        QuantumKapitzaDirac2D::rasterize(scene, config, 0.15, potential);
    if (!during.valid || !during.active || during.finished || !(potential.maximum() > 15.0f)) {
        error = "Kapitza-Dirac standing-wave pulse did not turn on";
        return false;
    }
    const QuantumKapitzaDiracPotentialState2D after =
        QuantumKapitzaDirac2D::rasterize(scene, config, 0.31, potential);
    if (!after.valid || after.active || !after.finished || std::fabs(potential.maximum()) > 1.0e-6f) {
        error = "Kapitza-Dirac standing-wave pulse did not turn off";
        return false;
    }
 
    const MomentumSpace2DResult momentum = MomentumSpace2D::compute(field, 1.0);
    const QuantumKapitzaDiracSample2D sample =
        QuantumKapitzaDirac2D::measure(momentum, config, 0.0, false);
    if (sample.momentumProbability.size() != field.width() ||
        sample.orderProbability.size() != 9u ||
        !(sample.norm > 0.999 && sample.norm < 1.001)) {
        error = "Kapitza-Dirac momentum analysis is inconsistent";
        return false;
    }
    const std::size_t zero = 4u;
    if (!(sample.orderProbability[zero] > 0.98f)) {
        error = "Kapitza-Dirac unpulsed beam is not concentrated in zero order";
        return false;
    }
 
    QuantumKapitzaDiracHistory2D history;
    if (!history.configure(config, 1.0)) {
        error = "Kapitza-Dirac history configuration failed";
        return false;
    }
    history.append(momentum, 0.0, false);
    if (history.empty() || history.samples().back().orderProbability.size() != 9u) {
        error = "Kapitza-Dirac history append failed";
        return false;
    }
 
    return true;
}
 
} // namespace tests
} // namespace qpe