ITHACAPOD Namespace Reference

namespace for the computation of the POD, it exploits bot the method of snapshots and SVD. More...

Functions
void	weightedGramSchmidt (Eigen::MatrixXd &matrix, Eigen::VectorXd &weights)
	Weighted Gram-Schmidt orthogonalization.
template<class Type, template< class > class PatchField, class GeoMesh>
void	getNestedSnapshotMatrix (PtrList< GeometricField< Type, PatchField, GeoMesh > > &snapshots, PtrList< GeometricField< Type, PatchField, GeoMesh > > &modes, word fieldName, label Npar, label NnestedOut)
	Nested-POD approach.
template void	getNestedSnapshotMatrix (PtrList< volScalarField > &snapshots, PtrList< volScalarField > &ModesGlobal, word fieldName, label Npar, label NnestedOut)
template void	getNestedSnapshotMatrix (PtrList< volVectorField > &snapshots, PtrList< volVectorField > &ModesGlobal, word fieldName, label Npar, label NnestedOut)
template<class Type, template< class > class PatchField, class GeoMesh>
void	getModes (PtrList< GeometricField< Type, PatchField, GeoMesh > > &snapshots, PtrList< GeometricField< Type, PatchField, GeoMesh > > &modes, word fieldName, bool podex, bool supex=0, bool sup=0, label nmodes=0, bool correctBC=true)
	Computes the bases or reads them for a field.
template void	getModes (PtrList< volVectorField > &snapshots, PtrList< volVectorField > &modes, word fieldName, bool podex, bool supex, bool sup, label nmodes, bool correctBC)
template void	getModes (PtrList< volScalarField > &snapshots, PtrList< volScalarField > &modes, word fieldName, bool podex, bool supex, bool sup, label nmodes, bool correctBC)
template void	getModes (PtrList< surfaceScalarField > &snapshots, PtrList< surfaceScalarField > &modes, word fieldName, bool podex, bool supex, bool sup, label nmodes, bool correctBC)
template<class Type, template< class > class PatchField, class GeoMesh>
void	getModesMemoryEfficient (GeometricField< Type, PatchField, GeoMesh > &templateField, word snapshotsPath, PtrList< GeometricField< Type, PatchField, GeoMesh > > &modes, word fieldName, bool podex=false, bool supex=false, bool sup=false, label nmodes=0, bool correctBC=true, autoPtr< GeometricField< Type, PatchField, GeoMesh > > meanField=NULL)
	Gets the modes in a memory-efficient manner.
template void	getModesMemoryEfficient (GeometricField< scalar, fvPatchField, volMesh > &, word, PtrList< GeometricField< scalar, fvPatchField, volMesh > > &, word, bool, bool, bool, label, bool, autoPtr< GeometricField< scalar, fvPatchField, volMesh > >)
template void	getModesMemoryEfficient (GeometricField< vector, fvPatchField, volMesh > &, word, PtrList< GeometricField< vector, fvPatchField, volMesh > > &, word, bool, bool, bool, label, bool, autoPtr< GeometricField< vector, fvPatchField, volMesh > >)
template<class Type, template< class > class PatchField, class GeoMesh>
void	getMeanMemoryEfficient (GeometricField< Type, PatchField, GeoMesh > &templateField, word snapshotsPath, autoPtr< GeometricField< Type, PatchField, GeoMesh > > &meanField, bool meanex=false)
	Gets the mean field in a memory-efficient manner.
template void	getMeanMemoryEfficient (GeometricField< scalar, fvPatchField, volMesh > &, word snapshotsPath, autoPtr< GeometricField< scalar, fvPatchField, volMesh > > &, bool)
template void	getMeanMemoryEfficient (GeometricField< vector, fvPatchField, volMesh > &, word snapshotsPath, autoPtr< GeometricField< vector, fvPatchField, volMesh > > &, bool)
template<class Type, template< class > class PatchField, class GeoMesh>
void	getWeightedModes (PtrList< GeometricField< Type, PatchField, GeoMesh > > &snapshots, PtrList< GeometricField< Type, PatchField, GeoMesh > > &modes, word fieldName, bool podex, bool supex=0, bool sup=0, label nmodes=0, bool correctBC=true)
	Computes the weighted bases (using the nested-pod approach) or read them for a vector field.
template void	getWeightedModes (PtrList< volScalarField > &snapshots, PtrList< volScalarField > &modes, word fieldName, bool podex, bool supex, bool sup, label nmodes, bool correctBC)
template void	getWeightedModes (PtrList< volVectorField > &snapshots, PtrList< volVectorField > &modes, word fieldName, bool podex, bool supex, bool sup, label nmodes, bool correctBC)
template<class Type, template< class > class PatchField, class GeoMesh>
void	getModesSVD (PtrList< GeometricField< Type, PatchField, GeoMesh > > &snapshots, PtrList< GeometricField< Type, PatchField, GeoMesh > > &modes, word fieldName, bool podex, bool supex=0, bool sup=0, label nmodes=0, bool correctBC=true)
	Gets the bases for a scalar field using SVD instead of the method of snapshots.
