Computes the volume or area of the element, as appropriate. This metric uses the average of the corner Jacobian determinants for the approximation to the volume of hex. More...

#include <CornerJacobianQualityMetric.hpp>

Inheritance diagram for CornerJacobianQualityMetric:

[legend]

Collaboration diagram for CornerJacobianQualityMetric:

[legend]

Public Member Functions
	CornerJacobianQualityMetric ()

virtual	~CornerJacobianQualityMetric ()
	virtual destructor ensures use of polymorphism during destruction More...

bool	evaluate_element (PatchData &pd, MsqMeshEntity *element, double &fval, MsqError &err)
	evaluate using mesquite objects More...

	CornerJacobianQualityMetric ()

virtual	~CornerJacobianQualityMetric ()
	virtual destructor ensures use of polymorphism during destruction More...

bool	evaluate_element (PatchData &pd, MsqMeshEntity *element, double &fval, MsqError &err)
	evaluate using mesquite objects More...

Public Member Functions inherited from VolumeQualityMetric
virtual	~VolumeQualityMetric ()

virtual	~VolumeQualityMetric ()

Public Member Functions inherited from QualityMetric
virtual	~QualityMetric ()

MetricType	get_metric_type ()

void	set_element_evaluation_mode (ElementEvaluationMode mode, MsqError &err)
	Sets the evaluation mode for the ELEMENT_BASED metrics. More...

ElementEvaluationMode	get_element_evaluation_mode ()
	Returns the evaluation mode for the metric. More...

void	set_averaging_method (AveragingMethod method, MsqError &err)

void	set_feasible_constraint (int alpha)

int	get_feasible_constraint ()
	Returns the feasible flag for this metric. More...

void	set_name (msq_std::string st)
	Sets the name of this metric. More...

msq_std::string	get_name ()
	Returns the name of this metric (as a string). More...

double	vertex_barrier_function (double det, double delta)
	Escobar Barrier Function for Shape and Other Metrics. More...

virtual bool	evaluate_vertex (PatchData &, MsqVertex *, double &, MsqError &err)
	Evaluate the metric for a vertex. More...

void	set_gradient_type (GRADIENT_TYPE grad)
	Sets gradType for this metric. More...

void	set_hessian_type (HESSIAN_TYPE ht)
	Sets hessianType for this metric. More...

bool	compute_vertex_gradient (PatchData &pd, MsqVertex &vertex, MsqVertex *vertices[], Vector3D grad_vec[], int num_vtx, double &metric_value, MsqError &err)
	Calls compute_vertex_numerical_gradient if gradType equals NUMERCIAL_GRADIENT. Calls compute_vertex_analytical_gradient if gradType equals ANALYTICAL_GRADIENT;. More...

bool	compute_element_gradient (PatchData &pd, MsqMeshEntity element, MsqVertex free_vtces[], Vector3D grad_vec[], int num_free_vtx, double &metric_value, MsqError &err)
	For MetricType == ELEMENT_BASED. Calls either compute_element_numerical_gradient() or compute_element_analytical_gradient() for gradType equal NUMERICAL_GRADIENT or ANALYTICAL_GRADIENT, respectively. More...

bool	compute_element_gradient_expanded (PatchData &pd, MsqMeshEntity element, MsqVertex free_vtces[], Vector3D grad_vec[], int num_free_vtx, double &metric_value, MsqError &err)

bool	compute_element_hessian (PatchData &pd, MsqMeshEntity element, MsqVertex free_vtces[], Vector3D grad_vec[], Matrix3D hessian[], int num_free_vtx, double &metric_value, MsqError &err)
	For MetricType == ELEMENT_BASED. Calls either compute_element_numerical_hessian() or compute_element_analytical_hessian() for hessianType equal NUMERICAL_HESSIAN or ANALYTICAL_HESSIAN, respectively. More...

void	set_negate_flag (int neg)

int	get_negate_flag ()
	Returns negateFlag. More...

