CalcKriging.cpp File Reference

#include "Db/DbGrid.hpp"
#include "Db/Db.hpp"
#include "Estimation/CalcKriging.hpp"
#include "Enum/EKrigOpt.hpp"
#include "Estimation/CalcKrigingSimpleCase.hpp"
#include "Estimation/KrigOpt.hpp"
#include "Estimation/KrigingSystem.hpp"
#include "Basic/OptDbg.hpp"
#include "Basic/OptCustom.hpp"
#include "Model/Model.hpp"
#include "Neigh/NeighBench.hpp"
#include "Neigh/NeighUnique.hpp"
#include <math.h>

Functions
int	kriging (Db *dbin, Db *dbout, ModelGeneric *model, ANeigh *neigh, bool flag_est, bool flag_std, bool flag_varz, const KrigOpt &krigopt, const NamingConvention &namconv)
int	krigcell (Db *dbin, Db *dbout, ModelGeneric *model, ANeigh *neigh, bool flag_est, bool flag_std, const KrigOpt &krigopt, const NamingConvention &namconv)
int	kribayes (Db *dbin, Db *dbout, ModelGeneric *model, ANeigh *neigh, const VectorDouble &prior_mean, const MatrixSymmetric &prior_cov, bool flag_est, bool flag_std, const NamingConvention &namconv)
Krigtest_Res	krigtest (Db *dbin, Db *dbout, ModelGeneric *model, ANeigh *neigh, int iech0, const KrigOpt &krigopt, bool verbose)
int	kriggam (Db *dbin, Db *dbout, ModelGeneric *model, ANeigh *neigh, AAnam *anam, const NamingConvention &namconv)
int	xvalid (Db *db, ModelGeneric *model, ANeigh *neigh, bool flag_kfold, int flag_xvalid_est, int flag_xvalid_std, int flag_xvalid_varz, const KrigOpt &krigopt, const NamingConvention &namconv)
int	test_neigh (Db *dbin, Db *dbout, ModelGeneric *model, ANeigh *neigh, const NamingConvention &namconv)

Function Documentation

kribayes()

int kribayes	(	Db *	dbin,
		Db *	dbout,
		ModelGeneric *	model,
		ANeigh *	neigh,
		const VectorDouble &	prior_mean,
		const MatrixSymmetric &	prior_cov,
		bool	flag_est,
		bool	flag_std,
		const NamingConvention &	namconv
	)

Estimation with Bayesian Drift

Returns

Error return code

Parameters

[in]	dbin	input Db structure
[in]	dbout	output Db structure
[in]	model	ModelGeneric structure
[in]	neigh	ANeigh structure
[in]	prior_mean	Array giving the prior means for the drift terms
[in]	prior_cov	Array containing the prior covariance matrix for the drift terms
[in]	flag_est	Pointer for the storing the estimation
[in]	flag_std	Pointer for the storing the standard deviation
[in]	namconv	Naming convention

krigcell()

int krigcell	(	Db *	dbin,
		Db *	dbout,
		ModelGeneric *	model,
		ANeigh *	neigh,
		bool	flag_est,
		bool	flag_std,
		const KrigOpt &	krigopt,
		const NamingConvention &	namconv
	)

Standard Block Kriging with variable cell dimension

Returns

Error return code

Parameters

[in]	dbin	Input Db structure
[in]	dbout	Output Db structure
[in]	model	ModelGeneric structure
[in]	neigh	ANeigh structure
[in]	flag_est	Option for the storing the estimation
[in]	flag_std	Option for the storing the standard deviation
[in]	krigopt	KrigOpt structure
[in]	namconv	Naming convention

kriggam()

int kriggam	(	Db *	dbin,
		Db *	dbout,
		ModelGeneric *	model,
		ANeigh *	neigh,
		AAnam *	anam,
		const NamingConvention &	namconv
	)

Punctual Kriging in the Anamorphosed Gaussian Model

Returns

Error return code

Parameters

[in]	dbin	input Db structure
[in]	dbout	output Db structure
[in]	model	ModelGeneric structure
[in]	neigh	ANeigh structure
[in]	anam	AAnam structure
[in]	namconv	Naming convention

kriging()

int kriging	(	Db *	dbin,
		Db *	dbout,
		ModelGeneric *	model,
		ANeigh *	neigh,
		bool	flag_est,
		bool	flag_std,
		bool	flag_varz,
		const KrigOpt &	krigopt,
		const NamingConvention &	namconv
	)

Standard Kriging

Returns

Error return code

Parameters

[in]	dbin	Input Db structure
[in]	dbout	Output Db structure
[in]	model	ModelGeneric structure
[in]	neigh	ANeigh structure
[in]	flag_est	Option for storing the estimation
[in]	flag_std	Option for storing the standard deviation
[in]	flag_varz	Option for storing the variance of the estimator (only available for stationary model)
[in]	krigopt	KrigOpt structure
[in]	namconv	Naming convention

krigtest()

Krigtest_Res krigtest	(	Db *	dbin,
		Db *	dbout,
		ModelGeneric *	model,
		ANeigh *	neigh,
		int	iech0,
		const KrigOpt &	krigopt,
		bool	verbose
	)

Perform kriging and return the calculation elements

Returns

A Krigtest_Res structure

Parameters

[in]	dbin	input Db structure
[in]	dbout	output Db structure
[in]	model	ModelGeneric structure
[in]	neigh	ANeigh structure
[in]	iech0	Rank of the target sample
[in]	krigopt	KrigOpt structure
[in]	verbose	When TRUE, the full debugging flag is switched ON (the current status is reset after the run)

test_neigh()

int test_neigh	(	Db *	dbin,
		Db *	dbout,
		ModelGeneric *	model,
		ANeigh *	neigh,
		const NamingConvention &	namconv
	)

Check the Neighborhood

Returns

Error return code (0: success, 1: error)

Parameters

[in]	dbin	input Db structure
[in]	dbout	output Db structure
[in]	model	ModelGeneric structure (optional)
[in]	neigh	ANeigh structure
[in]	namconv	Naming Convention

Remarks

This procedure creates the following arrays:

1 - The number of selected samples

2 - The maximum neighborhood distance

3 - The minimum neighborhood distance

4 - The number of non-empty sectors

5 - The number of consecutive empty sectors

xvalid()

int xvalid	(	Db *	db,
		ModelGeneric *	model,
		ANeigh *	neigh,
		bool	flag_kfold,
		int	flag_xvalid_est,
		int	flag_xvalid_std,
		int	flag_xvalid_varz,
		const KrigOpt &	krigopt,
		const NamingConvention &	namconv
	)

Standard Cross-Validation

Parameters

db	Db structure
model	ModelGeneric structure
neigh	ANeigh structure
flag_kfold	True if a code (K-FOLD) is used
flag_xvalid_est	Option for storing the estimation: 1 for Z*-Z; -1 for Z*; 0 not stored
flag_xvalid_std	Option for storing the standard deviation: 1:for (Z*-Z)/S; -1 for S; 0 not stored
flag_xvalid_varz	Option for storing the variance of the estimator: 1 to store and 0 not stored
krigopt	KrigOpt structure
namconv	Naming Convention

Returns

Error return code