#include "geoslib_old_f.h"
#include "geoslib_f.h"
#include "Enum/EAnam.hpp"
#include "Basic/AException.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
static void	st_modify_optvar_for_anam (Model *model, Option_VarioFit &optvar)

static void	st_name_range (int ivar)

static void	st_name_scale (int ivar)

static void	st_name_rotation (int rank)

static void	st_parid_decode (int parid, int imod, int icov, EConsElem icons, int ivar, int *jvar)

static int	st_parid_encode (int imod, int icov, const EConsElem &icons, int ivar, int jvar)

static int	st_parid_alloc (StrMod *strmod, int npar0)

static StrMod *	st_model_auto_strmod_free (StrMod *strmod)

static StrMod *	st_model_auto_strmod_alloc (Model model1, Model model2, int npar0, int norder, double hmax, VectorDouble &angles, const Option_VarioFit &optvar, int *npar)

static int	st_get_vario_dimension (Vario vario, int nbexp_ret, int *npadir_ret)

static int	st_get_vmap_dimension (const Db dbmap, int nvar, int nbexp_ret, int *npadir_ret)

static void	st_mauto_rescale (int nvar, VectorDouble &varchol, Option_AutoFit &mauto)

static void	st_goulard_verbose (int mode, Option_AutoFit &mauto)

static std::vector< StrExp >	st_strexp_manage (int nbexp, int ndim)

static void	st_compress_array (const Vario *vario, int npadir, VectorDouble &tabin, VectorDouble &tabout)

static double	st_get_c00 (const Vario *vario, int idir, int ivar, int jvar)

static void	st_load_gg (const Vario *vario, int npadir, std::vector< StrExp > &strexps, VectorDouble &gg)

static void	st_prepar_goulard_vario (int imod)

static void	st_load_ge (const Vario vario, Model model, int npadir, VectorDouble &dd, VectorDouble &ge)

static void	st_load_wt (const Vario *vario, int wmode, int npadir, VectorDouble &wt)

static void	st_goulard_score (const Option_AutoFit &mauto, int mode, int ncova, int iter, double crit)

static void	st_goulard_debug_title (int nvar, int ncova)

static void	st_goulard_debug_current (int nvar, int ncova, int iter, VectorDouble &sill, double crit)

static void	st_keypair_sill (int mode, Model *model)

static void	st_keypair_results (int mode, int icov, int nvar, double valpro, double vecpro)

static void	st_sill_reset (int nvar, int ncova, VectorDouble &sill)

static int	st_goulard_without_constraint (const Option_AutoFit &mauto, int nvar, int ncova, int npadir, VectorDouble &wt, VectorDouble &gg, VectorDouble &ge, VectorDouble &sill, double *crit_arg)

static void	st_affect (int rank, double def_val, double lower_val, double upper_val, VectorDouble &param, VectorDouble &lower, VectorDouble &upper)

static int	st_compress_parid (int n_init, VectorInt &parid, VectorDouble &param, VectorDouble &lower, VectorDouble &upper)

static void	st_print (const char *name, int flag_end, int rank, VectorDouble &param, VectorDouble &lower, VectorDouble &upper)

static void	st_model_auto_strmod_print (int flag_title, StrMod *strmod, const Constraints &constraints, const Option_AutoFit &mauto, VectorDouble &param, VectorDouble &lower, VectorDouble &upper, int npar, int nbexp)

static void	st_model_auto_scldef (StrMod *strmod, int npar, double hmax, VectorDouble &varchol, VectorDouble &scale)

static void	st_model_auto_constraints_apply (StrMod *strmod, int npar, const Constraints &constraints, VectorDouble &param, VectorDouble &lower, VectorDouble &upper)

static void	st_model_auto_pardef (StrMod *strmod, int npar, double hmax, VectorDouble &varchol, VectorDouble &angles, VectorDouble &param, VectorDouble &lower, VectorDouble &upper)

static void	st_model_auto_strmod_define (StrMod *strmod, int npar, VectorDouble &param)

static int	st_structure_reduce (StrMod *strmod, int imod, int icov, double hmax, double gmax, double tolsigma)

static void	st_evaluate_vario (int imod, int nbexp, std::vector< StrExp > &strexps, StrMod *strmod, VectorDouble &tabge)

static void	st_evaluate_vmap (int imod, StrMod *strmod, VectorDouble &tabge)

static int	st_parid_match (StrMod *strmod, int npar, int imod0, int icov0, const EConsElem &icons0, int ivar0, int jvar0)

static int	st_check_definite_positive (Model *model)

static int	st_truncate_negative_eigen (int nvar, int icov0, VectorDouble &matcor, VectorDouble &matcoru)

static double	st_sum_sills (int ivar0, int ncova, int nvar, VectorDouble &alpha)

static double	st_score (int nvar, int nvs2, int ncova, int npadir, VectorDouble &wt, VectorDouble &gg, VectorDouble &ge, VectorDouble &matcor)

static int	st_combineVariables (int ivar0, int jvar0)

static double	st_minimize_P4 (int icov0, int ivar0, int ncova, int nvar, int npadir, double xrmax, VectorDouble &xr, VectorDouble &alpha, VectorDouble &wt, VectorDouble &gg, VectorDouble &ge, const VectorDouble &consSill)

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

static void	st_updateOtherSills (int icov0, int ivar0, int ncova, int nvar, VectorDouble &alpha, VectorDouble &xr, VectorDouble &matcor)

static void	st_updateCurrentSillGoulard (int icov0, int ivar0, int npadir, int nvar, int ncova, VectorDouble &wt, VectorDouble &ge, VectorDouble &gg, const VectorDouble &consSill, VectorDouble &matcor)

static void	st_updateAlphaNoDiag (int icov0, int ivar0, int nvar, VectorDouble &xr, const VectorDouble &consSill, VectorDouble &matcor, VectorDouble &alpha)

static void	st_updateCurrentSillDiag (int icov0, int ivar0, int nvar, VectorDouble &alpha, VectorDouble &xr, VectorDouble &matcor)

