initial_conditions.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: © 2025 PRISMS Center at the University of Michigan
// SPDX-License-Identifier: GNU Lesser General Public Version 2.1
 
#pragma once
 
#include <deal.II/base/function.h>
#include <deal.II/base/point.h>
#include <deal.II/lac/vector.h>
 
#include <prismspf/core/type_enums.h>
 
#include <prismspf/config.h>
 
PRISMS_PF_BEGIN_NAMESPACE
 
template <unsigned int dim>
class UserInputParameters;
 
struct InitialConditionFile;
 
template <unsigned int dim>
struct SpatialDiscretization;
 
template <unsigned int dim, unsigned int degree, typename number>
class PDEOperatorBase;
 
template <unsigned int dim, typename number>
class ReadFieldBase;
 
template <unsigned int dim, unsigned int degree, typename number>
class InitialCondition : public dealii::Function<dim, number>
{
public:
  InitialCondition(const unsigned int                         &_index,
                   const TensorRank                           &_field_type,
                   const PDEOperatorBase<dim, degree, number> &_pde_operator);
 
  // NOLINTBEGIN(readability-identifier-length)
 
  void
  vector_value(const dealii::Point<dim> &p, dealii::Vector<number> &value) const override;
 
  // NOLINTEND(readability-identifier-length)
 
private:
  unsigned int index;
 
  TensorRank field_type;
 
  const PDEOperatorBase<dim, degree, number> *pde_operator;
};
 
template <unsigned int dim, typename number>
class ReadInitialCondition : public dealii::Function<dim, number>
{
public:
  ReadInitialCondition(std::string                       _field_name,
                       const TensorRank                 &_field_type,
                       const InitialConditionFile       &ic_file,
                       const SpatialDiscretization<dim> &spatial_discretization);
 
  // NOLINTBEGIN(readability-identifier-length)
 
  void
  vector_value(const dealii::Point<dim> &p, dealii::Vector<number> &value) const override;
 
  // NOLINTEND(readability-identifier-length)
 
private:
  std::string field_name;
 
  TensorRank field_type;
 
  std::shared_ptr<ReadFieldBase<dim, number>> reader;
};
 
PRISMS_PF_END_NAMESPACE