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 <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 void
      backward_euler_step(const StateVector<Number> &old_state_vector,
                          const Number old_t,
                          StageVectors stage_state_vectors,
                          StageWeights stage_weights,
                          StateVector<Number> &new_state_vector,
                          Number tau) const = 0;
 
      virtual void
      crank_nicolson_step(const StateVector<Number> &old_state_vector,
                          const Number old_t,
                          StateVector<Number> &new_state_vector,
                          Number tau) const = 0;
 
 
    protected:
      using Module<dim, Number>::Module;
    };
  } // namespace Interface
} // namespace ryujin