template void	getModesSVD (PtrList< volScalarField > &snapshots, PtrList< volScalarField > &modes, word fieldName, bool podex, bool supex, bool sup, label nmodes, bool correctBC)
template void	getModesSVD (PtrList< volVectorField > &snapshots, PtrList< volVectorField > &modes, word fieldName, bool podex, bool supex, bool sup, label nmodes, bool correctBC)
template<>
Eigen::MatrixXd	corMatrix (PtrList< volScalarField > &snapshots)
	Construct the Correlation Matrix for Scalar Field.
template<>
Eigen::MatrixXd	corMatrix (PtrList< volVectorField > &snapshots)
	Construct the Correlation Matrix for Vector Field.
template<>
Eigen::MatrixXd	corMatrix (List< Eigen::SparseMatrix< double > > &snapshots)
	Construct the Correlation Matrix for Vector Field.
template<>
Eigen::MatrixXd	corMatrix (List< Eigen::VectorXd > &snapshots)
	Construct the Correlation Matrix for Vector Field.
template<class Type, template< class > class PatchField, class GeoMesh>
void	exportBases (PtrList< GeometricField< Type, PatchField, GeoMesh > > &s, PtrList< GeometricField< Type, PatchField, GeoMesh > > &bases, word fieldName, bool sup)
	Export the Bases.
template void	exportBases (PtrList< volVectorField > &s, PtrList< volVectorField > &bases, word fieldName, bool sup)
template void	exportBases (PtrList< volScalarField > &s, PtrList< volScalarField > &bases, word fieldName, bool sup)
void	exportEigenvalues (scalarField Eigenvalues, fileName name, bool sup=0)
	Exports the eigenvalues as a txt file.
void	exportcumEigenvalues (scalarField cumEigenvalues, fileName name, bool sup=0)
	Export the cumulative eigenvalues as a txt file.
std::tuple< List< Eigen::SparseMatrix< double > >, List< Eigen::VectorXd > >	DEIMmodes (List< Eigen::SparseMatrix< double > > &A, List< Eigen::VectorXd > &b, label nmodesA, label nmodesB, word MatrixName="Matrix")
	Get the DEIM modes for a generic non a parametrized matrix coming from a differential operator function.
void	GrammSchmidt (Eigen::MatrixXd &Matrix)
	Performs GrammSchmidt orthonormalization on an Eigen Matrix.
template<class Type, template< class > class PatchField, class GeoMesh>
void	getModes (PtrList< GeometricField< Type, PatchField, GeoMesh > > &snapshots, PtrList< GeometricField< Type, PatchField, GeoMesh > > &modes, PtrList< volScalarField > &Volumes, word fieldName, bool podex, bool supex, bool sup, label nmodes, bool correctBC)
template void	getModes (PtrList< volScalarField > &snapshots, PtrList< volScalarField > &modes, PtrList< volScalarField > &Volumes, word fieldName, bool podex, bool supex, bool sup, label nmodes, bool correctBC)
template void	getModes (PtrList< volVectorField > &snapshots, PtrList< volVectorField > &modes, PtrList< volScalarField > &Volumes, word fieldName, bool podex, bool supex, bool sup, label nmodes, bool correctBC)
template<typename type_matrix>
std::tuple< List< Eigen::SparseMatrix< double > >, List< Eigen::VectorXd > >	DEIMmodes (PtrList< type_matrix > &MatrixList, label nmodesA, label nmodesB, word MatrixName="Matrix")
	Get the DEIM modes for a generic non-parametrized matrix coming from a differential operator function.
template std::tuple< List< Eigen::SparseMatrix< double > >, List< Eigen::VectorXd > >	DEIMmodes (PtrList< fvScalarMatrix > &MatrixList, label nmodesA, label nmodesB, word MatrixName)
template std::tuple< List< Eigen::SparseMatrix< double > >, List< Eigen::VectorXd > >	DEIMmodes (PtrList< fvVectorMatrix > &MatrixList, label nmodesA, label nmodesB, word MatrixName)
template<class Type, template< class > class PatchField, class GeoMesh>
PtrList< GeometricField< Type, PatchField, GeoMesh > >	DEIMmodes (PtrList< GeometricField< Type, PatchField, GeoMesh > > &SnapShotsMatrix, label nmodes, word FunctionName, word FieldName)
	Get the DEIM modes for a generic non linear function.
