#include <KrigingSystem.hpp>

Public Member Functions
	KrigingSystem (Db dbin, Db dbout, const Model model, ANeigh neigh)

	KrigingSystem (const KrigingSystem &m)=delete

KrigingSystem &	operator= (const KrigingSystem &m)=delete

virtual	~KrigingSystem ()

int	setKrigOptCalcul (const EKrigOpt &calcul, const VectorInt &ndiscs=VectorInt(), bool flag_per_cell=false)

int	setKrigOptXValid (bool flag_xvalid, bool flag_kfold, bool optionXValidEstim=false, bool optionXValidStdev=false, bool optionXValidVarZ=false)

int	setKrigOptColCok (const VectorInt &rank_colcok)

int	setKrigOptBayes (bool flag_bayes, const VectorDouble &prior_mean, const MatrixSquareSymmetric &prior_cov, int seed=414371)

int	setKrigOptImage (int seed=133271)

int	setKrigOptDataWeights (int iptrWeights, bool flagSet=true)

int	setKrigOptMatLC (const MatrixRectangular *matLC)

int	setKrigoptCode (bool flag_code)

int	setKrigOptFlagSimu (bool flagSimu, int nbsimu=0, int rankPGS=-1)

int	setKrigOptSaveWeights (bool flag_save)

int	setKrigOptDGM (bool flag_dgm, double eps=EPSILON6)

int	setKrigOptFlagGlobal (bool flag_global)

int	setKrigOptFlagLTerm (bool flag_lterm)

int	setKrigOptAnamophosis (AAnam *anam)

int	setKrigOptFactorKriging (bool flag_factor_kriging)

int	updKrigOptEstim (int iptrEst, int iptrStd, int iptrVarZ)

int	updKrigOptIclass (int index_class, int nclasses)

int	updKrigOptNeighOnly (int iptrNeigh)

bool	isReady ()

int	estimate (int iech_out)

void	conclusion ()

int	getNDim () const

int	getNech () const

int	getNeq () const

int	getNRed () const

VectorInt	getSampleIndices () const

VectorVectorDouble	getSampleCoordinates () const

VectorDouble	getSampleData () const

MatrixRectangular	getZam () const

MatrixSquareSymmetric	getLHSC () const

MatrixRectangular	getRHSC () const

MatrixRectangular	getWeights () const

MatrixSquareGeneral	getVariance () const

double	getLTerm () const

VectorDouble	getRHSC (int ivar) const

VectorDouble	getZamC () const

VectorDouble	getLHSInvC () const

Constructor & Destructor Documentation

◆ KrigingSystem() [1/2]

KrigingSystem::KrigingSystem	(	Db *	dbin,
		Db *	dbout,
		const Model *	model,
		ANeigh *	neigh
	)

◆ KrigingSystem() [2/2]

KrigingSystem::KrigingSystem ( const KrigingSystem & m )

delete

◆ ~KrigingSystem()

KrigingSystem::~KrigingSystem ( )

virtual

Member Function Documentation

◆ conclusion()

void KrigingSystem::conclusion ( )

This method closes the use of a KrigingSystem sequence

◆ estimate()

int KrigingSystem::estimate ( int iech_out )

Perform the Kriging of target

Parameters

iech_out Rank of the target

Returns

◆ getLHSC()

MatrixSquareSymmetric KrigingSystem::getLHSC ( ) const

inline

◆ getLHSInvC()

VectorDouble KrigingSystem::getLHSInvC ( ) const

inline

◆ getLTerm()

double KrigingSystem::getLTerm ( ) const

inline

◆ getNDim()

int KrigingSystem::getNDim ( ) const

◆ getNech()

int KrigingSystem::getNech ( ) const

Returns the number of samples in the neighborhood

Returns: Number of samples

Remarks: If collocated option is switch ON, samples are alreadycount (rank < 0)

◆ getNeq()

int KrigingSystem::getNeq ( ) const

◆ getNRed()

