parabolic_module.h

Go to the documentation of this file.

//
// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
// Copyright (C) 2026 by the ryujin authors
//
 
#pragma once
 
#include <ryujin/base/compile_time_options.h>
#include <ryujin/interface/module.h>
#include <ryujin/linear_algebra/state_vector.h>
 
#include <limits>
#include <string>
 
namespace ryujin
{
  namespace Interface
  {
    template <int dim, typename Number = double>
    class ParabolicModule : public Module<dim, Number>
    {
    public:
 
      using typename Module<dim, Number>::StageVectors;
 
      using typename Module<dim, Number>::StageWeights;
 
 
 
      virtual Number backward_euler_step(
          const StateVector<Number> &old_state_vector,
          StageVectors stage_state_vectors,
          StageWeights stage_weights,
          StateVector<Number> &new_state_vector,
          Number t,
          Number tau,
          Number tau_max = std::numeric_limits<Number>::max()) const = 0;
 
      virtual Number crank_nicolson_step(
          const StateVector<Number> &old_state_vector,
          StateVector<Number> &new_state_vector,
          Number t,
          Number tau,
          Number tau_max = std::numeric_limits<Number>::max()) const = 0;
 
 
    protected:
      using Module<dim, Number>::Module;
    };
  } // namespace Interface
} // namespace ryujin