Efficient concrete implementation subclass for SPMD vectors. More...

#include <Thyra_DefaultSpmdVector_decl.hpp>

Inheritance diagram for Thyra::DefaultSpmdVector< Scalar >:

Constructors/initializers
	DefaultSpmdVector ()
	Construct to uninitialized. More...

	DefaultSpmdVector (const RCP< const SpmdVectorSpaceBase< Scalar > > &spmdSpace, const ArrayRCP< Scalar > &localValues, const Ordinal stride)
	Calls `initialize()`. More...

void	initialize (const RCP< const SpmdVectorSpaceBase< Scalar > > &spmdSpace, const ArrayRCP< Scalar > &localValues, const Ordinal stride)
	Initialize. More...

void	uninitialize (RCP< const SpmdVectorSpaceBase< Scalar > > spmdSpace=NULL, ArrayRCP< Scalar > localValues=NULL, Ordinal *stride=NULL)
	Set to an uninitialized state. More...

Accessors (inlined for minimal overhead)
ArrayRCP< Scalar >	getRCPtr ()

ArrayRCP< const Scalar >	getRCPtr () const

Scalar *	getPtr ()

const Scalar *	getPtr () const

Ordinal	getStride () const

Overridden from SpmdMultiVectorBase
RCP< const SpmdVectorSpaceBase < Scalar > >	spmdSpaceImpl () const

Overridden from SpmdVectorBase
void	getNonconstLocalVectorDataImpl (const Ptr< ArrayRCP< Scalar > > &localValues)

void	getLocalVectorDataImpl (const Ptr< ArrayRCP< const Scalar > > &localValues) const

