MixedBC Class Reference

Mimetic Mixed Boundary Condition operator. More...

#include <mixedbc.h>

Inheritance diagram for MixedBC:

Public Member Functions
	MixedBC (u16 k, u32 m, Real dx, const std::string &left, const std::vector< Real > &coeffs_left, const std::string &right, const std::vector< Real > &coeffs_right)
	1-D Constructor
	MixedBC (u16 k, u32 m, Real dx, u32 n, Real dy, const std::string &left, const std::vector< Real > &coeffs_left, const std::string &right, const std::vector< Real > &coeffs_right, const std::string &bottom, const std::vector< Real > &coeffs_bottom, const std::string &top, const std::vector< Real > &coeffs_top)
	2-D Constructor
	MixedBC (u16 k, u32 m, Real dx, u32 n, Real dy, u32 o, Real dz, const std::string &left, const std::vector< Real > &coeffs_left, const std::string &right, const std::vector< Real > &coeffs_right, const std::string &bottom, const std::vector< Real > &coeffs_bottom, const std::string &top, const std::vector< Real > &coeffs_top, const std::string &front, const std::vector< Real > &coeffs_front, const std::string &back, const std::vector< Real > &coeffs_back)
	3-D Constructor

Detailed Description

Mimetic Mixed Boundary Condition operator.

Definition at line 27 of file mixedbc.h.

Constructor & Destructor Documentation

MixedBC() [1/3]

MixedBC::MixedBC	(	u16	k,
		u32	m,
		Real	dx,
		const std::string &	left,
		const std::vector< Real > &	coeffs_left,
		const std::string &	right,
		const std::vector< Real > &	coeffs_right )

1-D Constructor

Parameters

k	Order of accuracy
m	Number of cells
dx	Spacing between cells
left	Type of boundary condition at the left boundary ('Dirichlet', 'Neumann', 'Robin')
coeffs_left	Coefficients for the left boundary condition
right	Type of boundary condition at the right boundary ('Dirichlet', 'Neumann', 'Robin')
coeffs_right	Coefficients for the right boundary condition

Definition at line 18 of file mixedbc.cpp.

                                                      {
  sp_mat A(m + 2, m + 2);
  sp_mat BG(m + 2, m + 2);
 
  Gradient *grad = nullptr;
 
  // Handle the left boundary condition
  if (left == "Dirichlet") {
    A.at(0, 0) = coeffs_left[0];
  } else if (left == "Neumann") {
    grad = new Gradient(k, m, dx);
    BG.row(0) = -coeffs_left[0] * grad->row(0);
  } else if (left == "Robin") {
    A.at(0, 0) = coeffs_left[0];
    grad = new Gradient(k, m, dx);
    BG.row(0) = -coeffs_left[1] * grad->row(0);
  } else {
    throw std::invalid_argument("Unknown boundary condition type");
  }
 
  // Handle the right boundary condition
  if (right == "Dirichlet") {
    A.at(m + 1, m + 1) = coeffs_right[0];
  } else if (right == "Neumann") {
    if (!grad)
      grad = new Gradient(k, m, dx);
    BG.row(m + 1) = coeffs_right[0] * grad->row(m);
  } else if (right == "Robin") {
    A.at(m + 1, m + 1) = coeffs_right[0];
    if (!grad)
      grad = new Gradient(k, m, dx);
    BG.row(m + 1) = coeffs_right[1] * grad->row(m);
  } else {
    throw std::invalid_argument("Unknown boundary condition type");
  }
 
  *this = A + BG;
 
  delete grad;
}

MixedBC() [2/3]

MixedBC::MixedBC	(	u16	k,
		u32	m,
		Real	dx,
		u32	n,
		Real	dy,
		const std::string &	left,
		const std::vector< Real > &	coeffs_left,
		const std::string &	right,
		const std::vector< Real > &	coeffs_right,
		const std::string &	bottom,
		const std::vector< Real > &	coeffs_bottom,
		const std::string &	top,
		const std::vector< Real > &	coeffs_top )

2-D Constructor

Parameters

