#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

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)

bool	isZero (double value, double eps)

bool	areEqual (double v1, double v2, double eps)

bool	isOne (double value, double eps)

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, const double tab, const double sel, const 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)

MatrixSquareGeneral	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)

operate_function	operate_Identify (int oper)

double	operate_Identity (double x)

double	operate_Inverse (double x)

double	operate_Square (double x)

double	operate_InverseSquare (double x)

double	operate_Sqrt (double x)

double	operate_InverseSqrt (double x)

double	modifyOperator (const EOperator &oper, double oldval, double value)

double	roundZero (double value, double eps)

double	truncateDecimals (double value, int ndec)

double	truncateDigits (double value, int ndigits)

void	setInternalDebug (bool status)

bool	isInternalDebug ()

Macro Definition Documentation

#define LSTACK 1000

#define MINI 10

Function Documentation

bool areEqual	(	double	v1,
		double	v2,
		double	eps
	)

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

int getClosestInteger ( double value )

int getITEST ( )

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

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)

double getMax	(	double	val1,
		double	val2
	)

double getMin	(	double	val1,
		double	val2
	)

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)

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

double getTEST ( )

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

bool isEven ( int number )

bool isInteger	(	double	value,
		double	eps
	)

bool isInternalDebug ( )

bool isMultiple	(	int	nbig,
		int	nsmall
	)

bool isOdd ( int number )

bool isOne	(	double	value,
		double	eps
	)

bool isZero	(	double	value,
		double	eps
	)

double modifyOperator	(	const EOperator &	oper,
		double	oldval,
		double	value
	)

Update an Old by a New value according to 'oper'

Parameters

oper	A keywork of EOperator enum
oldval	Old value
value	New value

operate_function operate_Identify ( int oper )

Identify the pointer to a function with following functionality: y = f(x)

Parameters

oper	Gives the type of operation to be performed 1: returns the value itslef (no change) -1: returns its inverse 2: returns the squared value -2: returns the inverse of the squared value 3: returns its square root -3: returns the inverse of the square root

Returns: Pointer to the specified function

double operate_Identity ( double x )

double operate_Inverse ( double x )

double operate_InverseSqrt ( double x )

double operate_InverseSquare ( double x )

double operate_Sqrt ( double x )

double operate_Square ( double x )

double roundZero	(	double	value,
		double	eps
	)

Round off the value if close enough to zero. This ensures that the printout of a very small value does not come out with a non-significant negative sign This trick should only serve to make printouts similar on different platforms.

Parameters

value	Input value
eps	Tolerance to check that the value is considered as small

Returns: The value itself or a very small positive value if the input value is too small.

void setInternalDebug ( bool status )

double truncateDecimals	(	double	value,
		int	ndec
	)

Rounding a double to a given number of decimals (from: https://stackoverflow.com/questions/304011/truncate-a-decimal-value-in-c/304013#304013)

Parameters

value	Value to be rounded up
ndec	Number of significant decimals

double truncateDigits	(	double	value,
		int	ndigits
	)

Rounding a double to a given number of decimals

Parameters

value	Value to be rounded up
ndigits	Number of significant digits

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

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)

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.

double ut_deg2rad ( double angle )

Translates from degree to radian

Parameters

[in] angle Angle in degrees

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

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

MatrixSquareGeneral 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

double ut_rad2deg ( double angle )

Translates from radian to degree

Parameters

[in] angle Angle in radian

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

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

StatResults ut_statistics	(	int	nech,
		const double *	tab,
		const double *	sel,
		const 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

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

Macros

Functions

Macro Definition Documentation

Function Documentation