solver_def.hpp

Go to the documentation of this file.

 
#ifndef ROOT_SOLVER_DEF_HPP
#define ROOT_SOLVER_DEF_HPP
 
#include <Eigen/Dense>
#include <Eigen/Sparse>
#include <functional>
#include <string>
 
#include "method.hpp"
#include "stepper_def.hpp"
 
template <typename T>
class StepperBase;
template <typename T>
class Solver {
  private:
    friend class SolverTester;  
    Method method;              
    int max_iterations;         
    double tolerance;           
    bool aitken_requirement;    
    bool verbose;               
    std::function<double(double)>
        derivative_or_function_g;  
    Eigen::MatrixX2d results;  
    std::function<double(double)>
        function;     
    T initial_guess;  
    void solver_step(int& iter, std::unique_ptr<StepperBase<T>>& stepper, double& err);
    void save_results(int iter, Eigen::Vector2d result_to_save);
    Eigen::Vector2d get_previous_result(int step_length);
    double calculate_error(double x_prev, double x_next);
    void save_starting_point();
    void convert_stepper(std::unique_ptr<StepperBase<T>>& stepper);
 
  public:
    Solver(std::function<double(double)> fun, T initial_guess, const Method method, int max_iterations,
           double tolerance, bool aitken_mode, bool verbose);
    Solver(std::function<double(double)> fun, T initial_guess, const Method method, int max_iterations,
           double tolerance, bool aitken_mode, bool verbose, std::function<double(double)> derivative_or_function_g);
    Eigen::MatrixX2d solve();
};
 
#endif  // ROOT_SOLVER_DEF_HPP