k	Order of accuracy
m	Number of cells along x-axis
dx	Spacing between cells along x-axis
n	Number of cells along y-axis
dy	Spacing between cells along y-axis
left	Type of boundary condition at the left boundary ('Dirichlet', 'Neumann', 'Robin')
coeffs_left	Coefficients for the left boundary condition
right	Type of boundary condition at the right boundary ('Dirichlet', 'Neumann', 'Robin')
coeffs_right	Coefficients for the right boundary condition
bottom	Type of boundary condition at the bottom boundary ('Dirichlet', 'Neumann', 'Robin')
coeffs_bottom	Coefficients for the bottom boundary condition
top	Type of boundary condition at the top boundary ('Dirichlet', 'Neumann', 'Robin')
coeffs_top	Coefficients for the top boundary condition

Definition at line 62 of file mixedbc.cpp.

                                                    {
  MixedBC Bm(k, m, dx, left, coeffs_left, right, coeffs_right);
  MixedBC Bn(k, n, dy, bottom, coeffs_bottom, top, coeffs_top);
 
  sp_mat Im = speye(m + 2, m + 2);
  sp_mat In = speye(n + 2, n + 2);
 
  In.at(0, 0) = 0;
  In.at(n + 1, n + 1) = 0;
 
  sp_mat BC1 = Utils::spkron(In, Bm);
  sp_mat BC2 = Utils::spkron(Bn, Im);
 
  *this = BC1 + BC2;
}

MixedBC() [3/3]

MixedBC::MixedBC	(	u16	k,
		u32	m,
		Real	dx,
		u32	n,
		Real	dy,
		u32	o,
		Real	dz,
		const std::string &	left,
		const std::vector< Real > &	coeffs_left,
		const std::string &	right,
		const std::vector< Real > &	coeffs_right,
		const std::string &	bottom,
		const std::vector< Real > &	coeffs_bottom,
		const std::string &	top,
		const std::vector< Real > &	coeffs_top,
		const std::string &	front,
		const std::vector< Real > &	coeffs_front,
		const std::string &	back,
		const std::vector< Real > &	coeffs_back )

3-D Constructor

Parameters

k	Order of accuracy
m	Number of cells along x-axis
dx	Spacing between cells along x-axis
n	Number of cells along y-axis
dy	Spacing between cells along y-axis
o	Number of cells along z-axis
dz	Spacing between cells along z-axis
left	Type of boundary condition at the left boundary ('Dirichlet', 'Neumann', 'Robin')
coeffs_left	Coefficients for the left boundary condition
right	Type of boundary condition at the right boundary ('Dirichlet', 'Neumann', 'Robin')
coeffs_right	Coefficients for the right boundary condition
bottom	Type of boundary condition at the bottom boundary ('Dirichlet', 'Neumann', 'Robin')
coeffs_bottom	Coefficients for the bottom boundary condition
top	Type of boundary condition at the top boundary ('Dirichlet', 'Neumann', 'Robin')
coeffs_top	Coefficients for the top boundary condition
front	Type of boundary condition at the front boundary ('Dirichlet', 'Neumann', 'Robin')
coeffs_front	Coefficients for the front boundary condition
back	Type of boundary condition at the back boundary ('Dirichlet', 'Neumann', 'Robin')
coeffs_back	Coefficients for the back boundary condition

Definition at line 84 of file mixedbc.cpp.

                                                     {
  MixedBC Bm(k, m, dx, left, coeffs_left, right, coeffs_right);
  MixedBC Bn(k, n, dy, bottom, coeffs_bottom, top, coeffs_top);
  MixedBC Bo(k, o, dz, front, coeffs_front, back, coeffs_back);
 
  sp_mat Im = speye(m + 2, m + 2);
  sp_mat In = speye(n + 2, n + 2);
  sp_mat Io = speye(o + 2, o + 2);
 
  Io.at(0, 0) = 0;
  Io.at(o + 1, o + 1) = 0;
 
  sp_mat In2 = In;
  In2.at(0, 0) = 0;
  In2.at(n + 1, n + 1) = 0;
 
  sp_mat BC1 = Utils::spkron(Utils::spkron(Io, In2), Bm);
  sp_mat BC2 = Utils::spkron(Utils::spkron(Io, Bn), Im);
  sp_mat BC3 = Utils::spkron(Utils::spkron(Bo, In), Im);
 
  *this = BC1 + BC2 + BC3;
}

---

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

• src/cpp/mixedbc.h
• src/cpp/mixedbc.cpp