static int	st_makeDefinitePositive (int icov0, int nvar, const VectorDouble &consSill, VectorDouble &matcor)

static int	st_optimize_under_constraints (int nvar, int nvs2, int ncova, int npadir, const Constraints &constraints, const Option_AutoFit &mauto, VectorDouble &wt, VectorDouble &gg, VectorDouble &ge, VectorDouble &matcor, double *score)

static int	st_initialize_goulard (int nvar, int nvs2, int ncova, int npadir, VectorDouble &wt, VectorDouble &gg, VectorDouble &ge, const VectorDouble &consSill, VectorDouble &matcor)

static void	st_goulard_sill_to_model (int nvar, int ncova, VectorDouble &sill, Model *model)

static int	st_goulard_with_constraints (const Constraints &constraints, const Option_AutoFit &mauto, int nvar, int ncova, int npadir, VectorDouble &wt, VectorDouble &gg, VectorDouble &ge, VectorDouble &sill)

static int	st_sill_fitting_int (Model *model, const Option_AutoFit &mauto, int npadir, VectorDouble &wt, VectorDouble &gg, VectorDouble &ge, VectorDouble &wt2, VectorDouble &ge1, VectorDouble &ge2, VectorDouble &gg2, VectorDouble &alphau, VectorDouble &sill1)

static int	st_goulard_fitting (int flag_reset, int flag_title, Model *model, const Constraints &constraints, const Option_AutoFit &mauto)

static int	st_model_has_intrinsic (Model model, const int filter)

static int	st_model_auto_strmod_reduce (StrMod strmod, int npar, double hmax, double gmax, VectorDouble &param, VectorDouble &lower, VectorDouble &upper, Constraints &constraints, Option_AutoFit &mauto)

static int	st_model_define (Model *model, const Option_VarioFit &optvar)

static int	st_alter_model_optvar (const Vario vario, Model model, Constraints &constraints, Option_VarioFit &optvar)

static int	st_alter_vmap_optvar (const Db dbmap, Model model, Constraints &constraints, Option_VarioFit &optvar)

static int	st_model_auto_count (const Vario vario, Model model1, Model *model2, Constraints &constraints, Option_VarioFit &optvar, VectorDouble &param, VectorDouble &lower, VectorDouble &upper)

static void	st_strmod_vario_evaluate (int nbexp, int npar, VectorDouble &param, VectorDouble &tabge)

static void	st_prepar_goulard_vmap (int imod)

static void	st_strmod_vmap_evaluate (int, int npar, VectorDouble &param, VectorDouble &tabge)

static void	st_vario_varchol_manage (const Vario vario, Model model, VectorDouble &varchol)

static void	st_vmap_varchol_manage (const Db *dbmap, VectorDouble &varchol)

static void	st_model_post_update (StrMod *strmod, const Option_VarioFit &optvar)

static int	st_manage_recint (const Option_AutoFit &mauto, int flag_exp, int ndim, int nvar, int nbexp, int ncova, int npadir)

static void	st_regularize_init ()

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)

static int	st_vmap_auto_count (const Db dbmap, Model model, Constraints &constraints, Option_VarioFit &optvar, VectorDouble &param, VectorDouble &lower, VectorDouble &upper)

static void	st_load_vmap (int npadir, VectorDouble &gg, VectorDouble &wt)

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)

Variables
static char	string [STRING_LENGTH]

static char	cov_name [STRING_LENGTH]

static int	CONGRUENCY = 50

static double	EpsFit = 1.e-12

static Regularize	REGULARIZE

static std::vector< StrExp >	STREXPS

static StrMod *	STRMOD = nullptr

static Option_AutoFit	MAUTO

static Constraints	CONSTRAINTS

static VectorInt	INDG1

static VectorInt	INDG2

static const DbGrid *	DBMAP

static void(*	ST_PREPAR_GOULARD )(int imod)

static Recint	RECINT

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_affect()

static void st_affect	(	int	rank,
		double	def_val,
		double	lower_val,
		double	upper_val,
		VectorDouble &	param,
		VectorDouble &	lower,
		VectorDouble &	upper
	)

static

Affect values for default, lower and upper parameters

Parameters

[in]	rank	Rank of the parameter
[in]	def_val	Default distance
[in]	lower_val	Lowest possible value
[in]	upper_val	Highest possible value
[in]	param	Array of current values
[in]	lower	Array of lower values
[in]	upper	Array of upper values

◆ st_alter_model_optvar()

static int st_alter_model_optvar	(	const Vario *	vario,
		Model *	model,
		Constraints &	constraints,
		Option_VarioFit &	optvar
	)

static

Define the options of the structure Opt_Vario

Returns: Error return code

Parameters

[in]	vario	Vario structure containing the exp. variogram
[in]	model	Model structure
[in]	constraints	Constraints structure
[out]	optvar	Opt_Vario structure

◆ st_alter_vmap_optvar()

static int st_alter_vmap_optvar	(	const Db *	dbmap,
		Model *	model,
		Constraints &	constraints,
		Option_VarioFit &	optvar
	)

static

Define the options of the structure Opt_Vario

Parameters

[in]	dbmap	Db Grid structure containing the Vmap
[in]	model	Model structure
[in]	constraints	Constraints structure
[out]	optvar	Opt_Vario structure

◆ st_check_definite_positive()

static int st_check_definite_positive ( Model * model )

static

Check if the resulting Model is definite positive

Returns: 1 if the model is not definite positive

Parameters

[in] model Model structure

◆ st_combineVariables()

static int st_combineVariables	(	int	ivar0,
		int	jvar0
	)

static

Return the rank of the variable pair, given the two variable ranks

Parameters

[in]	ivar0	Index of the first variable
[in]	jvar0	Index of the second variable

◆ st_compress_array()

static void st_compress_array	(	const Vario *	vario,
		int	npadir,
		VectorDouble &	tabin,
		VectorDouble &	tabout
	)

static

Compress the weights for the experimental variograms

Parameters

[in]	vario	Vario structure
[in]	npadir	Total number of lags
[in]	tabin	Uncompressed array
[out]	tabout	Compressed array

