#include "geoslib_old_f.h"
#include "Polygon/Polygons.hpp"
#include "Basic/Utilities.hpp"
#include "Basic/GlobalEnvironment.hpp"
#include "Db/Db.hpp"
#include "Db/DbGrid.hpp"
#include <math.h>

Functions
int	compat_NDIM (Db db1, Db db2)

void	db_grid_print (Db *db)

int	get_LOCATOR_NITEM (const Db *db, const ELoc &locatorType)

double *	db_sample_free (double *tab)

double *	db_sample_alloc (const Db *db, const ELoc &locatorType)

double	distance_inter (const Db db1, const Db db2, int iech1, int iech2, double *dist_vect)

double	distance_intra (const Db db, int iech1, int iech2, double dist_vect)

double	distance_grid (DbGrid db, int flag_moins1, int iech1, int iech2, double dist_vect)

void	db_sample_print (Db *db, int iech, int flag_ndim, int flag_nvar, int flag_nerr)

int	db_center (Db db, double center)

int	db_attribute_range (const Db db, int iatt, double mini, double maxi, double delta)

int	db_grid_define_coordinates (DbGrid *db)

void	db_attribute_copy (Db *db, int iatt_in, int iatt_out)

int	db_gradient_update (Db *db)

int	db_selref (int ndim, const int nx, const int ref, const double tabin, double tabout)

int	db_locator_attribute_add (Db db, const ELoc &locatorType, int number, int r_tem, double valinit, int iptr)

int	db_grid_copy (DbGrid db1, DbGrid db2, const int ind1, const int ind2, int ncol, int *cols)

int	db_grid_copy_dilate (DbGrid db1, int iatt1, DbGrid db2, int iatt2, int mode, const int *nshift)

int	point_to_point (Db db, const double coor)

int	point_to_grid (const DbGrid db, const double coor, int flag_outside, int *indg)

int	point_to_bench (const DbGrid db, double coor, int flag_outside, int *indb)

int	index_point_to_grid (const Db dbin, int iech, int flag_outside, const DbGrid dbout, double *coor)

int	point_inside_grid (Db db, int iech, const DbGrid dbgrid)

void	db_monostat (Db db, int iatt, double wtot, double mean, double var, double mini, double maxi)

int	db_proportion (Db db, DbGrid dbgrid, int nfac1max, int nfac2max, int *nclout)

int	db_merge (Db db, int ncol, int cols)

int	db_count_defined (Db *db, int icol)

void	db_locators_correct (VectorString &strings, const VectorInt &current, int flag_locnew)

int	db_prop_read (DbGrid db, int ix, int iy, double props)

int	db_prop_write (DbGrid db, int ix, int iy, double props)

double *	db_distances_general (Db db1, Db db2, int niso, int mode, int flag_same, int n1, int n2, double dmin, double dmax)

int	db_is_isotropic (const Db db, int iech, double data)

int	is_grid_multiple (DbGrid db1, DbGrid db2)

DbGrid *	db_create_grid_multiple (DbGrid *dbin, const VectorInt &nmult, bool flagAddSampleRank)

DbGrid *	db_create_grid_divider (DbGrid *dbin, const VectorInt &nmult, bool flagAddSampleRank)

DbGrid *	db_create_grid_dilate (DbGrid *dbin, int mode, const VectorInt &nshift, bool flagAddSampleRank)

int	db_gradient_modang_to_component (Db *db, int ang_conv, int iad_mod, int iad_ang, int iad_gx, int iad_gy)

int	db_gradient_component_to_modang (Db *db, int verbose, int iad_gx, int iad_gy, int iad_mod, int iad_ang, double scale, double ve)

int	db_get_rank_absolute_to_relative (Db *db, int iech0)

int	db_get_rank_relative_to_absolute (Db *db, int iech0)

double	get_grid_value (DbGrid *dbgrid, int iptr, VectorInt &indg, int ix, int iy, int iz)