template<class Field_type, class Field_type_2>
void	getModes (PtrList< Field_type > &snapshots, PtrList< Field_type > &modes, PtrList< Field_type_2 > &fields2, word fieldName, bool podex, bool supex, bool sup, label nmodes, bool correctBC)
template void	getModes (PtrList< surfaceScalarField > &snapshots, PtrList< surfaceScalarField > &modes, PtrList< volVectorField > &fields2, word fieldName, bool podex, bool supex, bool sup, label nmodes, bool correctBC)
template void	getModes (PtrList< volScalarField > &snapshots, PtrList< volScalarField > &modes, PtrList< volVectorField > &fields2, word fieldName, bool podex, bool supex, bool sup, label nmodes, bool correctBC)
template PtrList< volScalarField >	DEIMmodes (PtrList< volScalarField > &SnapShotsMatrix, label nmodes, word FunctionName, word FieldName)
template PtrList< volVectorField >	DEIMmodes (PtrList< volVectorField > &SnapShotsMatrix, label nmodes, word FunctionName, word FieldName)
template<>
scalar	computeInnerProduct (const GeometricField< scalar, fvPatchField, volMesh > &field1, const GeometricField< scalar, fvPatchField, volMesh > &field2)
template<>
scalar	computeInnerProduct (const GeometricField< vector, fvPatchField, volMesh > &field1, const GeometricField< vector, fvPatchField, volMesh > &field2)
template<>
scalar	computeFrobeniusInnerProduct (const GeometricField< scalar, fvPatchField, volMesh > &field1, const GeometricField< scalar, fvPatchField, volMesh > &field2)
template<>
scalar	computeFrobeniusInnerProduct (const GeometricField< vector, fvPatchField, volMesh > &field1, const GeometricField< vector, fvPatchField, volMesh > &field2)
template<class Field_type>
Eigen::MatrixXd	corMatrix (Field_type &snapshots)
	Computes the correlation matrix given a vector field snapshot Matrix using different norms depending on the input snapshots.
template<class Type, template< class > class PatchField, class GeoMesh>
void	getModes (PtrList< GeometricField< Type, PatchField, GeoMesh > > &snapshots, PtrList< GeometricField< Type, PatchField, GeoMesh > > &modes, PtrList< volScalarField > &Volumes, word fieldName, bool podex, bool supex=0, bool sup=0, label nmodes=0)
	Gets the modes.
template<class Field_type, class Field_type_2>
void	getModes (PtrList< Field_type > &snapshots, PtrList< Field_type > &modes, PtrList< Field_type_2 > &fields2, word fieldName, bool podex, bool supex=0, bool sup=0, label nmodes=0)
	Gets Modes based on eigenvalue decomposition of another Field.
template<class Type, template< class > class PatchField, class GeoMesh>
scalar	computeInnerProduct (const GeometricField< Type, PatchField, GeoMesh > &field1, const GeometricField< Type, PatchField, GeoMesh > &field2)
	Helper function to compute inner product value for different field types.
template<class Type, template< class > class PatchField, class GeoMesh>
scalar	computeFrobeniusInnerProduct (const GeometricField< Type, PatchField, GeoMesh > &field1, const GeometricField< Type, PatchField, GeoMesh > &field2)
	Helper function to compute Frobenius norm inner product.

Detailed Description

namespace for the computation of the POD, it exploits bot the method of snapshots and SVD.

Function Documentation

computeFrobeniusInnerProduct() [1/3]

template<>

scalar ITHACAPOD::computeFrobeniusInnerProduct	(	const GeometricField< scalar, fvPatchField, volMesh > &	field1,
		const GeometricField< scalar, fvPatchField, volMesh > &	field2 )

Definition at line 2047 of file ITHACAPOD.C.

computeFrobeniusInnerProduct() [2/3]

template<class Type, template< class > class PatchField, class GeoMesh>

scalar ITHACAPOD::computeFrobeniusInnerProduct	(	const GeometricField< Type, PatchField, GeoMesh > &	field1,
		const GeometricField< Type, PatchField, GeoMesh > &	field2 )

Helper function to compute Frobenius norm inner product.

Parameters

[in]	field1	First field
[in]	field2	Second field

Returns

Inner product value as scalar

computeFrobeniusInnerProduct() [3/3]

template<>

scalar ITHACAPOD::computeFrobeniusInnerProduct	(	const GeometricField< vector, fvPatchField, volMesh > &	field1,
		const GeometricField< vector, fvPatchField, volMesh > &	field2 )

Definition at line 2055 of file ITHACAPOD.C.

computeInnerProduct() [1/3]

template<>

scalar ITHACAPOD::computeInnerProduct	(	const GeometricField< scalar, fvPatchField, volMesh > &	field1,
		const GeometricField< scalar, fvPatchField, volMesh > &	field2 )

Definition at line 2031 of file ITHACAPOD.C.

computeInnerProduct() [2/3]

template<class Type, template< class > class PatchField, class GeoMesh>

scalar ITHACAPOD::computeInnerProduct	(	const GeometricField< Type, PatchField, GeoMesh > &	field1,
		const GeometricField< Type, PatchField, GeoMesh > &	field2 )

Helper function to compute inner product value for different field types.

Parameters

[in]	field1	First field
[in]	field2	Second field

Returns

Inner product value as scalar

computeInnerProduct() [3/3]

template<>

