StaticVectorSolver::Gamma< dim, stage > Class Template Reference

Implements the coefficient \(\gamma\) on the left-hand side of the Robin boundary condition (iii) of the static vector boundary value problem. More...

#include <static_vector_input.hpp>

Inheritance diagram for StaticVectorSolver::Gamma< dim, stage >:

Collaboration diagram for StaticVectorSolver::Gamma< dim, stage >:

Public Member Functions
void	value_list (const std::vector< Point< dim >> &r, const std::vector< Tensor< 1, dim >> &n, types::boundary_id bid, types::material_id mid, unsigned int cuid, unsigned int fuid, std::vector< double > &values) const
	Computes values of the parameter \(\gamma\) of the Robin boundary condition at quadrature points. More...
void	value_list (const std::vector< Point< 3 >> &r, const std::vector< Tensor< 1, 3 >> &n, types::boundary_id bid, types::material_id mid, unsigned int cuid, unsigned int fuid, std::vector< double > &values) const
void	value_list (const std::vector< Point< 3 >> &r, const std::vector< Tensor< 1, 3 >> &n, types::boundary_id bid, types::material_id mid, unsigned int cuid, unsigned int fuid, std::vector< double > &values) const
void	value_list (const std::vector< Point< 3 >> &r, const std::vector< Tensor< 1, 3 >> &n, types::boundary_id bid, types::material_id mid, unsigned int cuid, unsigned int fuid, std::vector< double > &values) const
void	value_list (const std::vector< Point< 2 >> &r, const std::vector< Tensor< 1, 2 >> &n, types::boundary_id bid, types::material_id mid, unsigned int cuid, unsigned int fuid, std::vector< double > &values) const
void	value_list (const std::vector< Point< 3 >> &r, const std::vector< Tensor< 1, 3 >> &n, types::boundary_id bid, types::material_id mid, unsigned int cuid, unsigned int fuid, std::vector< double > &values) const
void	value_list (const std::vector< Point< 3 >> &r, const std::vector< Tensor< 1, 3 >> &n, types::boundary_id bid, types::material_id mid, unsigned int cuid, unsigned int fuid, std::vector< double > &values) const
void	value_list (const std::vector< Point< 3 >> &r, const std::vector< Tensor< 1, 3 >> &n, types::boundary_id bid, types::material_id mid, unsigned int cuid, unsigned int fuid, std::vector< double > &values) const
void	value_list (const std::vector< Point< 3 >> &r, const std::vector< Tensor< 1, 3 >> &n, types::boundary_id bid, types::material_id mid, unsigned int cuid, unsigned int fuid, std::vector< double > &values) const

Detailed Description

template<int dim, int stage = 1>
class StaticVectorSolver::Gamma< dim, stage >

Implements the coefficient \(\gamma\) on the left-hand side of the Robin boundary condition (iii) of the static vector boundary value problem.

This class is declared in shared/include/static_vector_input.hpp but must be implemented in xyz/src/static_vector_input.cpp, where xyz is the directory of the current numerical experiment. That is, the declaration is shared between all numerical experiments while implementation is specific for each individual numerical experiment. See the structure of the code for more details.

The dim template parameter is, as per usual, the amount of spatial dimensions. The purpose of the stage template parameter is discussed in here.

The user is supposed to implement

void value_list(...)

member function of this class.

Definition at line 220 of file static_vector_input.hpp.

Member Function Documentation

value_list()

template<int dim, int stage = 1>

void StaticVectorSolver::Gamma< dim, stage >::value_list	(	const std::vector< Point< dim >> &	r,
		const std::vector< Tensor< 1, dim >> &	n,
		types::boundary_id	bid,
		types::material_id	mid,
		unsigned int	cuid,
		unsigned int	fuid,
		std::vector< double > &	values
	)		const

Computes values of the parameter \(\gamma\) of the Robin boundary condition at quadrature points.

This function is called by the StaticVectorSolver::Solver1 and StaticVectorSolver::Solver2 classes during the assembly of the system matrix. It is called ones per each boundary face if the face ID is an even integer greater than zero, see [boundary ID convention](veis_bnd_convention). This function must fill the vector values with the values of \(\gamma\). The values[i] is interpreted as \(\gamma\) at quadrature point r[i]. The code snippet below provides an example in which \(\gamma\) equals

\[ \gamma = (1 / \mu) ( \sqrt{x^2+y^2}+2 ) \]

,

where \(\mu\) is the permeability.

#pragma GCC diagnostic push
#pragma GCC diagnostic ignored "-Wunused-parameter"
template<>
void Gamma<3>::value_list(
  const std::vector<Point<3>>& r,
  const std::vector<Tensor<1, 3>>& n,
  types::boundary_id bid,
  types::material_id mid,
  unsigned int cuid,
  unsigned int fuid,
  std::vector<double>& values) const
{
  Assert(r.size() == values.size(),
    ExcDimensionMismatch(r.size(), values.size()));
 
  auto v = values.begin();
  for (auto p: r)
  {
    double mu = MU0 * (pow(p[0],2) + pow(p[1],2) + 1);
    *v = (1 / mu) * (sqrt(pow(p[0],2) + pow(p[1],2)) + 2);
    v++;
  }
}
#pragma GCC diagnostic pop

Parameters

[in]	r	- A vector that contains the quadrature points of the cell being processed.
[in]	n	- A vector that contains the normal vectors at quadrature points.
[in]	bid	- The boundary ID.
[in]	mid	- The material ID.
[in]	cuid	- The cell user ID.
[in]	fuid	- The face user ID.
[out]	values	- The output data.

---

The documentation for this class was generated from the following file:

• static/shared/include/static_vector_input.hpp