#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 "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 |
| 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 |