void	set_grid_value (DbGrid *dbgrid, int iptr, VectorInt &indg, int ix, int iy, int iz, double value)

DbGrid *	db_grid_reduce (DbGrid db_grid, int iptr, const int margin, const int *limmin, int flag_sel, int flag_copy, int verbose, double vmin, double vmax)

int	db_grid_patch (DbGrid ss_grid, DbGrid db_grid, int iptr_ss, int iptr_db, int iptr_rank, int new_rank, int oper, int verbose)

int	db_name_identify (Db *db, const String &string)

static void	st_rotate (int ndim, double *rotmat, VectorDouble &coor, VectorDouble &mini, VectorDouble &maxi)

void	db_extension_rotated (Db db, double rotmat, VectorDouble &mini, VectorDouble &maxi)

VectorDouble	db_get_grid_axis (DbGrid *dbgrid, int idim)

VectorDouble	db_get_attribute (Db *db, int iatt, bool verbose)

VectorInt	db_identify_variables_by_name (Db *db, const String &pattern)

void	grid_iterator_init (Grid *grid, const VectorInt &order)

VectorInt	grid_iterator_next (Grid *grid)

Function Documentation

◆ compat_NDIM()

int compat_NDIM	(	Db *	db1,
		Db *	db2
	)

Checks if Space Dimension of the first Db is at least equal to the Space dimension of the second Db

Returns: Error return code (1 for correct and 0 otherwise).; If an error is found, a message is issued

Parameters

[in]	db1	first Db descriptor (usually dbin)
[in]	db2	second Db descriptor (usually dbout)

Remarks: This function checks that the we can perform an operation; (i.e. migration) of a 2D data set onto a 3D one.; The opposite is not correct.

◆ db_attribute_copy()

void db_attribute_copy	(	Db *	db,
		int	iatt_in,
		int	iatt_out
	)

Copy an attribute into another one

Parameters

[in]	db	Db structure
[in]	iatt_in	Input attribute
[in]	iatt_out	Output attribute

Remarks: The output attribute must already be allocated

◆ db_attribute_range()

int db_attribute_range	(	const Db *	db,
		int	iatt,
		double *	mini,
		double *	maxi,
		double *	delta
	)

Returns the range for a variable

Returns: Error return code

Parameters

[in]	db	Db structure
[in]	iatt	Rank of the target attribute
[out]	mini	Minimum
[out]	maxi	Maximum
[out]	delta	Extension

◆ db_center()

int db_center	(	Db *	db,
		double *	center
	)

Returns the center of a data set

Parameters

[in]	db	Db structure
[out]	center	Array containing the coordinates of the center (Dimension = get_NDIM(db))

◆ db_count_defined()

int db_count_defined	(	Db *	db,
		int	icol
	)

Count the number of defined values for a Column

Returns: Number of defined values

Parameters

[in]	db	Db structure
[in]	icol	Rank of the Column

◆ db_create_grid_dilate()

DbGrid* db_create_grid_dilate	(	DbGrid *	dbin,
		int	mode,
		const VectorInt &	nshift,
		bool	flagAddSampleRank
	)

Create a Grid Db extended (or compressed) from another Grid Db

Returns: Pointer to the newly created Db grid structure

Parameters

[in]	dbin	Initial Db Grid
[in]	mode	1 for extending; -1 for compressing
[in]	nshift	Array of shifts
[in]	flagAddSampleRank	true to add the 'rank' as first column

◆ db_create_grid_divider()

DbGrid* db_create_grid_divider	(	DbGrid *	dbin,
		const VectorInt &	nmult,
		bool	flagAddSampleRank
	)

Create a Grid Db as a divider of another Grid Db

Returns: Pointer to the newly created Db grid structure

Parameters

[in]	dbin	Initial Db Grid
[in]	nmult	Array of subdivision coefficients
[in]	flagAddSampleRank	true to add the 'rank' as first column

◆ db_create_grid_multiple()

