kinsol.hh

Go to the documentation of this file.

#pragma once
 
// standard library
#include <functional>
 
// external libraries
#include <deal.II/sundials/kinsol.h>
 
namespace DiFfRG
{
  using namespace dealii;
  template <typename VectorType_> class KINSOL
  {
  public:
    using VectorType = VectorType_;
 
    KINSOL(double abstol_ = 1e-10, double reltol_ = 1e-6, double assemble_threshold_ = 0., unsigned int max_steps_ = 21,
           unsigned int n_stepsize_iterations_ = 21);
 
    void reinit(const VectorType &u);
 
    void operator()(VectorType &iterate);
 
    double threshold(const double thr);
 
    const VectorType &get_residual() const;
    const VectorType &get_initial() const;
    const VectorType &get_iterate() const;
 
    double average_time_newton() const;
    unsigned int num_newton_calls() const;
 
    void set_ignore_nonconv(bool x);
 
    std::function<void(VectorType &, const VectorType &)> residual;
    std::function<void(VectorType &, const VectorType &)> lin_solve;
    std::function<void(const VectorType &)> update_jacobian;
 
    double get_error();
 
    unsigned int get_step();
 
    unsigned int get_jacobians();
 
  private:
    VectorType residual_vector;
    VectorType initial_vector;
    VectorType iterate_vector;
    VectorType eest_tmp;
    double abstol, reltol, assemble_threshold;
    unsigned int max_steps, n_stepsize_iterations;
 
    bool converged;
    bool ignore_nonconv;
    double resnorm;
    unsigned int step, jacobians;
 
    std::vector<double> timings_newton;
 
    double get_EEst();
 
    bool check();
 
    std::shared_ptr<SUNDIALS::KINSOL<VectorType>> kinsol;
  };
} // namespace DiFfRG