A section evaluator using quadratic interpolation over a 1D domain. More...

#include <quadratic_1d.h>

Inheritance diagram for fiber_bundle::quadratic_1d:

QUADRATIC_1D FACET
	quadratic_1d ()
	Default constructor. More...

	quadratic_1d (const quadratic_1d &xother)
	Copy constructor. More...

virtual	~quadratic_1d ()
	Destructor. More...

LINEAR_FCN_SPACE FACET
Todo: Fix!
virtual int	dl () const
	The dimension of this function space. More...

virtual void	basis_at_coord (const dof_type xlocal_coord[], size_type xlocal_coord_ub)
	Computes the value of each basis function at local coordinates xlocal_coord. More...

virtual void	basis_derivs_at_coord (const dof_type xlocal_coords[], size_type xlocal_coords_ub)
	Computes the value of the derivatives of each basis function at local coordinates xlocal_coords. More...

INTEGRABLE_SECTION_EVALUATOR FACET
virtual value_type	jacobian_determinant (const dof_type xcoord_dofs[], size_type xcoord_dofs_ub, size_type xdf, const coord_type xlocal_coords[], size_type xlocal_coords_ub)
	Computes the the determinant of the jacobian matrix at local coordinates xlocal_coords with coordinate dofs xcoord_dofs. More...

virtual value_type	volume (const dof_type xcoord_dofs[], size_type xcoord_dofs_ub, size_type xdf)
	Volume for specified coordinate dofs xcoord_dofs and fiber space dimension xdf. More...

virtual void	integrate (const dof_type xcoord_dofs[], size_type xcoord_dofs_ub, size_type xdf, const dof_type xintegrands[], size_type xintegrands_ub, value_type xresult_integrals[], size_type xresult_integrals_ub)
	Computes the value of the integral of the integrand array... More...

virtual void	integrate (const dof_type xcoord_dofs[], size_type xcoord_dofs_ub, size_type xdf, const dof_type xintegrand, value_type xresult_integrals[], size_type xresult_integrals_ub)
	Computes the value of the integral of the integrand... More...

virtual void	gauss_point (pod_index_type xindex, coord_type xresult[], size_type xresult_ub)
	The local coordinates of the gauss point with index xindex. More...

DIFFERENTIABLE_SECTION_EVALUATOR FACET
virtual void	dxi_local (size_type xlocal_coord_index, const dof_type xsource_dofs[], size_type xsource_dofs_ub, dof_type xresult_dofs[], size_type xresult_dofs_ub) const
	First partial derivative of this with respect to local coordinate xlocal_coord_index. More...

virtual void	jacobian (const dof_type xcoord_dofs[], size_type xcoord_dofs_ub, size_type xdf, const dof_type xlocal_coords[], size_type xlocal_coords_ub)
	Computes the the jacobian matrix at local coordinates xlocal_coords with coordinate dofs xcoord_dofs. More...

DOMAIN FACET
virtual int	db () const
	The base dimension; the dimension of the local coordinates (independent variable). More...

virtual void	local_coordinates (pod_index_type xindex, coord_type xresult[], size_type xresult_ub) const
	The local coordinates of the dof with local index xindex. More...

virtual bool	in_standard_domain (const dof_type xlocal_coords[], size_type xlocal_coords_ub) const
	Return true if the specified local coordinates are in the "standard" domain; otherwise return false. More...

EVALUATION FACET
virtual void	coord_at_value (const dof_type xdofs[], size_type xdofs_ub, const dof_type xglobal_coords[], size_type xglobal_coord_ub, dof_type xlocal_coords[], size_type xlocal_coords_ub) const
	The local coordinates of a point at which the field has the value xvalue. More...

ANY FACET
virtual quadratic_1d *	clone () const
	Virtual constructor, creates a new instance of the same type as this. More...

virtual quadratic_1d &	operator= (const section_evaluator &xother)
	Assignment operator. More...

quadratic_1d &	operator= (const quadratic_1d &xother)
	Assignment operator. More...

virtual bool	invariant () const
	Class invariant. More...

