quantumhalledge2d_tests.cpp

Go to the documentation of this file.

/**
 * @file quantumhalledge2d_tests.cpp
 * @brief Implements regression and invariant tests for quantumhalledge2d tests.
 * @ingroup qpe_tests
 *
 * Source: @c tests/quantumhalledge2d_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 "quantumhalledge2d_tests.h"
#include "core/math/complexfield2d.h"
#include "core/quantum/potentialfield2d.h"
#include "core/quantum/quantumhalledge2d.h"
#include "core/quantum/splitstepfouriersolver2d.h"
#include <cmath>
namespace qpe { namespace tests {
namespace {
/**
 * @brief Runs one.
 * @param b The b value used by the operation.
 * @param drift The drift value used by the operation.
 * @param edgeProbability The edge probability value used by the operation.
 * @param error Output diagnostic text on failure.
 */
bool runOne(double b, double &drift, double &edgeProbability, std::string &error)
{
    QuantumHallEdge2DConfig c; c.enabled=true; c.fieldStrength=b; c.edgeY=-8.0; c.wallHeight=120.0;
    c.launchY=-5.3; c.orbitRadius=4.2; c.packetSigma=1.1; c.edgeBand=4.0; c.analysisEndTime=8.0;
    ComplexField2D f(256u,256u,0.16f,0.16f); PotentialField2D v; U1GaugeField2D g;
    if (!QuantumHallEdge2D::buildFields(f,c,v,g,&error) ||
        !QuantumHallEdge2D::prepareState(f,c,1.0,1.0,&error)) return false;
    QuantumHallEdgeHistory2D h; if(!h.configure(f,c)){error="Hall history configuration failed";return false;} h.append(0.0,f);
    SplitStepFourierSolver2D s; SplitStepFourier2DParams p; p.mass=1.0f;p.hbar=1.0f;p.dt=0.02f;
    s.setParameters(p);s.prepare(f);s.preparePotential(v);
    for(int i=0;i<400;++i)s.stepGaugeCovariant(f,v,g);
    h.append(s.time(),f); const auto &a=h.samples().back(); drift=a.driftFromStart; edgeProbability=a.edgeProbability;
    if(std::fabs(a.norm-1.0)>3.0e-3){error="Hall edge propagation lost norm";return false;} return true;
}
}
/**
 * @brief Runs quantum hall edge2 d tests.
 * @param error Output diagnostic text on failure.
 */
bool runQuantumHallEdge2DTests(std::string &error)
{
    double dp=0.0,dm=0.0,ep=0.0,em=0.0;
    if(!runOne(+0.45,dp,ep,error)||!runOne(-0.45,dm,em,error))return false;
    if(dp<9.0||dm>-9.0||ep<0.80||em<0.80){error="Hall edge packet did not show robust chiral edge transport";return false;}
    if(std::fabs(dp+dm)>0.35){error="Hall edge chirality is not symmetric under B reversal";return false;}
    error.clear();return true;
}
} }