virtual void	change_metric_type (MetricType t, MsqError &err)

virtual	~QualityMetric ()

MetricType	get_metric_type ()

void	set_element_evaluation_mode (ElementEvaluationMode mode, MsqError &err)
	Sets the evaluation mode for the ELEMENT_BASED metrics. More...

ElementEvaluationMode	get_element_evaluation_mode ()
	Returns the evaluation mode for the metric. More...

void	set_averaging_method (AveragingMethod method, MsqError &err)

void	set_feasible_constraint (int alpha)

int	get_feasible_constraint ()
	Returns the feasible flag for this metric. More...

void	set_name (msq_std::string st)
	Sets the name of this metric. More...

msq_std::string	get_name ()
	Returns the name of this metric (as a string). More...

double	vertex_barrier_function (double det, double delta)
	Escobar Barrier Function for Shape and Other Metrics. More...

virtual bool	evaluate_vertex (PatchData &, MsqVertex *, double &, MsqError &err)
	Evaluate the metric for a vertex. More...

void	set_gradient_type (GRADIENT_TYPE grad)
	Sets gradType for this metric. More...

void	set_hessian_type (HESSIAN_TYPE ht)
	Sets hessianType for this metric. More...

bool	compute_vertex_gradient (PatchData &pd, MsqVertex &vertex, MsqVertex *vertices[], Vector3D grad_vec[], int num_vtx, double &metric_value, MsqError &err)

bool	compute_element_gradient (PatchData &pd, MsqMeshEntity element, MsqVertex free_vtces[], Vector3D grad_vec[], int num_free_vtx, double &metric_value, MsqError &err)
	For MetricType == ELEMENT_BASED. Calls either compute_element_numerical_gradient() or compute_element_analytical_gradient() for gradType equal NUMERICAL_GRADIENT or ANALYTICAL_GRADIENT, respectively. More...

bool	compute_element_gradient_expanded (PatchData &pd, MsqMeshEntity element, MsqVertex free_vtces[], Vector3D grad_vec[], int num_free_vtx, double &metric_value, MsqError &err)

bool	compute_element_hessian (PatchData &pd, MsqMeshEntity element, MsqVertex free_vtces[], Vector3D grad_vec[], Matrix3D hessian[], int num_free_vtx, double &metric_value, MsqError &err)
	For MetricType == ELEMENT_BASED. Calls either compute_element_numerical_hessian() or compute_element_analytical_hessian() for hessianType equal NUMERICAL_HESSIAN or ANALYTICAL_HESSIAN, respectively. More...

void	set_negate_flag (int neg)

int	get_negate_flag ()
	Returns negateFlag. More...

virtual void	change_metric_type (MetricType t, MsqError &err)

Additional Inherited Members
Public Types inherited from QualityMetric
enum	MetricType { MT_UNDEFINED, VERTEX_BASED, ELEMENT_BASED, VERTEX_BASED_FREE_ONLY, MT_UNDEFINED, VERTEX_BASED, ELEMENT_BASED, VERTEX_BASED_FREE_ONLY }

enum	ElementEvaluationMode { EEM_UNDEFINED, ELEMENT_VERTICES, LINEAR_GAUSS_POINTS, QUADRATIC_GAUSS_POINTS, CUBIC_GAUSS_POINTS, EEM_UNDEFINED, ELEMENT_VERTICES, LINEAR_GAUSS_POINTS, QUADRATIC_GAUSS_POINTS, CUBIC_GAUSS_POINTS }

enum	AveragingMethod { NONE, LINEAR, RMS, HMS, MINIMUM, MAXIMUM, HARMONIC, GEOMETRIC, SUM, SUM_SQUARED, GENERALIZED_MEAN, STANDARD_DEVIATION, MAX_OVER_MIN, MAX_MINUS_MIN, SUM_OF_RATIOS_SQUARED, NONE, LINEAR, RMS, HMS, MINIMUM, MAXIMUM, HARMONIC, GEOMETRIC, SUM, SUM_SQUARED, GENERALIZED_MEAN, STANDARD_DEVIATION, MAX_OVER_MIN, MAX_MINUS_MIN, SUM_OF_RATIOS_SQUARED }