DbGrid* db_create_grid_multiple	(	DbGrid *	dbin,
		const VectorInt &	nmult,
		bool	flagAddSampleRank
	)

Create a Grid Db as a multiple of another Grid Db

Returns: Pointer to the newly created Db grid structure

Parameters

[in]	dbin	Initial Db Grid
[in]	nmult	Array of multiplicity coefficients
[in]	flagAddSampleRank	true to add the 'rank' as first column **

◆ db_distances_general()

double* db_distances_general	(	Db *	db1,
		Db *	db2,
		int	niso,
		int	mode,
		int	flag_same,
		int *	n1,
		int *	n2,
		double *	dmin,
		double *	dmax
	)

Evaluate the array of distances between samples of two data sets Takes the selection into account (if any)

Returns: Array of distances

Parameters

[in]	db1	Db first descriptor
[in]	db2	Db second descriptor
[in]	niso	Number of variables tested for isotopy
[in]	mode	Type of array returned 0 : extreme distances 1 : the distance to the closest sample 2 : all point-to-point distances
[in]	flag_same	1 if both Db coincide
[out]	n1	First dimension of the returned array
[out]	n2	Second dimension of the returned array
[out]	dmin	Minimum distance
[out]	dmax	Maximum distance

Remarks: The returned array must be freed by calling routine; When the two Dbs coincide, the distance calculation excludes; the comparison between one sample and itself

◆ db_extension_rotated()

void db_extension_rotated	(	Db *	db,
		double *	rotmat,
		VectorDouble &	mini,
		VectorDouble &	maxi
	)

Returns the extension of the field (after rotation) along each axis

Parameters

[in]	db	Db structure
[in]	rotmat	Rotation matrix (optional)
[out]	mini	Array containing the minimum (Dimension = ndim)
[out]	maxi	Array containing the maximum (Dimension = ndim)

Remarks: This function does nothing if:; - no rotation matrix is defined; - space dimension is 1 or larger than 3

◆ db_get_attribute()

VectorDouble db_get_attribute	(	Db *	db,
		int	iatt,
		bool	verbose
	)

Returns a vector containing an attribute from the Db

Parameters

[in]	db	Db structure
[in]	iatt	Attribute rank
[in]	verbose	Verbose flag

◆ db_get_grid_axis()

VectorDouble db_get_grid_axis	(	DbGrid *	dbgrid,
		int	idim
	)

Returns a vector containing the coordinates of a Grid along one space dimension

Parameters

[in]	dbgrid	Db structure
[in]	idim	Rank of the space dimension

◆ db_get_rank_absolute_to_relative()

int db_get_rank_absolute_to_relative	(	Db *	db,
		int	iech0
	)

Returns the relative rank of a sample from its absolute rank These are different due to the presence of a selection Returns -1 if not found

Returns: Relative rank of a sample

Parameters

[in]	db	Db structure
[in]	iech0	Absolute sample rank

◆ db_get_rank_relative_to_absolute()

int db_get_rank_relative_to_absolute	(	Db *	db,
		int	iech0
	)

Returns the absolute rank of a sample from its relative rank These are different due to the presence of a selection Returns -1 if not found

Returns: Relative rank of a sample

Parameters

[in]	db	Db structure
[in]	iech0	Relative sample rank

◆ db_gradient_component_to_modang()

int db_gradient_component_to_modang	(	Db *	db,
		int	verbose,
		int	iad_gx,
		int	iad_gy,
		int	iad_mod,
		int	iad_ang,
		double	scale,
		double	ve
	)

Transform a set of gradients defined by (gx,gy) into (module,angle) Only defined in the 2-D case

Returns: Error return code

Parameters

[in]	db	Initial Db
[in]	verbose	1 for the verbose option
[in]	iad_gx	Rank of the 'gx' attribute (defined in degrees)
[in]	iad_gy	Rank of the 'gy' attribute (defined in degrees)
[in]	iad_mod	Rank of the 'modulus' attribute
[in]	iad_ang	Rank of the 'angle' attribute (defined in degrees)
[in]	scale	Scaling factor for the modulus
[in]	ve	Moderation factor