◆ st_compress_parid()

static int st_compress_parid	(	int	n_init,
		VectorInt &	parid,
		VectorDouble &	param,
		VectorDouble &	lower,
		VectorDouble &	upper
	)

static

Compress the default, lower and upper arrays

Returns: nb_tot : Number of parameters (input)

Parameters

[in]	n_init	Initial number of parameters
[in,out]	parid	Array of parameters identificators
[in,out]	param	Array of current values
[in,out]	lower	Array of lower values
[in,out]	upper	Array of upper values

◆ st_evaluate_vario()

static void st_evaluate_vario	(	int	imod,
		int	nbexp,
		std::vector< StrExp > &	strexps,
		StrMod *	strmod,
		VectorDouble &	tabge
	)

static

Calculates the values of the model corresponding to the experimentations stored in strexp

Parameters

[in]	nbexp	Number of experimental conditions
[in]	imod	Rank of the model
[in]	strexps	Array of StrExp structures to be freed
[in]	strmod	StrMod structure
[out]	tabge	Array of generic covariance values

◆ st_evaluate_vmap()

static void st_evaluate_vmap	(	int	imod,
		StrMod *	strmod,
		VectorDouble &	tabge
	)

static

Calculates the values of the model corresponding to the experimentations stored in vmap

Parameters

[in]	imod	Rank of the model
[in]	strmod	StrMod structure
[out]	tabge	Array of generic covariance values

◆ st_get_c00()

static double st_get_c00	(	const Vario *	vario,
		int	idir,
		int	ivar,
		int	jvar
	)

static

Find a possible value for C00 (for asymetric case only)

Returns: The value of C00

Parameters

[in]	vario	Vario structure
[in]	idir	Direction structure
[in]	ivar	First variable
[in]	jvar	Second variable

Remarks: This function is meant to provide a reasonable value for C00.; By default it returns the value of the central point; If this value if zero (with sw=0), it returns the closest non-zero; value

◆ st_get_vario_dimension()

static int st_get_vario_dimension	(	Vario *	vario,
		int *	nbexp_ret,
		int *	npadir_ret
	)

static

Return the number of experimental conditions

Returns: Error return code if the returned argument are null

Parameters

[in]	vario	Vario structure
[out]	nbexp_ret	Total number of experimental values (number of valid directions/lags/variables)
[out]	npadir_ret	Total number of lags for all directions

◆ st_get_vmap_dimension()

static int st_get_vmap_dimension	(	const Db *	dbmap,
		int	nvar,
		int *	nbexp_ret,
		int *	npadir_ret
	)

static

Return the number of experimental conditions

Returns: Error return code if the returned argument are null

Parameters

[in]	dbmap	Db structure
[in]	nvar	Number of variables
[out]	nbexp_ret	Total number of experimental values (number of valid directions/lags/variables)
[out]	npadir_ret	Maximum number of lags for all directions

◆ st_goulard_debug_current()

static void st_goulard_debug_current	(	int	nvar,
		int	ncova,
		int	iter,
		VectorDouble &	sill,
		double	crit
	)

static

Display the current status for the Goulard algorithm

Parameters

[in]	nvar	Number of variables
[in]	ncova	Number of basic structures
[in]	iter	Rank of the iteration
[in]	sill	Matrix of sills
[in]	crit	Valeu for the convergence criterion

◆ st_goulard_debug_title()

static void st_goulard_debug_title	(	int	nvar,
		int	ncova
	)

static

Display the title for the Goulard algorithm

Parameters

[in]	nvar	Number of variables
[in]	ncova	Number of basic structures

◆ st_goulard_fitting()

static int st_goulard_fitting	(	int	flag_reset,
		int	flag_title,
		Model *	model,
		const Constraints &	constraints,
		const Option_AutoFit &	mauto
	)

static

General Routine for fitting a model using an experimental variogram

Returns: Error return code

Parameters

[in]	flag_reset	1 to reset the array of sill before usage
[in]	flag_title	1 to print the title
[in]	model	Model to be fitted
[in]	constraints	Constraints structure
[in]	mauto	Option_AutoFit structure

◆ st_goulard_score()

static void st_goulard_score	(	const Option_AutoFit &	mauto,
		int	mode,
		int	ncova,
		int	iter,
		double	crit
	)

static

Display the Goulard final score

Parameters

[in]	mode	0 without; 1 with constraints
[in]	mauto	Option_AutoFit structure
[in]	ncova	Number of basic structures
[in]	iter	Number of iterations
[in]	crit	Convergence criterion

◆ st_goulard_sill_to_model()

static void st_goulard_sill_to_model	(	int	nvar,
		int	ncova,
		VectorDouble &	sill,
		Model *	model
	)

static

Copy the resulting sill into the Model

Parameters

[in]	nvar	Number of variables
[in]	ncova	Number of covariances
[in]	sill	Array of sills
[out]	model	Fitted Model structure

◆ st_goulard_verbose()

static void st_goulard_verbose	(	int	mode,
		Option_AutoFit &	mauto
	)

static

Manage the verbose option for iterative Goulard when called from Foxleg

Parameters

[in]	mode	0 (before) or 1 (after) Goulard
[in]	mauto	Option_AutoFit

◆ st_goulard_with_constraints()

static int st_goulard_with_constraints	(	const Constraints &	constraints,
		const Option_AutoFit &	mauto,
		int	nvar,
		int	ncova,
		int	npadir,
		VectorDouble &	wt,
		VectorDouble &	gg,
		VectorDouble &	ge,
		VectorDouble &	sill
	)

static

Internal function for Goulard under constraints

Returns: Error returned code

Parameters

[in]	constraints	Constraints structure
[in]	mauto	Option_AutoFit structure
[in]	nvar	Number of variables
[in]	ncova	Number of covariances
[in]	npadir	Maximum number of lags for all directions
[in]	wt	Array of weights (Dimension: npadir)
[in]	gg	Array of experimental values (Dimension: npadir)
[in]	ge	Array of model values (Dimension: npadir)
[out]	sill	Array of resulting sills

