ShiftOpCs Class Reference

Shift Operator for performing the basic tasks of SPDE. More...

#include <ShiftOpCs.hpp>

Inheritance diagram for ShiftOpCs:

Detailed Description

Shift Operator for performing the basic tasks of SPDE.

Public Member Functions
	ShiftOpCs (const CGParam params=CGParam())
	ShiftOpCs (const AMesh *amesh, Model *model, const Db *dbout=nullptr, int igrf=0, int icov=0, const CGParam params=CGParam(), bool verbose=false)
	ShiftOpCs (const MatrixSparse *S, const VectorDouble &TildeC, const VectorDouble &Lambda, Model *model, const CGParam params=CGParam(), bool verbose=false)
	ShiftOpCs (const ShiftOpCs &shift)
ShiftOpCs &	operator= (const ShiftOpCs &shift)
virtual	~ShiftOpCs ()
void	_evalDirect (const VectorDouble &inv, VectorDouble &outv) const override
int	initFromMesh (const AMesh *amesh, Model *model, const Db *dbout=nullptr, int igrf=0, int icov=0, bool flagAdvection=false, bool verbose=false)
int	initGradFromMesh (const AMesh *amesh, Model *model, int igrf=0, int icov=0, bool verbose=false, double tol=EPSILON10)
int	initFromCS (const MatrixSparse *S, const VectorDouble &TildeC, const VectorDouble &Lambda, Model *model, bool verbose=false)
int	getSize () const override
int	getNDim () const
int	getNModelGradParam () const
void	prodTildeC (const VectorDouble &x, VectorDouble &y, const EPowerPT &power) const
void	prodLambda (const VectorDouble &x, VectorDouble &y, const EPowerPT &power) const
void	prodLambdaOnSqrtTildeC (const VectorDouble &inv, VectorDouble &outv, double puis=2) const
double	getMaxEigenValue () const
MatrixSparse *	getS () const
MatrixSparse *	getTildeCGrad (int iapex, int igparam) const
MatrixSparse *	getSGrad (int iapex, int igparam) const
const VectorDouble &	getTildeC () const
const VectorDouble &	getLambdas () const
double	getLambda (int iapex) const
const VectorDouble &	getLambdaGrads (int idim) const
double	getLambdaGrad (int idim, int iapex) const
int	getSGradAddress (int iapex, int igparam) const
int	getLambdaGradSize () const
Public Member Functions inherited from ALinearOp
	ALinearOp (const CGParam params=CGParam())
	ALinearOp (const ALinearOp &m)
ALinearOp &	operator= (const ALinearOp &m)
virtual	~ALinearOp ()
virtual void	evalInverse (const VectorDouble &inv, VectorDouble &outv) const
void	evalDirect (const VectorDouble &inv, VectorDouble &outv) const
void	setX0 (const VectorDouble &x0)
void	mustShowStats (bool status)
const LogStats &	getLogStats () const

Static Public Member Functions
static ShiftOpCs *	create (const AMesh *amesh, Model *model, const Db *dbout=nullptr, int igrf=0, int icov=0, const CGParam params=CGParam(), bool verbose=false)
static ShiftOpCs *	createFromSparse (const MatrixSparse *S, const VectorDouble &TildeC, const VectorDouble &Lambda, Model *model, const CGParam params=CGParam(), bool verbose=false)

Constructor & Destructor Documentation

ShiftOpCs() [1/4]

ShiftOpCs::ShiftOpCs

(

const CGParam

params = CGParam()

)

ShiftOpCs() [2/4]

ShiftOpCs::ShiftOpCs	(	const AMesh *	amesh,
		Model *	model,
		const Db *	dbout = nullptr,
		int	igrf = 0,
		int	icov = 0,
		const CGParam	params = CGParam(),
		bool	verbose = false
	)

ShiftOpCs() [3/4]

ShiftOpCs::ShiftOpCs	(	const MatrixSparse *	S,
		const VectorDouble &	TildeC,
		const VectorDouble &	Lambda,
		Model *	model,
		const CGParam	params = CGParam(),
		bool	verbose = false
	)

ShiftOpCs() [4/4]

ShiftOpCs::ShiftOpCs

(

const ShiftOpCs &

shift

)

~ShiftOpCs()

ShiftOpCs::~ShiftOpCs ( )

virtual

Member Function Documentation

_evalDirect()

void ShiftOpCs::_evalDirect ( const VectorDouble &  x, VectorDouble &  y  ) const

overridevirtual

Perform the operation: y = S * x

Parameters

[in]	x	Input vector
[in]	y	Output vector

Remarks

'S' is a member that stands as a sparse matrix

Implements ALinearOp.

create()

