SAMRAI: SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE

The patch hierarchy and set of levels within that hierarcy over which the operations will be performed are set in the constructor. However, note that the constructor accepts default arguments for the coarsest and finest level numbers. If the level numbers are not specified when calling the constructor the levels which exist in the hierarchy will be assumed in all operations. The hierarchy and levels may be changed at any time using the proper member functions.

Note that, when it makes sense, an operation accept a boolean argument which indicates whether the operation should be performed on all of the data or just those data elements corresponding to the patch interiors. If no boolean argument is provided, the default behavior is to treat only the patch interiors. Also, a similar set of operations for complex and integer node-centered data is provided in the classes HierarchyNodeDataOpsComplex<DIM> and HierarchyNodeDataOpsInteger<DIM>, respectively.

Constructor & Destructor Documentation

template<int DIM, class TYPE>

SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::HierarchyNodeDataOpsReal	(	tbox::Pointer< hier::PatchHierarchy< DIM > >	hierarchy,
		const int	coarsest_level = `-1`,
		const int	finest_level = `-1`
	)

The constructor for the HierarchyNodeDataOpsReal<DIM> class sets the default patch hierarchy and coarsest and finest patch levels in that hierarchy over which operations will be performed. The hierarchy and operations may be reset using the member fuctions setPatchHierarchy() and resetLevels() below. If no level number arguments are given here, the levels over which the operations will be performed are those already existing in the hierarchy. If the hierarchy level configuration changes, the operations must be explicitly reset by calling the resetLevels() function.

template<int DIM, class TYPE>

SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::~HierarchyNodeDataOpsReal ( ) [virtual]

Virtual destructor for the HierarchyNodeDataOpsReal<DIM> class.

Member Function Documentation

template<int DIM, class TYPE>

void SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::setPatchHierarchy ( tbox::Pointer< hier::PatchHierarchy< DIM > > hierarchy ) [virtual]

Reset patch hierarchy over which operations occur.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

void SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::resetLevels	(	const int	coarsest_level,
		const int	finest_level
	)			`[virtual]`

Reset range of patch levels over which operations occur. Levels must exist in hierarchy or an assertion will result.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

const tbox::Pointer< hier::PatchHierarchy< DIM > > SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::getPatchHierarchy ( ) const [virtual]

Return const pointer to patch hierarchy associated with operations.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

void SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::copyData	(	const int	dst_id,
		const int	src_id,
		const bool	interior_only = `true`
	)			const `[virtual]`

Copy source data to destination data.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

void SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::swapData	(	const int	data1_id,
		const int	data2_id
	)			const `[virtual]`

Swap data pointers (i.e., storage) between two data components.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

void SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::printData	(	const int	data_id,
		std::ostream &	s,
		const bool	interior_only = `true`
	)			const `[virtual]`

Print data over multiple levels to specified output stream.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

void SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::setToScalar	(	const int	data_id,
		const TYPE &	alpha,
		const bool	interior_only = `true`
	)			const `[virtual]`

Set data component to given scalar.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

void SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::scale	(	const int	dst_id,
		const TYPE &	alpha,
		const int	src_id,
		const bool	interior_only = `true`
	)			const `[virtual]`

Set destination to source multiplied by given scalar, pointwise.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

void SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::addScalar	(	const int	dst_id,
		const int	src_id,
		const TYPE &	alpha,
		const bool	interior_only = `true`
	)			const `[virtual]`

Add scalar to each entry in source data and set destination to result.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

void SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::add	(	const int	dst_id,
		const int	src1_id,
		const int	src2_id,
		const bool	interior_only = `true`
	)			const `[virtual]`

Set destination to sum of two source components, pointwise.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

void SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::subtract	(	const int	dst_id,
		const int	src1_id,
		const int	src2_id,
		const bool	interior_only = `true`
	)			const `[virtual]`

Subtract second source component from first source component pointwise and set destination data component to result.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

void SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::multiply	(	const int	dst_id,
		const int	src1_id,
		const int	src2_id,
		const bool	interior_only = `true`
	)			const `[virtual]`

Set destination component to product of two source components, pointwise.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

void SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::divide	(	const int	dst_id,
		const int	src1_id,
		const int	src2_id,
		const bool	interior_only = `true`
	)			const `[virtual]`

Divide first data component by second source component pointwise and set destination data component to result.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

void SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::reciprocal	(	const int	dst_id,
		const int	src_id,
		const bool	interior_only = `true`
	)			const `[virtual]`

Set each entry of destination component to reciprocal of corresponding source data component entry.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

void SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::linearSum	(	const int	dst_id,
		const TYPE &	alpha,
		const int	src1_id,
		const TYPE &	beta,
		const int	src2_id,
		const bool	interior_only = `true`
	)			const `[virtual]`

Set $d = \alpha s_1 + \beta s_2$ , where $d$ is the destination patch data component and $s_1, s_2$ are the first and second source components, respectively. Here $\alpha, \beta$ are scalar values.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

void SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::axpy	(	const int	dst_id,
		const TYPE &	alpha,
		const int	src1_id,
		const int	src2_id,
		const bool	interior_only = `true`
	)			const `[virtual]`

Set $d = \alpha s_1 + s_2$ , where $d$ is the destination patch data component and $s_1, s_2$ are the first and second source components, respectively. Here $\alpha$ is a scalar.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

void SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::axmy	(	const int	dst_id,
		const TYPE &	alpha,
		const int	src1_id,
		const int	src2_id,
		const bool	interior_only = `true`
	)			const `[virtual]`

Set $d = \alpha s_1 - s_2$ , where $d$ is the destination patch data component and $s_1, s_2$ are the first and second source components, respectively. Here $\alpha$ is a scalar.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

void SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::abs	(	const int	dst_id,
		const int	src_id,
		const bool	interior_only = `true`
	)			const `[virtual]`

Set destination data to absolute value of source data, pointwise.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

TYPE SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::min	(	const int	data_id,
		const bool	interior_only = `true`
	)			const `[virtual]`

Return minimum data value over all patches in the collection of levels.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

TYPE SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::max	(	const int	data_id,
		const bool	interior_only = `true`
	)			const `[virtual]`

Return maximum data value over all patches in the collection of levels.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

void SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::setRandomValues	(	const int	data_id,
		const TYPE &	width,
		const TYPE &	low,
		const bool	interior_only = `true`
	)			const `[virtual]`

Set data entries to random values. See the operations in the array data operation classes for each data type for details on the generation of the random values.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

int SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::numberOfEntries	(	const int	data_id,
		const bool	interior_only = `true`
	)			const `[virtual]`

Return the total number of data values for the component on the set of hierarchy levels. If the boolean argument is true, the number of elements will be summed over patch interiors in a unique way which avoids multiple counting of redundant values (recall the definition of node points on a patch interior). If the boolean argument is false, all elements will be counted (including ghost values) over all patches.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

double SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::sumControlVolumes	(	const int	data_id,
		const int	vol_id
	)			const `[virtual]`

Return sum of the control volumes associated with the data component. Note that if the ontrol volumes are set propery, this is equivalent to integrating a data component containing all ones over the collection of hierarchy levels.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

double SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::L1Norm	(	const int	data_id,
		const int	vol_id = `-1`,
		bool	local_only = `false`
	)			const `[virtual]`

Return discrete $L_1$ -norm of the data using the control volume to weight the contribution of each data entry to the sum. That is, the return value is the sum $\sum_i ( \| data_i \| cvol_i )$ . If the control volume is not defined (vol_id < 0), the return value is $\sum_i ( \| data_i \| )$ . If local_only is true, the global reduction is not performed (thus each process will get only local results).

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

double SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::L2Norm	(	const int	data_id,
		const int	vol_id = `-1`,
		bool	local_only = `false`
	)			const `[virtual]`

Return discrete $L_2$ -norm of the data using the control volume to weight the contribution of each data entry to the sum. That is, the return value is the sum $\sqrt{ \sum_i ( (data_i)^2 cvol_i ) }$ . If the control volume is not defined (vol_id < 0), the return value is $\sqrt{ \sum_i ( (data_i)^2 cvol_i ) }$ . If local_only is true, the global reduction is not performed (thus each process will get only local results).

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

double SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::weightedL2Norm	(	const int	data_id,
		const int	weight_id,
		const int	vol_id = `-1`
	)			const `[virtual]`

Return discrete weighted $L_2$ -norm of the data using the control volume to weight the contribution of the data and weight entries to the sum. That is, the return value is the sum $\sqrt{ \sum_i ( (data_i * weight_i)^2 cvol_i ) }$ . If the control volume is not defined (vol_id < 0), the return value is $\sqrt{ \sum_i ( (data_i * weight_i)^2 ) }$ .

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

double SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::RMSNorm	(	const int	data_id,
		const int	vol_id = `-1`
	)			const `[virtual]`

Return discrete root mean squared norm of the data. If the control volume is defined, the return value is the $L_2$ -norm divided by the square root of the sum of the control volumes. Otherwise, the return value is the $L_2$ -norm divided by the square root of the number of data entries.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

double SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::weightedRMSNorm	(	const int	data_id,
		const int	weight_id,
		const int	vol_id = `-1`
	)			const `[virtual]`

Return discrete weighted root mean squared norm of the data. If the control volume is define, the return value is the weighted $L_2$ -norm divided by the square root of the sum of the control volumes. Otherwise, the return value is the weighted $L_2$ -norm divided by the square root of the number of data entries.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

double SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::maxNorm	(	const int	data_id,
		const int	vol_id = `-1`,
		bool	local_only = `false`
	)			const `[virtual]`

Return the $\max$ -norm of the data using the control volume to weight the contribution of each data entry to the maximum. That is, the return value is $\max_i ( \| data_i \| )$ , where the max is over the data elements where $cvol_i > 0$ . If the control volume is undefined (vol_id < 0), it is ignored during the computation of the maximum. If local_only is true, the global reduction is not performed (thus each process will get only local results).

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

TYPE SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::dot	(	const int	data1_id,
		const int	data2_id,
		const int	vol_id = `-1`,
		bool	local_only = `false`
	)			const `[virtual]`

Return the dot product of the two data arrays using the control volume to weight the contribution of each product to the sum. That is, the return value is the sum $\sum_i ( data1_i * data2_i * cvol_i )$ . If the control volume is undefined (vol_id < 0), it is ignored during the summation. If local_only is true, the global reduction is not performed (thus each process will get only local results).

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

TYPE SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::integral	(	const int	data_id,
		const int	vol_id
	)			const

Return the integral of the function represented by the data array. The return value is the sum $\sum_i ( data_i * vol_i )$ .

template<int DIM, class TYPE>

int SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::computeConstrProdPos	(	const int	data1_id,
		const int	data2_id,
		const int	vol_id = `-1`
	)			const `[virtual]`

Return 1 if $\|data2_i\| > 0$ and $data1_i * data2_i \leq 0$ , for any $i$ in the set of patch data indices, where $cvol_i > 0$ . Otherwise, return 0. If the control volume is undefined (vol_id < 0), all values on the patch interiors are considered.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

void SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::compareToScalar	(	const int	dst_id,
		const int	src_id,
		const TYPE &	alpha,
		const int	vol_id = `-1`
	)			const `[virtual]`

Wherever $cvol_i > 0$ in the set of patch data indices, set $dst_i = 1$ if $\|src_i\| > \alpha$ , and $dst_i = 0$ otherwise. If the control volume is undefined (vol_id < 0), all values on the patch interiors are considered.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

int SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::testReciprocal	(	const int	dst_id,
		const int	src_id,
		const int	vol_id = `-1`
	)			const `[virtual]`

Wherever $cvol_i > 0$ in the set of patch data indices, set $dst_i = 1/src_i$ if $src_i \neq 0$ , and $dst_i = 0$ otherwise. If anywhere, 0 is the return value. Otherwise 1 is returned. If the control volume is undefined (vol_id < 0), all values on the patch interiors are considered.

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

TYPE SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::maxPointwiseDivide	(	const int	numer_id,
		const int	denom_id,
		bool	local_only = `false`
	)			const `[virtual]`

Compute max of "conditional" quotients of two arrays.

Return the maximum of pointwise "conditional" quotients of the numerator and denominator.

The "conditional" quotient is defined as |numerator/denominator| where the denominator is nonzero. Otherwise, it is defined as |numerator|.

Note: This method is currently intended to support the PETSc-2.1.6 vector wrapper only. Please do not use it!

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

template<int DIM, class TYPE>

TYPE SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >::minPointwiseDivide	(	const int	numer_id,
		const int	denom_id,
		bool	local_only = `false`
	)			const `[virtual]`

Compute min of quotients of two arrays.

Return the minimum of pointwise quotients of the numerator and denominator.

The quotient is defined as (numerator/denominator) where the denominator is nonzero. When the denominator is zero, the entry is skipped. If the denominator is always zero, the value of tbox::IEEE::getFLT_MAX() is returned (see SAMRAI::tbox::IEEE).