enum	GRADIENT_TYPE { NUMERICAL_GRADIENT, ANALYTICAL_GRADIENT, NUMERICAL_GRADIENT, ANALYTICAL_GRADIENT }

enum	HESSIAN_TYPE { NUMERICAL_HESSIAN, ANALYTICAL_HESSIAN, NUMERICAL_HESSIAN, ANALYTICAL_HESSIAN }

enum	MetricType { MT_UNDEFINED, VERTEX_BASED, ELEMENT_BASED, VERTEX_BASED_FREE_ONLY, MT_UNDEFINED, VERTEX_BASED, ELEMENT_BASED, VERTEX_BASED_FREE_ONLY }

enum	ElementEvaluationMode { EEM_UNDEFINED, ELEMENT_VERTICES, LINEAR_GAUSS_POINTS, QUADRATIC_GAUSS_POINTS, CUBIC_GAUSS_POINTS, EEM_UNDEFINED, ELEMENT_VERTICES, LINEAR_GAUSS_POINTS, QUADRATIC_GAUSS_POINTS, CUBIC_GAUSS_POINTS }

enum	AveragingMethod { NONE, LINEAR, RMS, HMS, MINIMUM, MAXIMUM, HARMONIC, GEOMETRIC, SUM, SUM_SQUARED, GENERALIZED_MEAN, STANDARD_DEVIATION, MAX_OVER_MIN, MAX_MINUS_MIN, SUM_OF_RATIOS_SQUARED, NONE, LINEAR, RMS, HMS, MINIMUM, MAXIMUM, HARMONIC, GEOMETRIC, SUM, SUM_SQUARED, GENERALIZED_MEAN, STANDARD_DEVIATION, MAX_OVER_MIN, MAX_MINUS_MIN, SUM_OF_RATIOS_SQUARED }

enum	GRADIENT_TYPE { NUMERICAL_GRADIENT, ANALYTICAL_GRADIENT, NUMERICAL_GRADIENT, ANALYTICAL_GRADIENT }

enum	HESSIAN_TYPE { NUMERICAL_HESSIAN, ANALYTICAL_HESSIAN, NUMERICAL_HESSIAN, ANALYTICAL_HESSIAN }

Protected Member Functions inherited from QualityMetric
	QualityMetric ()

void	set_metric_type (MetricType t)
	This function should be used in the constructor of every concrete quality metric. More...

double	average_metrics (const double metric_values[], const int &num_values, MsqError &err)
	average_metrics takes an array of length num_values and averages the contents using averaging method data member avgMethod . More...

double	average_metric_and_weights (double metric_values[], int num_metric_values, MsqError &err)
	Given a list of metric values, calculate the average metric valude according to the current avgMethod and write into the passed metric_values array the the value weight/count to use when averaging gradient vectors for the metric. More...

double	weighted_average_metrics (const double coef[], const double metric_values[], const int &num_values, MsqError &err)
	takes an array of coefficients and an array of metrics (both of length num_value) and averages the contents using averaging method 'method'. More...

bool	compute_vertex_numerical_gradient (PatchData &pd, MsqVertex &vertex, MsqVertex *vertices[], Vector3D grad_vec[], int num_vtx, double &metric_value, MsqError &err)

bool	compute_element_numerical_gradient (PatchData &pd, MsqMeshEntity element, MsqVertex free_vtces[], Vector3D grad_vec[], int num_free_vtx, double &metric_value, MsqError &err)
	Non-virtual function which numerically computes the gradient of a QualityMetric of a given element for a given set of free vertices on that element. This is used by metric which mType is ELEMENT_BASED. For parameters, see compute_element_gradient() . More...

