#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 "Basic/MathFunc.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 "Db/Db.hpp"
#include "Db/DbGrid.hpp"
#include "Variogram/Vario.hpp"
#include "Geometry/GeometryHelper.hpp"
#include "Matrix/MatrixSquareGeneral.hpp"
#include "LinearOp/CholeskyDense.hpp"
#include "Core/Keypair.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 ¶m, VectorDouble &lower, VectorDouble &upper) |
static int | st_compress_parid (int n_init, VectorInt &parid, VectorDouble ¶m, VectorDouble &lower, VectorDouble &upper) |
static void | st_print (const char *name, int flag_end, int rank, VectorDouble ¶m, VectorDouble &lower, VectorDouble &upper) |
static void | st_model_auto_strmod_print (int flag_title, StrMod *strmod, const Constraints &constraints, const Option_AutoFit &mauto, VectorDouble ¶m, 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 ¶m, VectorDouble &lower, VectorDouble &upper) |
static void | st_model_auto_pardef (StrMod *strmod, int npar, double hmax, VectorDouble &varchol, VectorDouble &angles, VectorDouble ¶m, VectorDouble &lower, VectorDouble &upper) |
static void | st_model_auto_strmod_define (StrMod *strmod, int npar, VectorDouble ¶m) |
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 ¶m, 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 ¶m, VectorDouble &lower, VectorDouble &upper) |
static void | st_strmod_vario_evaluate (int nbexp, int npar, VectorDouble ¶m, VectorDouble &tabge) |
static void | st_prepar_goulard_vmap (int imod) |
static void | st_strmod_vmap_evaluate (int, int npar, VectorDouble ¶m, 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 ¶m, 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) |
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 |
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
[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
[in] | vario | Vario structure |
[in,out] | model | Model to be fitted |
[in] | constraints | Constraints structure |
[in] | mauto | Option_AutoFit structure |
|
static |
|
static |
Define the options of the structure Opt_Vario
[in] | vario | Vario structure containing the exp. variogram |
[in] | model | Model structure |
[in] | constraints | Constraints structure |
[out] | optvar | Opt_Vario structure |
|
static |
Define the options of the structure Opt_Vario
[in] | dbmap | Db Grid structure containing the Vmap |
[in] | model | Model structure |
[in] | constraints | Constraints structure |
[out] | optvar | Opt_Vario structure |
|
static |
|
static |
Return the rank of the variable pair, given the two variable ranks
[in] | ivar0 | Index of the first variable |
[in] | jvar0 | Index of the second variable |
|
static |
Compress the weights for the experimental variograms
[in] | vario | Vario structure |
[in] | npadir | Total number of lags |
[in] | tabin | Uncompressed array |
[out] | tabout | Compressed array |
|
static |
|
static |
|
static |
|
static |
Find a possible value for C00 (for asymetric case only)
[in] | vario | Vario structure |
[in] | idir | Direction structure |
[in] | ivar | First variable |
[in] | jvar | Second variable |
|
static |
Return the number of experimental conditions
[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 |
|
static |
Return the number of experimental conditions
[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 |
|
static |
Display the current status for the Goulard algorithm
[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 |
|
static |
Display the title for the Goulard algorithm
[in] | nvar | Number of variables |
[in] | ncova | Number of basic structures |
|
static |
General Routine for fitting a model using an experimental variogram
[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 |
|
static |
Display the Goulard final score
[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 |
|
static |
|
static |
Manage the verbose option for iterative Goulard when called from Foxleg
[in] | mode | 0 (before) or 1 (after) Goulard |
[in] | mauto | Option_AutoFit |
|
static |
Internal function for Goulard under constraints
[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 |
|
static |
Routine for fitting a model using an experimental variogram
[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 |
|
static |
Initialize the system for Goulard algorithm
[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) |
|
static |
|
static |
Save the sill matrices in the keypair mechanism
[in] | mode | 1 for setting; -1 for deleting |
[in] | model | Model structure containing the basic structures |
|
static |
|
static |
Fill the array of pointers on the experimental conditions
[in] | vario | Vario structure |
[in] | npadir | Total number of lags |
[out] | strexps | Allocated array of StrExp (optional) |
[out] | gg | Allocated array of experimental values |
|
static |
Fill the array of pointers on the experimental conditions
[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 |
|
static |
|
static |
Make sure the current matrix of sills is definite positive (diagonal unchanged)
[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 |
|
static |
Manage memory for variogram fitting
[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 |
|
static |
Update the parameters of Option_AutoFit structure according to the parameters of the Experimental Variogram
[in] | nvar | Number of variables |
[in] | varchol | Cholesky decomposition of the variance matrix |
[in] | mauto | Option_AutoFit |
|
static |
Minimization of the order-4 polynomial
[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) |
|
static |
Update default values, lower and upper bounds
[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 |
|
static |
Count the number of 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 |
|
static |
Manage the default values of the 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 |
|
static |
|
static |
Allocate the pointers on the StrMod structure
[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 |
|
static |
|
static |
Print the resulting parameters of the Model
[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 |
|
static |
Reduce the Model by discarding unnecessary basic structure(s)
[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 |
|
static |
Define the options of the model
[in] | model | Model structure containing the basic structures |
[in] | optvar | Opt_Vario structure |
TODO [Cova] : To be restored? Why this?
|
static |
|
static |
Update of the resulting model after the automatic fit
[in,out] | strmod | StrMod structure |
[in] | optvar | Opt_Vario structure |
|
static |
|
static |
Compute the name of the range
[in] | ivar | Rank of the variable |
|
static |
Compute the name of the rotation
[in] | rank | Rank of the angle |
|
static |
Compute the name of the scale factor
[in] | ivar | Rank of the variable |
|
static |
Optimization under constraints
[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 |
|
static |
Define the list of parameters identificators
[in] | strmod | StrMod structure |
[in] | npar0 | Initial number of parameters to be inferred |
|
static |
Decode the parameter identificator
[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 |
|
static |
Encode the parameter identificator
[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 |
|
static |
Find the rank of the 'parid' which matches several criteria
[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) |
|
static |
Prepare the array for Goulard's algorithm in the case of Variogram calculation
[in] | imod | Rank of the model |
|
static |
Prepare the array for Goulard's algorithm in the case of VarioMap calculation
[in] | imod | Rank of the model |
|
static |
|
static |
Initialize the internal structure for regularization
|
static |
Evaluate the score
[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 |
|
static |
Routine for fitting a model using an experimental variogram for the Intrinsic case
[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 |
|
static |
Reset the array of sills
[in] | nvar | Number of variables |
[in] | ncova | Number of covariances |
[out] | sill | Array of resulting sills |
|
static |
Manage the array of StrExp structures
[in] | nbexp | Number of items |
[in] | ndim | Space dimension |
|
static |
Evaluate the model for an experiment
[in] | nbexp | Number of experimental conditions |
[in] | npar | Number of parameters |
[in] | param | Current values for parameters |
[out] | tabge | Array of resulting values |
|
static |
Evaluate the model for an experiment
[in] | npar | Number of parameters |
[in] | param | Current values for parameters |
[out] | tabge | Array of resulting values |
|
static |
|
static |
Sum the sill of all the structures
[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 ) |
|
static |
Cut negative eigen values
[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) |
void st_updateAlphaDiag | ( | int | icov0, |
int | ivar0, | ||
int | ncova, | ||
int | nvar, | ||
VectorDouble & | xr, | ||
VectorDouble & | alpha, | ||
const VectorDouble & | consSill | ||
) |
Update alpha from xr (diagonal only)
[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) |
|
static |
Update 'sills' for the structures other than the current one
[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 |
|
static |
Update the sill matrix for the current structure 'icov0' (diagonal only)
[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 |
|
static |
Update the sill matrix for the current structure (except diagonal in the constrained case)
[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 |
|
static |
Update 'sills' for the structures other than the current one
[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 |
|
static |
Define Cholesky decomposition of the variance-covariance matrix
|
static |
Count the number of 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 |
|
static |
Define Cholesky decomposition of the variance-covariance matrix
[in] | dbmap | Db structure |
[out] | varchol | Cholesky array |
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
[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 |
|
static |
|
static |
|
static |
|
static |
|
static |
|
static |
|
static |
|
static |
|
static |
|
static |
|
static |
|
static |
|
static |
|
static |