dae_time_integrator.h

Go to the documentation of this file.

#ifndef dae_time_integrator_h
#define dae_time_integrator_h

#include <deal.II/base/logstream.h>
#include <deal.II/base/exceptions.h>
#include <deal.II/base/parameter_handler.h>

// For time integration.
#include <ida/ida.h>
#include <ida/ida_dense.h>
#include <ida/ida_lapack.h>
#include <ida/ida_spgmr.h>
#include <nvector/nvector_serial.h>
#include <sundials/sundials_math.h>
#include <sundials/sundials_types.h>


#include "ode_argument.h"


class DAETimeIntegrator
{
public:
  DAETimeIntegrator(OdeArgument &solver);

  DAETimeIntegrator(OdeArgument &solver,
                    double initial_step_size,
                    double min_step_size,
                    double initial_time,
                    double final_time,
                    double abs_tol,
                    double rel_tol,
                    unsigned int max_n_steps,
                    double outputs_period,
                    unsigned int ic_type,
                    bool use_iterative,
                    bool provide_jac,
                    bool provide_jac_prec);

  ~DAETimeIntegrator();

  static void declare_parameters(ParameterHandler &prm);

  void parse_parameters(ParameterHandler &prm);

  unsigned int start_ode(Vector<double> &solution,
                         Vector<double> &solution_dot,
                         const unsigned int max_steps);

  void reset_ode(Vector<double> &y, Vector<double> &yp,
                 double t, double h, unsigned int max_steps);


  double final_time;

  double initial_time;

private:
  OdeArgument &solver;

  double initial_step_size;

  double min_step_size;

  double abs_tol;

  double rel_tol;

  unsigned int max_n_steps;

  double outputs_period;

  unsigned int ic_type;

  bool is_initialized;

  bool use_iterative;

  bool provide_jac;
  bool provide_jac_prec;


  void *ida_mem;

  N_Vector yy;
  N_Vector yp;
  N_Vector abs_tolls;
  N_Vector diff_id;


};

#endif