scalar ITHACAPOD::computeInnerProduct	(	const GeometricField< vector, fvPatchField, volMesh > &	field1,
		const GeometricField< vector, fvPatchField, volMesh > &	field2 )

Definition at line 2039 of file ITHACAPOD.C.

corMatrix() [1/5]

template<class Field_type>

Eigen::MatrixXd ITHACAPOD::corMatrix

(

Field_type &

snapshots

)

Computes the correlation matrix given a vector field snapshot Matrix using different norms depending on the input snapshots.

Parameters

[in]

snapshots

List of snapshots.

Template Parameters

Field_type

Type of the input list you are passing, it can be: PtrList<volVectorField>, PtrList<volScalarField>, List<Eigen::SparseMatrix<double>>, List<Eigen::VectorXd>.

Returns

the Eigen::MatrixXd correlation matrix.

corMatrix() [2/5]

template<>

Eigen::MatrixXd ITHACAPOD::corMatrix

(

List< Eigen::SparseMatrix< double > > &

snapshots

)

Construct the Correlation Matrix for Vector Field.

Definition at line 978 of file ITHACAPOD.C.

corMatrix() [3/5]

template<>

Eigen::MatrixXd ITHACAPOD::corMatrix

(

List< Eigen::VectorXd > &

snapshots

)

Construct the Correlation Matrix for Vector Field.

Definition at line 1011 of file ITHACAPOD.C.

corMatrix() [4/5]

template<>

Eigen::MatrixXd ITHACAPOD::corMatrix

(

PtrList< volScalarField > &

snapshots

)

Construct the Correlation Matrix for Scalar Field.

Definition at line 923 of file ITHACAPOD.C.

corMatrix() [5/5]

template<>

Eigen::MatrixXd ITHACAPOD::corMatrix

(

PtrList< volVectorField > &

snapshots

)

Construct the Correlation Matrix for Vector Field.

Definition at line 951 of file ITHACAPOD.C.

DEIMmodes() [1/7]

std::tuple< List< Eigen::SparseMatrix< double > >, List< Eigen::VectorXd > > ITHACAPOD::DEIMmodes	(	List< Eigen::SparseMatrix< double > > &	A,
		List< Eigen::VectorXd > &	b,
		label	nmodesA,
		label	nmodesB,
		word	MatrixName = "Matrix" )

Get the DEIM modes for a generic non a parametrized matrix coming from a differential operator function.

Parameters

	A	List of matrices for the lhs
	b	List of vectors for the rhs
[in]	nmodesA	The number of modes
[in]	nmodesB	The nmodes b
[in]	MatrixName	The matrix name, used to save the matrix
[in]	MatrixList	The matrix list as a list of pointers

Template Parameters

The	type of matrix, can be fvScalarMatrix or fvVectorMatrix

Returns

It returns a tuple containing the list of modes for the matrix and the list of modes for the source term

Examples

08DEIM.C.

Definition at line 1138 of file ITHACAPOD.C.

DEIMmodes() [2/7]

template std::tuple< List< Eigen::SparseMatrix< double > >, List< Eigen::VectorXd > > ITHACAPOD::DEIMmodes	(	PtrList< fvScalarMatrix > &	MatrixList,
		label	nmodesA,
		label	nmodesB,
		word	MatrixName )

DEIMmodes() [3/7]

template std::tuple< List< Eigen::SparseMatrix< double > >, List< Eigen::VectorXd > > ITHACAPOD::DEIMmodes	(	PtrList< fvVectorMatrix > &	MatrixList,
		label	nmodesA,
		label	nmodesB,
		word	MatrixName )

DEIMmodes() [4/7]

template<class Type, template< class > class PatchField, class GeoMesh>

PtrList< GeometricField< Type, PatchField, GeoMesh > > ITHACAPOD::DEIMmodes	(	PtrList< GeometricField< Type, PatchField, GeoMesh > > &	SnapShotsMatrix,
		label	nmodes,
		word	FunctionName,
		word	FieldName )

Get the DEIM modes for a generic non linear function.

Parameters

[in]	SnapShotsMatrix	The snapshots matrix
[in]	nmodes	The number of modes
[in]	FunctionName	The function name

Template Parameters

Type	templated object for the snapshots matrix type

Returns

Type the POD modes

Definition at line 1681 of file ITHACAPOD.C.

DEIMmodes() [5/7]

template<typename type_matrix>

std::tuple< List< Eigen::SparseMatrix< double > >, List< Eigen::VectorXd > > ITHACAPOD::DEIMmodes	(	PtrList< type_matrix > &	MatrixList,
		label	nmodesA,
		label	nmodesB,
		word	MatrixName = "Matrix" )

Get the DEIM modes for a generic non-parametrized matrix coming from a differential operator function.

Parameters

[in]	MatrixList	The matrix list as a list of pointers
[in]	nmodesA	The number of modes for A
[in]	nmodesB	The number of modes for B
[in]	MatrixName	The matrix name, used to save the matrix

Template Parameters

type_matrix

