QuantumDynamics.hpp

Go to the documentation of this file.

#ifndef QUANTUMDYNAMICS_QUANTUMDYNAMICS_HPP
#define QUANTUMDYNAMICS_QUANTUMDYNAMICS_HPP

#include <armadillo>
#include <string>

#include "../QuantumSystem/QuantumSystem.hpp"
#include "../QuantumState/QuantumState.hpp"
#include "../Parameters/Parameters.hpp"
#include "../Hamiltonian/TermsEnabled.hpp"
#include "DynamicsSchedule.hpp"
#include "Solvers/QuantumDynamicSolver.hpp"
#include "Solvers/SolverSwitcher.hpp"

namespace solid
{

template <template <typename> class T1, typename T2, typename T3>
class QuantumDynamics
{

public:
    void Create(
        QuantumSystem<T1, T2> &qSystem,
        QuantumState<T3> &initQuantumState,
        DynamicsSchedule<arma::SpMat<T2>> &dynSchedule,
        MeasurementSchedule<T1, T2, T3> &mesSchedule,
        const std::string label = "rk4")
    {
        terms = qSystem.hamiltonian.termsEnabled;
        quantumSystem = qSystem;
        quantumState = initQuantumState;
        param = qSystem.parameters;
        dynamicsSchedule = dynSchedule;
        measurementSchedule = mesSchedule;
        solver = SolverSwitcher<T1, T2, T3>::Switch(label);
    }

    void Run();

private:
    void LoadParameters();

    void Propagate();

    void Measure();

    // TODO reorganize friend class, add getter setter methods

    template <template <typename> class G1, typename G2, typename G3>
    friend class RK4;

    template <template <typename> class G1, typename G2, typename G3>
    friend class AdiabaticSolver;

    template <template <typename> class G1, typename G2, typename G3>
    friend class IQuantumDynamicSolver;

    double time;                                         
    TermsEnabled terms;                                  
    QuantumState<T3> quantumState;                       
    QuantumSystem<T1, T2> quantumSystem;                 
    Parameters<T2> param;                                
    DynamicsSchedule<arma::SpMat<T2>> dynamicsSchedule;  
    MeasurementSchedule<T1, T2, T3> measurementSchedule; 
    IQuantumDynamicSolver<T1, T2, T3> *solver;           
};

} // namespace solid
#endif