Additional Inherited Members
Public Member Functions inherited from Thyra::SpmdVectorDefaultBase< Scalar >
virtual void	applyOpImplWithComm (const Ptr< const Teuchos::Comm< Ordinal > > &comm, 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
	Implementation of applyOpImpl(...) that uses an input Comm. More...

	SpmdVectorDefaultBase ()

std::string	description () const

Teuchos::RCP< const VectorSpaceBase< Scalar > >	space () const
	Returns `this->spmdSpace()`. More...

Public Member Functions inherited from Thyra::SpmdVectorBase< Scalar >
RTOpPack::SubVectorView< Scalar >	getNonconstLocalSubVector ()
	Get a non-const generalized view of local vector data. More...

RTOpPack::ConstSubVectorView < Scalar >	getLocalSubVector () const
	Get a const generalized view of local vector data. More...

void	getNonconstLocalData (const Ptr< ArrayRCP< Scalar > > &localValues)
	Returns a non-`const` pointer to the beginning of the local vector data. More...

void	getLocalData (const Ptr< ArrayRCP< const Scalar > > &localValues) const
	Returns a `const` pointer to the beginning of the local vector data. More...

Public Member Functions inherited from Thyra::VectorBase< Scalar >
void	assign (const VectorBase< Scalar > &x)
	Vector assignment: More...

void	randomize (Scalar l, Scalar u)
	Random vector generation: More...

void	update (Scalar alpha, const VectorBase< Scalar > &x)
	AXPY: More...

void	linear_combination (const ArrayView< const Scalar > &alpha, const ArrayView< const Ptr< const VectorBase< Scalar > > > &x, const Scalar &beta)
	Linear combination: More...

Scalar	dot (const VectorBase< Scalar > &x) const
	Euclidean dot product: `result = x^H * this`. More...

Teuchos::ScalarTraits< Scalar > ::magnitudeType	norm_1 () const
	One (1) norm: `result = \|\|v\|\|1`. More...

Teuchos::ScalarTraits< Scalar > ::magnitudeType	norm_2 () const
	Euclidean (2) norm: `result = \|\|v\|\|2`. More...

Teuchos::ScalarTraits< Scalar > ::magnitudeType	norm_2 (const VectorBase< Scalar > &x) const
	Weighted Euclidean (2) norm: `result = \|\|v\|\|2`. More...

Teuchos::ScalarTraits< Scalar > ::magnitudeType	norm_inf () const
	Infinity norm: `result = \|\|v\|\|inf`. More...

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(). More...

void	acquireDetachedView (const Range1D &rng, RTOpPack::ConstSubVectorView< Scalar > *sub_vec) const
	Calls acquireDetachedVectorViewImpl(). More...

void	releaseDetachedView (RTOpPack::ConstSubVectorView< Scalar > *sub_vec) const
	Calls releaseDetachedVectorViewImpl(). More...

void	acquireDetachedView (const Range1D &rng, RTOpPack::SubVectorView< Scalar > *sub_vec)
	Calls acquireNonconstDetachedVectorViewImpl(). More...

void	commitDetachedView (RTOpPack::SubVectorView< Scalar > *sub_vec)
	Calls commitDetachedView(). More...

void	setSubVector (const RTOpPack::SparseSubVectorT< Scalar > &sub_vec)
	Calls setSubVectorImpl(). More...

Public Member Functions inherited from Thyra::MultiVectorBase< Scalar >
void	assign (Scalar alpha)
	V = alpha. More...

void	assign (const MultiVectorBase< Scalar > &mv)
	V = mv. More...

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),` More...

void	dots (const MultiVectorBase< Scalar > &mv, const ArrayView< Scalar > &prods) const
	Column-wise Euclidean dot product. More...

void	norms_1 (const ArrayView< typename ScalarTraits< Scalar >::magnitudeType > &norms) const
	Column-wise 1-norms. More...

void	norms_2 (const ArrayView< typename ScalarTraits< Scalar >::magnitudeType > &norms) const
	Column-wise 2-norms. More...

void	norms_inf (const ArrayView< typename ScalarTraits< Scalar >::magnitudeType > &norms) const
	Column-wise infinity-norms. More...

RCP< const VectorBase< Scalar > >	col (Ordinal j) const
	Calls colImpl(). More...

RCP< VectorBase< Scalar > >	col (Ordinal j)
	Calls nonconstColImpl(). More...

RCP< const MultiVectorBase < Scalar > >	subView (const Range1D &colRng) const
	Calls contigSubViewImpl(). More...

RCP< MultiVectorBase< Scalar > >	subView (const Range1D &colRng)
	Calls nonconstContigSubViewImpl(). More...

RCP< const MultiVectorBase < Scalar > >	subView (const ArrayView< const int > &cols) const
	nonContigSubViewImpl(). More...

RCP< MultiVectorBase< Scalar > >	subView (const ArrayView< const int > &cols)
	nonconstNonContigSubViewImpl(). More...

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(). More...

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(). More...

void	acquireDetachedView (const Range1D &rowRng, const Range1D &colRng, RTOpPack::ConstSubMultiVectorView< Scalar > *sub_mv) const
	Calls acquireDetachedMultiVectorViewImpl(). More...

void	releaseDetachedView (RTOpPack::ConstSubMultiVectorView< Scalar > *sub_mv) const
	Calls releaseDetachedMultiVectorViewImpl(). More...

void	acquireDetachedView (const Range1D &rowRng, const Range1D &colRng, RTOpPack::SubMultiVectorView< Scalar > *sub_mv)
	Calls acquireNonconstDetachedMultiVectorViewImpl(). More...

void	commitDetachedView (RTOpPack::SubMultiVectorView< Scalar > *sub_mv)
	Calls commitNonconstDetachedMultiVectorViewImpl(). More...

RCP< const LinearOpBase< Scalar > >	clone () const
	This function is simply overridden to return `this->clone_mv()`. More...

Public Member Functions inherited from Thyra::LinearOpBase< Scalar >
bool	opSupported (EOpTransp M_trans) const
	Return if the `M_trans` operation of `apply()` is supported or not. More...

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`. More...

