thresh.cpp File Reference

#include "geoslib_old_f.h"
#include "Basic/Utilities.hpp"
#include "Basic/String.hpp"
#include "Basic/OptDbg.hpp"
#include "Basic/VectorHelper.hpp"
#include "LithoRule/RuleProp.hpp"
#include "LithoRule/PropDef.hpp"
#include "LithoRule/Rule.hpp"
#include "LithoRule/RuleShift.hpp"
#include "LithoRule/RuleShadow.hpp"
#include "Db/Db.hpp"
#include "Db/DbGrid.hpp"
#include "Variogram/Vario.hpp"
#include "Model/Model.hpp"

Functions
Rule *	rule_free (const Rule *rule)
void	proportion_rule_process (PropDef *propdef, const EProcessOper &mode)
void	proportion_print (PropDef *propdef)
int	rule_thresh_define_shadow (PropDef *propdef, Db *db, const RuleShadow *rule, int facies, int iech, int isimu, int nbsimu, double *t1min, double *t1max, double *t2min, double *t2max, double *sh_dsup, double *sh_down)
int	rule_thresh_define (PropDef *propdef, Db *db, const Rule *rule, int facies, int iech, int isimu, int nbsimu, int flag_check, double *t1min, double *t1max, double *t2min, double *t2max)
int	db_rule_shadow (Db *db, Db *dbprop, RuleShadow *rule, Model *model, const VectorDouble &props, int flag_stat, int nfacies)
int	_db_rule (Db *db, const RuleProp *ruleprop, Model *model, const NamingConvention &namconv)
int	db_bounds_shadow (Db *db, Db *dbprop, RuleShadow *rule, Model *model, const VectorDouble &props, int flag_stat, int nfacies)
int	_db_bounds (Db *db, const RuleProp *ruleprop, Model *model, const NamingConvention &namconv)
void	propdef_reset (PropDef *propdef)
PropDef *	proportion_manage (int mode, int flag_facies, int flag_stat, int ngrf1, int ngrf2, int nfac1, int nfac2, Db *db, const Db *dbprop, const VectorDouble &propcst, PropDef *proploc)
int	_db_threshold (Db *db, const RuleProp *ruleprop, Model *model, const NamingConvention &namconv)
Model *	model_rule_combine (const Model *model1, const Model *model2, const Rule *rule)

Function Documentation

int _db_bounds	(	Db *	db,
		const RuleProp *	ruleprop,
		Model *	model,
		const NamingConvention &	namconv
	)

Apply the Rule transformation to derive the bounds variables for a Db

Returns

Error return code

Parameters

[in]	db	Db structure
[in]	ruleprop	RuleProp structure
[in]	model	First Model structure (only for SHIFT)
[in]	namconv	Naming convention

int _db_rule	(	Db *	db,
		const RuleProp *	ruleprop,
		Model *	model,
		const NamingConvention &	namconv
	)

Apply the Rule transformation to convert a set of Gaussian vectors into the corresponding Facies in a Db

Returns

Error return code

Parameters

[in]	db	Output Db structure
[in]	ruleprop	RuleProp structure
[in]	model	First Model structure (only for SHIFT)
[in]	namconv	Naming convention

Remarks

The input variable must be locatorized as Z or ELoc::SIMU

int _db_threshold	(	Db *	db,
		const RuleProp *	ruleprop,
		Model *	model,
		const NamingConvention &	namconv
	)

Calculate all the thresholds at each sample of a Db

Returns

Error return code

Parameters

[in]	db	Db structure
[in]	ruleprop	RuleProp structure
[in]	model	First Model structure (only for SHIFT)
[in]	namconv	Naming Convention

int db_bounds_shadow	(	Db *	db,
		Db *	dbprop,
		RuleShadow *	rule,
		Model *	model,
		const VectorDouble &	props,
		int	flag_stat,
		int	nfacies
	)

Apply the Rule transformation to derive the bounds variables for a Db (Shadow case)

Returns

Error return code

Parameters

[in]	db	Db structure
[in]	dbprop	Db structure used for proportions (non-stationary case)
[in]	rule	Lithotype Rule definition
[in]	model	First Model structure (only for SHIFT)
[in]	props	Array of proportions for the facies
[in]	flag_stat	1 for stationary; 0 otherwise
[in]	nfacies	Number of facies

int db_rule_shadow	(	Db *	db,
		Db *	dbprop,
		RuleShadow *	rule,
		Model *	model,
		const VectorDouble &	props,
		int	flag_stat,
		int	nfacies
	)

Apply the Rule transformation to the GRFs of a Db (Shadow case)

Returns

Error return code

Parameters

