Interpol Class Reference

Mimetic Interpolator operator. More...

#include <interpol.h>

Inheritance diagram for Interpol:

Public Member Functions
	Interpol (u32 m, Real c)
	1-D Mimetic Interpolator Constructor
	Interpol (u32 m, u32 n, Real c1, Real c2)
	2-D Mimetic Interpolator Constructor
	Interpol (u32 m, u32 n, u32 o, Real c1, Real c2, Real c3)
	3-D Mimetic Interpolator Constructor
	Interpol (bool type, u32 m, Real c)
	1-D Mimetic Interpolator Constructor
	Interpol (bool type, u32 m, u32 n, Real c1, Real c2)
	2-D Mimetic Interpolator Constructor
	Interpol (bool type, u32 m, u32 n, u32 o, Real c1, Real c2, Real c3)
	3-D Mimetic Interpolator Constructor

Detailed Description

Mimetic Interpolator operator.

Definition at line 25 of file interpol.h.

Constructor & Destructor Documentation

Interpol() [1/6]

Interpol::Interpol	(	u32	m,
		Real	c )

1-D Mimetic Interpolator Constructor

Parameters

m	Number of cells
c	Weight for ends, can be any value from 0.0<=c<=1.0

Definition at line 19 of file interpol.cpp.

                                : sp_mat(m + 1, m + 2) {
  assert(m >= 4);
  assert(c >= 0 && c <= 1);
 
  at(0, 0) = 1;
  at(m, m + 1) = 1;
 
  for (u32 i = 1; i < m; i++) {
    at(i, i) = c;
    at(i, i + 1) = 1 - c;
  }
}

Interpol() [2/6]

Interpol::Interpol	(	u32	m,
		u32	n,
		Real	c1,
		Real	c2 )

2-D Mimetic Interpolator Constructor

Parameters

m	Number of cells in x-direction
n	Number of cells in y-direction
c1	Weight for ends in x-direction, can be any value from 0.0<=c<=1.0
c2	Weight for ends in y-direction, can be any value from 0.0<=c<=1.0

Definition at line 33 of file interpol.cpp.

                                                 {
  Interpol Ix(m, c1);
  Interpol Iy(n, c2);
 
  sp_mat Im = speye(m + 2, m + 2);
  sp_mat In = speye(n + 2, n + 2);
 
  Im.shed_row(0);
  Im.shed_row(m);
  In.shed_row(0);
  In.shed_row(n);
 
  sp_mat I1 = Utils::spkron(In, Ix);
  sp_mat I2 = Utils::spkron(Iy, Im);
 
  // Dimensions = 2*m*n+m+n, (m+2)*(n+2)
  if (m != n)
    *this = Utils::spjoin_cols(I1, I2);
  else {
    sp_mat A1(2, 1);
    sp_mat A2(2, 1);
    A1(0, 0) = A2(1, 0) = 1.0;
    *this = Utils::spkron(A1, I1) + Utils::spkron(A2, I2);
  }
}

Interpol() [3/6]

Interpol::Interpol	(	u32	m,
		u32	n,
		u32	o,
		Real	c1,
		Real	c2,
		Real	c3 )

3-D Mimetic Interpolator Constructor

Parameters

m	Number of cells in x-direction
n	Number of cells in y-direction
o	Number of cells in z-direction
c1	Weight for ends in x-direction, can be any value from 0.0<=c<=1.0
c2	Weight for ends in y-direction, can be any value from 0.0<=c<=1.0
c3	Weight for ends in z-direction, can be any value from 0.0<=c<=1.0

Definition at line 60 of file interpol.cpp.

                                                                 {
  Interpol Ix(m, c1);
  Interpol Iy(n, c2);
  Interpol Iz(o, c3);
 
  sp_mat Im = speye(m + 2, m + 2);
  sp_mat In = speye(n + 2, n + 2);
  sp_mat Io = speye(o + 2, o + 2);
 
  Im.shed_row(0);
  Im.shed_row(m);
  In.shed_row(0);
  In.shed_row(n);
  Io.shed_row(0);
  Io.shed_row(o);
 
  sp_mat I1 = Utils::spkron(Utils::spkron(Io, In), Ix);
  sp_mat I2 = Utils::spkron(Utils::spkron(Io, Iy), Im);
  sp_mat I3 = Utils::spkron(Utils::spkron(Iz, In), Im);
 
  // Dimensions = 3*m*n*o+m*n+m*o+n*o, (m+2)*(n+2)*(o+2)
  if ((m != n) || (n != o))
    *this = Utils::spjoin_cols(Utils::spjoin_cols(I1, I2), I3);
  else {
    sp_mat A1(3, 1);
    sp_mat A2(3, 1);
    sp_mat A3(3, 1);
    A1(0, 0) = A2(1, 0) = A3(2, 0) = 1.0;
    *this =
        Utils::spkron(A1, I1) + Utils::spkron(A2, I2) + Utils::spkron(A3, I3);
  }
}

