gstlearn  1.0.0
CCC
Utilities.cpp File Reference
#include "geoslib_d.h"
#include "geoslib_old_f.h"
#include "Enum/EDbg.hpp"
#include "Basic/Utilities.hpp"
#include "Basic/VectorHelper.hpp"
#include "Basic/Law.hpp"
#include <cmath>
#include <map>

Macros

#define LSTACK   1000
 
#define MINI   10
 
#define M(j, i)   (m[(i) * ndim + (j)])
 

Functions

bool isInteger (double value, double eps)
 
int getClosestInteger (double value)
 
bool isMultiple (int nbig, int nsmall)
 
bool isOdd (int number)
 
bool isEven (int number)
 
double getMin (double val1, double val2)
 
double getMax (double val1, double val2)
 
double getTEST ()
 
int getITEST ()
 
int FFFF (double value)
 
int IFFFF (int value)
 
double ut_deg2rad (double angle)
 
double ut_rad2deg (double angle)
 
void ut_sort_double (int safe, int nech, int *ind, double *value)
 
StatResults ut_statistics (int nech, double *tab, double *sel, double *wgt)
 
void ut_stats_mima_print (const char *title, int nech, double *tab, double *sel)
 
void ut_facies_statistics (int nech, double *tab, double *sel, int *nval, int *mini, int *maxi)
 
void ut_classify (int nech, double *tab, double *sel, int nclass, double start, double pas, int *nmask, int *ntest, int *nout, int *classe)
 
double ut_median (double *tab, int ntab)
 
double ut_cnp (int n, int k)
 
double * ut_pascal (int ndim)
 
int * ut_combinations (int n, int maxk, int *ncomb)
 
void ut_shuffle_array (int nrow, int ncol, double *tab)
 
VectorInt getListActiveToAbsolute (const VectorDouble &sel)
 
std::map< int, int > getMapAbsoluteToRelative (const VectorDouble &sel, bool verbose)
 
int getRankMapAbsoluteToRelative (const std::map< int, int > &map, int iabs)
 
int getRankMapRelativeToAbsolute (const std::map< int, int > &map, int irel)
 

Macro Definition Documentation

◆ LSTACK

#define LSTACK   1000

◆ M

#define M (   j,
 
)    (m[(i) * ndim + (j)])

◆ MINI

#define MINI   10

Function Documentation

◆ FFFF()

int FFFF ( double  value)

Checks if a double value is TEST

Returns
1 if a TEST value is encountered; 0 otherwise
Parameters
[in]valueValue to be tested

◆ getClosestInteger()

int getClosestInteger ( double  value)

◆ getITEST()

int getITEST ( )

◆ getListActiveToAbsolute()

VectorInt getListActiveToAbsolute ( const VectorDouble sel)

Returns the list of absolute indices for the only active samples A sample is active if its 'sel' value is equal to 1

Parameters
selVector giving the status of all samples (Dimension: absolute)
Returns

◆ getMapAbsoluteToRelative()

std::map<int, int> getMapAbsoluteToRelative ( const VectorDouble sel,
bool  verbose 
)

Returns the map such that MAP[iabs] = iact. A sample is active if its 'sel' value is equal to 1

Parameters
selVector giving the status of all samples (Dimension: absolute)
verboseVerbose flag
Returns
The map (dimension: nrel)

◆ getMax()

double getMax ( double  val1,
double  val2 
)

◆ getMin()

double getMin ( double  val1,
double  val2 
)

◆ getRankMapAbsoluteToRelative()

int getRankMapAbsoluteToRelative ( const std::map< int, int > &  map,
int  iabs 
)

Returns the rank of the relative grid node from its absolute index using the Map

Parameters
mapThe <int,int> map
iabsAbsolute rank of the grid node
Returns
Rank of the corresponding active (relative) grid node (or -1 is not found)

◆ getRankMapRelativeToAbsolute()

int getRankMapRelativeToAbsolute ( const std::map< int, int > &  map,
int  irel 
)

◆ getTEST()

double getTEST ( )

◆ IFFFF()

int IFFFF ( int  value)

Checks if an integer value is TEST

Returns
1 if a ITEST value is encountered; 0 otherwise
Parameters
[in]valueValue to be tested

