#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

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

void constraints_print ( const Constraints & constraints )

Print the Auto Fitting Constraints Structure

Parameters

[in] constraints Constraints structure

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

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

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

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)

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

Functions

Function Documentation