virtual bool	compute_vertex_analytical_gradient (PatchData &pd, MsqVertex &vertex, MsqVertex *vertices[], Vector3D grad_vec[], int num_vtx, double &metric_value, MsqError &err)
	Virtual function that computes the gradient of the QualityMetric analytically. The base class implementation of this function simply prints a warning and calls compute_numerical_gradient to calculate the gradient. This is used by metric which mType is VERTEX_BASED. More...

virtual bool	compute_element_analytical_gradient (PatchData &pd, MsqMeshEntity element, MsqVertex free_vtces[], Vector3D grad_vec[], int num_free_vtx, double &metric_value, MsqError &err)
	Virtual function that computes the gradient of the QualityMetric analytically. The base class implementation of this function simply prints a warning and calls compute_numerical_gradient to calculate the gradient. This is used by metric which mType is ELEMENT_BASED. For parameters, see compute_element_gradient() . More...

bool	compute_element_numerical_hessian (PatchData &pd, MsqMeshEntity element, MsqVertex free_vtces[], Vector3D grad_vec[], Matrix3D hessian[], int num_free_vtx, double &metric_value, MsqError &err)

virtual bool	compute_element_analytical_hessian (PatchData &pd, MsqMeshEntity element, MsqVertex free_vtces[], Vector3D grad_vec[], Matrix3D hessian[], int num_free_vtx, double &metric_value, MsqError &err)

	QualityMetric ()

void	set_metric_type (MetricType t)
	This function should be used in the constructor of every concrete quality metric. More...

double	average_metrics (const double metric_values[], const int &num_values, MsqError &err)
	average_metrics takes an array of length num_values and averages the contents using averaging method data member avgMethod . More...

double	average_metric_and_weights (double metric_values[], int num_metric_values, MsqError &err)
	Given a list of metric values, calculate the average metric valude according to the current avgMethod and write into the passed metric_values array the the value weight/count to use when averaging gradient vectors for the metric. More...

double	weighted_average_metrics (const double coef[], const double metric_values[], const int &num_values, MsqError &err)
	takes an array of coefficients and an array of metrics (both of length num_value) and averages the contents using averaging method 'method'. More...

bool	compute_vertex_numerical_gradient (PatchData &pd, MsqVertex &vertex, MsqVertex *vertices[], Vector3D grad_vec[], int num_vtx, double &metric_value, MsqError &err)

bool	compute_element_numerical_gradient (PatchData &pd, MsqMeshEntity element, MsqVertex free_vtces[], Vector3D grad_vec[], int num_free_vtx, double &metric_value, MsqError &err)
	Non-virtual function which numerically computes the gradient of a QualityMetric of a given element for a given set of free vertices on that element. This is used by metric which mType is ELEMENT_BASED. For parameters, see compute_element_gradient() . More...

virtual bool	compute_vertex_analytical_gradient (PatchData &pd, MsqVertex &vertex, MsqVertex *vertices[], Vector3D grad_vec[], int num_vtx, double &metric_value, MsqError &err)
	Virtual function that computes the gradient of the QualityMetric analytically. The base class implementation of this function simply prints a warning and calls compute_numerical_gradient to calculate the gradient. This is used by metric which mType is VERTEX_BASED. More...

virtual bool	compute_element_analytical_gradient (PatchData &pd, MsqMeshEntity element, MsqVertex free_vtces[], Vector3D grad_vec[], int num_free_vtx, double &metric_value, MsqError &err)
	Virtual function that computes the gradient of the QualityMetric analytically. The base class implementation of this function simply prints a warning and calls compute_numerical_gradient to calculate the gradient. This is used by metric which mType is ELEMENT_BASED. For parameters, see compute_element_gradient() . More...

bool	compute_element_numerical_hessian (PatchData &pd, MsqMeshEntity element, MsqVertex free_vtces[], Vector3D grad_vec[], Matrix3D hessian[], int num_free_vtx, double &metric_value, MsqError &err)

