Concrete implementation of a clustered Spmd-based product vector. More...

#include <Thyra_DefaultClusteredSpmdProductVector_decl.hpp>

Inheritance diagram for Thyra::DefaultClusteredSpmdProductVector< Scalar >:

[legend]

Overridden protected members from VectorBase
void	applyOpImpl (const RTOpPack::RTOpT< Scalar > &op, const ArrayView< const Ptr< const VectorBase< Scalar > > > &vecs, const ArrayView< const Ptr< VectorBase< Scalar > > > &targ_vecs, const Ptr< RTOpPack::ReductTarget > &reduct_obj, const Ordinal global_offset) const

Constructors/initializers/accessors
	DefaultClusteredSpmdProductVector ()
	Constructs to uninitialized (see postconditions for uninitialize()).
	DefaultClusteredSpmdProductVector (const Teuchos::RCP< const DefaultClusteredSpmdProductVectorSpace< Scalar > > &productSpace, const Teuchos::RCP< VectorBase< Scalar > > vecs[])
	Constructs to initialized (calls initialize()).
void	initialize (const Teuchos::RCP< const DefaultClusteredSpmdProductVectorSpace< Scalar > > &productSpace, const Teuchos::RCP< VectorBase< Scalar > > vecs[])
	Initialize.
void	uninitialize (Teuchos::RCP< const DefaultClusteredSpmdProductVectorSpace< Scalar > > *productSpace=NULL, Teuchos::RCP< VectorBase< Scalar > > vecs[]=NULL)
	Uninitialize.

Overridden from ProductVectorBase
Teuchos::RCP< VectorBase< Scalar > >	getNonconstVectorBlock (const int k)
Teuchos::RCP< const VectorBase< Scalar > >	getVectorBlock (const int k) const

Overridden from ProductMultiVectorBase
Teuchos::RCP< const ProductVectorSpaceBase< Scalar > >	productSpace () const
bool	blockIsConst (const int k) const
Teuchos::RCP< MultiVectorBase< Scalar > >	getNonconstMultiVectorBlock (const int k)
Teuchos::RCP< const MultiVectorBase< Scalar > >	getMultiVectorBlock (const int k) const

Overridden from VectorBase
Teuchos::RCP< const VectorSpaceBase< Scalar > >	space () const

Additional Inherited Members
void	assign (const VectorBase< Scalar > &x)
	Vector assignment:
void	randomize (Scalar l, Scalar u)
	Random vector generation:
void	update (Scalar alpha, const VectorBase< Scalar > &x)
	AXPY:
void	linear_combination (const ArrayView< const Scalar > &alpha, const ArrayView< const Ptr< const VectorBase< Scalar > > > &x, const Scalar &beta)
	Linear combination:
Scalar	dot (const VectorBase< Scalar > &x) const
	Euclidean dot product: result = x^H * this.
Teuchos::ScalarTraits< Scalar >::magnitudeType	norm_1 () const
	One (1) norm: result = \|\|v\|\|1.
Teuchos::ScalarTraits< Scalar >::magnitudeType	norm_2 () const
	Euclidean (2) norm: result = \|\|v\|\|2.
Teuchos::ScalarTraits< Scalar >::magnitudeType	norm_2 (const VectorBase< Scalar > &x) const
	Weighted Euclidean (2) norm: result = \|\|v\|\|2.
Teuchos::ScalarTraits< Scalar >::magnitudeType	norm_inf () const
	Infinity norm: result = \|\|v\|\|inf.
void	applyOp (const RTOpPack::RTOpT< Scalar > &op, const ArrayView< const Ptr< const VectorBase< Scalar > > > &vecs, const ArrayView< const Ptr< VectorBase< Scalar > > > &targ_vecs, const Ptr< RTOpPack::ReductTarget > &reduct_obj, const Ordinal global_offset) const
	Calls applyOpImpl().
void	acquireDetachedView (const Range1D &rng, RTOpPack::ConstSubVectorView< Scalar > *sub_vec) const
	Calls acquireDetachedVectorViewImpl().
void	releaseDetachedView (RTOpPack::ConstSubVectorView< Scalar > *sub_vec) const
	Calls releaseDetachedVectorViewImpl().
void	acquireDetachedView (const Range1D &rng, RTOpPack::SubVectorView< Scalar > *sub_vec)
	Calls acquireNonconstDetachedVectorViewImpl().
void	commitDetachedView (RTOpPack::SubVectorView< Scalar > *sub_vec)
	Calls commitDetachedView().
void	setSubVector (const RTOpPack::SparseSubVectorT< Scalar > &sub_vec)
	Calls setSubVectorImpl().
void	assign (Scalar alpha)
	V = alpha.
void	assign (const MultiVectorBase< Scalar > &mv)
	V = mv.