virtual bool	is_ancestor_of (const any *xother) const
	Conformance test; true if other conforms to this. More...

Additional Inherited Members
Public Types inherited from fiber_bundle::section_evaluator
typedef chart_point_coord_type	coord_type
	The type of local coordinate; the scalar type for the local coordinate vector space. More...

typedef vd_value_type	value_type
	The type of component in the value; the scalar type in the range vector space. More...

typedef sec_vd_dof_type	dof_type
	The type of degree of freedom. More...

Public Member Functions inherited from fiber_bundle::linear_fcn_space
	linear_fcn_space (const linear_fcn_space &xother)
	Copy constructor. More...

virtual	~linear_fcn_space ()
	Destructor. More...

const value_type *	basis_values () const
	The result of the preceding call to basis_at_coord. More...

const value_type *	basis_deriv_values () const
	The result of the preceding call to basis_derivs_at_coord. More...

virtual size_type	dof_ct () const
	The number of dofs required for each component of the dependent variable. More...

virtual void	value_at_coord (const dof_type xdofs[], size_type xdofs_ub, const dof_type xlocal_coords[], size_type xlocal_coords_ub, dof_type xresult[], size_type xresult_ub) const
	Computes the value of the function at local coordinates xlocal_coords using the degrees of freedom xdofs. More...

linear_fcn_space &	operator= (const linear_fcn_space &xother)
	Assignment operator. More...

Public Member Functions inherited from fiber_bundle::integrable_section_evaluator
	integrable_section_evaluator (const integrable_section_evaluator &xother)
	Copy constructor. More...

virtual	~integrable_section_evaluator ()
	Destructor. More...

virtual integrable_section_evaluator &	operator= (const integrable_section_evaluator &xother)
	Assignment operator. More...

Public Member Functions inherited from fiber_bundle::differentiable_section_evaluator
	differentiable_section_evaluator (const differentiable_section_evaluator &xother)
	Copy constructor. More...

virtual	~differentiable_section_evaluator ()
	Destructor. More...

void	dxi_local (size_type xlocal_coord_index, const block< dof_type > &xsource_dofs, block< dof_type > &xresult_dofs) const
	First partial derivative of this with respect to local coordinate xlocal_coord_index. More...

const value_type *	jacobian_values () const
	The result of the preceding call to jacobian. More...

virtual differentiable_section_evaluator &	operator= (const differentiable_section_evaluator &xother)
	Assignment operator. More...

Public Member Functions inherited from fiber_bundle::section_evaluator
	section_evaluator (const section_evaluator &xother)
	Copy constructor. More...

virtual	~section_evaluator ()
	Destructor. More...

void	local_coordinates (pod_index_type xindex, block< coord_type > &xresult) const
	The local coordinates of the dof with local index xindex. More...

virtual void	center (coord_type xresult[], size_type xresult_ub) const
	The local coordinates at the center of the evaluator. More...

void	center (block< coord_type > &xresult) const
	The local coordinates at the center of the evaluator. More...

bool	in_standard_domain (const block< dof_type > &xlocal_coords) const
	True if the specified local coordinates are in the "standard" domain. More...

bool	df_is_variable () const
	True if and only if the fiber dimension (number of components of the dependent variable) is variable. More...

virtual int	df () const
	The fiber dimension (number of components of dependent variable). More...

virtual dof_type	value_at_coord (const dof_type xdofs[], size_type xdofs_ub, const dof_type xlocal_coords[], size_type xlocal_coords_ub) const
	Value at a specified local_coordinate. Single component version. More...

dof_type	value_at_coord (const block< dof_type > &xdofs, const block< dof_type > &xlocal_coords) const
	Value at a specified local_coordinate. Single component version. More...

void	value_at_coord (const block< dof_type > &xdofs, const block< dof_type > &xlocal_coords, block< dof_type > &xresult) const
	Values at a specified local_coordinate. Multiple component version. More...

void	coord_at_value (const block< dof_type > &xdofs, const block< dof_type > &xvalue, block< dof_type > &xlocal_coords) const
	The local coordinates of a point at which the field has the value xvalue. The dofs are assumed to be interleaved. More...

