model_auto.cpp File Reference

#include "geoslib_f.h"
#include "geoslib_old_f.h"
#include "Enum/EAnam.hpp"
#include "Enum/EJustify.hpp"
#include "Basic/AException.hpp"
#include "Basic/File.hpp"
#include "Basic/Utilities.hpp"
#include "Basic/String.hpp"
#include "Basic/OptDbg.hpp"
#include "Basic/VectorHelper.hpp"
#include "Covariances/CovAniso.hpp"
#include "Covariances/CovLMCTapering.hpp"
#include "Covariances/CovLMCConvolution.hpp"
#include "Covariances/CovLMCAnamorphosis.hpp"
#include "Model/Option_AutoFit.hpp"
#include "Model/Model.hpp"
#include "Model/Constraints.hpp"
#include "Model/ConsItem.hpp"
#include "Db/Db.hpp"
#include "Db/DbGrid.hpp"
#include "Variogram/Vario.hpp"
#include "Geometry/GeometryHelper.hpp"
#include "Matrix/MatrixSquareGeneral.hpp"
#include <math.h>

Functions
void	st_updateAlphaDiag (int icov0, int ivar0, int ncova, int nvar, VectorDouble &xr, VectorDouble &alpha, const VectorDouble &consSill)
int	model_auto_fit (Vario *vario, Model *model, bool verbose, const Option_AutoFit &mauto_arg, const Constraints &cons_arg, const Option_VarioFit &optvar_arg)
int	model_fitting_sills (Vario *vario, Model *model, const Constraints &constraints, const Option_AutoFit &mauto)
int	vmap_auto_fit (const DbGrid *dbmap, Model *model, bool verbose, const Option_AutoFit &mauto_arg, const Constraints &cons_arg, const Option_VarioFit &optvar_arg)
void	constraints_print (const Constraints &constraints)
int	modify_constraints_on_sill (Constraints &constraints)
double	constraints_get (const Constraints &constraints, const EConsType &icase, int igrf, int icov, const EConsElem &icons, int iv1, int iv2)

Function Documentation

constraints_get()

double constraints_get	(	const Constraints &	constraints,
		const EConsType &	icase,
		int	igrf,
		int	icov,
		const EConsElem &	icons,
		int	iv1,
		int	iv2
	)

Return the constraint value (if defined) or TEST

Returns

Returned value or TEST

Parameters

[in,out]	constraints	Constraints structure
[in]	icase	Parameter type (EConsType)
[in]	igrf	Rank of the Gaussian Random Function
[in]	icov	Rank of the structure (starting from 0)
[in]	icons	Type of the constraint (EConsElem)
[in]	iv1	Rank of the first variable
[in]	iv2	Rank of the second variable

constraints_print()

void constraints_print

(

const Constraints &

constraints

)

Print the Auto Fitting Constraints Structure

Parameters

[in]

constraints

Constraints structure

model_auto_fit()

int model_auto_fit	(	Vario *	vario,
		Model *	model,
		bool	verbose,
		const Option_AutoFit &	mauto_arg,
		const Constraints &	cons_arg,
		const Option_VarioFit &	optvar_arg
	)

Automatic model fitting

Returns

Error returned code

Parameters

[in]	vario	Vario structure containing the exp. variogram
[in]	model	Model structure containing the basic structures
[in]	verbose	Verbose flag
[in]	mauto_arg	Option_AutoFit structure
[in]	cons_arg	Constraints structure
[in]	optvar_arg	Opt_Vario structure

model_fitting_sills()

int model_fitting_sills	(	Vario *	vario,
		Model *	model,
		const Constraints &	constraints,
		const Option_AutoFit &	mauto
	)

Fitting a model using an experimental variogram

Returns

Error return code

Parameters

[in]	vario	Vario structure
[in,out]	model	Model to be fitted
[in]	constraints	Constraints structure
[in]	mauto	Option_AutoFit structure

modify_constraints_on_sill()

int modify_constraints_on_sill

(

Constraints &

constraints

)

If a constraint concerns a sill, take its square root as it corresponds to a constraints on AIC (not on a sill directly) due to the fact that it will be processed in FOXLEG (not in GOULARD) This transform only makes sense for MONOVARIATE case (the test should have been performed beforehand)

Returns

Error code (if the sill constraint is negative)

Parameters

[in]

constraints

Constraints structure

st_updateAlphaDiag()

void st_updateAlphaDiag	(	int	icov0,
		int	ivar0,
		int	ncova,
		int	nvar,
		VectorDouble &	xr,
		VectorDouble &	alpha,
		const VectorDouble &	consSill
	)

Update alpha from xr (diagonal only)

Parameters

[in]	icov0	Target basic structure
[in]	ivar0	Target variable
[in]	ncova	Number of basic structures
[in]	nvar	Number of variables
[in]	xr	Current vector of sqrt(constraint/(sum of the sills))
[in,out]	alpha	Current auxiliary matrices alpha
[in]	consSill	Vector of constant Sill (optional)

vmap_auto_fit()

int vmap_auto_fit	(	const DbGrid *	dbmap,
		Model *	model,
		bool	verbose,
		const Option_AutoFit &	mauto_arg,
		const Constraints &	cons_arg,
		const Option_VarioFit &	optvar_arg
	)

Automatic model fitting

Returns

Error returned code

Parameters

[in]	dbmap	Db Grid structure containing the Vmap
[in]	model	Model structure containing the basic structures
[in]	verbose	Verbose flag
[in]	mauto_arg	Option_AutoFit structure
[in]	cons_arg	Constraints structure
[in]	optvar_arg	Opt_Vario structure