#include "Db/DbGrid.hpp"
#include "Db/Db.hpp"
#include "Estimation/CalcKriging.hpp"
#include "Estimation/KrigingSystem.hpp"
#include "Basic/OptDbg.hpp"
#include <math.h>

Functions
int	kriging (Db dbin, Db dbout, Model model, ANeigh neigh, const EKrigOpt &calcul, bool flag_est, bool flag_std, bool flag_varz, const VectorInt &ndiscs, const VectorInt &rank_colcok, const MatrixRectangular *matLC, const NamingConvention &namconv)

int	krigcell (Db dbin, Db dbout, Model model, ANeigh neigh, bool flag_est, bool flag_std, const VectorInt &ndiscs, const VectorInt &rank_colcok, const NamingConvention &namconv)

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

int	krigprof (Db dbin, Db dbout, Model model, ANeigh neigh, bool flag_est, bool flag_std, const NamingConvention &namconv)

Krigtest_Res	krigtest (Db dbin, Db dbout, Model model, ANeigh neigh, int iech0, const EKrigOpt &calcul, const VectorInt &ndiscs, bool flagPerCell, bool verbose)

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

int	xvalid (Db db, Model model, ANeigh *neigh, bool flag_kfold, int flag_xvalid_est, int flag_xvalid_std, int flag_xvalid_varz, const VectorInt &rank_colcok, const NamingConvention &namconv)

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

Function Documentation

int kribayes	(	Db *	dbin,
		Db *	dbout,
		Model *	model,
		ANeigh *	neigh,
		const VectorDouble &	prior_mean,
		const MatrixSquareSymmetric &	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	Model 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

int krigcell	(	Db *	dbin,
		Db *	dbout,
		Model *	model,
		ANeigh *	neigh,
		bool	flag_est,
		bool	flag_std,
		const VectorInt &	ndiscs,
		const VectorInt &	rank_colcok,
		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	Model structure
[in]	neigh	ANeigh structure
[in]	ndiscs	Array giving the discretization counts
[in]	flag_est	Option for the storing the estimation
[in]	flag_std	Option for the storing the standard deviation
[in]	rank_colcok	Option for running Collocated Cokriging
[in]	namconv	Naming convention

int kriggam	(	Db *	dbin,
		Db *	dbout,
		Model *	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	Model structure
[in]	neigh	ANeigh structure
[in]	anam	AAnam structure
[in]	namconv	Naming convention

int kriging	(	Db *	dbin,
		Db *	dbout,
		Model *	model,
		ANeigh *	neigh,
		const EKrigOpt &	calcul,
		bool	flag_est,
		bool	flag_std,
		bool	flag_varz,
		const VectorInt &	ndiscs,
		const VectorInt &	rank_colcok,
		const MatrixRectangular *	matLC,
		const NamingConvention &	namconv
	)

Standard Kriging

Returns: Error return code

Parameters

[in]	dbin	Input Db structure
[in]	dbout	Output Db structure
[in]	model	Model structure
[in]	neigh	ANeigh structure
[in]	calcul	Kriging calculation option (EKrigOpt)
[in]	ndiscs	Array giving the discretization counts
[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]	rank_colcok	Option for running Collocated Cokriging
[in]	matLC	Matrix of linear combination (or NULL) (Dimension: nvarLC * model->getNVar())
[in]	namconv	Naming convention

int krigprof	(	Db *	dbin,
		Db *	dbout,
		Model *	model,
		ANeigh *	neigh,
		bool	flag_est,
		bool	flag_std,
		const NamingConvention &	namconv
	)

Punctual Kriging based on profiles

Returns: Error return code

Parameters

[in]	dbin	input Db structure
[in]	dbout	output Db structure
[in]	model	Model 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]	namconv	Naming convention

Krigtest_Res krigtest	(	Db *	dbin,
		Db *	dbout,
		Model *	model,
		ANeigh *	neigh,
		int	iech0,
		const EKrigOpt &	calcul,
		const VectorInt &	ndiscs,
		bool	flagPerCell,
		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	Model structure
[in]	neigh	ANeigh structure
[in]	iech0	Rank of the target sample
[in]	calcul	Kriging calculation option (EKrigOpt)
[in]	ndiscs	Array giving the discretization counts
[in]	flagPerCell	Use local block extensions (when defined)
[in]	verbose	When TRUE, the full debugging flag is switched ON (the current status is reset after the run)

int test_neigh	(	Db *	dbin,
		Db *	dbout,
		Model *	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	Model 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

int xvalid	(	Db *	db,
		Model *	model,
		ANeigh *	neigh,
		bool	flag_kfold,
		int	flag_xvalid_est,
		int	flag_xvalid_std,
		int	flag_xvalid_varz,
		const VectorInt &	rank_colcok,
		const NamingConvention &	namconv
	)

Standard Cross-Validation

Parameters

db	Db structure
model	Model 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
flag_xvalid_std	Option for storing the standard deviation: 1:for (Z*-Z)/S; -1 for S
flag_xvalid_varz	Option for storing the variance of the estimator
rank_colcok	Option for running Collocated Cokriging
namconv	Naming Convention

Returns: Error return code

Functions

Function Documentation