Interpol() [4/6]

Interpol::Interpol	(	bool	type,
		u32	m,
		Real	c )

1-D Mimetic Interpolator Constructor

Parameters

type	Dummy holder to trigger overloaded function
m	Number of cells
c	Weight for ends, can be any value from 0.0<=c<=1.0

Definition at line 94 of file interpol.cpp.

                                           : sp_mat(m + 2, m + 1) {
  assert(m >= 4 && "m >= 4");
  assert(c >= 0 && c <= 1 && "0 <= c <= 1");
 
  at(0, 0) = 1;
  at(m + 2 - 1, m + 1 - 1) = 1;
 
  vec avg = {c, 1 - c};
 
  int j = 0;
  for (int i = 1; i < m + 1; ++i) {
    at(i, j) = avg(0);
    at(i, j + 1) = avg(1);
    j++;
  }
}

Interpol() [5/6]

Interpol::Interpol	(	bool	type,
		u32	m,
		u32	n,
		Real	c1,
		Real	c2 )

2-D Mimetic Interpolator Constructor

Parameters

type	Dummy holder to trigger overloaded function
m	Number of cells in x-direction
n	Number of cells in y-direction
c1	Weight for ends in x-direction, can be any value from 0.0<=c<=1.0
c2	Weight for ends in y-direction, can be any value from 0.0<=c<=1.0

Definition at line 112 of file interpol.cpp.

                                                            {
  Interpol Ix(true, m, c1);
  Interpol Iy(true, n, c2);
 
  sp_mat Im(m + 2, m);
  Im.submat(1, 0, m, m - 1) = speye(m, m);
 
  sp_mat In(n + 2, n);
  In.submat(1, 0, n, n - 1) = speye(n, n);
 
  sp_mat Sx = Utils::spkron(In, Ix);
  sp_mat Sy = Utils::spkron(Iy, Im);
 
  *this = Utils::spjoin_rows(Sx, Sy);
}

Interpol() [6/6]

Interpol::Interpol	(	bool	type,
		u32	m,
		u32	n,
		u32	o,
		Real	c1,
		Real	c2,
		Real	c3 )

3-D Mimetic Interpolator Constructor

Parameters

type	Dummy holder to trigger overloaded function
m	Number of cells in x-direction
n	Number of cells in y-direction
o	Number of cells in z-direction
c1	Weight for ends in x-direction, can be any value from 0.0<=c<=1.0
c2	Weight for ends in y-direction, can be any value from 0.0<=c<=1.0
c3	Weight for ends in z-direction, can be any value from 0.0<=c<=1.0

Definition at line 129 of file interpol.cpp.

                                                                            {
  Interpol Ix(true, m, c1);
  Interpol Iy(true, n, c2);
  Interpol Iz(true, o, c3);
 
  sp_mat Im(m + 2, m);
  Im.submat(1, 0, m, m - 1) = speye(m, m);
 
  sp_mat In(n + 2, n);
  In.submat(1, 0, n, n - 1) = speye(n, n);
 
  sp_mat Io(o + 2, o);
  Io.submat(1, 0, o, o - 1) = speye(o, o);
 
  sp_mat Sx = Utils::spkron(Utils::spkron(Io, In), Ix);
  sp_mat Sy = Utils::spkron(Utils::spkron(Io, Iy), Im);
  sp_mat Sz = Utils::spkron(Utils::spkron(Iz, In), Im);
 
  *this = Utils::spjoin_rows(Utils::spjoin_rows(Sx, Sy), Sz);
}

---

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

• src/cpp/interpol.h
• src/cpp/interpol.cpp