SolverSwitcher.hpp

Go to the documentation of this file.

#ifndef QUANTUMDYNAMICS_SOLVERS_SOLVERSWITCHER_HPP
#define QUANTUMDYNAMICS_SOLVERS_SOLVERSWITCHER_HPP

#include "QuantumDynamicSolver.hpp"
#include "RungeKutta4.hpp"
#include "AdiabaticSolver.hpp"
#include "../../Misc/StrToInt.hpp"

#include <string>
#include <cassert>

namespace solid
{



template <template <typename> class T1, typename T2, typename T3>
class SolverSwitcher
{
public:
    static IQuantumDynamicSolver<T1, T2, T3> *Switch(const std::string label)
    {
        // default solver label
        const std::string defaultSolver = RK4<T1, T2, T3>::label;

        switch (StrToInt(label.c_str()))
        {

        // Runge Kutta 4 order Solver
        case StrToInt(RK4<T1, T2, T3>::label):
            if constexpr (std::is_same<T3, arma::cx_double>::value)
                return new RK4<T1, T2, T3>();
            else
                assert(!"selected solver requires complex type! (T3 = arma::cx_double)");
            break;

        // Adiabatic Solver
        case StrToInt(AdiabaticSolver<T1, T2, T3>::label):
            return new AdiabaticSolver<T1, T2, T3>();

        // default Solver
        default:
            // TODO move to info
            std::cout << "Unknown solver: " << label << ". Running with default solver " << defaultSolver << std::endl;
            return Switch(defaultSolver);
            break;
        }
    }
};

} // namespace solid

#endif