virtual dof_type	min (const dof_type xdofs[], size_type xdofs_ub) const
	The minimum value of the scalar or component section defined by xdofs. More...

dof_type	min (const block< dof_type > &xdofs) const
	The minimum value of the scalar or component section defined by xdofs. More...

virtual void	min (const dof_type xdofs[], size_type xdofs_ub, dof_type xresult[], size_type xresult_ub) const
	The minimum value of the vector section defined by xdofs. More...

void	min (const block< dof_type > &xdofs, block< dof_type > &xresult) const
	The minimum value of the vector section defined by xdofs. More...

virtual dof_type	max (const dof_type xdofs[], size_type xdofs_ub) const
	The maximum value of the scalar or component section defined by xdofs. More...

dof_type	max (const block< dof_type > &xdofs) const
	The maximum value of the scalar or component section defined by xdofs. More...

virtual void	max (const dof_type xdofs[], size_type xdofs_ub, dof_type xresult[], size_type xresult_ub) const
	The maximum value of the vector section defined by xdofs. More...

void	max (const block< dof_type > &xdofs, block< dof_type > &xresult) const
	The maximum value of the vector section defined by xdofs. More...

Public Member Functions inherited from sheaf::any
bool	is_same_type (const any *other) const
	True if other is the same type as this. More...

virtual	~any ()
	Destructor. More...

bool	invariant_check () const
	True if invariant checking is enabled. More...

void	enable_invariant_check () const
	Enable invariant checking. More...

void	disable_invariant_check () const
	Disable invariant check. Intended for preventing recursive calls to invariant and for suppressing invariant checking during multi-phase initialization. More...

int	disable_invariance_check_request_depth () const
	Number of times disable_invariant_check has been called without matching call to enable_invariant_check. More...

Protected Member Functions inherited from fiber_bundle::linear_fcn_space
	linear_fcn_space ()
	Default constructor. More...

Protected Member Functions inherited from fiber_bundle::integrable_section_evaluator
	integrable_section_evaluator ()
	Default constructor. More...

Protected Member Functions inherited from fiber_bundle::differentiable_section_evaluator
	differentiable_section_evaluator ()
	Default constructor. More...

Protected Member Functions inherited from fiber_bundle::section_evaluator
	section_evaluator ()
	Default constructor. More...

Protected Member Functions inherited from sheaf::any
	any ()
	default constructor More...

Protected Attributes inherited from fiber_bundle::linear_fcn_space
value_type *	_basis_values
	The result of the preceding call to basis_at_coord. More...

value_type *	_basis_deriv_values
	The result of the preceding call to basis_derivs_at_coord. More...

Protected Attributes inherited from fiber_bundle::differentiable_section_evaluator
value_type *	_jacobian_values
	The result of the preceding call to "jacobian". More...

Detailed Description

A section evaluator using quadratic interpolation over a 1D domain.

Definition at line 38 of file quadratic_1d.h.

Constructor & Destructor Documentation

◆ quadratic_1d() [1/2]

fiber_bundle::quadratic_1d::quadratic_1d ( )

Default constructor.

Postcondition

invariant()

Definition at line 35 of file quadratic_1d.cc.

◆ quadratic_1d() [2/2]

fiber_bundle::quadratic_1d::quadratic_1d ( const quadratic_1d & xother )

Copy constructor.

Postcondition

invariant()

Definition at line 55 of file quadratic_1d.cc.

References ~quadratic_1d().

◆ ~quadratic_1d()

fiber_bundle::quadratic_1d::~quadratic_1d ( )

virtual

Destructor.

Postcondition

invariant()

Definition at line 77 of file quadratic_1d.cc.

References dl().

Referenced by quadratic_1d().

Member Function Documentation

◆ basis_at_coord()

void fiber_bundle::quadratic_1d::basis_at_coord	(	const dof_type	xlocal_coord[],
		size_type	xlocal_coord_ub
	)

virtual

