CppNoddy  0.92
Loading...
Searching...
No Matches
Public Member Functions | List of all members
CppNoddy::Example::OS_evp_equation Class Reference

Define the OS equation for the global QZ EVP. More...

Inheritance diagram for CppNoddy::Example::OS_evp_equation:
CppNoddy::Equation_2matrix< D_complex > CppNoddy::Equation_2matrix< D_complex > CppNoddy::Equation_1matrix< _Type, _Xtype > CppNoddy::Equation_1matrix< _Type, _Xtype > CppNoddy::Residual_with_coords< _Type, _Xtype > CppNoddy::Residual_with_coords< _Type, _Xtype > CppNoddy::Residual< _Type > CppNoddy::Residual< _Type >

Public Member Functions

 OS_evp_equation ()
 The OS equation is a 4th order complex ODE. More...
 
void residual_fn (const DenseVector< D_complex > &z, DenseVector< D_complex > &g) const
 The OS equation. More...
 
void matrix0 (const DenseVector< D_complex > &z, DenseMatrix< D_complex > &m) const
 matrix to multiply the BVP coordinate More...
 
void matrix1 (const DenseVector< D_complex > &z, DenseMatrix< D_complex > &m) const
 Define the unsteady terms by providing the mass matrix Here we define the eigenvalue contribution to the g[ psid ] equation. More...
 
 OS_evp_equation ()
 The OS equation is a 4th order complex ODE. More...
 
void residual_fn (const DenseVector< D_complex > &z, DenseVector< D_complex > &g) const
 The OS equation. More...
 
void matrix0 (const DenseVector< D_complex > &z, DenseMatrix< D_complex > &m) const
 matrix to multiply the BVP coordinate More...
 
void matrix1 (const DenseVector< D_complex > &z, DenseMatrix< D_complex > &m) const
 Define the unsteady terms by providing the mass matrix Here we define the eigenvalue contribution to the g[ psid ] equation. More...
 
- Public Member Functions inherited from CppNoddy::Equation_2matrix< D_complex >
 Equation_2matrix (const unsigned &order)
 Constructor for equation class. More...
 
virtual ~Equation_2matrix ()
 An empty destructor, virtual since we have virtual methods. More...
 
void update (const DenseVector< D_complex > &state)
 Update the Equation object for the current set of state variables. More...
 
const DenseMatrix< D_complex > & matrix1 () const
 Return a handle to the matrix member data. More...
 
virtual void get_jacobian_of_matrix1_mult_vector (const DenseVector< D_complex > &state, const DenseVector< D_complex > &vec, DenseMatrix< D_complex > &h) const
 Return the product of the Jacobian-of-the-matrix and a vector 'vec' when the equation has a given 'state'. More...
 
- Public Member Functions inherited from CppNoddy::Equation_1matrix< _Type, _Xtype >
 Equation_1matrix (const unsigned &order)
 Constructor for equation class. More...
 
virtual ~Equation_1matrix ()
 An empty destructor, virtual since we have virtual methods. More...
 
void update (const DenseVector< _Type > &state)
 Update the Equation object for the current set of state variables. More...
 
const DenseMatrix< _Type > & matrix0 () const
 Return a handle to the matrix. More...
 
virtual void get_jacobian_of_matrix0_mult_vector (const DenseVector< _Type > &state, const DenseVector< _Type > &vec, DenseMatrix< _Type > &h) const
 Return the product of the Jacobian-of-the-matrix and a vector 'vec' when the equation has a given 'state'. More...
 
- Public Member Functions inherited from CppNoddy::Residual_with_coords< _Type, _Xtype >
 Residual_with_coords (const unsigned &order, const unsigned &ncoords)
 Constructor for a 'square' residual object that is, N residuals for N unknowns. More...
 
 Residual_with_coords (const unsigned &order, const unsigned &nvars, const unsigned &ncoords)
 Constructor for a 'non-square' residual object that is, there are less residual constraints than unknowns. More...
 
virtual ~Residual_with_coords ()
 An empty destructor. More...
 
_Xtype & coord (const unsigned &i)
 General handle access to the coordinates. More...
 
const _Xtype & coord (const unsigned &i) const
 General handle access to the coordinates. More...
 
- Public Member Functions inherited from CppNoddy::Residual< _Type >
 Residual (const unsigned &order)
 Constructor for a 'square' residual object that is, N residuals for N unknowns. More...
 
 Residual (const unsigned &order, const unsigned &nvars)
 Constructor for a 'non-square' residual object that is, there are less residual constraints than unknowns. More...
 
virtual ~Residual ()
 An empty destructor, virtual since we have virtual methods. More...
 
void update (const DenseVector< _Type > &state)
 Update the Residual object for the current set of state variables. More...
 