Public Member Functions inherited from Thyra::RowStatLinearOpBase< Scalar >
bool	rowStatIsSupported (const RowStatLinearOpBaseUtils::ERowStat rowStat) const
	Determine if a given row stat is supported. More...

void	getRowStat (const RowStatLinearOpBaseUtils::ERowStat rowStat, const Ptr< VectorBase< Scalar > > &rowStatVec) const
	Get some statistics about a supported row. More...

Public Member Functions inherited from Thyra::ScaledLinearOpBase< Scalar >
bool	supportsScaleLeft () const
	Determines if this objects supports left scaling. More...

bool	supportsScaleRight () const
	Determines if this objects supports right scaling. More...

void	scaleLeft (const VectorBase< Scalar > &row_scaling)
	Left scales operator with diagonal scaling operator. More...

void	scaleRight (const VectorBase< Scalar > &col_scaling)
	Right scales operator with diagonal scaling operator. More...

Public Member Functions inherited from Thyra::SpmdMultiVectorBase< Scalar >
RCP< const SpmdVectorSpaceBase < Scalar > >	spmdSpace () const
	Returns the SPMD vector space object for the range of `*this` multi-vector. More...

RTOpPack::SubMultiVectorView < Scalar >	getNonconstLocalSubMultiVector ()
	Get a non-const generalized view of local multi-vector data. More...

RTOpPack::ConstSubMultiVectorView < Scalar >	getLocalSubMultiVector () const
	Get a const generalized view of local multi-vector data. More...

void	getNonconstLocalData (const Ptr< ArrayRCP< Scalar > > &localValues, const Ptr< Ordinal > &leadingDim)
	Returns a non-`const` pointer to a Fortran-style view of the local multi-vector data. More...

void	getLocalData (const Ptr< ArrayRCP< const Scalar > > &localValues, const Ptr< Ordinal > &leadingDim) const
	Returns a `const` pointer to a Fortran-style view of the local multi-vector data. More...

Public Member Functions inherited from Thyra::VectorDefaultBase< Scalar >
virtual void	describe (Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel) const
	Generates a default outputting for all vectors. More...

virtual RCP< const VectorSpaceBase< Scalar > >	range () const
	Returns `this->space()`. More...

virtual RCP< const VectorSpaceBase< Scalar > >	domain () const
	Returns a `DefaultSerialVectorSpace` object with dimension 1. More...

RCP< MultiVectorBase< Scalar > >	clone_mv () const
	Returns `this->clone_v()`. More...

RCP< VectorBase< Scalar > >	clone_v () const
	Simply creates a new vector and copies the contents from `*this`. More...

Protected Member Functions inherited from Thyra::SpmdVectorDefaultBase< Scalar >
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
	Calls applyOpImplWithComm(null,op,...). More...

void	acquireDetachedVectorViewImpl (const Range1D &rng, RTOpPack::ConstSubVectorView< Scalar > *sub_vec) const
	Implemented through `this->getLocalData()` More...

void	releaseDetachedVectorViewImpl (RTOpPack::ConstSubVectorView< Scalar > *sub_vec) const
	Implemented through `this->freeLocalData()` More...

void	acquireNonconstDetachedVectorViewImpl (const Range1D &rng, RTOpPack::SubVectorView< Scalar > *sub_vec)
	Implemented through `this->getLocalData()` More...

void	commitNonconstDetachedVectorViewImpl (RTOpPack::SubVectorView< Scalar > *sub_vec)
	Implemented through `this->commitLocalData()` More...

RTOpPack::SubMultiVectorView < Scalar >	getNonconstLocalSubMultiVectorImpl ()

RTOpPack::ConstSubMultiVectorView < Scalar >	getLocalSubMultiVectorImpl () const

void	getNonconstLocalMultiVectorDataImpl (const Ptr< ArrayRCP< Scalar > > &localValues, const Ptr< Ordinal > &leadingDim)