Remarks: Attributes 'gx' and 'gy' may coincide with 'modulus' and 'angle'; Nothing is done (more than the conversion if scale=1 and ve=0)

◆ db_gradient_modang_to_component()

int db_gradient_modang_to_component	(	Db *	db,
		int	ang_conv,
		int	iad_mod,
		int	iad_ang,
		int	iad_gx,
		int	iad_gy
	)

Transform a set of gradients defined by (modulus,angle) into (gx,gy) Only defined in the 2-D case

Returns: Error return code

Parameters

[in]	db	Initial Db
[in]	ang_conv	Convention on the angles 1: Trigonometry (From East counter-clockwise) 2: From North clockwise
[in]	iad_mod	Rank of the 'modulus' attribute
[in]	iad_ang	Rank of the 'angle' attribute (defined in degrees)
[in]	iad_gx	Rank of the 'gx' attribute (defined in degrees)
[in]	iad_gy	Rank of the 'gy' attribute (defined in degrees)

Remarks: Attributes 'gx' and 'gy' may coincide with 'modulus' and 'angle'

◆ db_gradient_update()

int db_gradient_update ( Db * db )

Update the Data Base for Kriging with Gradient components This update is limited to the 2-D Monovariate case

Returns: Error return code

Parameters

[in] db Input/output Db

Remarks: In the case of Kriging with Gradient, the gradient components; are transformed into additional variables

◆ db_grid_copy()

int db_grid_copy	(	DbGrid *	db1,
		DbGrid *	db2,
		const int *	ind1,
		const int *	ind2,
		int	ncol,
		int *	cols
	)

Copy a set of variables from a grid Db to another grid Db

Returns: Error return code

Parameters

[in]	db1	Input Grid Db structure
[in]	db2	Output Grid Db structure
[in]	ind1	Vector of indices in the first Db
[in]	ind2	Vector of indices in the second Db
[in]	ncol	Number of variables to be copied
[in]	cols	Array of input variable columns

◆ db_grid_copy_dilate()

int db_grid_copy_dilate	(	DbGrid *	db1,
		int	iatt1,
		DbGrid *	db2,
		int	iatt2,
		int	mode,
		const int *	nshift
	)

Copy one variable from a grid Db to another grid Db

Returns: Error return code

Parameters

[in]	db1	Input Grid Db structure
[in]	iatt1	Rank of the attribute in Db1
[in]	db2	Output Grid Db structure
[in]	iatt2	Rank of the attribute in Db2
[in]	mode	1 for dilation; -1 for compression
[in]	nshift	Vector of dilation/compression defined in cell number along each space direction

◆ db_grid_define_coordinates()

int db_grid_define_coordinates ( DbGrid * db )

Define the coordinates in a Grid structure

Returns: Error return code

Parameters

[in] db The Db structure

Remarks: This function considers the grid characteristics and updates; the locators dedicated to coordinates; This makes sense when a new grid is generated or when the grid; characteristics have changed

◆ db_grid_patch()

int db_grid_patch	(	DbGrid *	ss_grid,
		DbGrid *	db_grid,
		int	iptr_ss,
		int	iptr_db,
		int	iptr_rank,
		int	new_rank,
		int	oper,
		int	verbose
	)

Patch a sub-grid within a main grid

Returns: Error return code

Parameters

[in]	ss_grid	Db sub-grid structure
[in]	db_grid	Db main grid structure
[in]	iptr_ss	Rank of the attribute in the sub-grid
[in]	iptr_db	Rank of the attribute in the main grid
[in]	iptr_rank	Rank of the attribute storing object rank If <0, no check is performed: always patch
[in]	new_rank	Rank of the current object to patch in main grid
[in]	oper	>0 for larger; <0 for smaller
[in]	verbose	Verbose flag