const DenseVector< _Type > & residual () const
 Return a handle to the residuals corresponding to the last update state. More...
 
const DenseMatrix< _Type > & jacobian () const
 Retrun a handle to the Jacobian of the residual corresponding to the last update state. More...
 
_Type & delta ()
 
const _Type & delta () const
 
unsigned get_order () const
 Get the order of the residual vector. More...
 
unsigned get_number_of_vars () const
 Get the number of variables that this residual condition is defined for. More...
 
virtual void residual_fn (const DenseVector< _Type > &state, DenseVector< _Type > &f) const
 A blank virtual residual function method. More...
 

Additional Inherited Members

virtual void matrix1 (const DenseVector< D_complex > &state, DenseMatrix< D_complex > &m) const
 Define the matrix in terms of the current state vector. More...
 
- Protected Member Functions inherited from CppNoddy::Equation_1matrix< _Type, _Xtype >
virtual void matrix0 (const DenseVector< _Type > &x, DenseMatrix< _Type > &m) const
 Define the matrix in terms of the current state vector. More...
 
- Protected Member Functions inherited from CppNoddy::Residual< _Type >
virtual void jacobian (const DenseVector< _Type > &state, DenseMatrix< _Type > &jac) const
 Because the residual evaluation at the current state is assumed to have already been done by the 'update' method, this routine is protected. More...
 
- Protected Attributes inherited from CppNoddy::Residual_with_coords< _Type, _Xtype >
std::vector< _Xtype > coords
 The coordinates stored for this residual. More...
 
- Protected Attributes inherited from CppNoddy::Residual< _Type >
DenseMatrix< _Type > JAC_AT_LAST_STATE
 Jacobian for the last state vector. More...
 
DenseVector< _Type > FN_AT_LAST_STATE
 Residual for the last state vector. More...
 
DenseVector< _Type > LAST_STATE
 The last state vector. More...
 
_Type DELTA
 A default step for FD computation of the Jacobian. More...
 
unsigned ORDER_OF_SYSTEM
 The order of the system of equations. More...
 
unsigned NUMBER_OF_VARS
 The number of elements in the state vector. More...
 

Detailed Description

Define the OS equation for the global QZ EVP.

Definition at line 47 of file EVPOrrSommerfeldEasy_lapack.cpp.

Constructor & Destructor Documentation

◆ OS_evp_equation() [1/2]

CppNoddy::Example::OS_evp_equation::OS_evp_equation ( )
inline

The OS equation is a 4th order complex ODE.

Definition at line 51 of file EVPOrrSommerfeldEasy_lapack.cpp.

51: Equation_2matrix<D_complex>( 4 ) {}
CppNoddy::Equation_2matrix
An equation object base class used in the PDE_double_IBVP class.
Definition: Equation_2matrix.h:21

◆ OS_evp_equation() [2/2]

CppNoddy::Example::OS_evp_equation::OS_evp_equation ( )
inline

The OS equation is a 4th order complex ODE.

Definition at line 53 of file EVPOrrSommerfeldNeutralCurve_lapack.cpp.

53: Equation_2matrix<D_complex>( 4 ) {}

Member Function Documentation

◆ matrix0() [1/2]

void CppNoddy::Example::OS_evp_equation::matrix0 ( const DenseVector< D_complex > &  z,
DenseMatrix< D_complex > &  m 
) const
inline

matrix to multiply the BVP coordinate

Definition at line 65 of file EVPOrrSommerfeldEasy_lapack.cpp.

66 {
67 Utils_Fill::fill_identity(m);
68 }
Utils_Fill::fill_identity
void fill_identity(CppNoddy::Sequential_Matrix_base< _Type > &A)
Fill diagonal with unit values.
Definition: Utils_Fill.h:22

References Utils_Fill::fill_identity(), and m.

◆ matrix0() [2/2]

void CppNoddy::Example::OS_evp_equation::matrix0 ( const DenseVector< D_complex > &  z,
DenseMatrix< D_complex > &  m 
) const
inline

matrix to multiply the BVP coordinate

Definition at line 67 of file EVPOrrSommerfeldNeutralCurve_lapack.cpp.

68 {
69 Utils_Fill::fill_identity(m);
70 }

References Utils_Fill::fill_identity(), and m.

◆ matrix1() [1/2]

void CppNoddy::Example::OS_evp_equation::matrix1 ( const DenseVector< D_complex > &  z,
DenseMatrix< D_complex > &  m 
) const
inlinevirtual

Define the unsteady terms by providing the mass matrix Here we define the eigenvalue contribution to the g[ psid ] equation.

  • D_complex( 0.0, 1.0 ) * alpha * Re * c * z[ psi ]

Reimplemented from CppNoddy::Equation_2matrix< D_complex >.