Computes the value of each basis function at local coordinates xlocal_coord.

Precondition

xlocal_coord != 0

xlocal_coord_ub >= db()

Postcondition

invariant()

Implements fiber_bundle::linear_fcn_space.

Definition at line 120 of file quadratic_1d.cc.

References basis_derivs_at_coord().

Referenced by dl().

◆ basis_derivs_at_coord()

void fiber_bundle::quadratic_1d::basis_derivs_at_coord	(	const dof_type	xlocal_coords[],
		size_type	xlocal_coords_ub
	)

virtual

Computes the value of the derivatives of each basis function at local coordinates xlocal_coords.

Precondition

xlocal_coord != 0

xlocal_coord_ub >= db()

basis_deriv_values() != 0

Postcondition

invariant()

Implements fiber_bundle::linear_fcn_space.

Definition at line 146 of file quadratic_1d.cc.

References jacobian_determinant().

Referenced by basis_at_coord().

◆ clone()

fiber_bundle::quadratic_1d * fiber_bundle::quadratic_1d::clone ( ) const

virtual

Virtual constructor, creates a new instance of the same type as this.

Postcondition

result != 0

is_same_type(result)

result->invariant()

Implements fiber_bundle::linear_fcn_space.

Definition at line 602 of file quadratic_1d.cc.

References invariant(), and operator=().

Referenced by coord_at_value().

◆ coord_at_value()

void fiber_bundle::quadratic_1d::coord_at_value	(	const dof_type	xdofs[],
		size_type	xdofs_ub,
		const dof_type	xglobal_coords[],
		size_type	xglobal_coord_ub,
		dof_type	xlocal_coords[],
		size_type	xlocal_coords_ub
	)		const

virtual

The local coordinates of a point at which the field has the value xvalue.

Precondition

xdofs != 0

xdofs_ub >= 3

xglobal_coords != 0

xglobal_coord_ub >= 1

xlocal_coords != 0

xlocal_coords_ub >= 1

Postcondition

invariant()

Issue:: Made unexecutable because it is true only in exact arithmetic.

Implements fiber_bundle::section_evaluator.

Definition at line 526 of file quadratic_1d.cc.

References clone().

Referenced by in_standard_domain().

◆ db()

int fiber_bundle::quadratic_1d::db ( ) const

virtual

The base dimension; the dimension of the local coordinates (independent variable).

Postcondition

result == 1

Implements fiber_bundle::section_evaluator.

Definition at line 440 of file quadratic_1d.cc.

References local_coordinates().

Referenced by jacobian().

◆ dl()

int fiber_bundle::quadratic_1d::dl ( ) const

virtual

The dimension of this function space.

Postcondition

result == 3

Implements fiber_bundle::linear_fcn_space.

Definition at line 97 of file quadratic_1d.cc.

References basis_at_coord().

Referenced by ~quadratic_1d().

◆ dxi_local()

void fiber_bundle::quadratic_1d::dxi_local	(	size_type	xlocal_coord_index,
		const dof_type	xsource_dofs[],
		size_type	xsource_dofs_ub,
		dof_type	xresult_dofs[],
		size_type	xresult_dofs_ub
	)		const

virtual

First partial derivative of this with respect to local coordinate xlocal_coord_index.

Precondition

xlocal_coord_index < db()

xsource_dofs != 0

xsource_dofs_ub >= 3

xresult_dofs != 0

xresult_dofs_ub >= 3

Postcondition

invariant()

Implements fiber_bundle::differentiable_section_evaluator.

Definition at line 368 of file quadratic_1d.cc.

References jacobian().

Referenced by gauss_point().

◆ gauss_point()

void fiber_bundle::quadratic_1d::gauss_point	(	pod_index_type	xindex,
		coord_type	xresult[],
		size_type	xresult_ub
	)

virtual

The local coordinates of the gauss point with index xindex.

Precondition

(0 <= xindex) && (xindex < dof_ct())

xresult_ub >= db()

Not Implemented.

Todo:: Implement.

Postcondition

in_standard_domain(xresult, xresult_ub)

invariant()