void	getLocalMultiVectorDataImpl (const Ptr< ArrayRCP< const Scalar > > &localValues, const Ptr< Ordinal > &leadingDim) const

RTOpPack::SubVectorView< Scalar >	getNonconstLocalSubVectorImpl ()
	Virtual implementation for getNonconstLocalSubVector(). More...

RTOpPack::ConstSubVectorView < Scalar >	getLocalSubVectorImpl () const
	Virtual implementation for getLocalSubVector(). More...

virtual void	updateSpmdSpace ()
	Subclasses must call this function whenever the structure of the VectorSpaceBase changes. More...

Protected Member Functions inherited from Thyra::SpmdVectorBase< Scalar >
Protected Member Functions inherited from Thyra::VectorBase< Scalar >
Protected Member Functions inherited from Thyra::MultiVectorBase< Scalar >
void	absRowSum (const Teuchos::Ptr< Thyra::VectorBase< Scalar > > &output) const

void	absColSum (const Teuchos::Ptr< Thyra::VectorBase< Scalar > > &output) const

virtual RCP< const VectorBase < Scalar > >	colImpl (Ordinal j) const
	Return a non-changeable view of a constituent column vector. More...

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)

Protected Member Functions inherited from Thyra::LinearOpBase< Scalar >
Protected Member Functions inherited from Thyra::RowStatLinearOpBase< Scalar >
Protected Member Functions inherited from Thyra::ScaledLinearOpBase< Scalar >
Protected Member Functions inherited from Thyra::SpmdMultiVectorBase< Scalar >
Protected Member Functions inherited from Thyra::VectorDefaultBase< Scalar >
virtual void	assignVecImpl (const VectorBase< Scalar > &x)
	Default implementation of assign(vector) using RTOps. More...

virtual void	randomizeImpl (Scalar l, Scalar u)
	Default implementation of randomize using RTOps. More...

virtual void	absImpl (const VectorBase< Scalar > &x)
	Default implementation of abs using RTOps. More...

virtual void	reciprocalImpl (const VectorBase< Scalar > &x)
	Default implementation of reciprocal using RTOps. More...

virtual void	eleWiseScaleImpl (const VectorBase< Scalar > &x)
	Default implementation of ele_wise_scale using RTOps. More...

virtual void	updateVecImpl (Scalar alpha, const VectorBase< Scalar > &x)
	Default implementation of update using RTOps. More...

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. More...

virtual Scalar	dotImpl (const VectorBase< Scalar > &x) const
	Default implementation of dot using RTOps. More...

virtual Teuchos::ScalarTraits < Scalar >::magnitudeType	norm1Impl () const
	Default implementation of norm_1 using RTOps. More...

virtual Teuchos::ScalarTraits < Scalar >::magnitudeType	norm2Impl () const
	Default implementation of norm_2 using RTOps. More...

virtual Teuchos::ScalarTraits < Scalar >::magnitudeType	norm2WeightedImpl (const VectorBase< Scalar > &x) const
	Default implementation of norm_2 (weighted) using RTOps. More...

virtual Teuchos::ScalarTraits < Scalar >::magnitudeType	normInfImpl () const
	Default implementation of norm_inf using RTOps. More...

virtual RCP< VectorBase< Scalar > >	nonconstColImpl (Ordinal j)
	Returns `Teuchos::rcp(this,false)`. More...

virtual RCP< const MultiVectorBase< Scalar > >	contigSubViewImpl (const Range1D &col_rng) const
	Returns `Teuchos::rcp(this,false)`. More...

virtual RCP< MultiVectorBase < Scalar > >	nonconstContigSubViewImpl (const Range1D &col_rng)
	Returns `Teuchos::rcp(this,false)`. More...

virtual RCP< const MultiVectorBase< Scalar > >	nonContigSubViewImpl (const ArrayView< const int > &cols) const
	Returns `Teuchos::rcp(this,false)`. More...