◆ st_goulard_without_constraint()

static int st_goulard_without_constraint	(	const Option_AutoFit &	mauto,
		int	nvar,
		int	ncova,
		int	npadir,
		VectorDouble &	wt,
		VectorDouble &	gg,
		VectorDouble &	ge,
		VectorDouble &	sill,
		double *	crit_arg
	)

static

Routine for fitting a model using an experimental variogram

Returns: Error return code

Parameters

[in]	mauto	Option_AutoFit structure
[in]	nvar	Number of variables
[in]	ncova	Number of covariances
[in]	npadir	Maximum number of lags for all directions
[in]	wt	Array of weights (Dimension: npadir)
[in]	gg	Array of experimental values (Dimension: npadir)
[in]	ge	Array of model values (Dimension: npadir)
[out]	sill	Array of resulting sills
[out]	crit_arg	Convergence criterion

Remarks: Internal arrays:; MP : Contains the current Model (ijvar,ipadir)

◆ st_initialize_goulard()

static int st_initialize_goulard	(	int	nvar,
		int	nvs2,
		int	ncova,
		int	npadir,
		VectorDouble &	wt,
		VectorDouble &	gg,
		VectorDouble &	ge,
		const VectorDouble &	consSill,
		VectorDouble &	matcor
	)

static

Initialize the system for Goulard algorithm

Returns: Error return code

Parameters

[in]	nvar	Number of variables
[in]	nvs2	Dimension of the triangular matrix (Dimension: nvar)
[in]	ncova	Number of basic structures
[in]	npadir	Total number of lags
[in]	wt	Array of weights attached to variogram lags
[in]	gg	Array of experimental values
[in]	ge	Array of generic covariance values
[in]	consSill	Array containing the constraints (optional)
[out]	matcor	Matrix of sills (Dimension: ncova * nvar * nvar)

◆ st_keypair_results()

static void st_keypair_results	(	int	mode,
		int	icov,
		int	nvar,
		double *	valpro,
		double *	vecpro
	)

static

Store the results using the keypair mechanism

Parameters

[in]	mode	1 for setting; -1 for deleting
[in]	icov	Rank of the covariance
[in]	nvar	Number of variables
[in]	valpro	Array of eigen values
[in]	vecpro	Array of eigen vectors

◆ st_keypair_sill()

static void st_keypair_sill	(	int	mode,
		Model *	model
	)

static

Save the sill matrices in the keypair mechanism

Parameters

[in]	mode	1 for setting; -1 for deleting
[in]	model	Model structure containing the basic structures

◆ st_load_ge()

static void st_load_ge	(	const Vario *	vario,
		Model *	model,
		int	npadir,
		VectorDouble &	dd,
		VectorDouble &	ge
	)

static

Calculates the values of a generic covariance model corresponding to the lags of an experimental variogram

Parameters

[in]	vario	Vario structure
[in]	model	Model structure
[in]	npadir	Total number of lags
[out]	dd	Array of distances (optional)
[out]	ge	Array of generic covariance values (optional)

◆ st_load_gg()

static void st_load_gg	(	const Vario *	vario,
		int	npadir,
		std::vector< StrExp > &	strexps,
		VectorDouble &	gg
	)

static

Fill the array of pointers on the experimental conditions

Parameters

[in]	vario	Vario structure
[in]	npadir	Total number of lags
[out]	strexps	Allocated array of StrExp (optional)
[out]	gg	Allocated array of experimental values

◆ st_load_vmap()

static void st_load_vmap	(	int	npadir,
		VectorDouble &	gg,
		VectorDouble &	wt
	)

static

Fill the array of pointers on the experimental conditions

Parameters

[in]	npadir	Maximum number of lags for all directions
[out]	gg	Allocated array of experimental uncompressed values
[out]	wt	Allocated array of experimental uncompressed weights

◆ st_load_wt()

static void st_load_wt	(	const Vario *	vario,
		int	wmode,
		int	npadir,
		VectorDouble &	wt
	)

static

Calculates the weighting factors for each experimental variogram value of each directional variogram

Parameters

[in]	vario	Vario structure
[in]	wmode	Type of the weighting procedure
[in]	npadir	Total number of lags
[out]	wt	Array of weights attached to variogram lags

◆ st_makeDefinitePositive()

static int st_makeDefinitePositive	(	int	icov0,
		int	nvar,
		const VectorDouble &	consSill,
		VectorDouble &	matcor
	)

static

Make sure the current matrix of sills is definite positive (diagonal unchanged)

Parameters

[in]	icov0	Index of the target basic structure
[in]	nvar	Number of variables
[in]	consSill	Vector of optional constant sills
[in,out]	matcor	Matrices of sills

◆ st_manage_recint()

static int st_manage_recint	(	const Option_AutoFit &	mauto,
		int	flag_exp,
		int	ndim,
		int	nvar,
		int	nbexp,
		int	ncova,
		int	npadir
	)

static

Manage memory for variogram fitting

Returns: Error returned code

Parameters

[in]	mauto	Option_AutoFit structure
[in]	flag_exp	1 for experimental variogram
[in]	ndim	Space dimension
[in]	nvar	Number of variables
[in]	nbexp	Number of experimental variogram values
[in]	ncova	Number of covariances
[in]	npadir	Total number of lags

◆ st_mauto_rescale()

static void st_mauto_rescale	(	int	nvar,
		VectorDouble &	varchol,
		Option_AutoFit &	mauto
	)

static

Update the parameters of Option_AutoFit structure according to the parameters of the Experimental Variogram

Parameters

[in]	nvar	Number of variables
[in]	varchol	Cholesky decomposition of the variance matrix
[in]	mauto	Option_AutoFit

◆ st_minimize_P4()

static double st_minimize_P4	(	int	icov0,
		int	ivar0,
		int	ncova,
		int	nvar,
		int	npadir,
		double	xrmax,
		VectorDouble &	xr,
		VectorDouble &	alpha,
		VectorDouble &	wt,
		VectorDouble &	gg,
		VectorDouble &	ge,
		const VectorDouble &	consSill
	)

