time_integrators.h

// ---------------------------------------------------------------------
//
// Copyright (C) 2020 by the hyper.deal authors
//
// This file is part of the hyper.deal library.
//
// The hyper.deal library is free software; you can use it, redistribute
// it, and/or modify it under the terms of the GNU Lesser General
// Public License as published by the Free Software Foundation; either
// version 3.0 of the License, or (at your option) any later version.
// The full text of the license can be found in the file LICENSE.MD at
// the top level directory of hyper.deal.
//
// ---------------------------------------------------------------------
 
#ifndef HYPERDEAL_FUNCTIONALITIES_TIME_INGEGRATORS
#define HYPERDEAL_FUNCTIONALITIES_TIME_INGEGRATORS
 
#include <hyper.deal/base/config.h>
 
#include <deal.II/base/config.h>
 
#include <deal.II/base/exceptions.h>
 
#include <deal.II/lac/petsc_block_vector.h>
#include <deal.II/lac/petsc_vector.h>
#include <deal.II/lac/trilinos_parallel_block_vector.h>
#include <deal.II/lac/trilinos_vector.h>
 
#include <functional>
#include <string>
#include <type_traits>
#include <vector>
 
namespace hyperdeal
{
  template <typename Number, typename VectorType>
  class LowStorageRungeKuttaIntegrator
  {
  public:
    LowStorageRungeKuttaIntegrator(VectorType &      vec_Ki,
                                   VectorType &      vec_Ti,
                                   const std::string type,
                                   const bool        only_Ti_is_ghosted = true);
 
    void
    perform_time_step(
      VectorType &  solution,
      const Number &current_time,
      const Number &time_step,
      const std::function<void(const VectorType &, VectorType &, const Number)>
        &op);
 
    unsigned int
    n_stages() const;
 
  private:
    VectorType &vec_Ki;
 
    VectorType &vec_Ti;
 
    const bool only_Ti_is_ghosted;
 
    std::vector<Number> ai, bi;
 
    static constexpr bool manual_compress =
      false
#ifdef DEAL_II_WITH_PETSC
      || std::is_same_v<VectorType, dealii::PETScWrappers::MPI::Vector> ||
      std::is_same_v<VectorType, dealii::PETScWrappers::MPI::BlockVector>
#endif
#ifdef DEAL_II_WITH_TRILINOS
      || std::is_same_v<VectorType, dealii::TrilinosWrappers::MPI::Vector> ||
      std::is_same_v<VectorType, dealii::TrilinosWrappers::MPI::BlockVector>
#endif
      ;
  };
 
 
} // namespace hyperdeal
 
#endif