virtual bool	compute_element_analytical_hessian (PatchData &pd, MsqMeshEntity element, MsqVertex free_vtces[], Vector3D grad_vec[], Matrix3D hessian[], int num_free_vtx, double &metric_value, MsqError &err)

Protected Attributes inherited from QualityMetric
AveragingMethod	avgMethod

int	feasible

msq_std::string	metricName

Detailed Description

Computes the volume or area of the element, as appropriate. This metric uses the average of the corner Jacobian determinants for the approximation to the volume of hex.

Todo:: MB: CornerJacobianQualityMetric is currently not being evaluated using the corner jacobian method (instead it is using the area functions from MsqMeshEntity... needs to be modified to do so.

Todo:: MB: CornerJacobianQualityMetric is currently not being evaluated using the corner jacobian method (instead it is using the area functions from MsqMeshEntity... needs to be modified to do so.

The metric does not use the sample point functionality or the compute_weighted_jacobian (except for possibly indirectly when evaluating the metric value for a hex). It evaluates the signed area of surface elements and the signed volume of volume elements. It does require a feasible region, and (in general) the metric needs to be minimized.

Definition at line 63 of file includeLinks/CornerJacobianQualityMetric.hpp.

Constructor & Destructor Documentation

CornerJacobianQualityMetric ( )

Definition at line 48 of file QualityMetric/Volume/CornerJacobianQualityMetric.cpp.

References QualityMetric::avgMethod, QualityMetric::ELEMENT_BASED, QualityMetric::ELEMENT_VERTICES, QualityMetric::feasible, MSQ_CHKERR, QualityMetric::NONE, QualityMetric::set_element_evaluation_mode(), QualityMetric::set_metric_type(), and QualityMetric::set_name().

 {
   MsqError err;
   set_metric_type(ELEMENT_BASED);
   set_element_evaluation_mode(ELEMENT_VERTICES, err); MSQ_CHKERR(err);
   avgMethod=QualityMetric::NONE;
   feasible=0;
   set_name("Corner Jacobian Volume (or Area)");
 }

Here is the call graph for this function:

virtual ~CornerJacobianQualityMetric ( )

inlinevirtual

virtual destructor ensures use of polymorphism during destruction

Definition at line 70 of file includeLinks/CornerJacobianQualityMetric.hpp.

71 {}

CornerJacobianQualityMetric ( )

virtual ~CornerJacobianQualityMetric ( )

inlinevirtual

virtual destructor ensures use of polymorphism during destruction

Definition at line 70 of file src/QualityMetric/Volume/CornerJacobianQualityMetric.hpp.

71 {}

Member Function Documentation

bool evaluate_element	(	PatchData &	pd,
		MsqMeshEntity *	element,
		double &	fval,
		MsqError &	err
	)

virtual

evaluate using mesquite objects

Reimplemented from QualityMetric.

Definition at line 58 of file QualityMetric/Volume/CornerJacobianQualityMetric.cpp.

References MsqMeshEntity::compute_unsigned_area(), MsqMeshEntity::compute_unsigned_volume(), MsqMeshEntity::get_element_type(), Mesquite::HEXAHEDRON, MSQ_ERRZERO, Mesquite::MSQ_MAX_CAP, Mesquite::QUADRILATERAL, Mesquite::TETRAHEDRON, and Mesquite::TRIANGLE.

 {
   switch(element->get_element_type()){
     case TRIANGLE:
     case QUADRILATERAL:
       fval=element->compute_unsigned_area(pd,err);  MSQ_ERRZERO(err);
       break;
     case TETRAHEDRON:
     case HEXAHEDRON:
       fval=element->compute_unsigned_volume(pd,err);  MSQ_ERRZERO(err);
       break;
     default:
       fval=MSQ_MAX_CAP;
   }// end switch over element type
   return true;
 }