static

Minimization of the order-4 polynomial

Returns: Value for which the fourth order polynomial is minimum.

Parameters

[in]	icov0	Index of the target basic structure
[in]	ivar0	Index of the target variable
[in]	ncova	Number of basic structures
[in]	nvar	Number of variables
[in]	npadir	Maximum number of lags for all directions
[in]	xrmax	Maximum value for the solution
[in]	xr	Current vector of sqrt(constraint/(sum of the sills))
[in]	alpha	Current auxiliary matrices alpha
[in]	wt	Array of weights attached to variogram lags
[in]	gg	Array of experimental values
[in]	ge	Array of generic covariance values
[in]	consSill	Vector of constant required sills (optional)

◆ st_model_auto_constraints_apply()

static void st_model_auto_constraints_apply	(	StrMod *	strmod,
		int	npar,
		const Constraints &	constraints,
		VectorDouble &	param,
		VectorDouble &	lower,
		VectorDouble &	upper
	)

static

Update default values, lower and upper bounds

Parameters

[in]	strmod	StrMod structure
[in]	npar	Number of parameters to be inferred
[in]	constraints	Constraints structure
[out]	param	Current values for parameters
[out]	lower	Lower values for parameters
[out]	upper	Upper values for parameters

◆ st_model_auto_count()

static int st_model_auto_count	(	const Vario *	vario,
		Model *	model1,
		Model *	model2,
		Constraints &	constraints,
		Option_VarioFit &	optvar,
		VectorDouble &	param,
		VectorDouble &	lower,
		VectorDouble &	upper
	)

static

Count the number of parameters

Returns: Number of parameters; -1 if one of the model components is not valid or; -2 if Goulard is not used althrough in multivariate

Parameters

[in]	vario	Vario structure containing the exp. variogram
[in]	model1	Model first structure
[in]	model2	Model second structure
[in]	constraints	Constraints structure
[in]	optvar	Opt_Vario structure
[out]	param	Array giving the default parameter value
[out]	lower	Array giving the minimum parameter value
[out]	upper	Array giving the maximum parameter value

Remarks: The arrays param, lower and upper are allocated by this function; (Dimension: Number of Parameters); They should be freed by the user

◆ st_model_auto_pardef()

static void st_model_auto_pardef	(	StrMod *	strmod,
		int	npar,
		double	hmax,
		VectorDouble &	varchol,
		VectorDouble &	angles,
		VectorDouble &	param,
		VectorDouble &	lower,
		VectorDouble &	upper
	)

static

Manage the default values of the parameters

Parameters

[in]	strmod	StrMod structure
[in]	npar	Number of parameters to be inferred
[in]	hmax	Maximum experimental distance
[in]	varchol	Cholesky decomposition of the variance matrix
[in]	angles	Reference angles
[out]	param	Current values for parameters
[out]	lower	Lower values for parameters
[out]	upper	Upper values for parameters

◆ st_model_auto_scldef()

static void st_model_auto_scldef	(	StrMod *	strmod,
		int	npar,
		double	hmax,
		VectorDouble &	varchol,
		VectorDouble &	scale
	)

static

Manage the scaling factors of the parameters

Parameters

[in]	strmod	StrMod structure
[in]	npar	Number of parameters to be inferred
[in]	hmax	Maximum experimental distance
[in]	varchol	Cholesky decomposition of the variance matrix
[out]	scale	Array of scales

◆ st_model_auto_strmod_alloc()

static StrMod* st_model_auto_strmod_alloc	(	Model *	model1,
		Model *	model2,
		int	npar0,
		int	norder,
		double	hmax,
		VectorDouble &	angles,
		const Option_VarioFit &	optvar,
		int *	npar
	)

static

Allocate the pointers on the StrMod structure

Returns: Pointer to the newly allocated StrMod structure

Parameters

[in]	model1	Model structure
[in]	model2	Model structure
[in]	npar0	Initial number of parameters
[in]	norder	Order of the Generalized Variogram
[in]	hmax	Maximum experimental distance
[in]	angles	Reference angles (coming from the Variogram)
[in]	optvar	Opt_Vario structure
[out]	npar	Final number of parameters

◆ st_model_auto_strmod_define()

static void st_model_auto_strmod_define	(	StrMod *	strmod,
		int	npar,
		VectorDouble &	param
	)

static

Load the parameters in the Model structure

Parameters

[in]	strmod	StrMod structure
[in]	npar	Number of parameters
[in]	param	Current values for parameters

◆ st_model_auto_strmod_free()

static StrMod* st_model_auto_strmod_free ( StrMod * strmod )

static

Free the pointer on the StrMod structure

Returns: Newly freed pointer

Parameters

[in] strmod Pointer to the StrMod to be freed

◆ st_model_auto_strmod_print()

static void st_model_auto_strmod_print	(	int	flag_title,
		StrMod *	strmod,
		const Constraints &	constraints,
		const Option_AutoFit &	mauto,
		VectorDouble &	param,
		VectorDouble &	lower,
		VectorDouble &	upper,
		int	npar,
		int	nbexp
	)

static

Print the resulting parameters of the Model

Parameters

[in]	flag_title	1 if the conditions must be printed
[in]	strmod	StrMod structure
[in]	constraints	Constraints structure
[in]	mauto	Option_AutoFit structure
[in]	param	Current values for parameters
[in]	lower	Array of lower values
[in]	upper	Array of upper values
[in]	npar	Number of parameters to be inferred
[in]	nbexp	Number of experiments

◆ st_model_auto_strmod_reduce()

static int st_model_auto_strmod_reduce	(	StrMod *	strmod,
		int *	npar,
		double	hmax,
		double	gmax,
		VectorDouble &	param,
		VectorDouble &	lower,
		VectorDouble &	upper,
		Constraints &	constraints,
		Option_AutoFit &	mauto
	)

static

Reduce the Model by discarding unnecessary basic structure(s)

Returns: 1 if the Model has been reduced; 0 otherwise

Parameters

