algorithms.hpp

Go to the documentation of this file.

/// @file solve/algorithms.hpp
/// @brief Algorithm tags: carry the numerics, not the mathematics.
#pragma once
 
#include "core/types.hpp"
#include "linalg/solvers/linear_solver.hpp"
 
namespace num {
 
struct Euler {
    double h = 1e-3;
};
 
struct RK4 {
    double h = 1e-3;
};
 
struct RK45 {
    double h = 1e-3;
    double rtol = 1e-6;
    double atol = 1e-9;
    idx max_steps = 1000000;
};
 
struct BackwardEuler {
    LinearSolver solver;
    double dt;
    int nstep;
};
 
struct Metropolis {
    int equilibration = 1000;
    int measurements = 500;
};
 
} // namespace num