virtual RCP< MultiVectorBase < Scalar > >	nonconstNonContigSubViewImpl (const ArrayView< const int > &cols)
	Returns `Teuchos::rcp(this,false)`. More...

virtual void	acquireDetachedMultiVectorViewImpl (const Range1D &rowRng, const Range1D &colRng, RTOpPack::ConstSubMultiVectorView< Scalar > *sub_mv) const
	Implemented in terms of `this->acquireDetachedView()`. More...

virtual void	releaseDetachedMultiVectorViewImpl (RTOpPack::ConstSubMultiVectorView< Scalar > *sub_mv) const
	Implemented in terms of `this->releaseDetachedView()`. More...

virtual void	acquireNonconstDetachedMultiVectorViewImpl (const Range1D &rowRng, const Range1D &colRng, RTOpPack::SubMultiVectorView< Scalar > *sub_mv)
	Implemented in terms of `this->acquireDetachedView()`. More...

virtual void	commitNonconstDetachedMultiVectorViewImpl (RTOpPack::SubMultiVectorView< Scalar > *sub_mv)
	Implemented in terms of `this->commitDetachedView()`. More...

virtual void	setSubVectorImpl (const RTOpPack::SparseSubVectorT< Scalar > &sub_vec)

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`. More...

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. More...

Protected Member Functions inherited from Thyra::MultiVectorDefaultBase< Scalar >
virtual void	assignImpl (Scalar alpha)
	Default implementation of assign(scalar) using RTOps. More...

virtual void	assignMultiVecImpl (const MultiVectorBase< Scalar > &mv)
	Default implementation of assign(MV) using RTOps. More...

virtual void	scaleImpl (Scalar alpha)
	Default implementation of scale using RTOps. More...

virtual void	updateImpl (Scalar alpha, const MultiVectorBase< Scalar > &mv)
	Default implementation of update using RTOps. More...

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. More...

virtual void	dotsImpl (const MultiVectorBase< Scalar > &mv, const ArrayView< Scalar > &prods) const
	Default implementation of dots using RTOps. More...

virtual void	norms1Impl (const ArrayView< typename ScalarTraits< Scalar >::magnitudeType > &norms) const
	Default implementation of norms_1 using RTOps. More...

virtual void	norms2Impl (const ArrayView< typename ScalarTraits< Scalar >::magnitudeType > &norms) const
	Default implementation of norms_2 using RTOps. More...

virtual void	normsInfImpl (const ArrayView< typename ScalarTraits< Scalar >::magnitudeType > &norms) const
	Default implementation of norms_inf using RTOps. More...

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

Protected Member Functions inherited from Thyra::LinearOpDefaultBase< Scalar >
Related Functions 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)`. More...

template<class Scalar >
Scalar	sum (const VectorBase< Scalar > &v)
	Sum of vector elements: `result = sum( v(i), i = 0...v.space()->dim()-1 )`. More...

template<class Scalar >
Scalar	scalarProd (const VectorBase< Scalar > &x, const VectorBase< Scalar > &y)
	Scalar product `result = <x,y>`. More...

template<class Scalar >
Scalar	inner (const VectorBase< Scalar > &x, const VectorBase< Scalar > &y)
	Inner/Scalar product `result = <x,y>`. More...

template<class Scalar >
Teuchos::ScalarTraits< Scalar > ::magnitudeType	norm (const VectorBase< Scalar > &v)
	Natural norm: `result = sqrt(<v,v>)`. More...

template<class Scalar >
Teuchos::ScalarTraits< Scalar > ::magnitudeType	norm_1 (const VectorBase< Scalar > &v)
	One (1) norm: `result = \|\|v\|\|1`. More...

template<class Scalar >
Teuchos::ScalarTraits< Scalar > ::magnitudeType	norm_2 (const VectorBase< Scalar > &v)
	Euclidean (2) norm: `result = \|\|v\|\|2`. More...

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)) )`. More...

template<class Scalar >
Teuchos::ScalarTraits< Scalar > ::magnitudeType	norm_inf (const VectorBase< Scalar > &v_rhs)
	Infinity norm: `result = \|\|v\|\|inf`. More...

template<class Scalar >
Scalar	dot (const VectorBase< Scalar > &x, const VectorBase< Scalar > &y)
	Dot product: `result = conj(x)'*y`. More...

