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,
			i
	)		(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]

value

Value 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

sel	Vector 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

sel	Vector giving the status of all samples (Dimension: absolute)
verbose	Verbose 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

map	The <int,int> map
iabs	Absolute 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]

value

Value 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]	nech	Number of samples
[in]	tab	Array of values
[in]	sel	Array containing the Selection or NULL
[in]	nclass	Number of sieve classes
[in]	start	Starting sieve value
[in]	pas	Width of the sieve
[out]	nmask	Number of masked values
[out]	ntest	Number of undefined values
[out]	nout	Number of values outside the classes
[out]	classe	Array 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]	n	Total number of objects (>= 1)
[in]	k	Selected 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]	n	Total number of objects (>1)
[in]	maxk	Selected number of objects (1<=maxk<n)
[out]	ncomb	Number 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]

angle

Angle 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]	nech	Number of samples
[in]	tab	Array of values
[in]	sel	Array containing the Selection or NULL
[out]	nval	Number of active values
[out]	mini	Minimum value
[out]	maxi	Maximum value

ut_median()

double ut_median	(	double *	tab,
		int	ntab
	)

Calculate the median from a table of values

Returns

The median value

Parameters

[in]	tab	Array of values
[in]	ntab	Number 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]

ndim

Size 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]

angle

Angle in radian

ut_shuffle_array()

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

Shuffle an array (by line)

Parameters

[in]	nrow	Number of rows
[in]	ncol	Number of columns
[in,out]	tab	Array 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]	safe	1 if the value array if preserved 0 if the value array is also sorted
[in]	nech	number of samples
[out]	ind	output int array
[out]	value	input 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]	nech	Number of samples
[in]	tab	Array of values
[in]	sel	Array containing the Selection or NULL
[in]	wgt	Array 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]	title	Title
[in]	nech	Number of samples
[in]	tab	Array of values
[in]	sel	Array containing the Selection or NULL