The type of matrix, can be fvScalarMatrix or fvVectorMatrix

Returns

It returns a tuple containing the list of modes for the matrix and the list of modes for the source term

Definition at line 1506 of file ITHACAPOD.C.

DEIMmodes() [6/7]

template PtrList< volScalarField > ITHACAPOD::DEIMmodes	(	PtrList< volScalarField > &	SnapShotsMatrix,
		label	nmodes,
		word	FunctionName,
		word	FieldName )

DEIMmodes() [7/7]

template PtrList< volVectorField > ITHACAPOD::DEIMmodes	(	PtrList< volVectorField > &	SnapShotsMatrix,
		label	nmodes,
		word	FunctionName,
		word	FieldName )

exportBases() [1/3]

template<class Type, template< class > class PatchField, class GeoMesh>

void ITHACAPOD::exportBases	(	PtrList< GeometricField< Type, PatchField, GeoMesh > > &	snapshots,
		PtrList< GeometricField< Type, PatchField, GeoMesh > > &	bases,
		word	fieldName,
		bool	sup = 0 )

Export the Bases.

Exports the basis for an OpenFOAM GeometricField into the ITHACAOutput/POD or ITHACAOutput/supremizer.

Parameters

	snapshots	The snapshots
[in]	bases	a PtrList of GeometricField where the bases are stored.
[in]	fieldName	The field name
[in]	sup	a boolean variable 1 if you want to export the supremizer bases (in ITHACAOutput/supremizer) 0 elsewhere (Default is 0).
[in]	s	a PtrList of GeometricField where the snapshots associated with the bases are stored.

Template Parameters

Type	vector or scalar.
PatchField	fvPatchField or fvsPatchField.
GeoMesh	volMesh or surfaceMesh.

Definition at line 1040 of file ITHACAPOD.C.

exportBases() [2/3]

template void ITHACAPOD::exportBases	(	PtrList< volScalarField > &	s,
		PtrList< volScalarField > &	bases,
		word	fieldName,
		bool	sup )

exportBases() [3/3]

template void ITHACAPOD::exportBases	(	PtrList< volVectorField > &	s,
		PtrList< volVectorField > &	bases,
		word	fieldName,
		bool	sup )

exportcumEigenvalues()

void ITHACAPOD::exportcumEigenvalues	(	scalarField	cumEigenvalues,
		fileName	name,
		bool	sup = 0 )

Export the cumulative eigenvalues as a txt file.

Parameters

[in]	cumEigenvalues	a scalarField of cumEigenvalues.
[in]	name	a fileName to indicate the name of the txt file.
[in]	sup	If 1, it exports the cumEigenValues in "ITHACAOutput/supremizer". If 0, no actions (Default is 0).

Definition at line 1107 of file ITHACAPOD.C.

exportEigenvalues()

void ITHACAPOD::exportEigenvalues	(	scalarField	Eigenvalues,
		fileName	name,
		bool	sup = 0 )

Exports the eigenvalues as a txt file.

Parameters

[in]	Eigenvalues	a scalarField of eigenvalues.
[in]	name	a fileName to indicate the name of the txt file.
[in]	sup	If 1, it exports the eigenValues in "ITHACAOutput/supremizer". If 0, no actions (Default is 0).

Definition at line 1078 of file ITHACAPOD.C.

getMeanMemoryEfficient() [1/3]

template void ITHACAPOD::getMeanMemoryEfficient	(	GeometricField< scalar, fvPatchField, volMesh > &	,
		word	snapshotsPath,
		autoPtr< GeometricField< scalar, fvPatchField, volMesh > > &	,
		bool	)

getMeanMemoryEfficient() [2/3]

template<class Type, template< class > class PatchField, class GeoMesh>

void ITHACAPOD::getMeanMemoryEfficient	(	GeometricField< Type, PatchField, GeoMesh > &	templateField,
		word	snapshotsPath,
		autoPtr< GeometricField< Type, PatchField, GeoMesh > > &	meanField,
		bool	meanex )

Gets the mean field in a memory-efficient manner.

Definition at line 633 of file ITHACAPOD.C.

getMeanMemoryEfficient() [3/3]

template void ITHACAPOD::getMeanMemoryEfficient	(	GeometricField< vector, fvPatchField, volMesh > &	,
		word	snapshotsPath,
		autoPtr< GeometricField< vector, fvPatchField, volMesh > > &	,
		bool	)

getModes() [1/12]

template<class Field_type, class Field_type_2>

void ITHACAPOD::getModes	(	PtrList< Field_type > &	snapshots,
		PtrList< Field_type > &	modes,
		PtrList< Field_type_2 > &	fields2,
		word	fieldName,
		bool	podex,
		bool	supex,
		bool	sup,
		label	nmodes,
		bool	correctBC )

Definition at line 1872 of file ITHACAPOD.C.

getModes() [2/12]

template<class Field_type, class Field_type_2>