ShiftOpCs * ShiftOpCs::create ( const AMesh *  amesh, Model *  model, const Db *  dbout = nullptr, int  igrf = 0, int  icov = 0, const CGParam  params = CGParam(), bool  verbose = false  )

static

createFromSparse()

ShiftOpCs * ShiftOpCs::createFromSparse ( const MatrixSparse *  S, const VectorDouble &  TildeC, const VectorDouble &  Lambda, Model *  model, const CGParam  params = CGParam(), bool  verbose = false  )

static

getLambda()

double ShiftOpCs::getLambda ( int  iapex ) const

inline

getLambdaGrad()

double ShiftOpCs::getLambdaGrad ( int  idim, int  iapex  ) const

inline

getLambdaGrads()

const VectorDouble& ShiftOpCs::getLambdaGrads ( int  idim ) const

inline

getLambdaGradSize()

int ShiftOpCs::getLambdaGradSize

(

)

const

getLambdas()

const VectorDouble& ShiftOpCs::getLambdas ( ) const

inline

getMaxEigenValue()

double ShiftOpCs::getMaxEigenValue

(

)

const

getNDim()

int ShiftOpCs::getNDim ( ) const

inline

getNModelGradParam()

int ShiftOpCs::getNModelGradParam ( ) const

inline

getS()

MatrixSparse* ShiftOpCs::getS ( ) const

inline

getSGrad()

MatrixSparse * ShiftOpCs::getSGrad	(	int	iapex,
		int	igparam
	)		const

getSGradAddress()

int ShiftOpCs::getSGradAddress	(	int	iapex,
		int	igparam
	)		const

Returns the internal address for a given vertex and a given parameter It returns -1 if the address is ivalid

Parameters

iapex	Rank of the target apex
igparam	Rank of the target parameter

Returns

getSize()

int ShiftOpCs::getSize ( ) const

inlineoverridevirtual

Implements ALinearOp.

getTildeC()

const VectorDouble& ShiftOpCs::getTildeC ( ) const

inline

getTildeCGrad()

MatrixSparse * ShiftOpCs::getTildeCGrad	(	int	iapex,
		int	igparam
	)		const

initFromCS()

int ShiftOpCs::initFromCS	(	const MatrixSparse *	S,
		const VectorDouble &	TildeC,
		const VectorDouble &	Lambda,
		Model *	model,
		bool	verbose = false
	)

Parameters

S	Sparse matrix describing the S information
TildeC	Diagonal array containing TildeC
Lambda	Normalization vector
model	Pointer to the Model structure
verbose	Verbose flag

Returns

initFromMesh()

int ShiftOpCs::initFromMesh	(	const AMesh *	amesh,
		Model *	model,
		const Db *	dbout = nullptr,
		int	igrf = 0,
		int	icov = 0,
		bool	flagAdvection = false,
		bool	verbose = false
	)

Parameters

amesh	Meshing description (New format)
model	Pointer to the Model structure
dbout	Pointer to the Db structure
igrf	Rank of the GRF
icov	Rank of the Covariance within the Model
flagAdvection	When TRUE, S is replaced by G
verbose	Verbose flag

Returns

Error return code

initGradFromMesh()

int ShiftOpCs::initGradFromMesh	(	const AMesh *	amesh,
		Model *	model,
		int	igrf = 0,
		int	icov = 0,
		bool	verbose = false,
		double	tol = EPSILON10
	)

Initialize the environment for calculation of derivatives of S

Parameters

amesh	Meshing description (New format)
model	Pointer to the Model structure
igrf	Rank of the GRF
icov	Rank of the Covariance within the Model
verbose	Verbose flag
tol	Smallest value below which the value is not stored in sparse matrix

Returns

Error return code

operator=()

ShiftOpCs & ShiftOpCs::operator=

(

const ShiftOpCs &

shift

)

prodLambda()

void ShiftOpCs::prodLambda	(	const VectorDouble &	x,
		VectorDouble &	y,
		const EPowerPT &	power
	)		const

prodLambdaOnSqrtTildeC()

void ShiftOpCs::prodLambdaOnSqrtTildeC	(	const VectorDouble &	inv,
		VectorDouble &	outv,
		double	puis = 2
	)		const

prodTildeC()

void ShiftOpCs::prodTildeC	(	const VectorDouble &	x,
		VectorDouble &	y,
		const EPowerPT &	power
	)		const

Perform the operation: y = x * C^power

Parameters

[in]	x	Input vector
[in]	y	Output vector
[in]	power	Value of the exponent

Remarks

'C' is a member (_TildeC) that stands as a vector

Specific coding has been realized for the cases

where 'power' is equal to 1, -1, 0.5 and -0.5

---

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

• include/LinearOp/ShiftOpCs.hpp
• src/LinearOp/ShiftOpCs.cpp