CppNoddy: src/ODE_EVP.cpp File Reference

A specification of a class for an $n^{th}$ -order ODE LINEAR EVP defined by. More...

#include <string>
#include <Types.h>
#include <Equation.h>
#include <Utility.h>
#include <ODE_EVP.h>
#include <Exceptions.h>
#include <DenseLinearEigenSystem.h>
#include <Equation_2matrix.h>

Go to the source code of this file.

Namespaces
namespace	CppNoddy
	A collection of OO numerical routines aimed at simple (typical) applied problems in continuum mechanics.

Detailed Description

A specification of a class for an $n^{th}$ -order ODE LINEAR EVP defined by.

${\underline f}^\prime (x) = {\underline R}( {\underline f}(x), x )\,,$

subject to zero Dirichlet conditions defined at $x = x_{left}$ or $x_{right}$ for some components of ${\underline f}(x)$ . The routine constructs a banded 2nd order finite-difference (banded) representation of the EVP in the form

$A_{MxM} {\underline x} = \lambda B_{MxM} {\underline x}$

which can then be solved via the LinearEigenSystem class where for an order O equation and N (uniformly spaced) mesh points. The default configuration uses the DenseLinearEigenSystem class for solution via LAPACK although ARPACK can be used if you really know the resulting system has the correct mass matrix properties (typically it wont!).

Definition in file ODE_EVP.cpp.

Namespaces

Detailed Description