SPDE.cpp File Reference

#include "Enum/ECov.hpp"
#include "API/SPDE.hpp"
#include "Model/ANoStat.hpp"
#include "Covariances/CovAniso.hpp"
#include "Mesh/MeshETurbo.hpp"
#include "Basic/AException.hpp"
#include "Basic/Law.hpp"
#include "Basic/VectorHelper.hpp"
#include "Basic/NamingConvention.hpp"
#include "Model/Model.hpp"
#include "LinearOp/ShiftOpCs.hpp"
#include "LinearOp/PrecisionOp.hpp"
#include "LinearOp/PrecisionOpMultiConditional.hpp"
#include "LinearOp/PrecisionOpMultiConditionalCs.hpp"
#include "LinearOp/ProjMatrix.hpp"
#include "Db/Db.hpp"
#include "Db/DbGrid.hpp"
#include <iostream>
#include <math.h>

Functions
int	krigingSPDE (Db *dbin, Db *dbout, Model *model, bool flag_est, bool flag_std, bool flag_varz, const AMesh *mesh, int useCholesky, SPDEParam params, bool verbose, bool showStats, const NamingConvention &namconv)
int	simulateSPDE (Db *dbin, Db *dbout, Model *model, int nbsimu, const AMesh *mesh, int useCholesky, SPDEParam params, int seed, bool verbose, bool showStats, const NamingConvention &namconv)

Function Documentation

krigingSPDE()

int krigingSPDE	(	Db *	dbin,
		Db *	dbout,
		Model *	model,
		bool	flag_est,
		bool	flag_std,
		bool	flag_varz,
		const AMesh *	mesh,
		int	useCholesky,
		SPDEParam	params,
		bool	verbose,
		bool	showStats,
		const NamingConvention &	namconv
	)

Perform the estimation by KRIGING under the SPDE framework

Parameters

dbin	Input Db (must contain the variable to be estimated)
dbout	Output Db where the estimation must be performed
model	Model definition
flag_est	True for the estimation
flag_std	True for the standard deviation of estimation error
flag_varz	True for the variance of the estimator
mesh	Mesh description (optional)
useCholesky	Define the choice regarding Cholesky
params	Set of parameters
verbose	Verbose flag
showStats	Show statistics for Linear Operations
namconv	Naming convention

Returns

Error return code

Remarks

You can provide an already existing mesh. Otherwise an optimal mesh will be created

internally: one per structure constituting the Model for Kriging.

Each mesh is created using the Turbo Meshing facility... based on an internal grid.

This internal grid is rotated according to the rotation of the structure. Its mesh size

is derived from the range (per direction) by dividing it by the refinement factor.

simulateSPDE()

int simulateSPDE	(	Db *	dbin,
		Db *	dbout,
		Model *	model,
		int	nbsimu,
		const AMesh *	mesh,
		int	useCholesky,
		SPDEParam	params,
		int	seed,
		bool	verbose,
		bool	showStats,
		const NamingConvention &	namconv
	)

Perform simulations under the SPDE framework

Parameters

dbin	Input Db. If defined, the simulations are conditional; non conditional otherwise
dbout	Output Db where the simulations must be performed
model	Model definition
nbsimu	Number of simulations
mesh	Mesh description (optional)
useCholesky	Define the choice regarding Cholesky
params	Set of parametes
seed	Seed used for the Random Number generator
verbose	Verbose flag
showStats	Show statistics for Linear Operations
namconv	Naming convention

Returns

Error return code

Remarks

You can provide an already existing mesh. Otherwise an optimal mesh will be created

internally: one per structure constituting the Model for Kriging and one for Simulating

Each mesh is created using the Turbo Meshing facility... based on an internal grid.

This internal grid is rotated according to the rotation of the structure. Its mesh size

is derived from the range (per direction) by dividing it by the refinement factor.