Implements fiber_bundle::integrable_section_evaluator.

Definition at line 336 of file quadratic_1d.cc.

References dxi_local().

Referenced by integrate().

◆ in_standard_domain()

bool fiber_bundle::quadratic_1d::in_standard_domain	(	const dof_type	xlocal_coords[],
		size_type	xlocal_coords_ub
	)		const

virtual

Return true if the specified local coordinates are in the "standard" domain; otherwise return false.

Precondition

xlocal_coords != 0

xlocal_coords_ub >= 1

Implements fiber_bundle::section_evaluator.

Definition at line 492 of file quadratic_1d.cc.

References coord_at_value().

Referenced by local_coordinates().

◆ integrate() [1/2]

void fiber_bundle::quadratic_1d::integrate	(	const dof_type	xcoord_dofs[],
		size_type	xcoord_dofs_ub,
		size_type	xdf,
		const dof_type	xintegrands[],
		size_type	xintegrands_ub,
		value_type	xresult_integrals[],
		size_type	xresult_integrals_ub
	)

virtual

Computes the value of the integral of the integrand array...

Todo:: Get these correct.

Precondition

xcoord_dofs != 0

xcoord_dofs_ub >= dl()*db()

xintegrands != 0

xintegrands_ub >= dl()

xresult_integrals != 0

xresult_integrals_ub > 0

Not Implemented.

Todo:: Implement.

Postcondition

invariant()

Implements fiber_bundle::integrable_section_evaluator.

Definition at line 265 of file quadratic_1d.cc.

Referenced by volume().

◆ integrate() [2/2]

void fiber_bundle::quadratic_1d::integrate	(	const dof_type	xcoord_dofs[],
		size_type	xcoord_dofs_ub,
		size_type	xdf,
		const dof_type	xintegrand,
		value_type	xresult_integrals[],
		size_type	xresult_integrals_ub
	)

virtual

Computes the value of the integral of the integrand...

Todo:: Get these correct.

Precondition

xcoord_dofs != 0

xcoord_dofs_ub >= dl()*db()

xresult_integrals != 0

xresult_integrals_ub >= dl()

Not Implemented.

Todo:: Implement.

Postcondition

invariant()

Reimplemented from fiber_bundle::integrable_section_evaluator.

Definition at line 302 of file quadratic_1d.cc.

References gauss_point().

◆ invariant()

bool fiber_bundle::quadratic_1d::invariant ( ) const

virtual

Class invariant.

Invariant

basis_values() != 0

Reimplemented from fiber_bundle::linear_fcn_space.

Definition at line 669 of file quadratic_1d.cc.

References is_ancestor_of().

Referenced by clone(), and operator=().

◆ is_ancestor_of()

bool fiber_bundle::quadratic_1d::is_ancestor_of ( const any * xother ) const

virtual

Conformance test; true if other conforms to this.

Precondition

xother != 0

Reimplemented from fiber_bundle::linear_fcn_space.

Definition at line 702 of file quadratic_1d.cc.

Referenced by invariant().

◆ jacobian()

void fiber_bundle::quadratic_1d::jacobian	(	const dof_type	xcoord_dofs[],
		size_type	xcoord_dofs_ub,
		size_type	xdf,
		const dof_type	xlocal_coords[],
		size_type	xlocal_coords_ub
	)

virtual

Computes the the jacobian matrix at local coordinates xlocal_coords with coordinate dofs xcoord_dofs.

Precondition

xcoord_dofs != 0

xcoord_dofs_ub >= db()*dl()

xlocal_coords != 0

xlocal_coords_ub >= db()

jacobian_values() != 0

Not Implemented.

Todo:: Implement.

Postcondition

invariant()

Reimplemented from fiber_bundle::differentiable_section_evaluator.

Definition at line 404 of file quadratic_1d.cc.

References db().

Referenced by dxi_local().

◆ jacobian_determinant()

fiber_bundle::section_evaluator::value_type fiber_bundle::quadratic_1d::jacobian_determinant	(	const dof_type	xcoord_dofs[],
		size_type	xcoord_dofs_ub,
		size_type	xdf,
		const coord_type	xlocal_coords[],
		size_type	xlocal_coords_ub
	)