void	scale (Scalar alpha)
void	update (Scalar alpha, const MultiVectorBase< Scalar > &mv)
void	linear_combination (const ArrayView< const Scalar > &alpha, const ArrayView< const Ptr< const MultiVectorBase< Scalar > > > &mv, const Scalar &beta)
	Y.col(j)(i) = betaY.col(j)(i) + sum( alpha[k]X[k].col(j)(i),
void	dots (const MultiVectorBase< Scalar > &mv, const ArrayView< Scalar > &prods) const
	Column-wise Euclidean dot product.
void	norms_1 (const ArrayView< typename ScalarTraits< Scalar >::magnitudeType > &norms) const
	Column-wise 1-norms.
void	norms_2 (const ArrayView< typename ScalarTraits< Scalar >::magnitudeType > &norms) const
	Column-wise 2-norms.
void	norms_inf (const ArrayView< typename ScalarTraits< Scalar >::magnitudeType > &norms) const
	Column-wise infinity-norms.
RCP< const VectorBase< Scalar > >	col (Ordinal j) const
	Calls colImpl().
RCP< VectorBase< Scalar > >	col (Ordinal j)
	Calls nonconstColImpl().
RCP< const MultiVectorBase< Scalar > >	subView (const Range1D &colRng) const
	Calls contigSubViewImpl().
RCP< MultiVectorBase< Scalar > >	subView (const Range1D &colRng)
	Calls nonconstContigSubViewImpl().
RCP< const MultiVectorBase< Scalar > >	subView (const ArrayView< const int > &cols) const
	nonContigSubViewImpl().
RCP< MultiVectorBase< Scalar > >	subView (const ArrayView< const int > &cols)
	nonconstNonContigSubViewImpl().
void	applyOp (const RTOpPack::RTOpT< Scalar > &primary_op, const ArrayView< const Ptr< const MultiVectorBase< Scalar > > > &multi_vecs, const ArrayView< const Ptr< MultiVectorBase< Scalar > > > &targ_multi_vecs, const ArrayView< const Ptr< RTOpPack::ReductTarget > > &reduct_objs, const Ordinal primary_global_offset) const
	Calls mvMultiReductApplyOpImpl().
void	applyOp (const RTOpPack::RTOpT< Scalar > &primary_op, const RTOpPack::RTOpT< Scalar > &secondary_op, const ArrayView< const Ptr< const MultiVectorBase< Scalar > > > &multi_vecs, const ArrayView< const Ptr< MultiVectorBase< Scalar > > > &targ_multi_vecs, const Ptr< RTOpPack::ReductTarget > &reduct_obj, const Ordinal primary_global_offset) const
	mvSingleReductApplyOpImpl().
void	acquireDetachedView (const Range1D &rowRng, const Range1D &colRng, RTOpPack::ConstSubMultiVectorView< Scalar > *sub_mv) const
	Calls acquireDetachedMultiVectorViewImpl().
void	releaseDetachedView (RTOpPack::ConstSubMultiVectorView< Scalar > *sub_mv) const
	Calls releaseDetachedMultiVectorViewImpl().
void	acquireDetachedView (const Range1D &rowRng, const Range1D &colRng, RTOpPack::SubMultiVectorView< Scalar > *sub_mv)
	Calls acquireNonconstDetachedMultiVectorViewImpl().
void	commitDetachedView (RTOpPack::SubMultiVectorView< Scalar > *sub_mv)
	Calls commitNonconstDetachedMultiVectorViewImpl().
RCP< const LinearOpBase< Scalar > >	clone () const
	This function is simply overridden to return this->clone_mv().
bool	opSupported (EOpTransp M_trans) const
	Return if the M_trans operation of apply() is supported or not.
void	apply (const EOpTransp M_trans, const MultiVectorBase< Scalar > &X, const Ptr< MultiVectorBase< Scalar > > &Y, const Scalar alpha, const Scalar beta) const
	Apply the linear operator to a multi-vector : Y = alphaop(M)X + beta*Y.
bool	rowStatIsSupported (const RowStatLinearOpBaseUtils::ERowStat rowStat) const
	Determine if a given row stat is supported.
void	getRowStat (const RowStatLinearOpBaseUtils::ERowStat rowStat, const Ptr< VectorBase< Scalar > > &rowStatVec) const
	Get some statistics about a supported row.
bool	supportsScaleLeft () const
	Determines if this objects supports left scaling.
bool	supportsScaleRight () const
	Determines if this objects supports right scaling.
void	scaleLeft (const VectorBase< Scalar > &row_scaling)
	Left scales operator with diagonal scaling operator.
void	scaleRight (const VectorBase< Scalar > &col_scaling)
	Right scales operator with diagonal scaling operator.
virtual RCP< const VectorBase< Scalar > >	colImpl (Ordinal j) const
	Return a non-changeable view of a constituent column vector.
virtual bool	rowStatIsSupportedImpl (const RowStatLinearOpBaseUtils::ERowStat rowStat) const
virtual void	getRowStatImpl (const RowStatLinearOpBaseUtils::ERowStat rowStat, const Ptr< VectorBase< Scalar > > &rowStatVec) const
virtual bool	supportsScaleLeftImpl () const
virtual bool	supportsScaleRightImpl () const
virtual void	scaleLeftImpl (const VectorBase< Scalar > &row_scaling)
virtual void	scaleRightImpl (const VectorBase< Scalar > &col_scaling)
void	absRowSum (const Teuchos::Ptr< Thyra::VectorBase< Scalar > > &output) const
void	absColSum (const Teuchos::Ptr< Thyra::VectorBase< Scalar > > &output) const
bool	opSupportedImpl (EOpTransp M_trans) const
	For complex Scalar types returns true for NOTRANS and CONJTRANS and for real types returns true for all values of M_trans.
void	applyImpl (const EOpTransp M_trans, const MultiVectorBase< Scalar > &X, const Ptr< MultiVectorBase< Scalar > > &Y, const Scalar alpha, const Scalar beta) const
	. Applies vector or its adjoint (transpose) as a linear operator.
virtual void	assignVecImpl (const VectorBase< Scalar > &x)
	Default implementation of assign(vector) using RTOps.
virtual void	randomizeImpl (Scalar l, Scalar u)
	Default implementation of randomize using RTOps.
virtual void	absImpl (const VectorBase< Scalar > &x)
	Default implementation of abs using RTOps.
virtual void	reciprocalImpl (const VectorBase< Scalar > &x)
	Default implementation of reciprocal using RTOps.
virtual void	eleWiseScaleImpl (const VectorBase< Scalar > &x)
	Default implementation of ele_wise_scale using RTOps.
virtual void	updateVecImpl (Scalar alpha, const VectorBase< Scalar > &x)
	Default implementation of update using RTOps.
virtual void	linearCombinationVecImpl (const ArrayView< const Scalar > &alpha, const ArrayView< const Ptr< const VectorBase< Scalar > > > &x, const Scalar &beta)
	Default implementation of linear_combination using RTOps.
virtual Scalar	dotImpl (const VectorBase< Scalar > &x) const
	Default implementation of dot using RTOps.
virtual Teuchos::ScalarTraits< Scalar >::magnitudeType	norm1Impl () const
	Default implementation of norm_1 using RTOps.
virtual Teuchos::ScalarTraits< Scalar >::magnitudeType	norm2Impl () const
	Default implementation of norm_2 using RTOps.
virtual Teuchos::ScalarTraits< Scalar >::magnitudeType	norm2WeightedImpl (const VectorBase< Scalar > &x) const
	Default implementation of norm_2 (weighted) using RTOps.
virtual Teuchos::ScalarTraits< Scalar >::magnitudeType	normInfImpl () const
	Default implementation of norm_inf using RTOps.
virtual RCP< VectorBase< Scalar > >	nonconstColImpl (Ordinal j)
	Returns Teuchos::rcp(this,false).
virtual RCP< const MultiVectorBase< Scalar > >	contigSubViewImpl (const Range1D &col_rng) const
	Returns Teuchos::rcp(this,false).
virtual RCP< MultiVectorBase< Scalar > >	nonconstContigSubViewImpl (const Range1D &col_rng)
	Returns Teuchos::rcp(this,false).
virtual RCP< const MultiVectorBase< Scalar > >	nonContigSubViewImpl (const ArrayView< const int > &cols) const
	Returns Teuchos::rcp(this,false).
virtual RCP< MultiVectorBase< Scalar > >	nonconstNonContigSubViewImpl (const ArrayView< const int > &cols)
	Returns Teuchos::rcp(this,false).
virtual void	acquireDetachedMultiVectorViewImpl (const Range1D &rowRng, const Range1D &colRng, RTOpPack::ConstSubMultiVectorView< Scalar > *sub_mv) const
	Implemented in terms of this->acquireDetachedView().
virtual void	releaseDetachedMultiVectorViewImpl (RTOpPack::ConstSubMultiVectorView< Scalar > *sub_mv) const
	Implemented in terms of this->releaseDetachedView().
virtual void	acquireNonconstDetachedMultiVectorViewImpl (const Range1D &rowRng, const Range1D &colRng, RTOpPack::SubMultiVectorView< Scalar > *sub_mv)
	Implemented in terms of this->acquireDetachedView().
virtual void	commitNonconstDetachedMultiVectorViewImpl (RTOpPack::SubMultiVectorView< Scalar > *sub_mv)
	Implemented in terms of this->commitDetachedView().
virtual void	acquireDetachedVectorViewImpl (const Range1D &rng, RTOpPack::ConstSubVectorView< Scalar > *sub_vec) const
virtual void	releaseDetachedVectorViewImpl (RTOpPack::ConstSubVectorView< Scalar > *sub_vec) const
virtual void	acquireNonconstDetachedVectorViewImpl (const Range1D &rng, RTOpPack::SubVectorView< Scalar > *sub_vec)
virtual void	commitNonconstDetachedVectorViewImpl (RTOpPack::SubVectorView< Scalar > *sub_vec)
virtual void	setSubVectorImpl (const RTOpPack::SparseSubVectorT< Scalar > &sub_vec)
virtual std::string	description () const
	Default description that gives the label, type, and dimenstion .
virtual void	describe (Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel) const
	Generates a default outputting for all vectors.
virtual RCP< const VectorSpaceBase< Scalar > >	range () const
	Returns this->space().
virtual RCP< const VectorSpaceBase< Scalar > >	domain () const
	Returns a DefaultSerialVectorSpace object with dimension 1.
RCP< MultiVectorBase< Scalar > >	clone_mv () const
	Returns this->clone_v().
RCP< VectorBase< Scalar > >	clone_v () const
	Simply creates a new vector and copies the contents from *this.
virtual void	assignImpl (Scalar alpha)
	Default implementation of assign(scalar) using RTOps.
virtual void	assignMultiVecImpl (const MultiVectorBase< Scalar > &mv)
	Default implementation of assign(MV) using RTOps.
virtual void	scaleImpl (Scalar alpha)
	Default implementation of scale using RTOps.
virtual void	updateImpl (Scalar alpha, const MultiVectorBase< Scalar > &mv)
	Default implementation of update using RTOps.
virtual void	linearCombinationImpl (const ArrayView< const Scalar > &alpha, const ArrayView< const Ptr< const MultiVectorBase< Scalar > > > &mv, const Scalar &beta)
	Default implementation of linear_combination using RTOps.
virtual void	dotsImpl (const MultiVectorBase< Scalar > &mv, const ArrayView< Scalar > &prods) const
	Default implementation of dots using RTOps.
virtual void	norms1Impl (const ArrayView< typename ScalarTraits< Scalar >::magnitudeType > &norms) const
	Default implementation of norms_1 using RTOps.
virtual void	norms2Impl (const ArrayView< typename ScalarTraits< Scalar >::magnitudeType > &norms) const
	Default implementation of norms_2 using RTOps.
virtual void	normsInfImpl (const ArrayView< typename ScalarTraits< Scalar >::magnitudeType > &norms) const
	Default implementation of norms_inf using RTOps.
virtual void	mvMultiReductApplyOpImpl (const RTOpPack::RTOpT< Scalar > &primary_op, const ArrayView< const Ptr< const MultiVectorBase< Scalar > > > &multi_vecs, const ArrayView< const Ptr< MultiVectorBase< Scalar > > > &targ_multi_vecs, const ArrayView< const Ptr< RTOpPack::ReductTarget > > &reduct_objs, const Ordinal primary_global_offset) const
virtual void	mvSingleReductApplyOpImpl (const RTOpPack::RTOpT< Scalar > &primary_op, const RTOpPack::RTOpT< Scalar > &secondary_op, const ArrayView< const Ptr< const MultiVectorBase< Scalar > > > &multi_vecs, const ArrayView< const Ptr< MultiVectorBase< Scalar > > > &targ_multi_vecs, const Ptr< RTOpPack::ReductTarget > &reduct_obj, const Ordinal primary_global_offset) const
Related Symbols inherited from Thyra::ProductVectorBase< Scalar >
template<class Scalar>
RCP< Thyra::ProductVectorBase< Scalar > >	nonconstProductVectorBase (const RCP< Thyra::VectorBase< Scalar > > &v)
	Dynamic cast from a VectorBase to a ProductVectorBase object and thow exception if this fails.
template<class Scalar>
RCP< const Thyra::ProductVectorBase< Scalar > >	productVectorBase (const RCP< const Thyra::VectorBase< Scalar > > &v)
	Dynamic cast from a const VectorBase to a const ProductVectorBase object and thow exception if this fails.
Related Symbols inherited from Thyra::VectorBase< Scalar >
template<class Scalar>
void	applyOp (const RTOpPack::RTOpT< Scalar > &op, const ArrayView< const Ptr< const VectorBase< Scalar > > > &vecs, const ArrayView< const Ptr< VectorBase< Scalar > > > &targ_vecs, const Ptr< RTOpPack::ReductTarget > &reduct_obj, const Ordinal global_offset=0)
	Apply a reduction/transformation operator over a set of vectors: op(op(v[0]...v[nv-1],z[0]...z[nz-1]),(reduct_obj)) -> z[0]...z[nz-1],(reduct_obj).
template<class Scalar>
Scalar	sum (const VectorBase< Scalar > &v)
	Sum of vector elements: result = sum( v(i), i = 0...v.space()->dim()-1 ).
template<class Scalar>
Scalar	scalarProd (const VectorBase< Scalar > &x, const VectorBase< Scalar > &y)
	Scalar product result = <x,y>.
template<class Scalar>
Scalar	inner (const VectorBase< Scalar > &x, const VectorBase< Scalar > &y)
	Inner/Scalar product result = <x,y>.
template<class Scalar>
Teuchos::ScalarTraits< Scalar >::magnitudeType	norm (const VectorBase< Scalar > &v)
	Natural norm: result = sqrt(<v,v>).
template<class Scalar>
Teuchos::ScalarTraits< Scalar >::magnitudeType	norm_1 (const VectorBase< Scalar > &v)
	One (1) norm: result = \|\|v\|\|1.
template<class Scalar>
Teuchos::ScalarTraits< Scalar >::magnitudeType	norm_2 (const VectorBase< Scalar > &v)
	Euclidean (2) norm: result = \|\|v\|\|2.
template<class Scalar>
Teuchos::ScalarTraits< Scalar >::magnitudeType	norm_2 (const VectorBase< Scalar > &w, const VectorBase< Scalar > &v)
	Weighted Euclidean (2) norm: result = sqrt( sum( w(i)conj(v(i))v(i)) ).
template<class Scalar>
Teuchos::ScalarTraits< Scalar >::magnitudeType	norm_inf (const VectorBase< Scalar > &v_rhs)
	Infinity norm: result = \|\|v\|\|inf.
template<class Scalar>
Scalar	dot (const VectorBase< Scalar > &x, const VectorBase< Scalar > &y)
	Dot product: result = conj(x)'*y.
template<class Scalar>
Scalar	get_ele (const VectorBase< Scalar > &v, Ordinal i)
	Get single element: result = v(i).
template<class Scalar>
void	set_ele (Ordinal i, Scalar alpha, const Ptr< VectorBase< Scalar > > &v)
	Set single element: v(i) = alpha.
template<class Scalar>
void	put_scalar (const Scalar &alpha, const Ptr< VectorBase< Scalar > > &y)
	Assign all elements to a scalar: y(i) = alpha, i = 0...y->space()->dim()-1.
template<class Scalar>
void	copy (const VectorBase< Scalar > &x, const Ptr< VectorBase< Scalar > > &y)
	Vector assignment: y(i) = x(i), i = 0...y->space()->dim()-1.
template<class Scalar>
void	add_scalar (const Scalar &alpha, const Ptr< VectorBase< Scalar > > &y)
	Add a scalar to all elements: y(i) += alpha, i = 0...y->space()->dim()-1.
template<class Scalar>
void	scale (const Scalar &alpha, const Ptr< VectorBase< Scalar > > &y)
	Scale all elements by a scalar: y(i) *= alpha, i = 0...y->space()->dim()-1.
template<class Scalar>
void	abs (const VectorBase< Scalar > &x, const Ptr< VectorBase< Scalar > > &y)
	Element-wise absolute value: y(i) = abs(x(i)), i = 0...y->space()->dim()-1.
template<class Scalar>
void	reciprocal (const VectorBase< Scalar > &x, const Ptr< VectorBase< Scalar > > &y)
	Element-wise reciprocal: y(i) = 1/x(i), i = 0...y->space()->dim()-1.
template<class Scalar>
void	ele_wise_prod (const Scalar &alpha, const VectorBase< Scalar > &x, const VectorBase< Scalar > &v, const Ptr< VectorBase< Scalar > > &y)
	Element-wise product update: y(i) += alpha * x(i) * v(i), i = 0...y->space()->dim()-1.
template<class Scalar>
void	pair_wise_max (const Scalar &alpha, const VectorBase< Scalar > &x, const VectorBase< Scalar > &v, const Ptr< VectorBase< Scalar > > &y)
	Element-wise maximum: y(i) = alpha * max(x(i), v(i)), i = 0...y->space()->dim()-1.
template<class Scalar>
void	ele_wise_conj_prod (const Scalar &alpha, const VectorBase< Scalar > &x, const VectorBase< Scalar > &v, const Ptr< VectorBase< Scalar > > &y)
	Element-wise conjugate product update: y(i) += alpha * conj(x(i)) * v(i), i = 0...y->space()->dim()-1.
template<class Scalar>
void	ele_wise_scale (const VectorBase< Scalar > &x, const Ptr< VectorBase< Scalar > > &y)
	Element-wise scaling: y(i) *= x(i), i = 0...y->space()->dim()-1.
template<class Scalar>
void	Vp_StVtV (const Ptr< VectorBase< Scalar > > &y, const Scalar &alpha, const VectorBase< Scalar > &x, const VectorBase< Scalar > &v)
	Element-wise product update: y(i) += alpha * x(i) * v(i), i = 0...y->space()->dim()-1.
template<class Scalar>
void	ele_wise_prod_update (const Scalar &alpha, const VectorBase< Scalar > &x, const Ptr< VectorBase< Scalar > > &y)
	Element-wise product update: y(i) = alpha x(i), i = 0...y->space()->dim()-1.
template<class Scalar>
void	pair_wise_max_update (const Scalar &alpha, const VectorBase< Scalar > &x, const Ptr< VectorBase< Scalar > > &y)
	Element-wise maximum update: y(i) = alpha * max(x(i), y(i)), i = 0...y->space()->dim()-1.
template<class Scalar>
void	Vt_StV (const Ptr< VectorBase< Scalar > > &y, const Scalar &alpha, const VectorBase< Scalar > &x)
	Element-wise product update: y(i) = alpha x(i), i = 0...y->space()->dim()-1.
template<class Scalar>
void	ele_wise_divide (const Scalar &alpha, const VectorBase< Scalar > &x, const VectorBase< Scalar > &v, const Ptr< VectorBase< Scalar > > &y)
	Element-wise division update: y(i) += alpha * x(i) / v(i), i = 0...y->space()->dim()-1.
template<class Scalar>
void	linear_combination (const ArrayView< const Scalar > &alpha, const ArrayView< const Ptr< const VectorBase< Scalar > > > &x, const Scalar &beta, const Ptr< VectorBase< Scalar > > &y)
	Linear combination: y(i) = betay(i) + sum( alpha[k]x[k](i), k=0...m-1 ), i = 0...y->space()->dim()-1.
template<class Scalar>
void	seed_randomize (unsigned int s)
	Seed the random number generator used in randomize().
template<class Scalar>
void	randomize (Scalar l, Scalar u, const Ptr< VectorBase< Scalar > > &v)
	Random vector generation: v(i) = rand(l,u), , i = 1...v->space()->dim().
template<class Scalar>
void	assign (const Ptr< VectorBase< Scalar > > &y, const Scalar &alpha)
	Assign all elements to a scalar: y(i) = alpha, i = 0...y->space()->dim()-1.
template<class Scalar>
void	assign (const Ptr< VectorBase< Scalar > > &y, const VectorBase< Scalar > &x)
	Vector assignment: y(i) = x(i), i = 0...y->space()->dim()-1.
template<class Scalar>
void	Vp_S (const Ptr< VectorBase< Scalar > > &y, const Scalar &alpha)
	Add a scalar to all elements: y(i) += alpha, i = 0...y->space()->dim()-1.
template<class Scalar>
void	Vt_S (const Ptr< VectorBase< Scalar > > &y, const Scalar &alpha)
	Scale all elements by a scalar: y(i) *= alpha, i = 0...y->space()->dim()-1.
template<class Scalar>
void	V_StV (const Ptr< VectorBase< Scalar > > &y, const Scalar &alpha, const VectorBase< Scalar > &x)
	Assign scaled vector: y(i) = alpha * x(i), i = 0...y->space()->dim()-1.
template<class Scalar>
void	Vp_StV (const Ptr< VectorBase< Scalar > > &y, const Scalar &alpha, const VectorBase< Scalar > &x)
	AXPY: y(i) = alpha * x(i) + y(i), i = 0...y->space()->dim()-1.
template<class Scalar>
void	Vp_V (const Ptr< VectorBase< Scalar > > &y, const VectorBase< Scalar > &x, const Scalar &beta=static_cast< Scalar >(1.0))
	y(i) = x(i) + beta*y(i), i = 0...y->space()->dim()-1.
template<class Scalar>
void	V_V (const Ptr< VectorBase< Scalar > > &y, const VectorBase< Scalar > &x)
	y(i) = x(i), i = 0...y->space()->dim()-1.
template<class Scalar>
void	V_S (const Ptr< VectorBase< Scalar > > &y, const Scalar &alpha)
	y(i) = alpha, i = 0...y->space()->dim()-1.
template<class Scalar>
void	V_VpV (const Ptr< VectorBase< Scalar > > &z, const VectorBase< Scalar > &x, const VectorBase< Scalar > &y)
	z(i) = x(i) + y(i), i = 0...z->space()->dim()-1.
template<class Scalar>
void	V_VmV (const Ptr< VectorBase< Scalar > > &z, const VectorBase< Scalar > &x, const VectorBase< Scalar > &y)
	z(i) = x(i) - y(i), i = 0...z->space()->dim()-1.
template<class Scalar>
void	V_StVpV (const Ptr< VectorBase< Scalar > > &z, const Scalar &alpha, const VectorBase< Scalar > &x, const VectorBase< Scalar > &y)
	z(i) = alpha*x(i) + y(i), i = 0...z->space()->dim()-1.
template<class Scalar>
void	V_VpStV (const Ptr< VectorBase< Scalar > > &z, const VectorBase< Scalar > &x, const Scalar &alpha, const VectorBase< Scalar > &y)
	z(i) = x(i) + alpha*y(i), i = 0...z->space()->dim()-1.
template<class Scalar>
void	V_StVpStV (const Ptr< VectorBase< Scalar > > &z, const Scalar &alpha, const VectorBase< Scalar > &x, const Scalar &beta, const VectorBase< Scalar > &y)
	z(i) = alphax(i) + betay(i), i = 0...z->space()->dim()-1.
template<class Scalar>
Scalar	min (const VectorBase< Scalar > &x)
	Min element: result = min{ x(i), i = 0...x.space()->dim()-1 } .
template<class Scalar>
void	min (const VectorBase< Scalar > &x, const Ptr< Scalar > &maxEle, const Ptr< Ordinal > &maxIndex)
	Min element and its index: Returns maxEle = x(k) and maxIndex = k such that x(k) <= x(i) for all i = 0...x.space()->dim()-1.
template<class Scalar>
void	minGreaterThanBound (const VectorBase< Scalar > &x, const Scalar &bound, const Ptr< Scalar > &minEle, const Ptr< Ordinal > &minIndex)
	Minimum element greater than some bound and its index: Returns minEle = x(k) and minIndex = k such that x(k) <= x(i) for all i where x(i) > bound.
template<class Scalar>
Scalar	max (const VectorBase< Scalar > &x)
	Max element: result = max{ x(i), i = 1...n } .
template<class Scalar>
void	max (const VectorBase< Scalar > &x, const Ptr< Scalar > &maxEle, const Ptr< Ordinal > &maxIndex)
	Max element and its index: Returns maxEle = x(k) and maxIndex = k such that x(k) >= x(i) for i = 0...x.space()->dim()-1.
template<class Scalar>
void	maxLessThanBound (const VectorBase< Scalar > &x, const Scalar &bound, const Ptr< Scalar > &maxEle, const Ptr< Ordinal > &maxIndex)
	Max element less than bound and its index: Returns maxEle = x(k) and maxIndex = k such that x(k) >= x(i) for all i where x(i) < bound.
Related Symbols inherited from Thyra::MultiVectorBase< Scalar >
template<class Scalar>
void	applyOp (const RTOpPack::RTOpT< Scalar > &primary_op, const ArrayView< const Ptr< const MultiVectorBase< Scalar > > > &multi_vecs, const ArrayView< const Ptr< MultiVectorBase< Scalar > > > &targ_multi_vecs, const ArrayView< const Ptr< RTOpPack::ReductTarget > > &reduct_objs, const Ordinal primary_global_offset=0)
	Apply a reduction/transformation operator column by column and return an array of the reduction objects.
template<class Scalar>
void	applyOp (const RTOpPack::RTOpT< Scalar > &primary_op, const RTOpPack::RTOpT< Scalar > &secondary_op, const ArrayView< const Ptr< const MultiVectorBase< Scalar > > > &multi_vecs, const ArrayView< const Ptr< MultiVectorBase< Scalar > > > &targ_multi_vecs, const Ptr< RTOpPack::ReductTarget > &reduct_obj, const Ordinal primary_global_offset=0)
	Apply a reduction/transformation operator column by column and reduce the intermediate reduction objects into one reduction object.
template<class Scalar>
void	norms (const MultiVectorBase< Scalar > &V, const ArrayView< typename ScalarTraits< Scalar >::magnitudeType > &norms)
	Column-wise multi-vector natural norm.
template<class Scalar, class NormOp>
void	reductions (const MultiVectorBase< Scalar > &V, const NormOp &op, const ArrayView< typename ScalarTraits< Scalar >::magnitudeType > &norms)
	Column-wise multi-vector reductions.
template<class Scalar>
void	norms_1 (const MultiVectorBase< Scalar > &V, const ArrayView< typename ScalarTraits< Scalar >::magnitudeType > &norms)
	Column-wise multi-vector one norm.
template<class Scalar>
void	norms_2 (const MultiVectorBase< Scalar > &V, const ArrayView< typename ScalarTraits< Scalar >::magnitudeType > &norms)
	Column-wise multi-vector 2 (Euclidean) norm.
template<class Scalar>
void	norms_inf (const MultiVectorBase< Scalar > &V, const ArrayView< typename ScalarTraits< Scalar >::magnitudeType > &norms)
	Column-wise multi-vector infinity norm.
template<class Scalar>
Array< typename ScalarTraits< Scalar >::magnitudeType >	norms_inf (const MultiVectorBase< Scalar > &V)
	Column-wise multi-vector infinity norm.
template<class Scalar>
void	dots (const MultiVectorBase< Scalar > &V1, const MultiVectorBase< Scalar > &V2, const ArrayView< Scalar > &dots)
	Multi-vector dot product.
template<class Scalar>
void	sums (const MultiVectorBase< Scalar > &V, const ArrayView< Scalar > &sums)
	Multi-vector column sum.
template<class Scalar>
ScalarTraits< Scalar >::magnitudeType	norm_1 (const MultiVectorBase< Scalar > &V)
	Take the induced matrix one norm of a multi-vector.
template<class Scalar>
void	scale (Scalar alpha, const Ptr< MultiVectorBase< Scalar > > &V)
	V = alpha*V.
template<class Scalar>
void	scaleUpdate (const VectorBase< Scalar > &a, const MultiVectorBase< Scalar > &U, const Ptr< MultiVectorBase< Scalar > > &V)
	A*U + V -> V (where A is a diagonal matrix with diagonal a).
template<class Scalar>
void	assign (const Ptr< MultiVectorBase< Scalar > > &V, Scalar alpha)
	V = alpha.
template<class Scalar>
void	assign (const Ptr< MultiVectorBase< Scalar > > &V, const MultiVectorBase< Scalar > &U)
	V = U.
template<class Scalar>
void	update (Scalar alpha, const MultiVectorBase< Scalar > &U, const Ptr< MultiVectorBase< Scalar > > &V)
	alpha*U + V -> V.
template<class Scalar>
void	update (const ArrayView< const Scalar > &alpha, Scalar beta, const MultiVectorBase< Scalar > &U, const Ptr< MultiVectorBase< Scalar > > &V)
	alpha[j]betaU(j) + V(j) - > V(j), for j = 0 ,,,
template<class Scalar>
void	update (const MultiVectorBase< Scalar > &U, const ArrayView< const Scalar > &alpha, Scalar beta, const Ptr< MultiVectorBase< Scalar > > &V)
	U(j) + alpha[j]betaV(j) - > V(j), for j = 0 ,,, U.domain()->dim()-1.
template<class Scalar>
void	linear_combination (const ArrayView< const Scalar > &alpha, const ArrayView< const Ptr< const MultiVectorBase< Scalar > > > &X, const Scalar &beta, const Ptr< MultiVectorBase< Scalar > > &Y)
	Y.col(j)(i) = betaY.col(j)(i) + sum( alpha[k]X[k].col(j)(i), k=0...m-1 ), for i = 0...Y->range()->dim()-1, j = 0...Y->domain()->dim()-1.
template<class Scalar>
void	randomize (Scalar l, Scalar u, const Ptr< MultiVectorBase< Scalar > > &V)
	Generate a random multi-vector with elements uniformly distributed elements.
template<class Scalar>
void	Vt_S (const Ptr< MultiVectorBase< Scalar > > &Z, const Scalar &alpha)
	Z(i,j) *= alpha, i = 0...Z->range()->dim()-1, j = 0...Z->domain()->dim()-1.
template<class Scalar>
void	Vp_S (const Ptr< MultiVectorBase< Scalar > > &Z, const Scalar &alpha)
	Z(i,j) += alpha, i = 0...Z->range()->dim()-1, j = 0...Z->domain()->dim()-1.
template<class Scalar>
void	Vp_V (const Ptr< MultiVectorBase< Scalar > > &Z, const MultiVectorBase< Scalar > &X)
	Z(i,j) += X(i,j), i = 0...Z->range()->dim()-1, j = 0...Z->domain()->dim()-1.
template<class Scalar>
void	V_VpV (const Ptr< MultiVectorBase< Scalar > > &Z, const MultiVectorBase< Scalar > &X, const MultiVectorBase< Scalar > &Y)
	Z(i,j) = X(i,j) + Y(i,j), i = 0...Z->range()->dim()-1, j = 0...Z->domain()->dim()-1.
template<class Scalar>
void	V_VmV (const Ptr< MultiVectorBase< Scalar > > &Z, const MultiVectorBase< Scalar > &X, const MultiVectorBase< Scalar > &Y)
	Z(i,j) = X(i,j) - Y(i,j), i = 0...Z->range()->dim()-1, j = 0...Z->domain()->dim()-1.
template<class Scalar>
void	V_StVpV (const Ptr< MultiVectorBase< Scalar > > &Z, const Scalar &alpha, const MultiVectorBase< Scalar > &X, const MultiVectorBase< Scalar > &Y)
	Z(i,j) = alpha*X(i,j) + Y(i), i = 0...z->space()->dim()-1, , j = 0...Z->domain()->dim()-1.
Related Symbols inherited from Thyra::LinearOpBase< Scalar >
template<class Scalar>
bool	isFullyUninitialized (const LinearOpBase< Scalar > &M)
	Determines if a linear operator is in the "Fully Uninitialized" state or not.
template<class Scalar>
bool	isPartiallyInitialized (const LinearOpBase< Scalar > &M)
	Determines if a linear operator is in the "Partially Initialized" state or not.
template<class Scalar>
bool	isFullyInitialized (const LinearOpBase< Scalar > &M)
	Determines if a linear operator is in the "Fully Initialized" state or not.
template<class Scalar>
bool	opSupported (const LinearOpBase< Scalar > &M, EOpTransp M_trans)
	Determines if an operation is supported for a single scalar type.
template<class Scalar>
void	apply (const LinearOpBase< Scalar > &M, const EOpTransp M_trans, const MultiVectorBase< Scalar > &X, const Ptr< MultiVectorBase< Scalar > > &Y, const Scalar alpha=static_cast< Scalar >(1.0), const Scalar beta=static_cast< Scalar >(0.0))
	Non-member function call for M.apply(...).
void	apply (const LinearOpBase< double > &M, const EOpTransp M_trans, const MultiVectorBase< double > &X, const Ptr< MultiVectorBase< double > > &Y, const double alpha=1.0, const double beta=0.0)
	Calls apply<double>(...).

Detailed Description

template<class Scalar>
class Thyra::DefaultClusteredSpmdProductVector< Scalar >

Concrete implementation of a clustered Spmd-based product vector.

ToDo: Finish documentation!

The default constructor is made private to avoid accidental default construction.

Definition at line 31 of file Thyra_DefaultClusteredSpmdProductVector_decl.hpp.

Constructor & Destructor Documentation

◆ DefaultClusteredSpmdProductVector() [1/2]

template<class Scalar>

Thyra::DefaultClusteredSpmdProductVector< Scalar >::DefaultClusteredSpmdProductVector ( )

Constructs to uninitialized (see postconditions for uninitialize()).

Definition at line 30 of file Thyra_DefaultClusteredSpmdProductVector_def.hpp.

◆ DefaultClusteredSpmdProductVector() [2/2]

template<class Scalar>

Thyra::DefaultClusteredSpmdProductVector< Scalar >::DefaultClusteredSpmdProductVector	(	const Teuchos::RCP< const DefaultClusteredSpmdProductVectorSpace< Scalar > > &	productSpace,
		const Teuchos::RCP< VectorBase< Scalar > >	vecs[] )

Constructs to initialized (calls initialize()).

Definition at line 37 of file Thyra_DefaultClusteredSpmdProductVector_def.hpp.

Member Function Documentation

◆ initialize()

template<class Scalar>

void Thyra::DefaultClusteredSpmdProductVector< Scalar >::initialize	(	const Teuchos::RCP< const DefaultClusteredSpmdProductVectorSpace< Scalar > > &	productSpace,
		const Teuchos::RCP< VectorBase< Scalar > >	vecs[] )

Initialize.

ToDo: Finish documentation.

Definition at line 47 of file Thyra_DefaultClusteredSpmdProductVector_def.hpp.

◆ uninitialize()

template<class Scalar>

void Thyra::DefaultClusteredSpmdProductVector< Scalar >::uninitialize	(	Teuchos::RCP< const DefaultClusteredSpmdProductVectorSpace< Scalar > > *	productSpace = NULL,
		Teuchos::RCP< VectorBase< Scalar > >	vecs[] = NULL )

Uninitialize.

ToDo: Finish documentation.

Definition at line 67 of file Thyra_DefaultClusteredSpmdProductVector_def.hpp.

◆ getNonconstVectorBlock()

template<class Scalar>

Teuchos::RCP< VectorBase< Scalar > > Thyra::DefaultClusteredSpmdProductVector< Scalar >::getNonconstVectorBlock ( const int k )

virtual

Implements Thyra::ProductVectorBase< Scalar >.

Definition at line 85 of file Thyra_DefaultClusteredSpmdProductVector_def.hpp.

◆ getVectorBlock()

template<class Scalar>

Teuchos::RCP< const VectorBase< Scalar > > Thyra::DefaultClusteredSpmdProductVector< Scalar >::getVectorBlock ( const int k ) const

virtual

Implements Thyra::ProductVectorBase< Scalar >.

Definition at line 95 of file Thyra_DefaultClusteredSpmdProductVector_def.hpp.

◆ productSpace()

template<class Scalar>

Teuchos::RCP< const ProductVectorSpaceBase< Scalar > > Thyra::DefaultClusteredSpmdProductVector< Scalar >::productSpace ( ) const

virtual

Implements Thyra::ProductMultiVectorBase< Scalar >.

Definition at line 108 of file Thyra_DefaultClusteredSpmdProductVector_def.hpp.

◆ blockIsConst()

template<class Scalar>

bool Thyra::DefaultClusteredSpmdProductVector< Scalar >::blockIsConst ( const int k ) const

virtual

Implements Thyra::ProductMultiVectorBase< Scalar >.

Definition at line 115 of file Thyra_DefaultClusteredSpmdProductVector_def.hpp.

◆ getNonconstMultiVectorBlock()

template<class Scalar>

Teuchos::RCP< MultiVectorBase< Scalar > > Thyra::DefaultClusteredSpmdProductVector< Scalar >::getNonconstMultiVectorBlock ( const int k )

virtual

Implements Thyra::ProductMultiVectorBase< Scalar >.

Definition at line 125 of file Thyra_DefaultClusteredSpmdProductVector_def.hpp.

◆ getMultiVectorBlock()

template<class Scalar>

Teuchos::RCP< const MultiVectorBase< Scalar > > Thyra::DefaultClusteredSpmdProductVector< Scalar >::getMultiVectorBlock ( const int k ) const

virtual

Implements Thyra::ProductMultiVectorBase< Scalar >.

Definition at line 133 of file Thyra_DefaultClusteredSpmdProductVector_def.hpp.

◆ space()

template<class Scalar>

Teuchos::RCP< const VectorSpaceBase< Scalar > > Thyra::DefaultClusteredSpmdProductVector< Scalar >::space ( ) const

virtual

Implements Thyra::VectorBase< Scalar >.

Definition at line 144 of file Thyra_DefaultClusteredSpmdProductVector_def.hpp.

◆ applyOpImpl()

template<class Scalar>

void Thyra::DefaultClusteredSpmdProductVector< Scalar >::applyOpImpl	(	const RTOpPack::RTOpT< Scalar > &	op,
		const ArrayView< const Ptr< const VectorBase< Scalar > > > &	vecs,
		const ArrayView< const Ptr< VectorBase< Scalar > > > &	targ_vecs,
		const Ptr< RTOpPack::ReductTarget > &	reduct_obj,
		const Ordinal	global_offset ) const

protectedvirtual

Implements Thyra::VectorBase< Scalar >.

Definition at line 154 of file Thyra_DefaultClusteredSpmdProductVector_def.hpp.

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

Overridden protected members from VectorBase

Constructors/initializers/accessors

Overridden from ProductVectorBase

Overridden from ProductMultiVectorBase

Overridden from VectorBase

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ DefaultClusteredSpmdProductVector() [1/2]

◆ DefaultClusteredSpmdProductVector() [2/2]

Member Function Documentation

◆ initialize()

◆ uninitialize()

◆ getNonconstVectorBlock()

◆ getVectorBlock()

◆ productSpace()

◆ blockIsConst()

◆ getNonconstMultiVectorBlock()

◆ getMultiVectorBlock()

◆ space()

◆ applyOpImpl()