void ITHACAPOD::getModes	(	PtrList< Field_type > &	snapshots,
		PtrList< Field_type > &	modes,
		PtrList< Field_type_2 > &	fields2,
		word	fieldName,
		bool	podex,
		bool	supex = 0,
		bool	sup = 0,
		label	nmodes = 0 )

Gets Modes based on eigenvalue decomposition of another Field.

Parameters

	snapshots	The snapshots of which we want to create the modes
	modes	The modes of the snapshots
	fields2	The fields that we want to use for the eigenvalue decomposition
[in]	fieldName	The field name
[in]	podex	boolean variable 1 if the modes have been already computed and stored (in this case the function is reading them) 0 elsewhere.
[in]	supex	boolean variable 1 if the supremizer modes have been already computed and stored (in this case the function is reading them) 0 elsewhere.
[in]	sup	boolean variable 1 if you want to compute the supremizer modes 0 elsewhere.
[in]	nmodes	label variable to set the number of modes to be stored, if set to 0 the maximum number of modes will computed.

Template Parameters

Field_type	Type of field snapshots
Field_type_2	Type of the other fields

getModes() [3/12]

template<class Type, template< class > class PatchField, class GeoMesh>

void ITHACAPOD::getModes	(	PtrList< GeometricField< Type, PatchField, GeoMesh > > &	snapshots,
		PtrList< GeometricField< Type, PatchField, GeoMesh > > &	modes,
		PtrList< volScalarField > &	Volumes,
		word	fieldName,
		bool	podex,
		bool	supex,
		bool	sup,
		label	nmodes,
		bool	correctBC )

Definition at line 1345 of file ITHACAPOD.C.

getModes() [4/12]

template<class Type, template< class > class PatchField, class GeoMesh>

void ITHACAPOD::getModes	(	PtrList< GeometricField< Type, PatchField, GeoMesh > > &	snapshots,
		PtrList< GeometricField< Type, PatchField, GeoMesh > > &	modes,
		PtrList< volScalarField > &	Volumes,
		word	fieldName,
		bool	podex,
		bool	supex = 0,
		bool	sup = 0,
		label	nmodes = 0 )

Gets the modes.

Parameters

[in]	snapshots	List of snapshots.
[out]	modes	A PtrList where modes are stored (it must be passed empty).
	Volumes	The volumes
[in]	fieldName	The field name
[in]	podex	If 1, the functions read the stored mode. If 0, the function computes the modes and stores them.
[in]	supex	If 1, the functions read the stored supremizer mode. If 0, the function computes and stores them.
[in]	sup	If 1 it computes the supremizer modes.
[in]	nmodes	Number of modes to be stored. If 0, the maximum number of modes will computed.

Template Parameters

Type	vector or scalar.
PatchField	fvPatchField or fvsPatchField.
GeoMesh	volMesh or surfaceMesh.

getModes() [5/12]

template<class Type, template< class > class PatchField, class GeoMesh>

void ITHACAPOD::getModes	(	PtrList< GeometricField< Type, PatchField, GeoMesh > > &	snapshots,
		PtrList< GeometricField< Type, PatchField, GeoMesh > > &	modes,
		word	fieldName,
		bool	podex,
		bool	supex = 0,
		bool	sup = 0,
		label	nmodes = 0,
		bool	correctBC = true )

Computes the bases or reads them for a field.

Parameters

[in]	snapshots	List of snapshots.
[out]	modes	A PtrList where modes are stored (it must be passed empty).
[in]	fieldName	The field name
[in]	podex	If 1, the functions read the stored mode. If 0, the function computes the modes and stores them.
[in]	supex	If 1, the functions read the stored supremizer mode. If 0, the function computes and stores them.
[in]	sup	If 1 it computes the supremizer modes.
[in]	nmodes	Number of modes to be stored. If 0, the maximum number of modes will computed.

Template Parameters

Type	vector or scalar.
PatchField	fvPatchField or fvsPatchField.
GeoMesh	volMesh or surfaceMesh.

Examples

02thermalBlock.C, 03steadyNS.C, 04unsteadyNS.C, 06POD_RBF.C, 09DEIM_ROM.C, 10UnsteadyBBEnclosed.C, 11UnsteadyBBOpen.C, 17YJunction.C, 19UnsteadyNSExplicit.C, and 20incrementalPOD.C.

Definition at line 90 of file ITHACAPOD.C.

getModes() [6/12]

template void ITHACAPOD::getModes	(	PtrList< surfaceScalarField > &	snapshots,
		PtrList< surfaceScalarField > &	modes,
		PtrList< volVectorField > &	fields2,
		word	fieldName,
		bool	podex,
		bool	supex,
		bool	sup,
		label	nmodes,
		bool	correctBC )

getModes() [7/12]

template void ITHACAPOD::getModes	(	PtrList< surfaceScalarField > &	snapshots,
		PtrList< surfaceScalarField > &	modes,
		word	fieldName,
		bool	podex,
		bool	supex,
		bool	sup,
		label	nmodes,
		bool	correctBC )