Remarks: When 'iptr_rank' is defined (>=0), defined pixels of the current; sub-grid overwrites the corresponding pixel within the main grid; only if its rank ('new_rank') is larger (oper>0) or smaller; (oper<0)) than the rank of the same pixel in the main grid; (attribute 'iptr_rank'); When 'iptr_rank' is undefined, arguments 'new_rank', 'oper' are; useless

◆ db_grid_print()

void db_grid_print ( Db * db )

Print the Summary of the Grid structure

Parameters

[in] db Pointer to the Db structure (organized as a grid)

◆ db_grid_reduce()

DbGrid* db_grid_reduce	(	DbGrid *	db_grid,
		int	iptr,
		const int *	margin,
		const int *	limmin,
		int	flag_sel,
		int	flag_copy,
		int	verbose,
		double	vmin,
		double	vmax
	)

Extract the subgrid (from a grid) which contains the only cells where the target variable lies wuthin the target interval Selection in the input grid is taken into account

Parameters

[in]	db_grid	Db structure
[in]	iptr	Rank of the column of the target variable
[in]	margin	Array of margins (or NULL)
[in]	limmin	Array of minimum dimensions (or NULL)
[in]	flag_sel	Create the selection
[in]	flag_copy	1 if the selection must be copied in sub-grid
[in]	verbose	Verbose flag
[in]	vmin	Lower bound (included)
[in]	vmax	Upper bound (excluded)

◆ db_identify_variables_by_name()

VectorInt db_identify_variables_by_name	(	Db *	db,
		const String &	pattern
	)

Identify the variables of a Db where names match a criterion

Parameters

[in]	db	Db structure
[in]	pattern	Matching pattern

◆ db_is_isotropic()

int db_is_isotropic	(	const Db *	db,
		int	iech,
		double *	data
	)

Check if a sample corresponding to a given rank is isotropic or not

Returns: 1 if all variables are defined (isotropic) or 0 otherwise

Parameters

[in]	db	Db descriptor
[in]	iech	Rank of the target sample
[out]	data	If the output vector is provided. Returns the data value at target sample as a vector

Remarks: If the sample is masked off, the function returns 0

◆ db_locator_attribute_add()

int db_locator_attribute_add	(	Db *	db,
		const ELoc &	locatorType,
		int	number,
		int	r_tem,
		double	valinit,
		int *	iptr
	)

Add and initiate several attributes corresponding to a given locator

Returns: Error return code

Parameters

[in]	db	Db structure
[in]	locatorType	Rank of the Pointer (ELoc)
[in]	number	Number of locators to be defined
[in]	r_tem	Rank of the first item in the pointer
[in]	valinit	Value to be used for initialization
[out]	iptr	Rank of the first new attribute

◆ db_locators_correct()

void db_locators_correct	(	VectorString &	strings,
		const VectorInt &	current,
		int	flag_locnew
	)

Define the array of locators

Parameters

[in]	strings	Array of locators
[in]	current	Array of ranks of the last modified locator
[in]	flag_locnew	Reset all locators

Remarks: The elements of the array current are numbered starting from 1

◆ db_merge()

int db_merge	(	Db *	db,
		int	ncol,
		int *	cols
	)

Merge a set of variables in a Db

Returns: Error return code

Parameters

[in]	db	Db structure
[in]	ncol	Number of variables to be merged
[in]	cols	Array of input variable columns

◆ db_monostat()

void db_monostat	(	Db *	db,
		int	iatt,
		double *	wtot,
		double *	mean,
		double *	var,
		double *	mini,
		double *	maxi
	)

Monovariate statistics

Parameters

[in]	db	Db structure
[in]	iatt	Rank of the attribute
[out]	wtot	Sum of the weights
[out]	mean	Mean of the variable
[out]	var	variance of the variable
[out]	mini	Minimum value
[out]	maxi	Maximum value

◆ db_name_identify()