template<class Scalar >
Scalar	get_ele (const VectorBase< Scalar > &v, Ordinal i)
	Get single element: `result = v(i)`. More...

template<class Scalar >
void	set_ele (Ordinal i, Scalar alpha, const Ptr< VectorBase< Scalar > > &v)
	Set single element: `v(i) = alpha`. More...

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`. More...

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`. More...

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`. More...

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`. More...

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`. More...

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`. More...

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`. More...

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`. More...

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`. More...

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`. More...

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`. More...

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`. More...

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`. More...

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`. More...

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`. More...

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`. More...

template<class Scalar >
void	seed_randomize (unsigned int s)
	Seed the random number generator used in `randomize()`. More...

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()`. More...

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`. More...

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`. More...

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`. More...

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`. More...

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`. More...

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`. More...

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`. More...

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`. More...

template<class Scalar >
void	V_S (const Ptr< VectorBase< Scalar > > &y, const Scalar &alpha)
	`y(i) = alpha, i = 0...y->space()->dim()-1`. More...

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`. More...

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`. More...

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`. More...

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`. More...

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`. More...

template<class Scalar >
Scalar	min (const VectorBase< Scalar > &x)
	Min element: `result = min{ x(i), i = 0...x.space()->dim()-1 }` . More...

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`. More...

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`. More...

template<class Scalar >
Scalar	max (const VectorBase< Scalar > &x)
	Max element: `result = max{ x(i), i = 1...n }` . More...

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`. More...

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`. More...

Related Functions 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. More...

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. More...

template<class Scalar >
void	norms (const MultiVectorBase< Scalar > &V, const ArrayView< typename ScalarTraits< Scalar >::magnitudeType > &norms)
	Column-wise multi-vector natural norm. More...

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. More...

template<class Scalar >
void	norms_1 (const MultiVectorBase< Scalar > &V, const ArrayView< typename ScalarTraits< Scalar >::magnitudeType > &norms)
	Column-wise multi-vector one norm. More...

template<class Scalar >
void	norms_2 (const MultiVectorBase< Scalar > &V, const ArrayView< typename ScalarTraits< Scalar >::magnitudeType > &norms)
	Column-wise multi-vector 2 (Euclidean) norm. More...

template<class Scalar >
void	norms_inf (const MultiVectorBase< Scalar > &V, const ArrayView< typename ScalarTraits< Scalar >::magnitudeType > &norms)
	Column-wise multi-vector infinity norm. More...

template<class Scalar >
Array< typename ScalarTraits < Scalar >::magnitudeType >	norms_inf (const MultiVectorBase< Scalar > &V)
	Column-wise multi-vector infinity norm. More...

template<class Scalar >
void	dots (const MultiVectorBase< Scalar > &V1, const MultiVectorBase< Scalar > &V2, const ArrayView< Scalar > &dots)
	Multi-vector dot product. More...

template<class Scalar >
void	sums (const MultiVectorBase< Scalar > &V, const ArrayView< Scalar > &sums)
	Multi-vector column sum. More...

template<class Scalar >
ScalarTraits< Scalar > ::magnitudeType	norm_1 (const MultiVectorBase< Scalar > &V)
	Take the induced matrix one norm of a multi-vector. More...

template<class Scalar >
void	scale (Scalar alpha, const Ptr< MultiVectorBase< Scalar > > &V)
	V = alpha*V. More...

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). More...

template<class Scalar >
void	assign (const Ptr< MultiVectorBase< Scalar > > &V, Scalar alpha)
	V = alpha. More...

