RobinBC Class Reference

Mimetic Robin Boundary Condition operator. More...

#include <robinbc.h>

Inheritance diagram for RobinBC:

Public Member Functions
	RobinBC (u16 k, u32 m, Real dx, Real a, Real b)
	1-D Robin boundary constructor
	RobinBC (u16 k, u32 m, Real dx, u32 n, Real dy, Real a, Real b)
	2-D Robin boundary constructor
	RobinBC (u16 k, u32 m, Real dx, u32 n, Real dy, u32 o, Real dz, Real a, Real b)
	3-D Robin boundary constructor

Detailed Description

Mimetic Robin Boundary Condition operator.

Definition at line 23 of file robinbc.h.

Constructor & Destructor Documentation

RobinBC() [1/3]

RobinBC::RobinBC	(	u16	k,
		u32	m,
		Real	dx,
		Real	a,
		Real	b )

1-D Robin boundary constructor

Parameters

k	mimetic order of accuracy
m	number of cells in x-dimension
dx	cell width in x-direction
a	Coefficient of the Dirichlet function
b	Coefficient of the Neumann function

Definition at line 16 of file robinbc.cpp.

                                                      {
  sp_mat A(m + 2, m + 2);
  sp_mat BG(m + 2, m + 2);
 
  A.at(0, 0) = a;
  A.at(m + 1, m + 1) = a;
 
  Gradient grad(k, m, dx);
 
  BG.row(0) = -b * grad.row(0);
  BG.row(m + 1) = b * grad.row(m);
 
  *this = A + BG;
}

RobinBC() [2/3]

RobinBC::RobinBC	(	u16	k,
		u32	m,
		Real	dx,
		u32	n,
		Real	dy,
		Real	a,
		Real	b )

2-D Robin boundary constructor

Parameters

k	mimetic order of accuracy
m	number of cells in x-dimension
dx	cell width in x-direction
n	number of cells in y-dimension
dy	cell width in y-direction
a	Coefficient of the Dirichlet function
b	Coefficient of the Neumann function

Note

Uses 1-D Robin to build the 2-D operator

Definition at line 32 of file robinbc.cpp.

                                                                      {
  RobinBC Bm(k, m, dx, a, b);
  RobinBC Bn(k, n, dy, a, b);
 
  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;
}

RobinBC() [3/3]

RobinBC::RobinBC	(	u16	k,
		u32	m,
		Real	dx,
		u32	n,
		Real	dy,
		u32	o,
		Real	dz,
		Real	a,
		Real	b )

3-D Robin boundary constructor

Parameters

k	mimetic order of accuracy
m	number of cells in x-dimension
dx	cell width in x-direction
n	number of cells in y-dimension
dy	cell width in y-direction
o	number of cells in z-dimension
dz	cell width in z-direction
a	Coefficient of the Dirichlet function
b	Coefficient of the Neumann function

Note

Uses 1-D Robin to build the 3-D operator

Definition at line 49 of file robinbc.cpp.

                         {
  RobinBC Bm(k, m, dx, a, b);
  RobinBC Bn(k, n, dy, a, b);
  RobinBC Bo(k, o, dz, a, b);
 
  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/robinbc.h
• src/cpp/robinbc.cpp