getModes() [8/12]

template void ITHACAPOD::getModes	(	PtrList< volScalarField > &	snapshots,
		PtrList< volScalarField > &	modes,
		PtrList< volScalarField > &	Volumes,
		word	fieldName,
		bool	podex,
		bool	supex,
		bool	sup,
		label	nmodes,
		bool	correctBC )

getModes() [9/12]

template void ITHACAPOD::getModes	(	PtrList< volScalarField > &	snapshots,
		PtrList< volScalarField > &	modes,
		PtrList< volVectorField > &	fields2,
		word	fieldName,
		bool	podex,
		bool	supex,
		bool	sup,
		label	nmodes,
		bool	correctBC )

getModes() [10/12]

template void ITHACAPOD::getModes	(	PtrList< volScalarField > &	snapshots,
		PtrList< volScalarField > &	modes,
		word	fieldName,
		bool	podex,
		bool	supex,
		bool	sup,
		label	nmodes,
		bool	correctBC )

getModes() [11/12]

template void ITHACAPOD::getModes	(	PtrList< volVectorField > &	snapshots,
		PtrList< volVectorField > &	modes,
		PtrList< volScalarField > &	Volumes,
		word	fieldName,
		bool	podex,
		bool	supex,
		bool	sup,
		label	nmodes,
		bool	correctBC )

getModes() [12/12]

template void ITHACAPOD::getModes	(	PtrList< volVectorField > &	snapshots,
		PtrList< volVectorField > &	modes,
		word	fieldName,
		bool	podex,
		bool	supex,
		bool	sup,
		label	nmodes,
		bool	correctBC )

getModesMemoryEfficient() [1/3]

template void ITHACAPOD::getModesMemoryEfficient	(	GeometricField< scalar, fvPatchField, volMesh > &	,
		word	,
		PtrList< GeometricField< scalar, fvPatchField, volMesh > > &	,
		word	,
		bool	,
		bool	,
		bool	,
		label	,
		bool	,
		autoPtr< GeometricField< scalar, fvPatchField, volMesh > >	)

getModesMemoryEfficient() [2/3]

template<class Type, template< class > class PatchField, class GeoMesh>

void ITHACAPOD::getModesMemoryEfficient	(	GeometricField< Type, PatchField, GeoMesh > &	templateField,
		word	snapshotsPath,
		PtrList< GeometricField< Type, PatchField, GeoMesh > > &	modes,
		word	fieldName,
		bool	podex = false,
		bool	supex = false,
		bool	sup = false,
		label	nmodes = 0,
		bool	correctBC = true,
		autoPtr< GeometricField< Type, PatchField, GeoMesh > >	meanField = NULL )

Gets the modes in a memory-efficient manner.

Parameters

[in]	templateField	The template field
[in]	snapshotsPath	The path to the snapshots
[out]	modes	The modes
[in]	fieldName	The field name
[in]	podex	If 1, the functions read the stored mode. If 0, the function computes the modes and stores them.
[in]	supex	If 1, the functions read the stored supremizer mode. If 0, the function computes and stores them.
[in]	sup	If 1, it computes the supremizer modes.
[in]	nmodes	Number of modes to be stored. If 0, the maximum number of modes will computed.
[in]	correctBC	If true, correct the boundary conditions.

Template Parameters

Type	The type of the field
PatchField	The patch field type
GeoMesh	The mesh type

Definition at line 317 of file ITHACAPOD.C.

getModesMemoryEfficient() [3/3]

template void ITHACAPOD::getModesMemoryEfficient	(	GeometricField< vector, fvPatchField, volMesh > &	,
		word	,
		PtrList< GeometricField< vector, fvPatchField, volMesh > > &	,
		word	,
		bool	,
		bool	,
		bool	,
		label	,
		bool	,
		autoPtr< GeometricField< vector, fvPatchField, volMesh > >	)

getModesSVD() [1/3]

template<class Type, template< class > class PatchField, class GeoMesh>

void ITHACAPOD::getModesSVD	(	PtrList< GeometricField< Type, PatchField, GeoMesh > > &	snapshots,
		PtrList< GeometricField< Type, PatchField, GeoMesh > > &	modes,
		word	fieldName,
		bool	podex,
		bool	supex = 0,
		bool	sup = 0,
		label	nmodes = 0,
		bool	correctBC = true )

Gets the bases for a scalar field using SVD instead of the method of snapshots.

Parameters

[in]	snapshots	List of snapshots.
[out]	modes	A PtrList where modes are stored (it must be passed empty).
[in]	fieldName	The file name
[in]	podex	If 1, the functions read the stored mode. If 0, the function computes the modes and stores them.
[in]	supex	If 1, the functions read the stored supremizer mode. If 0, the function computes and stores them.
[in]	sup	If 1 it computes the supremizer modes.
[in]	nmodes	Number of modes to be stored. If 0, the maximum number of modes will computed.

Template Parameters

