/* Authors Martin Schlather, schlather@math.uni-mannheim.de Copyright (C) 2017 -- 2017 Martin Schlather This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 3 of the License, or (at your option) any later version. This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. */ #ifndef RFuser_H #define RFuser_H 1 #include "AutoRandomFields.h" #define startmode normal #define NM ((int) startmode) /* normal mode */ // {careless, sloppy, easygoing, normal, precise, pedantic, neurotic} #define NAT_SCALE 1 #define MAX_CE_MEM 16777216 #define generalN 20 // IMPORTANT: all names of general must be at least 3 letters long !!! #define GENERAL_MODUS 0 #define GENERAL_STORING 1 #define GENERAL_EXACTNESS 7 #define GENERAL_CLOSE 10 #define GENERAL_REPORTCOORD 16 // extern const char *general[generalN]; typedef struct general_param { char pch; /* character shown after each simulation just for entertainment of the user except for "!", then the numbers are shown */ bool allowdist0, na_rm_lines, vdim_close_together, storing, /* true: intermediate results are stored: might be rather memory consuming, but the simulation can (depending on the method chosen) be much faster when done the second time with exactly the same parameters do not forget to call DeleteAllKeys when all the simulation are done. false: intermediate results are not stored if SimulateRF is called or stored only until DoGauss is called. minimum memory consumption, but might be rather slow if many simulations are performed with exactly the same parameters to be safe call DeleteAllKeys when all the simulations are done Warning! init_foo may depend on GENERAL_STORING! This may cause problems if GENERAL_STORING is changed between init_foo and do_foo from false to true and if do_foo is called more than once! (SimulateRF is safe in this case) indifferent behaviour for the simulation methods if parameters are changed after each simulation. In case, only MAXFIELDS different parameter sets are used, but after each step the parameter set is changed, use different keynr for each parametere set, and STORING==true, to get fast simulation */ sp_conform, /* should the simulation result be return in as an sp class or in the old form ? --- getting obsolete in future --- todo */ detailed_output, returncall; int mode, /* hightailing, fast, normal, save, pedantic */ output, /* output mode, #alternative to mode that changes several other parameters; vice versa, spConform sets 'output' */ reportcoord ; /* see convert.T PL_* and RF.h PL_* */ int naturalscaling; /* is called PracticalRange in R (see Chiles&Delfiner?!) Note: RFparameters() allows only 0 and 1 as values! has an effect only if cov (and not only cov.local, e.g. Brownian Motion) is defined 0 : using the covariance function as defined 1 : rescaling of cov fctn such that cov(1)~=0.05, if rescaling function does not exist then failure, e.g. for Bessel model 2 : exact or approximate value (e.g. expPLUScirc) 3 : MLE (special needs taken into account, long memory covariance functions have too long tails, which are shortened (so threshold 0.05 is changed to higher values in this case)) +10: if any of the above fails : numerical evaluation! else : using rescaling if rescaling function exists, otherwise without rescaling */ int expected_number_simu, every, // every 'every' line is announced if every>0 // bool aniso; // currerntly cannot be changed by the user !! // siehe InternalCov.cc line 350, simu.cc line 394 /* expand scale immediately to aniso; this allows higher flexibility in the choice of simulation methods (tbm2 does not like dimension reduction), but it is slower */ Ttriple, set, seed_incr, seed_sub_incr; double gridtolerance; usr_bool exactness; } general_param; #define general_START \ {pch[NM], \ allowdistance0[NM], false, false, false /*storing*/, true, /* 5 */ \ false, false, \ startmode/* mode */ , output_sp, reportcoord_warnings, \ false, \ 1, 0, \ NA_INTEGER, 0, 0, 101101, \ 1e-6, exactness[NM] \ } #define gaussN 6 // extern const char *gauss[gaussN]; #define GAUSS_BEST_DIRECT 3 #define GAUSS_BOXCOX_OPTION 5 typedef struct gauss_param{ usr_bool stationary_only; // logical + NA double approx_zero, boxcox[2 * MAXBOXCOXVDIM]; bool paired, loggauss; int direct_bestvariables; } gauss_param; #define gauss_START {Nan, 0.05, \ {RF_INF, 0, RF_INF, 0, RF_INF, 0, RF_INF, 0, RF_INF, 0, \ RF_INF, 0, RF_INF, 0, RF_INF, 0, RF_INF, 0, RF_INF, 0}, \ false, false, 1200} #define MINSPLITN(locmaxn) ((locmaxn) / 40) #define MEDSPLITN(locmaxn) ((locmaxn) / 8) #define MAXSPLITN(locmaxn) (((locmaxn) * 5) / 8) #define krigeN 5 #define KRIGE_SPLITN 2 // extern const char *krige[krigeN]; typedef struct krige_param { bool ret_variance, fillall; int locmaxn, locsplitn[3], // 0:when splitting is done; 1:min pts in a neighbourhood ; 2:max pts when still neighbourhoods are included locsplitfactor; } krige_param; #define krige_START {false, true, locmaxn[NM], \ {MINSPLITN(locmaxn[NM]), MEDSPLITN(locmaxn[NM]), MAXSPLITN(locmaxn[NM])},\ 2} #define TRIVIALSTRATEGY 0 #define STRATEGY_PARTIAL 1 #define LASTSTRATEGY STRATEGY_PARTIAL #define CEN 12 // extern const char *CE[CEN]; typedef struct ce_param { bool force, useprimes, dependent; char strategy; int trials, maxgridsize, maxmem; double maxGB, tol_re, tol_im, mmin[MAXCEDIM], approx_grid_step; } ce_param; #define ce_START \ {ce_force[NM], true, ce_dependent[NM], \ TRIVIALSTRATEGY, \ 3, MAX_CE_MEM, MAXINT, \ 1, ce_tolRe[NM], ce_tolIm[NM], {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, \ ce_approx_step[NM]} #define spectralN 4 #define SPECTRAL_PROPFACTOR 2 // extern const char * spectral[spectralN]; typedef struct spectral_param { bool grid; double prop_factor, sigma; int lines[MAXTBMSPDIM]; } spectral_param; #define spectral_START {sp_grid[NM], 50, 0.0, {2500, 2500, 2500, 2500}} #define pTBMN 9 // extern const char * pTBM[pTBMN]; typedef struct tbm_param { bool grid; int tbmdim, fulldim, points, lines[MAXTBMSPDIM]; // number of lines simulated usr_bool layers; // 0, 1, NA double linesimufactor, /*factor by which simulation on the line is finer than the grid */ linesimustep, /* grid lag defined absolutely */ center[MAXTBMSPDIM]; } tbm_param; #define tbm_START \ {sp_grid[NM], -2, 3, 0, {1, 60, 500}, Nan, tbm_linefactor[NM], 0.0, \ {RF_NA, RF_NA, RF_NA, RF_NA}} #define MAX_DIRECT_MAXVAR 30000 #define DIRECT_ORIG_MAXVAR 12000 #define directN 1 #define DIRECT_MAXVAR_PARAM 0 // extern const char *direct[directN]; typedef struct direct_param { int maxvariables; } direct_param; #define direct_START { DIRECT_ORIG_MAXVAR } #define sequN 2 // extern const char * sequ[sequN]; typedef struct sequ_param{ int back, initial; } sequ_param; #define sequ_START {10, -10} typedef struct ave_param { } ave_param; #define ave_START {} #define pnuggetN 1 // extern const char * pnugget[pnuggetN]; typedef struct nugget_param { double tol; } nugget_param; #define nugget_START {nugget_tol[NM]} #define mppN 6 // extern const char * mpp[mppN]; typedef struct mpp_param { int n_estim_E, scatter_max[MAXMPPDIM]; double intensity[MAXMPPDIM], // intensity factor for e.g. unif_initu about_zero, shape_power, scatter_step[MAXMPPDIM]; } mpp_param; #define mpp_START \ {50000, \ {NA_INTEGER, NA_INTEGER, NA_INTEGER, NA_INTEGER}, \ {mpp_intensity[NM], mpp_intensity[NM], mpp_intensity[NM],mpp_intensity[NM]},\ mpp_zero[NM], /* about zero */ \ 2.0, \ {RF_NA, RF_NA, RF_NA, RF_NA} \ } #define HYPER_UNIFORM 0 #define HYPER_FRECHET 1 #define HYPER_BERNOULLI 2 #define hyperN 4 // extern const char * hyper[hyperN]; typedef struct hyper_param { int superpos, maxlines, mar_distr; double mar_param; } hyper_param; #define hyper_START {700, 1000, HYPER_UNIFORM, RF_NA} #define extremeN 12 // extern const char * extreme[extremeN]; #define EXTREME_FLAT 5 typedef struct extremes_param { usr_bool flathull; int maxpoints, check_every, min_n_zhou, max_n_zhou, mcmc_zhou, scatter_method; double standardmax, GEV_xi, density_ratio, eps_zhou, min_shape_gumbel; } extremes_param; #define extreme_START \ {False, MAXINT, 30, \ 1000, 10000000, 20, POISSON_SCATTER_ANY , \ max_max_gauss[NM], \ 1.0, 0.0, 0.01, -1e15 /* == -inf */ \ } #define brN 7 // extern const char * br[brN]; typedef struct br_param { int BRmaxmem, BRvertnumber, BRoptim, deltaAM; double BRmeshsize, BRoptimtol, variobound; } br_param; #define br_START \ {10000000, 7, 1, 300, \ 0.1, 0.01, 8.0} \ #define distrN 9 // extern const char * distr[distrN]; typedef struct distr_param{ double safety, minsteplen, innermin, outermax; int maxsteps, parts, maxit, mcmc_n, repetitions; } distr_param; #define distr_START {0.08, 0, 1e-20, 1e5, \ DISTMAXSTEPS, 8, 20, 15, 1000} // distr (rectangular) // todo should be 500 and better algorithm for approximation! #define MINCLIQUE(locmaxn) ((locmaxn) / 25) #define MEDCLIQUE(locmaxn) ((locmaxn) / 5) #define MAXCLIQUE(locmaxn) (((locmaxn) * 3) / 5) #define FLAT_UNDETERMINED -1 #define fitN 41 #define FIT_BC_LB 10 #define FIT_BC_UB 11 #define FIT_MAXNEIGHBOUR 23 #define FIT_CLIQUE 24 // extern const char * fit[fitN]; typedef struct fit_param{ double bin_dist_factor, upperbound_scale_factor, lowerbound_scale_factor, lowerbound_scale_LS_factor, upperbound_var_factor, lowerbound_var_factor, //lowerbound_sill, scale_max_relative_factor, minbounddistance, minboundreldist, BC_lambdaLB[2 * MAXBOXCOXVDIM], BC_lambdaUB[2 * MAXBOXCOXVDIM], scale_ratio, min_diag, pgtol, pgtol_recall, factr, factr_recall;// logical + NA // 19 int approximate_functioncalls, bins, nphi, ntheta, ntime, critical, n_crit, locmaxn, locsplitn[3], // 0:when splitting is done; 1:min pts in a neighbourhood ; 2:max pts when still neighbourhoods are included locsplitfactor, smalldataset, optimiser, algorithm, likelihood, split; // 16 usr_bool estimate_variance; bool use_naturalscaling, onlyuser, reoptimise, ratiotest_approx, split_refined, cross_refit, use_pseudovariogram; // 8 char lengthshortname; // 1..255 } fit_param; #define fit_START\ {0.4, 3.0, 5.0, 3.0, 10.0, 1e4, /* 6 */ \ 1000.0, 0.001, 0.02, /* 11 */ \ {-10, -10.0, -10, -10.0, -10, -10.0, -10, -10.0, -10, -10.0, \ -10, -10.0, -10, -10.0, -10, -10.0, -10, -10.0, -10, -10.0}, \ {10, 10.0, 10, 10.0, 10, 10.0, 10, 10.0, 10, 10.0, \ 10, 10.0, 10, 10.0, 10, 10.0, 10, 10.0, 10, 10.0}, \ 0.1, 1e-7, fit_pgtol[NM], fit_pgtol_recall[NM], /* 17 */ \ fit_factr[NM], fit_factr_recall[NM], \ /* int: */ \ 50, 20, 1, 1, 20, 0 /* critical */, /* 6 */ \ 5 /* ncrit */ , maxclique[NM], \ {MINCLIQUE(maxclique[NM]), MEDCLIQUE(maxclique[NM]),\ MAXCLIQUE(maxclique[NM]) }, 2, 2000, /* 11 */ \ 0, -1 /* algorithm */, 0, fit_split[NM], /* 16 */ \ Nan, \ false, false, fit_reoptimise[NM], /* 5 */\ fit_ratiotest_approx[NM], true, fit_cross_refit[NM], true, \ 12} // fit #define empvarioN 6 // extern const char * empvario[empvarioN]; typedef struct empvario_param{ double phi0, theta0, tol; bool pseudovariogram, fft, halveangles; } empvario_param; #define empvario_START {0.0, 0.0, 1e-13, false, true, true} #define guiN 3 // extern const char * gui[guiN]; #define GUI_SIZE 2 typedef struct gui_param{ bool alwaysSimulate; int method, size[2]; } gui_param; #define gui_START {true, CircEmbed, {1024, 64}} #define graphicsN 13 // extern const char *graphics[graphicsN]; #define GRAPHICS_UPTO 3 typedef struct graphics_param{ usr_bool always_open, always_close; bool split_screen, close_screen, onefile, grDefault; double height, width, resolution; int PL, increase_upto[2], number; char filename[100]; } graphics_param; #define graphics_START { True, Nan, \ true, true, false, false, \ 6.0, RF_NA, 72, 1, {3, 4}, 0, ""} #define registersN 3 // extern const char *registers[registersN]; typedef struct registers_param { int keynr, predict, likelihood; } registers_param; #define register_START {0, MODEL_PREDICT, MODEL_USER} #define internalN 26 // extern const char * internals[internalN]; #define INTERNALS_NEWANISO 2 #define INTERNALS_STORED_INIT 6 #define INTERNALS_ONGRID 9 #define INTERNALS_COORD_CHANGE 12 #define INTERNALS_ZENIT 14 #define INTERNALS_DO_TESTS 15 #define INTERNALS_EX_RED 20 #define INTERNALS_PARALLEL 21 // extern const char *internals[internalsN]; typedef struct internal_param{ // if changed, CHANGE ALSO RestWarnings in 'userinterfaces.cc'; bool warn_oldstyle, warn_newstyle, warn_Aniso, warn_ambiguous, warn_normal_mode, warn_mode, stored_init, warn_scale, warn_coordinates, warn_on_grid, warn_new_definitions, warn_aspect_ratio, warn_coord_change, warn_color_palette, warn_zenit, // not used anymore do_tests, warn_constant, warn_negvar, warn_onlyvar, // warn_def, ex_red warn_parallel, allow_duplicated_loc, warn_seed, warn_modus_operandi; usr_bool warn_mathdef;// int examples_reduced, declare_PL; } internal_param; #define internal_START\ {true, true, true, false, true, \ true, false, true, true, true, \ true, true, true, true, true, \ DO_TESTS, true, true, true, true, \ false, true, true, \ Nan, \ 0, 0} #define coordsN 14 #define COORDS_XYZNOTATION 0 #define COORDS_DATAIDX 5 #define COORDS_DATANAMES 6 #define COORDS_XIDX 7 #define COORDS_XNAMES 8 #define ZENIT 10 // extern const char *coords[coordsN]; typedef struct coords_param { usr_bool xyz_notation; // 0; + Exception for RFcov double zenit[2]; // 8 coord_sys_enum coord_system, // 1 new_coord_system; // 7 char newunits[MAXCOOORDNAMES][MAXUNITSCHAR], // 2; only to read curunits[MAXCOOORDNAMES][MAXUNITSCHAR], // 3 varunits[MAXDATANAMES][MAXUNITSCHAR], // 4 // 2 user variables for data.frames data_names[MAXDATANAMES][MAXCHAR],// (5C) x_names[MAXCOOORDNAMES][MAXCHAR]; // (6 C) // auxiliary variables int data_nr_names, // 5; needed by data_names and x_names x_nr_names, // 5; needed by data_names and x_names data_idx[2], x_idx[2]; // (5B, 6B) alternatives for data_names and x_namwe bool polar_coord, // 9 allow_earth2cart; // 11 angle_modes anglemode; // 10 } coords_param; #define coords_START \ {Nan, {1, RF_NA}, coord_auto, coord_keep, {""}, {""}, {""}, {""}, {""}, \ 0, 0, {1, NA_INTEGER}, {1, NA_INTEGER}, \ false, false, radians} #define specialN 1 // extern const char * special[specialN]; typedef struct special_param { int multcopies; } special_param; #define special_START {20} #define obsoleteN 10 // extern const char * obsolete[obsoleteN]; typedef struct globalparam{ general_param general; gauss_param gauss; krige_param krige; ce_param ce; spectral_param spectral; tbm_param tbm; direct_param direct; sequ_param sequ; ave_param ave; nugget_param nugget; mpp_param mpp; hyper_param hyper; extremes_param extreme; br_param br; distr_param distr; fit_param fit; empvario_param empvario; gui_param gui; graphics_param graphics; registers_param registers; internal_param internal; coords_param coords; special_param special; } globalparam; #define prefixN 23 //extern const char * prefixlist[prefixN], **all[prefixN]; //extern int allN[prefixN]; void setparameter(int i, int j, SEXP el, char name[200], bool isList, int local); void getRFoptions(SEXP *sublist); void finalparameter(int local); #endif