int KrigingSystem::getNRed ( ) const

inline

◆ getRHSC() [1/2]

MatrixRectangular KrigingSystem::getRHSC ( ) const

inline

◆ getRHSC() [2/2]

VectorDouble KrigingSystem::getRHSC ( int ivar ) const

◆ getSampleCoordinates()

VectorVectorDouble KrigingSystem::getSampleCoordinates ( ) const

Returns the coordinates of the neighboring samples

Returns: Array organized by Coordinate (minor) then by Sample (major)

◆ getSampleData()

VectorDouble KrigingSystem::getSampleData ( ) const

◆ getSampleIndices()

VectorInt KrigingSystem::getSampleIndices ( ) const

inline

◆ getVariance()

MatrixSquareGeneral KrigingSystem::getVariance ( ) const

inline

◆ getWeights()

MatrixRectangular KrigingSystem::getWeights ( ) const

inline

◆ getZam()

MatrixRectangular KrigingSystem::getZam ( ) const

inline

◆ getZamC()

VectorDouble KrigingSystem::getZamC ( ) const

◆ isReady()

bool KrigingSystem::isReady ( )

Performs the last operations before launching the loop on Estimations

Returns

◆ operator=()

KrigingSystem& KrigingSystem::operator= ( const KrigingSystem & m )

delete

◆ setKrigOptAnamophosis()

int KrigingSystem::setKrigOptAnamophosis ( AAnam * anam )

Perform Gaussian Anamoprhosis kriging

Parameters

anam	Pointer to the AAnam structure

Returns

◆ setKrigOptBayes()

int KrigingSystem::setKrigOptBayes	(	bool	flag_bayes,
		const VectorDouble &	prior_mean,
		const MatrixSquareSymmetric &	prior_cov,
		int	seed = `414371`
	)

◆ setKrigOptCalcul()

int KrigingSystem::setKrigOptCalcul	(	const EKrigOpt &	calcul,
		const VectorInt &	ndiscs = `VectorInt()`,
		bool	flag_per_cell = `false`
	)

◆ setKrigoptCode()

int KrigingSystem::setKrigoptCode ( bool flag_code )

Define the option for Kriging By Profile

Parameters

flag_code True if the option is switched ON

Returns

Remarks: When the Option for Kriging By Profile is ON, the covariance term in LHS; is incremented by V(iech) when both samples 'iech' and 'jech' have; the same Code.

◆ setKrigOptColCok()

int KrigingSystem::setKrigOptColCok ( const VectorInt & rank_colcok )

Check the consistency of the Colocation specification

Returns: Error return code

Parameters

[in] rank_colcok Array of ranks of colocated variables

Remarks: The array 'rank_colcok' (if present) must be dimensioned; to the number of variables in Dbin.; Each element gives the rank of the colocated variable within Dbout; or -1 if not colocated; If the array 'rank_colcok' is absent, colocation option is OFF.; In input, the numbering in ; rank_colcok' starts from 1; In output, the numbering starts from 0

◆ setKrigOptDataWeights()

int KrigingSystem::setKrigOptDataWeights	(	int	iptrWeights,
		bool	flagSet = `true`
	)

◆ setKrigOptDGM()

int KrigingSystem::setKrigOptDGM	(	bool	flag_dgm,
		double	eps = `EPSILON6`
	)

◆ setKrigOptFactorKriging()

int KrigingSystem::setKrigOptFactorKriging ( bool flag_factor_kriging )

◆ setKrigOptFlagGlobal()

int KrigingSystem::setKrigOptFlagGlobal ( bool flag_global )

◆ setKrigOptFlagLTerm()

int KrigingSystem::setKrigOptFlagLTerm ( bool flag_lterm )

Ask for the specific calculation of Z * A-1 * Z

Parameters

flag_lterm Flag for asking this specific calculation

Returns

Remarks: The calculated value can be retrieved using _getLTerm() method

◆ setKrigOptFlagSimu()