Type	vector or scalar.
PatchField	fvPatchField or fvsPatchField.
GeoMesh	volMesh or surfaceMesh.

Definition at line 826 of file ITHACAPOD.C.

getModesSVD() [2/3]

template void ITHACAPOD::getModesSVD	(	PtrList< volScalarField > &	snapshots,
		PtrList< volScalarField > &	modes,
		word	fieldName,
		bool	podex,
		bool	supex,
		bool	sup,
		label	nmodes,
		bool	correctBC )

getModesSVD() [3/3]

template void ITHACAPOD::getModesSVD	(	PtrList< volVectorField > &	snapshots,
		PtrList< volVectorField > &	modes,
		word	fieldName,
		bool	podex,
		bool	supex,
		bool	sup,
		label	nmodes,
		bool	correctBC )

getNestedSnapshotMatrix() [1/3]

template<class Type, template< class > class PatchField, class GeoMesh>

void ITHACAPOD::getNestedSnapshotMatrix	(	PtrList< GeometricField< Type, PatchField, GeoMesh > > &	snapshots,
		PtrList< GeometricField< Type, PatchField, GeoMesh > > &	modes,
		word	fieldName,
		label	Npar,
		label	NnestedOut )

Nested-POD approach.

Computes the nested snapshot matrix and weighted bases for a vector field

Parameters

[in]	snapshots	List of snapshots.
	modes	The modes
[in]	fieldName	The field name
[in]	Npar	Number of parameters
[in]	NnestedOut	Number of weighted modes to be stored for the assembly of the nested matrix
[out]	UModesGlobal	List of snapshots where modes are stored (it must be passed empty).

Template Parameters

Type	vector or scalar.
PatchField	fvPatchField or fvsPatchField.
GeoMesh	volMesh or surfaceMesh.

Definition at line 41 of file ITHACAPOD.C.

getNestedSnapshotMatrix() [2/3]

template void ITHACAPOD::getNestedSnapshotMatrix	(	PtrList< volScalarField > &	snapshots,
		PtrList< volScalarField > &	ModesGlobal,
		word	fieldName,
		label	Npar,
		label	NnestedOut )

getNestedSnapshotMatrix() [3/3]

template void ITHACAPOD::getNestedSnapshotMatrix	(	PtrList< volVectorField > &	snapshots,
		PtrList< volVectorField > &	ModesGlobal,
		word	fieldName,
		label	Npar,
		label	NnestedOut )

getWeightedModes() [1/3]

template<class Type, template< class > class PatchField, class GeoMesh>

void ITHACAPOD::getWeightedModes	(	PtrList< GeometricField< Type, PatchField, GeoMesh > > &	snapshots,
		PtrList< GeometricField< Type, PatchField, GeoMesh > > &	modes,
		word	fieldName,
		bool	podex,
		bool	supex = 0,
		bool	sup = 0,
		label	nmodes = 0,
		bool	correctBC = true )

Computes the weighted bases (using the nested-pod approach) or read them for a vector field.

Parameters

[in]	snapshots	List of snapshots.
[out]	modes	A PtrList where modes are stored (it must be passed empty).
[in]	fieldName	The field name
[in]	podex	If 1, the functions read the stored mode. If 0, the function computes the modes and stores them.
[in]	supex	If 1, the functions read the stored supremizer mode. If 0, the function computes and stores them.
[in]	sup	If 1, it computes the supremizer modes.
[in]	nmodes	Number of modes to be stored. If 0, the maximum number of modes will computed.

Template Parameters

Type	vector or scalar.
PatchField	fvPatchField or fvsPatchField.
GeoMesh	volMesh or surfaceMesh.

Definition at line 687 of file ITHACAPOD.C.

getWeightedModes() [2/3]

template void ITHACAPOD::getWeightedModes	(	PtrList< volScalarField > &	snapshots,
		PtrList< volScalarField > &	modes,
		word	fieldName,
		bool	podex,
		bool	supex,
		bool	sup,
		label	nmodes,
		bool	correctBC )

getWeightedModes() [3/3]

template void ITHACAPOD::getWeightedModes	(	PtrList< volVectorField > &	snapshots,
		PtrList< volVectorField > &	modes,
		word	fieldName,
		bool	podex,
		bool	supex,
		bool	sup,
		label	nmodes,
		bool	correctBC )

GrammSchmidt()

void ITHACAPOD::GrammSchmidt

(

Eigen::MatrixXd &

Matrix

)

Performs GrammSchmidt orthonormalization on an Eigen Matrix.

Parameters

[out]

Matrix

The matrix

Definition at line 1324 of file ITHACAPOD.C.

weightedGramSchmidt()

void ITHACAPOD::weightedGramSchmidt	(	Eigen::MatrixXd &	matrix,
		Eigen::VectorXd &	weights )

Weighted Gram-Schmidt orthogonalization.

Parameters

[in]	matrix	Matrix to be orthogonalized
[in]	weights	Vector of weights

Definition at line 190 of file incrementalPOD.H.