virtual

Computes the the determinant of the jacobian matrix at local coordinates xlocal_coords with coordinate dofs xcoord_dofs.

Precondition

xcoord_dofs != 0

xcoord_dofs_ub >= db()*dl()

xlocal_coords != 0

xlocal_coords_ub >= db()

Not Implemented.

Todo:: Implement.

Postcondition

invariant()

Reimplemented from fiber_bundle::differentiable_section_evaluator.

Definition at line 178 of file quadratic_1d.cc.

References volume().

Referenced by basis_derivs_at_coord().

◆ local_coordinates()

void fiber_bundle::quadratic_1d::local_coordinates	(	pod_index_type	xindex,
		coord_type	xresult[],
		size_type	xresult_ub
	)		const

virtual

The local coordinates of the dof with local index xindex.

Precondition

(0 <= xindex) && (xindex < dof_ct())

xresult_ub >= db()

Postcondition

in_standard_domain(xresult, xresult_ub)

Implements fiber_bundle::section_evaluator.

Definition at line 464 of file quadratic_1d.cc.

References in_standard_domain().

Referenced by db().

◆ operator=() [1/2]

fiber_bundle::quadratic_1d & fiber_bundle::quadratic_1d::operator= ( const section_evaluator & xother )

virtual

Assignment operator.

Precondition

is_ancestor_of(&xother)

Not Implemented.

Postcondition

invariant()

Reimplemented from fiber_bundle::linear_fcn_space.

Definition at line 626 of file quadratic_1d.cc.

Referenced by clone().

◆ operator=() [2/2]

fiber_bundle::quadratic_1d & fiber_bundle::quadratic_1d::operator= ( const quadratic_1d & xother )

Assignment operator.

Precondition

is_ancestor_of(&xother)

Not Implemented.

Postcondition

invariant()

Definition at line 646 of file quadratic_1d.cc.

References invariant().

◆ volume()

fiber_bundle::section_evaluator::value_type fiber_bundle::quadratic_1d::volume	(	const dof_type	xcoord_dofs[],
		size_type	xcoord_dofs_ub,
		size_type	xdf
	)

virtual

Volume for specified coordinate dofs xcoord_dofs and fiber space dimension xdf.

Todo:: Fix.

Precondition

xcoord_dofs != 0

xcoord_dofs_ub >= dl()*db()

xdf == 1 || xdf == 2 || xdf == 3

Issue:: Note that this is "isoparametric" which means that for xdf == 2 or xdf == 3 the 1d segment may not be straight. For now we assume that it is.

Issue:: For xdf == 1 we can just take the absolute value of x1-x0 (no need to take the sqrt of the product).

Postcondition

invariant()

result >= 0.0

Implements fiber_bundle::integrable_section_evaluator.

Definition at line 212 of file quadratic_1d.cc.

References integrate(), and fiber_bundle::sec_at0_algebra::sqrt().

Referenced by jacobian_determinant().

The documentation for this class was generated from the following files:

QUADRATIC_1D FACET

LINEAR_FCN_SPACE FACET

INTEGRABLE_SECTION_EVALUATOR FACET

DIFFERENTIABLE_SECTION_EVALUATOR FACET

DOMAIN FACET

EVALUATION FACET

ANY FACET

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ quadratic_1d() [1/2]

◆ quadratic_1d() [2/2]

◆ ~quadratic_1d()

Member Function Documentation

◆ basis_at_coord()

◆ basis_derivs_at_coord()

◆ clone()

◆ coord_at_value()

◆ db()

◆ dl()

◆ dxi_local()

◆ gauss_point()

◆ in_standard_domain()

◆ integrate() [1/2]

◆ integrate() [2/2]

◆ invariant()

◆ is_ancestor_of()

◆ jacobian()

◆ jacobian_determinant()

◆ local_coordinates()

◆ operator=() [1/2]

◆ operator=() [2/2]

◆ volume()