nonuniform_dirichlet.h

// SPDX-FileCopyrightText: © 2025 PRISMS Center at the University of Michigan
// SPDX-License-Identifier: GNU Lesser General Public Version 2.1
 
#ifndef nonuniform_dirichlet_h
#define nonuniform_dirichlet_h
 
#include <deal.II/base/function.h>
#include <deal.II/base/point.h>
#include <deal.II/lac/vector.h>
 
#include <prismspf/config.h>
#include <prismspf/core/type_enums.h>
#include <prismspf/user_inputs/user_input_parameters.h>
 
PRISMS_PF_BEGIN_NAMESPACE
 
template <int dim>
class customNonuniformDirichlet;
 
template <int dim, fieldType field_type = fieldType::SCALAR>
class nonuniformDirichlet : public dealii::Function<dim, double>
{
public:
  nonuniformDirichlet(const unsigned int             &_index,
                      const unsigned int             &_boundary_id,
                      const userInputParameters<dim> &_user_inputs);
 
  double
  value(const dealii::Point<dim> &p, const unsigned int component = 0) const override;
 
  void
  vector_value(const dealii::Point<dim> &p, dealii::Vector<double> &value) const override;
 
private:
  const unsigned int index;
 
  const unsigned int boundary_id;
 
  const userInputParameters<dim> *user_inputs;
 
  customNonuniformDirichlet<dim> custom_nonuniform_dirichlet;
};
class nonuniformDirichlet : public dealii::Function<dim, double> {…};
 
template <int dim, fieldType field_type>
nonuniformDirichlet<dim, field_type>::nonuniformDirichlet(
  const unsigned int             &_index,
  const unsigned int             &_boundary_id,
  const userInputParameters<dim> &_user_inputs)
  : dealii::Function<dim>((field_type == fieldType::VECTOR) ? dim : 1)
  , index(_index)
  , boundary_id(_boundary_id)
  , user_inputs(&_user_inputs)
{}
nonuniformDirichlet<dim, field_type>::nonuniformDirichlet( {…}
 
template <int dim, fieldType field_type>
inline double
nonuniformDirichlet<dim, field_type>::value(
  const dealii::Point<dim>           &p,
  [[maybe_unused]] const unsigned int component) const
{
  // Initialize passed variables to zero
  double                 scalar_value = 0.0;
  dealii::Vector<double> vector_value(dim);
 
  // Pass variables to user-facing function to evaluate
  custom_nonuniform_dirichlet.set_nonuniform_dirichlet(index,
                                                       boundary_id,
                                                       0,
                                                       p,
                                                       scalar_value,
                                                       vector_value(0),
                                                       *user_inputs);
 
  return scalar_value;
}
nonuniformDirichlet<dim, field_type>::value( {…}
 
template <int dim, fieldType field_type>
inline void
nonuniformDirichlet<dim, field_type>::vector_value(const dealii::Point<dim> &p,
                                                   dealii::Vector<double>   &value) const
{
  // Initialize passed variables to zero
  double                 scalar_value = 0.0;
  dealii::Vector<double> vector_value(dim);
 
  // Pass variables to user-facing function to evaluate
  for (unsigned int i = 0; i < dim; i++)
    {
      custom_nonuniform_dirichlet.set_nonuniform_dirichlet(index,
                                                           boundary_id,
                                                           i,
                                                           p,
                                                           scalar_value,
                                                           vector_value(i),
                                                           *user_inputs);
    }
 
  value = vector_value;
}
nonuniformDirichlet<dim, field_type>::vector_value(const dealii::Point<dim> &p, {…}
 
template <int dim>
class customNonuniformDirichlet
{
public:
  customNonuniformDirichlet() = default;
 
  void
  set_nonuniform_dirichlet(const unsigned int             &index,
                           const unsigned int             &boundary_id,
                           const unsigned int             &component,
                           const dealii::Point<dim>       &point,
                           double                         &scalar_value,
                           double                         &vector_component_value,
                           const userInputParameters<dim> &user_inputs) const;
};
class customNonuniformDirichlet {…};
 
PRISMS_PF_END_NAMESPACE
 
#endif