Note: This method is currently intended to support the SUNDIALS vector wrapper only. Please do not use it!

Implements SAMRAI::math::HierarchyDataOpsReal< DIM, TYPE >.

SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE > Class Template Reference

Public Member Functions

Detailed Description

template<int DIM, class TYPE>
class SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >

Constructor & Destructor Documentation

Member Function Documentation


Public Member Functions
	HierarchyNodeDataOpsReal (tbox::Pointer< hier::PatchHierarchy< DIM > > hierarchy, const int coarsest_level=-1, const int finest_level=-1)
virtual	~HierarchyNodeDataOpsReal ()
void	setPatchHierarchy (tbox::Pointer< hier::PatchHierarchy< DIM > > hierarchy)
void	resetLevels (const int coarsest_level, const int finest_level)
const tbox::Pointer< hier::PatchHierarchy< DIM > >	getPatchHierarchy () const
void	copyData (const int dst_id, const int src_id, const bool interior_only=true) const
void	swapData (const int data1_id, const int data2_id) const
void	printData (const int data_id, std::ostream &s, const bool interior_only=true) const
void	setToScalar (const int data_id, const TYPE &alpha, const bool interior_only=true) const
void	scale (const int dst_id, const TYPE &alpha, const int src_id, const bool interior_only=true) const
void	addScalar (const int dst_id, const int src_id, const TYPE &alpha, const bool interior_only=true) const
void	add (const int dst_id, const int src1_id, const int src2_id, const bool interior_only=true) const
void	subtract (const int dst_id, const int src1_id, const int src2_id, const bool interior_only=true) const
void	multiply (const int dst_id, const int src1_id, const int src2_id, const bool interior_only=true) const
void	divide (const int dst_id, const int src1_id, const int src2_id, const bool interior_only=true) const
void	reciprocal (const int dst_id, const int src_id, const bool interior_only=true) const
void	linearSum (const int dst_id, const TYPE &alpha, const int src1_id, const TYPE &beta, const int src2_id, const bool interior_only=true) const
void	axpy (const int dst_id, const TYPE &alpha, const int src1_id, const int src2_id, const bool interior_only=true) const
void	axmy (const int dst_id, const TYPE &alpha, const int src1_id, const int src2_id, const bool interior_only=true) const
void	abs (const int dst_id, const int src_id, const bool interior_only=true) const
TYPE	min (const int data_id, const bool interior_only=true) const
TYPE	max (const int data_id, const bool interior_only=true) const
void	setRandomValues (const int data_id, const TYPE &width, const TYPE &low, const bool interior_only=true) const
int	numberOfEntries (const int data_id, const bool interior_only=true) const
double	sumControlVolumes (const int data_id, const int vol_id) const
double	L1Norm (const int data_id, const int vol_id=-1, bool local_only=false) const
double	L2Norm (const int data_id, const int vol_id=-1, bool local_only=false) const
double	weightedL2Norm (const int data_id, const int weight_id, const int vol_id=-1) const
double	RMSNorm (const int data_id, const int vol_id=-1) const
double	weightedRMSNorm (const int data_id, const int weight_id, const int vol_id=-1) const
double	maxNorm (const int data_id, const int vol_id=-1, bool local_only=false) const
TYPE	dot (const int data1_id, const int data2_id, const int vol_id=-1, bool local_only=false) const
TYPE	integral (const int data_id, const int vol_id) const
int	computeConstrProdPos (const int data1_id, const int data2_id, const int vol_id=-1) const
void	compareToScalar (const int dst_id, const int src_id, const TYPE &alpha, const int vol_id=-1) const
int	testReciprocal (const int dst_id, const int src_id, const int vol_id=-1) const
TYPE	maxPointwiseDivide (const int numer_id, const int denom_id, bool local_only=false) const
	Compute max of "conditional" quotients of two arrays.
TYPE	minPointwiseDivide (const int numer_id, const int denom_id, bool local_only=false) const
	Compute min of quotients of two arrays.

SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE > Class Template Reference

Public Member Functions

Detailed Description

template<int DIM, class TYPE> class SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >

Constructor & Destructor Documentation

Member Function Documentation

template<int DIM, class TYPE>
class SAMRAI::math::HierarchyNodeDataOpsReal< DIM, TYPE >