int KrigingSystem::setKrigOptFlagSimu	(	bool	flagSimu,
		int	nbsimu = `0`,
		int	rankPGS = `-1`
	)

◆ setKrigOptImage()

int KrigingSystem::setKrigOptImage ( int seed = 133271 )

◆ setKrigOptMatLC()

int KrigingSystem::setKrigOptMatLC ( const MatrixRectangular * matLC )

Define the output as Linear Combinations of the Input Variables

Parameters

matLC Vector of Vectors of weights (see remarks)

Returns

Remarks: The first dimension of 'matLC' is the number of Output variables; The second dimension is the number of input Variables.

◆ setKrigOptSaveWeights()

int KrigingSystem::setKrigOptSaveWeights ( bool flag_save )

Switch the option for saving the Kriging Weights using Keypair mechanism

Parameters

flag_save Value of the switch

◆ setKrigOptXValid()

int KrigingSystem::setKrigOptXValid	(	bool	flag_xvalid,
		bool	flag_kfold,
		bool	optionXValidEstim = `false`,
		bool	optionXValidStdev = `false`,
		bool	optionXValidVarZ = `false`
	)

Set the flag for performing Cross-Validation

Parameters

flag_xvalid	True if the Cross-Validation option is switched ON
flag_kfold	True if the KFold option is switch ON
optionXValidEstim	True for Z-Z; False for Z
optionXValidStdev	True for (Z*-Z)/S; False for S
optionXValidVarZ	True for Var(Z*)

Returns

Remarks: The KFold option requires a Code to be assigned to each Datum; In the Neighborhood search for a Target sample, some other data points; are discarded:; - either the Target sample itself (Leave-One-Point-Out) if KFold is False; - all samples with same code as Target if KFold is True

◆ updKrigOptEstim()

int KrigingSystem::updKrigOptEstim	(	int	iptrEst,
		int	iptrStd,
		int	iptrVarZ
	)

Set the calculation options

Parameters

iptrEst	UID for storing the estimation(s)
iptrStd	UID for storing the Standard deviations(s)
iptrVarZ	UID for storing the Variance(s) of estimator

Returns

Remarks: If a term must not be calculated, its UID must be negative

◆ updKrigOptIclass()

int KrigingSystem::updKrigOptIclass	(	int	index_class,
		int	nclasses
	)

◆ updKrigOptNeighOnly()

int KrigingSystem::updKrigOptNeighOnly ( int iptrNeigh )

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

include/Estimation/KrigingSystem.hpp
src/Estimation/KrigingSystem.cpp

Public Member Functions

Constructor & Destructor Documentation

◆ KrigingSystem() [1/2]

◆ KrigingSystem() [2/2]

◆ ~KrigingSystem()

Member Function Documentation

◆ conclusion()

◆ estimate()

◆ getLHSC()

◆ getLHSInvC()

◆ getLTerm()

◆ getNDim()

◆ getNech()

◆ getNeq()

◆ getNRed()

◆ getRHSC() [1/2]

◆ getRHSC() [2/2]

◆ getSampleCoordinates()

◆ getSampleData()

◆ getSampleIndices()

◆ getVariance()

◆ getWeights()

◆ getZam()

◆ getZamC()

◆ isReady()

◆ operator=()

◆ setKrigOptAnamophosis()

◆ setKrigOptBayes()

◆ setKrigOptCalcul()

◆ setKrigoptCode()

◆ setKrigOptColCok()

◆ setKrigOptDataWeights()

◆ setKrigOptDGM()

◆ setKrigOptFactorKriging()

◆ setKrigOptFlagGlobal()

◆ setKrigOptFlagLTerm()

◆ setKrigOptFlagSimu()

◆ setKrigOptImage()

◆ setKrigOptMatLC()

◆ setKrigOptSaveWeights()

◆ setKrigOptXValid()

◆ updKrigOptEstim()

◆ updKrigOptIclass()

◆ updKrigOptNeighOnly()