Definition at line 73 of file EVPOrrSommerfeldEasy_lapack.cpp.

74 {
75 // the eigenvalue is in equation 3, and multiplies unknown 2
76 m( 3, 2 ) = D_complex( 0.0, 1.0 ) * alpha * Re;
77 }
CppNoddy::Example::Re
double Re
Globally define the Reynolds number and wavenumber.
Definition: EVPOrrSommerfeldEasy_lapack.cpp:34
CppNoddy::D_complex
std::complex< double > D_complex
A complex double precision number using std::complex.
Definition: Types.h:98

References CppNoddy::Example::alpha, m, and CppNoddy::Example::Re.

◆ matrix1() [2/2]

void CppNoddy::Example::OS_evp_equation::matrix1 ( const DenseVector< D_complex > &  z,
DenseMatrix< D_complex > &  m 
) const
inlinevirtual

Define the unsteady terms by providing the mass matrix Here we define the eigenvalue contribution to the g[ psid ] equation.

  • D_complex( 0.0, 1.0 ) * alpha * Re * c * z[ psi ]

Reimplemented from CppNoddy::Equation_2matrix< D_complex >.

Definition at line 75 of file EVPOrrSommerfeldNeutralCurve_lapack.cpp.

76 {
77 // the eigenvalue is in equation 3, and multiplies unknown 2
78 m( 3, 2 ) = D_complex( 0.0, 1.0 ) * alpha * Re;
79 }

References CppNoddy::Example::alpha, m, and CppNoddy::Example::Re.

◆ residual_fn() [1/2]

void CppNoddy::Example::OS_evp_equation::residual_fn ( const DenseVector< D_complex > &  z,
DenseVector< D_complex > &  g 
) const
inline

The OS equation.

Definition at line 54 of file EVPOrrSommerfeldEasy_lapack.cpp.

55 {
56 // define the equation as 4 1st order equations
57 g[ phi ] = z[ phid ];
58 g[ phid ] = z[ psi ] + alpha * alpha * z[ phi ];
59 g[ psi ] = z[ psid ];
60 g[ psid ] = alpha * alpha * z[ psi ]
61 + D_complex( 0.0, 1.0 ) * alpha * Re * ( U( coord(0) ) * z[ psi ] - Udd( coord(0) ) * z[ phi ] );
62 }
U
@ U
Definition: BVPKarman.cpp:20
phid
@ phid
Definition: EVPOrrSommerfeldEasy_lapack.cpp:26
psid
@ psid
Definition: EVPOrrSommerfeldEasy_lapack.cpp:26
psi
@ psi
Definition: EVPOrrSommerfeldEasy_lapack.cpp:26
phi
@ phi
Definition: EVPOrrSommerfeldEasy_lapack.cpp:26
CppNoddy::Residual_with_coords::coord
_Xtype & coord(const unsigned &i)
General handle access to the coordinates.
Definition: Residual_with_coords.h:54
CppNoddy::Example::Udd
double Udd(double y)
Globally define the base flow curvature.
Definition: EVPOrrSommerfeldEasy_lapack.cpp:41
CppNoddy::Example::g
double g(1.0)
gravitational acceleration

References CppNoddy::Example::alpha, CppNoddy::Residual_with_coords< _Type, _Xtype >::coord(), CppNoddy::Example::g(), phi, phid, psi, psid, CppNoddy::Example::Re, U, CppNoddy::Example::Udd(), and CppNoddy::Example::z().

◆ residual_fn() [2/2]

void CppNoddy::Example::OS_evp_equation::residual_fn ( const DenseVector< D_complex > &  z,
DenseVector< D_complex > &  g 
) const
inline

The OS equation.

Definition at line 56 of file EVPOrrSommerfeldNeutralCurve_lapack.cpp.

57 {
58 // define the equation as four 1st order equations
59 g[ phi ] = z[ phid ];
60 g[ phid ] = z[ psi ] + alpha * alpha * z[ phi ];
61 g[ psi ] = z[ psid ];
62 g[ psid ] = alpha * alpha * z[ psi ]
63 + D_complex( 0.0, 1.0 ) * alpha * Re * ( U( coord(0) ) * z[ psi ] - Udd( coord(0) ) * z[ phi ] );
64 }

References CppNoddy::Example::alpha, CppNoddy::Residual_with_coords< _Type, _Xtype >::coord(), CppNoddy::Example::g(), phi, phid, psi, psid, CppNoddy::Example::Re, U, CppNoddy::Example::Udd(), and CppNoddy::Example::z().

The documentation for this class was generated from the following files:
Generated on Tue Jan 1 1980 00:00:00 for CppNoddy by   doxygen 1.9.5

© 2012

R.E. Hewitt