[in]	strmod	StrMod structure
[in]	npar	Number of parameters
[in]	hmax	Maximum distance value
[in]	gmax	Maximum variogram value
[in]	param	Current values for parameters
[in]	lower	Lower values for parameters
[in]	upper	Upper values for parameters
[in]	constraints	Constraints structure
[in]	mauto	Option_AutoFit structure

◆ st_model_define()

static int st_model_define	(	Model *	model,
		const Option_VarioFit &	optvar
	)

static

Define the options of the model

Returns: Error return code

Parameters

[in]	model	Model structure containing the basic structures
[in]	optvar	Opt_Vario structure

TODO [Cova] : To be restored? Why this?

◆ st_model_has_intrinsic()

static int st_model_has_intrinsic	(	Model *	model,
		const int *	filter
	)

static

Check if the model has (at least) one intrinsic structure

Returns: 1 if model contains at least one intrinsic structure; 0 otherwise

Parameters

[in]	model	Model structure containing the basic structures
[in]	filter	Array specifying if a basic structure is filtered (1) or not (0). This array is optional

◆ st_model_post_update()

static void st_model_post_update	(	StrMod *	strmod,
		const Option_VarioFit &	optvar
	)

static

Update of the resulting model after the automatic fit

Parameters

[in,out]	strmod	StrMod structure
[in]	optvar	Opt_Vario structure

◆ st_modify_optvar_for_anam()

static void st_modify_optvar_for_anam	(	Model *	model,
		Option_VarioFit &	optvar
	)

static

◆ st_name_range()

static void st_name_range ( int ivar )

static

Compute the name of the range

Parameters

[in] ivar Rank of the variable

◆ st_name_rotation()

static void st_name_rotation ( int rank )

static

Compute the name of the rotation

Parameters

[in] rank Rank of the angle

◆ st_name_scale()

static void st_name_scale ( int ivar )

static

Compute the name of the scale factor

Parameters

[in] ivar Rank of the variable

◆ st_optimize_under_constraints()

static int st_optimize_under_constraints	(	int	nvar,
		int	nvs2,
		int	ncova,
		int	npadir,
		const Constraints &	constraints,
		const Option_AutoFit &	mauto,
		VectorDouble &	wt,
		VectorDouble &	gg,
		VectorDouble &	ge,
		VectorDouble &	matcor,
		double *	score
	)

static

Optimization under constraints

Returns: The number of iteration required to reach convergence

Parameters

[in]	nvar	Number of variables
[in]	nvs2	Dimension of the triangular matrix (Dimension: nvar)
[in]	ncova	Number of basic structures
[in]	npadir	Total number of lags
[in]	constraints	Constraints structure
[in]	mauto	Option_AutoFit structure
[in]	wt	Array of weights attached to variogram lags
[in]	gg	Array of experimental values
[in]	ge	Array of generic covariance values
[in]	matcor	Matrix of sills
[out]	score	The convergence score

◆ st_parid_alloc()

static int st_parid_alloc	(	StrMod *	strmod,
		int	npar0
	)

static

Define the list of parameters identificators

Returns: Number of parameters

Parameters

[in]	strmod	StrMod structure
[in]	npar0	Initial number of parameters to be inferred

◆ st_parid_decode()

static void st_parid_decode	(	int	parid,
		int *	imod,
		int *	icov,
		EConsElem *	icons,
		int *	ivar,
		int *	jvar
	)

static

Decode the parameter identificator

Parameters

[in]	parid	Parameter identificator
[out]	imod	Rank of the model
[out]	icov	Rank of the covariance
[out]	icons	Type of the constraint (EConsElem)
[out]	ivar	Rank of the first index
[out]	jvar	Rank of the second index

◆ st_parid_encode()

static int st_parid_encode	(	int	imod,
		int	icov,
		const EConsElem &	icons,
		int	ivar,
		int	jvar
	)

static

Encode the parameter identificator

Returns: Encoded parameter identificator

Parameters

[in]	imod	Rank of the model
[in]	icov	Rank of the covariance
[in]	icons	Type of the constraint (EConsElem)
[in]	ivar	Rank of the first index
[in]	jvar	Rank of the second index

◆ st_parid_match()

static int st_parid_match	(	StrMod *	strmod,
		int	npar,
		int	imod0,
		int	icov0,
		const EConsElem &	icons0,
		int	ivar0,
		int	jvar0
	)

static

Find the rank of the 'parid' which matches several criteria

Returns: Rank of the 'parid' or -1 if not found

Parameters

[in]	strmod	StrMod structure
[in]	npar	Number of parid
[in]	imod0	Rank of the target model (or -1)
[in]	icov0	Rank of the target covariance (or -1)
[in]	icons0	Target type of the constraint (EConsElem)
[in]	ivar0	Target first variable (or -1)
[in]	jvar0	Target second variable (or -1)

◆ st_prepar_goulard_vario()

static void st_prepar_goulard_vario ( int imod )

static

Prepare the array for Goulard's algorithm in the case of Variogram calculation

Parameters

[in] imod Rank of the model

◆ st_prepar_goulard_vmap()

static void st_prepar_goulard_vmap ( int imod )

static

Prepare the array for Goulard's algorithm in the case of VarioMap calculation

Parameters

[in] imod Rank of the model

◆ st_print()

static void st_print	(	const char *	name,
		int	flag_end,
		int	rank,
		VectorDouble &	param,
		VectorDouble &	lower,
		VectorDouble &	upper
	)

static

Prints a parameter of the basic structure

Parameters

[in]	name	Name of the parameter
[in]	flag_end	1 to add the end_of_line symbol - 0 otherwise
[in]	rank	Rank of the parameter
[in]	param	Array of current values
[in]	lower	Array of lower values
[in]	upper	Array of upper values

◆ st_regularize_init()

static void st_regularize_init ( )

static

Initialize the internal structure for regularization

◆ st_score()

static double st_score	(	int	nvar,
		int	nvs2,
		int	ncova,
		int	npadir,
		VectorDouble &	wt,
		VectorDouble &	gg,
		VectorDouble &	ge,
		VectorDouble &	matcor
	)