int db_name_identify	(	Db *	db,
		const String &	string
	)

Identify the attribute by its name

Returns: Rank of the variable starting from 0 (or -1 if not found)

Parameters

[in]	db	Db descriptor
[in]	string	attribute name

◆ db_prop_read()

int db_prop_read	(	DbGrid *	db,
		int	ix,
		int	iy,
		double *	props
	)

Read the proportions per VPC

Returns: Error return code

Parameters

[in]	db	Db structure
[in]	ix	Rank of the cell along X
[in]	iy	Rank of the cell along Y
[out]	props	Array of proportions

Remarks: This procedure is meant for a 3-D grid file

◆ db_prop_write()

int db_prop_write	(	DbGrid *	db,
		int	ix,
		int	iy,
		double *	props
	)

Write the proportions (per VPC)

Returns: Error return code

Parameters

[in]	db	Db structure
[in]	ix	Rank of the cell along X
[in]	iy	Rank of the cell along Y
[in]	props	Array of proportions

Remarks: This procedure is meant for a 3-D grid file

◆ db_proportion()

int db_proportion	(	Db *	db,
		DbGrid *	dbgrid,
		int	nfac1max,
		int	nfac2max,
		int *	nclout
	)

Calculates the proportions of facies within a grid

Returns: Error return code

Parameters

[in]	db	Input Db structure
[in]	dbgrid	Output Grid Db structure
[in]	nfac1max	Maximum number of facies for the first variable
[in]	nfac2max	Maximum number of facies for the second variable
[out]	nclout	Total number of classes

Remarks: This procedure is designed for one or two Z-variables; When used for 2 variables, the proportions are for pairs of; variables

◆ db_sample_alloc()

double* db_sample_alloc	(	const Db *	db,
		const ELoc &	locatorType
	)

Allocates the array for storing a sample

Returns: A pointer to the allocated array

Parameters

[in]	db	Db descriptor
[in]	locatorType	vector type (ELoc)

Remarks: The allocated array must be freed using db_sample_free(); A fatal error occurs if the core allocation fails.

◆ db_sample_free()

double* db_sample_free ( double * tab )

Frees the array for storing a sample

Returns: A pointer to the array to be freed

Parameters

[in] tab Sample array to be freed

◆ db_sample_print()

void db_sample_print	(	Db *	db,
		int	iech,
		int	flag_ndim,
		int	flag_nvar,
		int	flag_nerr
	)

Print a sample

Parameters

[in]	db	Db structure
[in]	iech	Rank of the sample
[in]	flag_ndim	1 if the coordinates must be printed
[in]	flag_nvar	1 if the variables must be printed
[in]	flag_nerr	1 if the error measurement variance must be printed

◆ db_selref()

int db_selref	(	int	ndim,
		const int *	nx,
		const int *	ref,
		const double *	tabin,
		double *	tabout
	)

Select a 2-D grid from a Grid Db

Returns: Error returned code

Parameters

[in]	ndim	Space dimension of the Db
[in]	nx	Number of meshes of the grid
[in]	ref	Array giving the indices of the reference planes Grid indices are counted starting from 1 -1 for the directions of interest
[in]	tabin	Input array
[out]	tabout	Output array

◆ distance_grid()

double distance_grid	(	DbGrid *	db,
		int	flag_moins1,
		int	iech1,
		int	iech2,
		double *	dist_vect
	)

Calculates the distance between two points in the same grid Db

Returns: The calculated distance

Parameters

[in]	db	Db structure
[in]	flag_moins1	1 to use the cell center
[in]	iech1	rank of the first sample
[in]	iech2	rank of the second sample
[out]	dist_vect	If the output vector is provided. Returns the distance as a vector

◆ distance_inter()

double distance_inter	(	const Db *	db1,
		const Db *	db2,
		int	iech1,
		int	iech2,
		double *	dist_vect
	)

Calculates the distance between two points