◆ isEven()

bool isEven ( int  number)

◆ isInteger()

bool isInteger ( double  value,
double  eps 
)

◆ isMultiple()

bool isMultiple ( int  nbig,
int  nsmall 
)

◆ isOdd()

bool isOdd ( int  number)

◆ ut_classify()

void ut_classify ( int  nech,
double *  tab,
double *  sel,
int  nclass,
double  start,
double  pas,
int *  nmask,
int *  ntest,
int *  nout,
int *  classe 
)

Classify the samples into integer sieves

Parameters
[in]nechNumber of samples
[in]tabArray of values
[in]selArray containing the Selection or NULL
[in]nclassNumber of sieve classes
[in]startStarting sieve value
[in]pasWidth of the sieve
[out]nmaskNumber of masked values
[out]ntestNumber of undefined values
[out]noutNumber of values outside the classes
[out]classeArray for number of samples per sieve

◆ ut_cnp()

double ut_cnp ( int  n,
int  k 
)

Compute combinations(n,k)

Returns
Return the number of combinations of 'k' objects amongst 'n'
Parameters
[in]nTotal number of objects (>= 1)
[in]kSelected number of objects (>= 1)

◆ ut_combinations()

int* ut_combinations ( int  n,
int  maxk,
int *  ncomb 
)

Return all the combinations of k within n

Returns
Return all the combinations of 'k' objects amongst 'n'
Parameters
[in]nTotal number of objects (>1)
[in]maxkSelected number of objects (1<=maxk<n)
[out]ncombNumber of combinations
Remarks
The calling function must free the returned array.

◆ ut_deg2rad()

double ut_deg2rad ( double  angle)

Translates from degree to radian

Parameters
[in]angleAngle in degrees

◆ ut_facies_statistics()

void ut_facies_statistics ( int  nech,
double *  tab,
double *  sel,
int *  nval,
int *  mini,
int *  maxi 
)

Returns the statistics of an array containing the facies

Parameters
[in]nechNumber of samples
[in]tabArray of values
[in]selArray containing the Selection or NULL
[out]nvalNumber of active values
[out]miniMinimum value
[out]maxiMaximum value

◆ ut_median()

double ut_median ( double *  tab,
int  ntab 
)

Calculate the median from a table of values

Returns
The median value
Parameters
[in]tabArray of values
[in]ntabNumber of samples

◆ ut_pascal()

double* ut_pascal ( int  ndim)

Create the matrix containing the Pascal Triangle coefficients

Returns
A matrix (Dimension: ndim * ndim) containing the coefficients
or NULL if core allocation problem has been encountered
Parameters
[in]ndimSize of the matrix
Remarks
The calling function must free the returned matrix

◆ ut_rad2deg()

double ut_rad2deg ( double  angle)

Translates from radian to degree

Parameters
[in]angleAngle in radian

◆ ut_shuffle_array()

void ut_shuffle_array ( int  nrow,
int  ncol,
double *  tab 
)

Shuffle an array (by line)

Parameters
[in]nrowNumber of rows
[in]ncolNumber of columns
[in,out]tabArray to be suffled

◆ ut_sort_double()

void ut_sort_double ( int  safe,
int  nech,
int *  ind,
double *  value 
)

Sorts the (double) array value() and the array ind() in the ascending order of value

Parameters
[in]safe1 if the value array if preserved 0 if the value array is also sorted
[in]nechnumber of samples
[out]indoutput int array
[out]valueinput and output array
Remarks
If ind = NULL, ind is ignored

◆ ut_statistics()

StatResults ut_statistics ( int  nech,
double *  tab,
double *  sel,
double *  wgt 
)

Returns the statistics of an array in a StatResults structure

Parameters
[in]nechNumber of samples
[in]tabArray of values
[in]selArray containing the Selection or NULL
[in]wgtArray containing the Weights or NULL

◆ ut_stats_mima_print()

void ut_stats_mima_print ( const char *  title,
int  nech,
double *  tab,
double *  sel 
)

Print minimum and maximum of an array

Parameters
[in]titleTitle
[in]nechNumber of samples
[in]tabArray of values
[in]selArray containing the Selection or NULL