static

Evaluate the score

Returns: Value for the score

Parameters

[in]	nvar	Number of variables
[in]	nvs2	Dimension of the triangular matrix (Dimension: nvar)
[in]	ncova	Number of basic structures
[in]	npadir	Total number of lags
[in]	wt	Array of weights attached to variogram lags
[in]	gg	Array of experimental values
[in]	ge	Array of generic covariance values
[in]	matcor	Matrix of sills

◆ st_sill_fitting_int()

static int st_sill_fitting_int	(	Model *	model,
		const Option_AutoFit &	mauto,
		int	npadir,
		VectorDouble &	wt,
		VectorDouble &	gg,
		VectorDouble &	ge,
		VectorDouble &	wt2,
		VectorDouble &	ge1,
		VectorDouble &	ge2,
		VectorDouble &	gg2,
		VectorDouble &	alphau,
		VectorDouble &	sill1
	)

static

Routine for fitting a model using an experimental variogram for the Intrinsic case

Returns: Error return code

Parameters

[in]	model	Model to be fitted
[in]	mauto	Option_AutoFit structure
[in]	npadir	Maximum number of lags for all directions
[in]	wt	Array of weights (Dimension: npadir)
[in]	gg	Array of experimental values (Dimension: npadir)
[in]	ge	Array of model values (Dimension: npadir)
[in]	wt2	Array of weights (Dimension: npadir)
[in]	ge1	Array of model values (Dimension: npadir)
[in]	ge2	Array of model values (Dimension: npadir)
[in]	gg2	Array of experimental values (Dimension: npadir)
[out]	alphau	Array
[out]	sill1	Array of resulting sills

Remarks: Internal arrays:; MP : Contains the current Model (ijvar,ipadir)

◆ st_sill_reset()

static void st_sill_reset	(	int	nvar,
		int	ncova,
		VectorDouble &	sill
	)

static

Reset the array of sills

Parameters

[in]	nvar	Number of variables
[in]	ncova	Number of covariances
[out]	sill	Array of resulting sills

◆ st_strexp_manage()

static std::vector<StrExp> st_strexp_manage	(	int	nbexp,
		int	ndim
	)

static

Manage the array of StrExp structures

Returns: The vector of StrExp structures

Parameters

[in]	nbexp	Number of items
[in]	ndim	Space dimension

◆ st_strmod_vario_evaluate()

static void st_strmod_vario_evaluate	(	int	nbexp,
		int	npar,
		VectorDouble &	param,
		VectorDouble &	tabge
	)

static

Evaluate the model for an experiment

Returns: Value of the model

Parameters

[in]	nbexp	Number of experimental conditions
[in]	npar	Number of parameters
[in]	param	Current values for parameters
[out]	tabge	Array of resulting values

◆ st_strmod_vmap_evaluate()

static void st_strmod_vmap_evaluate	(	int	,
		int	npar,
		VectorDouble &	param,
		VectorDouble &	tabge
	)

static

Evaluate the model for an experiment

Returns: Value of the model

Parameters

[in]	npar	Number of parameters
[in]	param	Current values for parameters
[out]	tabge	Array of resulting values

◆ st_structure_reduce()

static int st_structure_reduce	(	StrMod *	strmod,
		int	imod,
		int	icov,
		double	hmax,
		double	gmax,
		double	tolsigma
	)

static

Check if the basic structure can be compressed

Returns: 1 if it can be removed; 0 otherwise

Parameters

[in]	strmod	StrMod structure
[in]	imod	Rank of the Model
[in]	icov	Rank of the basic structure
[in]	hmax	Maximum distance
[in]	gmax	Maximum variance
[in]	tolsigma	Percentage of the variance

◆ st_sum_sills()

static double st_sum_sills	(	int	ivar0,
		int	ncova,
		int	nvar,
		VectorDouble &	alpha
	)

static

Sum the sill of all the structures

Returns: The sum of the sills over all the sill matrices for a given variable

Parameters

[in]	ivar0	Index of the variable
[in]	ncova	Number of basic structures
[in]	nvar	Number of variables
[in]	alpha	The coregionalisation matrices (Dim: nvar^2 x ncova )

◆ st_truncate_negative_eigen()

static int st_truncate_negative_eigen	(	int	nvar,
		int	icov0,
		VectorDouble &	matcor,
		VectorDouble &	matcoru
	)

static

Cut negative eigen values

Returns: 1 if all eigen values are strictly positive

Parameters

[in]	nvar	Number of variables
[in]	icov0	Rank of the target basic structure
[in]	matcor	Matrix of sills
[out]	matcoru	Matrix of sills (can coincide with input matcor)

◆ 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)

◆ st_updateAlphaNoDiag()

static void st_updateAlphaNoDiag	(	int	icov0,
		int	ivar0,
		int	nvar,
		VectorDouble &	xr,
		const VectorDouble &	consSill,
		VectorDouble &	matcor,
		VectorDouble &	alpha
	)

static

Update 'sills' for the structures other than the current one

Parameters

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

◆ st_updateCurrentSillDiag()

static void st_updateCurrentSillDiag	(	int	icov0,
		int	ivar0,
		int	nvar,
		VectorDouble &	alpha,
		VectorDouble &	xr,
		VectorDouble &	matcor
	)

static

Update the sill matrix for the current structure 'icov0' (diagonal only)

Parameters

[in]	icov0	Target basic structure
[in]	ivar0	Index of the variable
[in]	nvar	Number of variables
[in]	alpha	Current auxiliary matrices alpha
[in]	xr	Current vector of sqrt(constraint/(sum of the sills))
[out]	matcor	Matrices of sills

◆ st_updateCurrentSillGoulard()

static void st_updateCurrentSillGoulard	(	int	icov0,
		int	ivar0,
		int	npadir,
		int	nvar,
		int	ncova,
		VectorDouble &	wt,
		VectorDouble &	ge,
		VectorDouble &	gg,
		const VectorDouble &	consSill,
		VectorDouble &	matcor
	)

static