Returns: The calculated distance or TEST if one coordinate not defined

Parameters

[in]	db1	Db structure for the first sample
[in]	db2	Db structure for the second sample
[in]	iech1	rank of the first sample
[in]	iech2	rank of the second sample
[out]	dist_vect	If the output vector is provided. Returns the distance as a vector

Remarks: : If both Data Bases do not share the same Space Dimension; : the test is performed on their minimum.; : The returned array 'vect' must be dimension to that value

◆ distance_intra()

double distance_intra	(	const Db *	db,
		int	iech1,
		int	iech2,
		double *	dist_vect
	)

Calculates the distance between two points in the same Db

Returns: The calculated distance or TEST if one coordinate not defined

Parameters

[in]	db	Db structure
[in]	iech1	rank of the first sample
[in]	iech2	rank of the second sample
[out]	dist_vect	If the output vector is provided. Returns the distance as a vector

◆ get_grid_value()

double get_grid_value	(	DbGrid *	dbgrid,
		int	iptr,
		VectorInt &	indg,
		int	ix,
		int	iy,
		int	iz
	)

Returns a value from the a 3-D (maximum) grid

Returns: Returned value

Parameters

[in]	dbgrid	Db Grid structure
[in]	iptr	Rank of the column
[in]	indg	Working index array (Dimension: get_NDIM(dbgrid))
[in]	ix	Rank of the node along first dimension
[in]	iy	Rank of the node along second dimension
[in]	iz	Rank of the node along third dimension

◆ get_LOCATOR_NITEM()

int get_LOCATOR_NITEM	(	const Db *	db,
		const ELoc &	locatorType
	)

Returns the number of items for a given locator in the Db

Returns: Number of items

Parameters

[in]	db	Db descriptor
[in]	locatorType	Rank of the pointer (ELoc)

◆ grid_iterator_init()

void grid_iterator_init	(	Grid *	grid,
		const VectorInt &	order
	)

Initialize the Grid iterator

◆ grid_iterator_next()

VectorInt grid_iterator_next ( Grid * grid )

Returns 1 when the last element of the iteration is reached

Increment the Grid iterator

◆ index_point_to_grid()

int index_point_to_grid	(	const Db *	dbin,
		int	iech,
		int	flag_outside,
		const DbGrid *	dbout,
		double *	coor
	)

Find the index of the output grid file which is the closest to the sample of the input file

Returns: Error return code; >= index of the grid node; -1 if the point is outside the grid

Parameters

[in]	dbin	descriptor of the input file
[in]	iech	Index of the data point
[in]	flag_outside	value returned for the point outside the grid 1 the index is set to the closest grid node 0 the index is set to -1 -1 do not correct the index
[in]	dbout	descriptor of the output grid file
[out]	coor	Working array (dimension: ndim)

◆ is_grid_multiple()

int is_grid_multiple	(	DbGrid *	db1,
		DbGrid *	db2
	)

Check if a grid is a multiple of the other grid

Returns: 1 if the two grid are multiple; 0 otherwise

Parameters

[in]	db1	Db1 small grid structure
[in]	db2	Db1 coarse grid structure

◆ point_inside_grid()

int point_inside_grid	(	Db *	db,
		int	iech,
		const DbGrid *	dbgrid
	)

Check if a sample of a Db lies within a Grid Db

Returns: 1 if the point lies within the grid; 0 otherwise

Parameters

[in]	db	Db structure
[in]	iech	Rank of the target sample
[in]	dbgrid	Grid Db structure

◆ point_to_bench()

int point_to_bench	(	const DbGrid *	db,
		double *	coor,
		int	flag_outside,
		int *	indb
	)

Converts from point coordinates to index of the bench to which it belongs

Returns: Error return code; 0 if the point is inside the grid; 1 if the point is outside the grid; -1 if one coordinate is undefined; -2 if the grid is defined in space less than 3D

Parameters