[in]	db	Output Db structure
[in]	dbprop	Db structure used for proportions (non-stationary case)
[in]	rule	Lithotype Rule definition
[in]	model	First Model structure (only for SHIFT)
[in]	props	Array of proportions for the facies
[in]	flag_stat	1 for stationary; 0 otherwise
[in]	nfacies	Number of facies

Remarks

The input variable must be locatorized as Z or ELoc::SIMU

It will be changed in this function to locator ELoc::SIMU

Model* model_rule_combine	(	const Model *	model1,
		const Model *	model2,
		const Rule *	rule
	)

Combine two basic models into a bivariate model (residuals model)

Returns

The newly Model structure

Parameters

[in]	model1	First input Model
[in]	model2	Second input Model
[in]	rule	Rule

Remarks

: The drift is not copied into the new model

void propdef_reset

(

PropDef *

propdef

)

Set memory proportion so as to provoke the update at first usage

Parameters

[in]

propdef

Pointer to Propdef structure

PropDef* proportion_manage	(	int	mode,
		int	flag_facies,
		int	flag_stat,
		int	ngrf1,
		int	ngrf2,
		int	nfac1,
		int	nfac2,
		Db *	db,
		const Db *	dbprop,
		const VectorDouble &	propcst,
		PropDef *	proploc
	)

Allocate or deallocate a proportion array

Returns

Pointer on the returned PropDef structure

Parameters

[in]	mode	1 for allocation; -1 for deallocation
[in]	flag_facies	1 if Gibbs is used for facies
[in]	flag_stat	1 if the proportions are stationary
[in]	ngrf1	Number of GRFs for the first PGS
[in]	ngrf2	Number of GRFs for the second PGS
[in]	nfac1	Number of facies for the first PGS
[in]	nfac2	Number of facies for the second PGS
[in]	db	Db structure containing the data
[in]	dbprop	Db structure containing the proportions (only used in the non-stationary case)
[in]	propcst	Constant set of proportions (used if flag_stat)
[in]	proploc	PropDef structure (used for mode<0)

void proportion_print

(

PropDef *

propdef

)

Print the (non-stationary) proportions

Parameters

[in]

propdef

PropDef structure

void proportion_rule_process	(	PropDef *	propdef,
		const EProcessOper &	mode
	)

Set the method to compute Proportions

Parameters

[in]	propdef	PropDef structure
[in]	mode	Type of operation (EProcessOper)

Rule* rule_free

(

const Rule *

rule

)

Free a Rule structure

Returns

Pointer to the newly freed Rule structure

Parameters

[in]

rule

Rule structure to be freed

int rule_thresh_define	(	PropDef *	propdef,
		Db *	db,
		const Rule *	rule,
		int	facies,
		int	iech,
		int	isimu,
		int	nbsimu,
		int	flag_check,
		double *	t1min,
		double *	t1max,
		double *	t2min,
		double *	t2max
	)

Set the (non-stationary) proportions and define thresholds

Returns

Error return code

Parameters

[in]	propdef	PropDef structure
[in]	db	Db input structure
[in]	rule	Rule structure
[in]	facies	Facies of interest (or ITEST) starting from 1
[in]	iech	Rank of the data in the input Db
[in]	isimu	Rank of the simulation (EProcessOper::CONDITIONAL)
[in]	nbsimu	Number of simulations
[in]	flag_check	1 if the consistency check with the actual proportion of the current facies must be done
[out]	t1min	Minimum threshold for Y1
[out]	t1max	Maximum threshold for Y1
[out]	t2min	Minimum threshold for Y2
[out]	t2max	Maximum threshold for Y2

int rule_thresh_define_shadow	(	PropDef *	propdef,
		Db *	db,
		const RuleShadow *	rule,
		int	facies,
		int	iech,
		int	isimu,
		int	nbsimu,
		double *	t1min,
		double *	t1max,
		double *	t2min,
		double *	t2max,
		double *	sh_dsup,
		double *	sh_down
	)

Set the (non-stationary) proportions and define thresholds (for shadow only)

Returns

Error return code

Parameters

[in]	propdef	PropDef structure
[in]	db	Db input structure
[in]	rule	Rule structure
[in]	facies	Facies of interest (or GV_ITEST)
[in]	iech	Rank of the data in the input Db
[in]	isimu	Rank of the simulation (EProcessOper::CONDITIONAL)
[in]	nbsimu	Number of simulations (EProcessOper::CONDITIONAL)
[out]	t1min	Minimum threshold for Y1
[out]	t1max	Maximum threshold for Y1
[out]	t2min	Minimum threshold for Y2
[out]	t2max	Maximum threshold for Y2
[out]	sh_dsup	Local or global upwards shift (shadow)
[out]	sh_down	Local or global downwards shift (shadow)