Update the sill matrix for the current structure (except diagonal in the constrained case)

Parameters

[in]	icov0	Target basic structure
[in]	ivar0	Target variable
[in]	npadir	Maximum number of lags for all directions
[in]	nvar	Number of variables
[in]	ncova	Number of basic structures
[in]	wt	Array of weights attached to variogram lags
[in]	ge	Array of generic covariance values
[in]	gg	Array of experimental values
[in]	consSill	Vector of optional constrained sills
[out]	matcor	Matrices of sills

◆ st_updateOtherSills()

static void st_updateOtherSills	(	int	icov0,
		int	ivar0,
		int	ncova,
		int	nvar,
		VectorDouble &	alpha,
		VectorDouble &	xr,
		VectorDouble &	matcor
	)

static

Update 'sills' for the structures other than the current one

Parameters

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

◆ st_vario_varchol_manage()

static void st_vario_varchol_manage	(	const Vario *	vario,
		Model *	model,
		VectorDouble &	varchol
	)

static

Define Cholesky decomposition of the variance-covariance matrix

Parameters

[in]	vario	Vario structure
[in]	model	Model structure
[out]	varchol	Cholesky array

Remarks: In the case of Automatic Model Fitting with Properties,; we consider that the variance (of the Variogram) cannot; serve as initial value for the sill of the Model

◆ st_vmap_auto_count()

static int st_vmap_auto_count	(	const Db *	dbmap,
		Model *	model,
		Constraints &	constraints,
		Option_VarioFit &	optvar,
		VectorDouble &	param,
		VectorDouble &	lower,
		VectorDouble &	upper
	)

static

Count the number of parameters

Returns: Number of parameters; -1 if one of the model components is not valid; -2 if Goulard is not used althrough in multivariate

Parameters

[in]	dbmap	Db grid structure containing the Vmap
[in]	model	Model structure
[in]	constraints	Constraint structure
[in]	optvar	Opt_Vario structure
[out]	param	Array giving the default parameter value
[out]	lower	Array giving the minimum parameter value
[out]	upper	Array giving the maximum parameter value

◆ st_vmap_varchol_manage()

static void st_vmap_varchol_manage	(	const Db *	dbmap,
		VectorDouble &	varchol
	)

static

Define Cholesky decomposition of the variance-covariance matrix

Parameters

[in]	dbmap	Db structure
[out]	varchol	Cholesky array

◆ 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

Variable Documentation

◆ CONGRUENCY

int CONGRUENCY = 50

static

◆ CONSTRAINTS

Constraints CONSTRAINTS

static

◆ cov_name

char cov_name[STRING_LENGTH]

static

◆ DBMAP

const DbGrid* DBMAP

static

◆ EpsFit

double EpsFit = 1.e-12

static

◆ INDG1

VectorInt INDG1

static

◆ INDG2

VectorInt INDG2

static

◆ MAUTO

Option_AutoFit MAUTO

static

◆ RECINT

Recint RECINT

static

◆ REGULARIZE

Regularize REGULARIZE

static

◆ ST_PREPAR_GOULARD

void(* ST_PREPAR_GOULARD) (int imod) ( int imod )

static

◆ STREXPS

std::vector<StrExp> STREXPS

static

◆ string

char string[STRING_LENGTH]

static

◆ STRMOD

StrMod* STRMOD = nullptr

static

Functions

Variables

Function Documentation

◆ constraints_get()

◆ constraints_print()

◆ model_auto_fit()

◆ model_fitting_sills()

◆ modify_constraints_on_sill()

◆ st_affect()

◆ st_alter_model_optvar()

◆ st_alter_vmap_optvar()

◆ st_check_definite_positive()

◆ st_combineVariables()

◆ st_compress_array()

◆ st_compress_parid()

◆ st_evaluate_vario()

◆ st_evaluate_vmap()

◆ st_get_c00()

◆ st_get_vario_dimension()

◆ st_get_vmap_dimension()

◆ st_goulard_debug_current()

◆ st_goulard_debug_title()

◆ st_goulard_fitting()

◆ st_goulard_score()

◆ st_goulard_sill_to_model()

◆ st_goulard_verbose()

◆ st_goulard_with_constraints()

◆ st_goulard_without_constraint()

◆ st_initialize_goulard()

◆ st_keypair_results()

◆ st_keypair_sill()

◆ st_load_ge()

◆ st_load_gg()

◆ st_load_vmap()

◆ st_load_wt()

◆ st_makeDefinitePositive()

◆ st_manage_recint()

◆ st_mauto_rescale()

◆ st_minimize_P4()

◆ st_model_auto_constraints_apply()

◆ st_model_auto_count()

◆ st_model_auto_pardef()

◆ st_model_auto_scldef()

◆ st_model_auto_strmod_alloc()

◆ st_model_auto_strmod_define()

◆ st_model_auto_strmod_free()

◆ st_model_auto_strmod_print()

◆ st_model_auto_strmod_reduce()

◆ st_model_define()

◆ st_model_has_intrinsic()

◆ st_model_post_update()

◆ st_modify_optvar_for_anam()

◆ st_name_range()

◆ st_name_rotation()

◆ st_name_scale()

◆ st_optimize_under_constraints()

◆ st_parid_alloc()

◆ st_parid_decode()

◆ st_parid_encode()

◆ st_parid_match()

◆ st_prepar_goulard_vario()

◆ st_prepar_goulard_vmap()

◆ st_print()

◆ st_regularize_init()

◆ st_score()

◆ st_sill_fitting_int()

◆ st_sill_reset()

◆ st_strexp_manage()

◆ st_strmod_vario_evaluate()

◆ st_strmod_vmap_evaluate()

◆ st_structure_reduce()

◆ st_sum_sills()

◆ st_truncate_negative_eigen()

◆ st_updateAlphaDiag()

◆ st_updateAlphaNoDiag()

◆ st_updateCurrentSillDiag()

◆ st_updateCurrentSillGoulard()

◆ st_updateOtherSills()

◆ st_vario_varchol_manage()

◆ st_vmap_auto_count()