[in]	db	descriptor of the grid parameters
[in]	coor	array of coordinates of the point
[in]	flag_outside	value returned for the point outside the grid 1 the index is set to the closest grid node 0 the index is set to -1 -1 do not correct the index
[out]	indb	index of the bench

Remarks: The bench corresponds to the third dimension of the grid provided; as reference

◆ point_to_grid()

int point_to_grid	(	const DbGrid *	db,
		const double *	coor,
		int	flag_outside,
		int *	indg
	)

Converts from point coordinates to nearest grid node indices

Returns: Error return code; 0 if the point is inside the grid; 1 if the point is outside the grid; -1 if one coordinate is undefined

Parameters

[in]	db	descriptor of the grid parameters
[in]	coor	array of coordinates of the point
[in]	flag_outside	value returned for the point outside the grid 1 the index is set to the closest grid node 0 the index is set to -1 -1 do not correct the index
[out]	indg	indices of the closest grid node

◆ point_to_point()

int point_to_point	(	Db *	db,
		const double *	coor
	)

Returns the rank of the closest isolated point

Returns: Output index

Parameters

[in]	db	descriptor of the Db
[in]	coor	array of coordinates of the point

◆ set_grid_value()

void set_grid_value	(	DbGrid *	dbgrid,
		int	iptr,
		VectorInt &	indg,
		int	ix,
		int	iy,
		int	iz,
		double	value
	)

Set the value in the a 3-D (maximum) grid

Parameters

[in]	dbgrid	Db Grid structure
[in]	iptr	Rank of the column
[in]	indg	Working index array (Dimension: get_NDIM(dbgrid))
[in]	ix	Rank of the node along first dimension
[in]	iy	Rank of the node along second dimension
[in]	iz	Rank of the node along third dimension
[in]	value	Assigned value

◆ st_rotate()

static void st_rotate	(	int	ndim,
		double *	rotmat,
		VectorDouble &	coor,
		VectorDouble &	mini,
		VectorDouble &	maxi
	)

static

Update the bounding box limits due to a point

Parameters

[in]	ndim	Space dimension
[in]	rotmat	Rotation matrix (optional)
[in]	coor	Array of sample coordinates
[out]	mini	Array containing the minimum
[out]	maxi	Array containing the maximum

Remarks: This function is limited to the 3-D case

Functions

Function Documentation

◆ compat_NDIM()

◆ db_attribute_copy()

◆ db_attribute_range()

◆ db_center()

◆ db_count_defined()

◆ db_create_grid_dilate()

◆ db_create_grid_divider()

◆ db_create_grid_multiple()

◆ db_distances_general()

◆ db_extension_rotated()

◆ db_get_attribute()

◆ db_get_grid_axis()

◆ db_get_rank_absolute_to_relative()

◆ db_get_rank_relative_to_absolute()

◆ db_gradient_component_to_modang()

◆ db_gradient_modang_to_component()

◆ db_gradient_update()

◆ db_grid_copy()

◆ db_grid_copy_dilate()

◆ db_grid_define_coordinates()

◆ db_grid_patch()

◆ db_grid_print()

◆ db_grid_reduce()

◆ db_identify_variables_by_name()

◆ db_is_isotropic()

◆ db_locator_attribute_add()

◆ db_locators_correct()

◆ db_merge()

◆ db_monostat()

◆ db_name_identify()

◆ db_prop_read()

◆ db_prop_write()

◆ db_proportion()

◆ db_sample_alloc()

◆ db_sample_free()

◆ db_sample_print()

◆ db_selref()

◆ distance_grid()

◆ distance_inter()

◆ distance_intra()

◆ get_grid_value()

◆ get_LOCATOR_NITEM()

◆ grid_iterator_init()

◆ grid_iterator_next()

◆ index_point_to_grid()

◆ is_grid_multiple()

◆ point_inside_grid()

◆ point_to_bench()

◆ point_to_grid()

◆ point_to_point()

◆ set_grid_value()

◆ st_rotate()