template<class Scalar >
void	assign (const Ptr< MultiVectorBase< Scalar > > &V, const MultiVectorBase< Scalar > &U)
	V = U. More...

template<class Scalar >
void	update (Scalar alpha, const MultiVectorBase< Scalar > &U, const Ptr< MultiVectorBase< Scalar > > &V)
	alpha*U + V -> V. More...

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 ,,, More...

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. More...

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`. More...

template<class Scalar >
void	randomize (Scalar l, Scalar u, const Ptr< MultiVectorBase< Scalar > > &V)
	Generate a random multi-vector with elements uniformly distributed elements. More...

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`. More...

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`. More...

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`. More...

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`. More...

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`. More...

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. More...

Related Functions 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. More...

template<class Scalar >
bool	isPartiallyInitialized (const LinearOpBase< Scalar > &M)
	Determines if a linear operator is in the "Partially Initialized" state or not. More...

template<class Scalar >
bool	isFullyInitialized (const LinearOpBase< Scalar > &M)
	Determines if a linear operator is in the "Fully Initialized" state or not. More...

template<class Scalar >
bool	opSupported (const LinearOpBase< Scalar > &M, EOpTransp M_trans)
	Determines if an operation is supported for a single scalar type. More...

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(...)`. More...

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>(...)`. More...

Detailed Description

template<class Scalar>
class Thyra::DefaultSpmdVector< Scalar >

Efficient concrete implementation subclass for SPMD vectors.

This subclass provides a very efficient and very general concrete implementation of a Thyra::VectorBase object for any SPMD platform.

Objects of this type generally should not be constructed directly by a client but instead by using the concrete vector space subclass Thyra::DefaultSpmdVectorSpace and using the function Thyra::createMember().

The storage type can be anything since an ArrayRCP is used to pass in the local values pointer into the constructor and initialize().

Definition at line 69 of file Thyra_DefaultSpmdVector_decl.hpp.

Constructor & Destructor Documentation

template<class Scalar >

Thyra::DefaultSpmdVector< Scalar >::DefaultSpmdVector ( )

Construct to uninitialized.

Definition at line 58 of file Thyra_DefaultSpmdVector_def.hpp.

template<class Scalar >

Thyra::DefaultSpmdVector< Scalar >::DefaultSpmdVector	(	const RCP< const SpmdVectorSpaceBase< Scalar > > &	spmdSpace,
		const ArrayRCP< Scalar > &	localValues,
		const Ordinal	stride
	)

Calls initialize().

Definition at line 64 of file Thyra_DefaultSpmdVector_def.hpp.

Member Function Documentation

template<class Scalar >

void Thyra::DefaultSpmdVector< Scalar >::initialize	(	const RCP< const SpmdVectorSpaceBase< Scalar > > &	spmdSpace,
		const ArrayRCP< Scalar > &	localValues,
		const Ordinal	stride
	)

Initialize.

Parameters

spmdSpace	[in] Smart pointer to `SpmdVectorSpaceBase` object that defines the data distribution for `spmdSpace()` and `space()`.
localValues	[in] Smart pointer to beginning of local strided vector data. This array must be at least of dimension `mpiRangeSpace->localDim()stride` and `(&localValues)[ i*stride ]` gives the local value of the zero-based entry `(i)` where `i=0...spmdSpace()->localSubDim()-1`.
stride	[in] Stride between local vector elements.

Preconditions:

spmdSpace.get()!=NULL
localValues.get()!=NULL
stride != 0

Postconditions:

this->getRCptr().get() == localValues.get()
this->getPtr() == &*localValues
this->getStride() == stride

Definition at line 75 of file Thyra_DefaultSpmdVector_def.hpp.

template<class Scalar >

void Thyra::DefaultSpmdVector< Scalar >::uninitialize	(	RCP< const SpmdVectorSpaceBase< Scalar > > *	spmdSpace = `NULL`,
		ArrayRCP< Scalar > *	localValues = `NULL`,
		Ordinal *	stride = `NULL`
	)