Thyra::VectorSpaceBase< Scalar > Class Template Reference
[C++ code for foundational Thyra operator/vector interfaces]

Abstract interface for objects that represent a space for vectors. More...

#include <Thyra_VectorSpaceBaseDecl.hpp>

Inheritance diagram for Thyra::VectorSpaceBase< Scalar >:

[legend]List of all members.

Public pure virtual functions that must be overridden
virtual Index	dim () const =0
	Return the dimension of the vector space.
virtual bool	isCompatible (const VectorSpaceBase< Scalar > &vecSpc) const =0
	Compare the compatibility of two vector spaces.
virtual RCP< const VectorSpaceFactoryBase< Scalar > >	smallVecSpcFcty () const =0
	Return a VectorSpaceFactoryBase object for the creation of (usually serial) vector spaces with a small dimension.
virtual Scalar	scalarProd (const VectorBase< Scalar > &x, const VectorBase< Scalar > &y) const =0
	Return the scalar product of two vectors in the vector space.
virtual void	scalarProds (const MultiVectorBase< Scalar > &X, const MultiVectorBase< Scalar > &Y, Scalar scalar_prods[]) const =0
	Return the scalar product of each column in two multi-vectors in the vector space.
Public virtual functions with default implementations
virtual bool	isEuclidean () const
	Return if this vector space has a Euclidean (identity) basis in which case the scalar product is the same as the dot product.
virtual bool	hasInCoreView (const Range1D &rng=Range1D(), const EViewType viewType=VIEW_TYPE_DETACHED, const EStrideType strideType=STRIDE_TYPE_NONUNIT) const
	Returns true if this->acquireDetachedView(rng,...) returns a direct view of the range of data requested.
virtual RCP< const VectorSpaceBase< Scalar > >	clone () const
	Clone this object (if supported).
Protected pure virtual functions that must be overridden
virtual RCP< VectorBase< Scalar > >	createMember () const =0
	Create a vector member from the vector space.
virtual RCP< MultiVectorBase< Scalar > >	createMembers (int numMembers) const =0
	Create a set of vector members (a MultiVectorBase) from the vector space.
virtual RCP< VectorBase< Scalar > >	createMemberView (const RTOpPack::SubVectorView< Scalar > &raw_v) const =0
	Create a vector member that is a non-const view of raw vector data.
virtual RCP< const VectorBase< Scalar > >	createMemberView (const RTOpPack::ConstSubVectorView< Scalar > &raw_v) const =0
	Create a vector member that is a const view of raw vector data.
virtual RCP< MultiVectorBase< Scalar > >	createMembersView (const RTOpPack::SubMultiVectorView< Scalar > &raw_mv) const =0
	Create a multi-vector member that is a non-const view of raw multi-vector data.
virtual RCP< const MultiVectorBase< Scalar > >	createMembersView (const RTOpPack::ConstSubMultiVectorView< Scalar > &raw_mv) const =0
	Create a multi-vector member that is a const view of raw multi-vector data.
Related Functions
(Note that these are not member functions.)
RCP< const VectorSpaceBase< Scalar > >	makeHaveOwnership (const RCP< const VectorSpaceBase< Scalar > > &vs)
	Helper function that clones a VectorSpaceBase object if the RCP does not have ownership.
RCP< VectorBase< Scalar > >	createMember (const RCP< const VectorSpaceBase< Scalar > > &vs, const std::string &label="")
	Create a vector member from the vector space.
RCP< VectorBase< Scalar > >	createMember (const VectorSpaceBase< Scalar > &vs, const std::string &label="")
	Calls createMember(Teuchos::rcp(&vs,false)).
RCP< MultiVectorBase< Scalar > >	createMembers (const RCP< const VectorSpaceBase< Scalar > > &vs, int numMembers, const std::string &label="")
	Create a set of vector members (a MultiVectorBase) from the vector space.
RCP< MultiVectorBase< Scalar > >	createMembers (const VectorSpaceBase< Scalar > &vs, int numMembers, const std::string &label="")
	Calls createMembers(Teuchos::rcp(&vs,false),numMembers).
RCP< VectorBase< Scalar > >	createMemberView (const RCP< const VectorSpaceBase< Scalar > > &vs, const RTOpPack::SubVectorView< Scalar > &raw_v, const std::string &label="")
	Create a vector member that is a non-const view of raw data.
RCP< VectorBase< Scalar > >	createMemberView (const VectorSpaceBase< Scalar > &vs, const RTOpPack::SubVectorView< Scalar > &raw_v, const std::string &label="")
	Calls createMemberView(Teuchos::rcp(&vs,false),raw_v).
RCP< const VectorBase< Scalar > >	createMemberView (const RCP< const VectorSpaceBase< Scalar > > &vs, const RTOpPack::ConstSubVectorView< Scalar > &raw_v, const std::string &label="")
	Create a vector member that is a const view of raw data.
RCP< const VectorBase< Scalar > >	createMemberView (const VectorSpaceBase< Scalar > &vs, const RTOpPack::ConstSubVectorView< Scalar > &raw_v, const std::string &label="")
	Calls createMemberView(Teuchos::rcp(&vs,false),raw_v).
RCP< MultiVectorBase< Scalar > >	createMembersView (const RCP< const VectorSpaceBase< Scalar > > &vs, const RTOpPack::SubMultiVectorView< Scalar > &raw_mv, const std::string &label="")
	Create a multi-vector member that is a non-const view of raw data.
RCP< MultiVectorBase< Scalar > >	createMembersView (const VectorSpaceBase< Scalar > &vs, const RTOpPack::SubMultiVectorView< Scalar > &raw_mv, const std::string &label="")
	Calls createMembersView(Teuchos::rcp(&vs,false),raw_mv).
RCP< const MultiVectorBase< Scalar > >	createMembersView (const RCP< const VectorSpaceBase< Scalar > > &vs, const RTOpPack::ConstSubMultiVectorView< Scalar > &raw_mv, const std::string &label="")
	Create a multi-vector member that is a const view of raw data.
RCP< const MultiVectorBase< Scalar > >	createMembersView (const VectorSpaceBase< Scalar > &vs, const RTOpPack::ConstSubMultiVectorView< Scalar > &raw_mv, const std::string &label="")
	Calls createMembersView(Teuchos::rcp(&vs,false),raw_mv).

---

Detailed Description

template<class Scalar>
class Thyra::VectorSpaceBase< Scalar >

Abstract interface for objects that represent a space for vectors.

This interface acts primarily as an "Abstract Factory" interface for creating VectorBase objects using the non-member function Thyra::createMember(). A VectorSpaceBase can also create MultiVectorBase objects which represent a compact collection of vectors. A secondary role for VectorSpaceBase objects is to test for compatibility of vector space objects using the isCompatible() method.

Clients can not directly create VectorBase and MultiVectorBase objects using the member functions createMember(), createMembers(), createMemberView(), and createMembersView() but instead must use the non-member Thyra_Op_Vec_createMember_grp.

Note that a VectorSpaceBase object must also be able to create MultiVectorBase objects with any number of column vectors using the Thyra::createMembers() function. Once created, the LinearOpBase::domain() function supported by a created MultiVectorBase object returns a vector space of dimension equal to the number of columns in the multi-vector. An interesting implication of this design is that the creation of a multi-vector provides a way for clients to create vector spaces of any arbitrary (although small usually) dimension. In order to give the client the same ability to create smaller vector spaces without having to create a dummy multi-vector object first, the method smallVecSpcFcty() is included. The method smallVecSpcFcty() returns a VectorSpaceFactoryBase object that can create (typically serial) vector spaces of any small dimension that are compatible with the domain spaces of MultiVectorBase objects created by the vector space object.

A vector space is also where the scalar product for the space is defined which is computed by the scalarProd() method. A scalar product allows the vector space to introduce scaling into many different types of numerical algorithms.

If the underlying object is not initialized, then dim()==0 will be true and none of the other methods should be called or exceptions will be thrown.

Notes for subclass developers

This is a fairly bare-bones interface class without much in the way of default function implementations. The subclass VectorSpaceDefaultBase provides a default multi-vector implementation and should the the first choice for subclass implementations.

---

Member Function Documentation

template<class Scalar> virtual Index Thyra::VectorSpaceBase< Scalar >::dim (   )  const [pure virtual]

Return the dimension of the vector space. If the underlying object is not initialized, then dim()==0 will be true and none of the other methods should be called.

template<class Scalar> virtual bool Thyra::VectorSpaceBase< Scalar >::isCompatible (  const VectorSpaceBase< Scalar > &  vecSpc  )  const [pure virtual]

Compare the compatibility of two vector spaces. If this function returns true, then vectors created from either of the vector spaces will be compatible and can be combined in vector operations. Preconditions: this->dim() > 0 Postconditions: [this->dim() != vecSpc.dim()] returnVal == false Invariants: [this->isCompatible(vecSpc) == true] vecSpc.isCompatible(*this) == true

template<class Scalar> virtual RCP< const VectorSpaceFactoryBase<Scalar> > Thyra::VectorSpaceBase< Scalar >::smallVecSpcFcty (   )  const [pure virtual]

Return a VectorSpaceFactoryBase object for the creation of (usually serial) vector spaces with a small dimension. Preconditions: this->dim() > 0

template<class Scalar> virtual Scalar Thyra::VectorSpaceBase< Scalar >::scalarProd (  const VectorBase< Scalar > &  x, const VectorBase< Scalar > &  y )  const [pure virtual]

Return the scalar product of two vectors in the vector space. Preconditions: x.space()->isCompatible(*this) (throw Exceptions::IncompatibleVectorSpaces) y.space()->isCompatible(*this) (throw Exceptions::IncompatibleVectorSpaces)

template<class Scalar> virtual void Thyra::VectorSpaceBase< Scalar >::scalarProds (  const MultiVectorBase< Scalar > &  X, const MultiVectorBase< Scalar > &  Y, Scalar  scalar_prods[] )  const [pure virtual]

Return the scalar product of each column in two multi-vectors in the vector space. Parameters: X  [in] Multi-vector. Y  [in] Multi-vector. scalar_prod  [out] Array (length X.domain()->dim()) containing the scalar products scalar_prod[j] = this->scalarProd(*X.col(j),*Y.col(j)), for j = 0 ... X.domain()->dim()-1. Preconditions: X.range()->isCompatible(*this) (throw Exceptions::IncompatibleVectorSpaces) Y.range()->isCompatible(*this) (throw Exceptions::IncompatibleVectorSpaces) X.domain()->isCompatible(*Y.domain()) (throw Exceptions::IncompatibleVectorSpaces)

template<class Scalar> bool Thyra::VectorSpaceBase< Scalar >::isEuclidean (   )  const [virtual]

Return if this vector space has a Euclidean (identity) basis in which case the scalar product is the same as the dot product. The default implementation returns false (even though on average the Euclidean scalar product is used).

template<class Scalar> bool Thyra::VectorSpaceBase< Scalar >::hasInCoreView (  const Range1D &  rng = Range1D(), const EViewType  viewType = VIEW_TYPE_DETACHED, const EStrideType  strideType = STRIDE_TYPE_NONUNIT )  const [virtual]

Returns true if this->acquireDetachedView(rng,...) returns a direct view of the range of data requested. Parameters: rng  [in] The range of elements for the view (see acquireDetachedView()). The default value is Range1D() (i.e. All of the elements in the vector). viewType  [in] The type of view allowed. strideType  [in] The type of stride the view is allowed to be. Preconditions: this->dim() > 0 full_range(rng,0,this->dim()-1).ubound < this->dim() There are three different questions about the behavior of the acquireDetachedView() that this query function can answer: The elements in rng are fairly cheaply accessble in local (i.e. in-core) memory if this->hasInCoreView(rng)==true. Note that this also allows for detached temporary copies of data. A direct view of the elements in rng is available in local (i.e. in-core) memory if this->hasInCoreView(rng,VIEW_TYPE_DIRECT)==true. No copy of data is allowed here. A direct view of the elements in rng with unit stride is available in local (i.e. in-core) memory if this->hasInCoreView(rng,VIEW_TYPE_DIRECT,STRIDE_TYPE_UNIT)==true No copy of data is allowed here. The default implementation returns false (i.e. by default we do not assume that any direct and/or contiguous views of any range of data are provided).

template<class Scalar> RCP< const VectorSpaceBase< Scalar > > Thyra::VectorSpaceBase< Scalar >::clone (   )  const [virtual]

Clone this object (if supported). It is allowed for returnVal.get()==NULL which means that this capability is optional. Preconditions: this->dim() > 0 The default implementation returns returnVal.get()==NULL.

template<class Scalar> virtual RCP< VectorBase<Scalar> > Thyra::VectorSpaceBase< Scalar >::createMember (   )  const [protected, pure virtual]

Create a vector member from the vector space. Preconditions: this->dim() > 0 Postconditions: returnVal.get() != NULL returnVal->space()->isCompatible(*this) == true Note: This function is not to be called directly since it is protected! See the Thyra_Op_Vec_createMember_grp. Returns: A smart reference counted pointer to a dynamically allocated vector object. After construction the values in the vector *returnVal are unspecified (uninitialized). This allows for faster execution times. Note that returnVal->space().get() == this need not be true.

template<class Scalar> virtual RCP< MultiVectorBase<Scalar> > Thyra::VectorSpaceBase< Scalar >::createMembers (  int  numMembers  )  const [protected, pure virtual]

Create a set of vector members (a MultiVectorBase) from the vector space. Preconditions: this->dim() > 0 num_vecs >= 1 Postconditions: returnVal->range()->isCompatible(*this) == true returnVal->domain()->dim() == numMembers Returns: A smart reference-counted pointer to a dynamically allocated multi-vector object. After construction, the values in *returnVal are unspecified (uninitialized). This allows for faster execution times. Note that returnVal->range().get()==this does not have to be true but will be in may cases.

template<class Scalar> virtual RCP<VectorBase<Scalar> > Thyra::VectorSpaceBase< Scalar >::createMemberView (  const RTOpPack::SubVectorView< Scalar > &  raw_v  )  const [protected, pure virtual]

Create a vector member that is a non-const view of raw vector data. Parameters: raw_v  [in] On input contains pointer (i.e. raw_v.values()) to array that the returned VectorBase will be a view of. The data pointed to by raw_v.values() must remain valid until the returned VectorBase object is destroyed. Preconditions: raw_v has been initialized to memory (i.e. raw_v.subDim()!=0 && raw_v.values()!=NULL). raw_v is consistent with the local storage of this vector spaces vector data. This precondition is purposefully vague since this function can be used an variety of specialized use-cases. Postconditions: See this->createMember() It is stated here that the client can not expect that the values pointed to by raw_v.values() to be changed until the smart pointer returnVal goes out of scope. This is to allow an implementation that temporarily copies data into and out of a VectorBase object using explicit vector access.

template<class Scalar> virtual RCP<const VectorBase<Scalar> > Thyra::VectorSpaceBase< Scalar >::createMemberView (  const RTOpPack::ConstSubVectorView< Scalar > &  raw_v  )  const [protected, pure virtual]

Create a vector member that is a const view of raw vector data. Parameters: raw_v  [in] On input contains pointer (i.e. raw_v.values()) to array that the returned VectorBase will be a view of. The data pointed to by raw_v.values() must remain valid until the returned VectorBase object is destroyed. This function works exactly the same as the previous version that takes a RTOpPack::SubVectorView object except that this version takes a RTOpPack::ConstSubVectorView and returns a smart pointer to a const VectorBase object. Preconditions: See the previous RTOpPack::SubVectorView version of this function. Postconditions: See this->createMember()

template<class Scalar> virtual RCP<MultiVectorBase<Scalar> > Thyra::VectorSpaceBase< Scalar >::createMembersView (  const RTOpPack::SubMultiVectorView< Scalar > &  raw_mv  )  const [protected, pure virtual]

Create a multi-vector member that is a non-const view of raw multi-vector data. Parameters: raw_mv  [in] On input contains pointer (i.e. raw_mv.values()) to array that the returned MultiVectorBase will be a view of. Preconditions: raw_mv has been initialized to memory (i.e. raw_mv.subDim()!=0 && raw_mv.values()!=NULL). raw_mv is consistent with the local storage of this spaces vector data. This precondition is purposefully vague since this function can be used an variety of specialized use-cases. Postconditions: See this->createMembers() where numMembers==raw_mv.numSubCols() It is stated here that the client can not expect that the values pointed to by raw_mv.values() to be changed until the smart pointer returnVal goes out of scope. This is to allow for an implementation that temporarily copies data into and out of a MultiVectorBase object using explicit vector access.

template<class Scalar> virtual RCP<const MultiVectorBase<Scalar> > Thyra::VectorSpaceBase< Scalar >::createMembersView (  const RTOpPack::ConstSubMultiVectorView< Scalar > &  raw_mv  )  const [protected, pure virtual]

Create a multi-vector member that is a const view of raw multi-vector data. Parameters: raw_mv  [in] On input contains pointer (i.e. raw_mv.values()) to array that the returned MultiVectorBase will be a view of. The data pointed to by raw_mv.values() must remain valid until the returned MultiVectorBase object is destroyed. This function works exactly the same as the previous version that takes a RTOpPack::SubMultiVectorView object except that this version takes a RTOpPack::ConstSubMultiVectorView object and returns a smart pointer to a const MultiVectorBase object. Preconditions: See the previous RTOpPack::SubMultiVectorView version of this function. Postconditions: See this->createMember()

---

Friends And Related Function Documentation

template<class Scalar> RCP< const VectorSpaceBase< Scalar > > makeHaveOwnership (  const RCP< const VectorSpaceBase< Scalar > > &  vs  )  [related]

Helper function that clones a VectorSpaceBase object if the RCP does not have ownership.

template<class Scalar> RCP< VectorBase< Scalar > > createMember (  const RCP< const VectorSpaceBase< Scalar > > &  vs, const std::string &  label = "" )  [related]

Create a vector member from the vector space. Calls VectorSpaceBase::createMember() on vs but the returned VectorBase object can live past vs.

template<class Scalar> RCP< VectorBase< Scalar > > createMember (  const VectorSpaceBase< Scalar > &  vs, const std::string &  label = "" )  [related]

Calls createMember(Teuchos::rcp(&vs,false)).

template<class Scalar> RCP< MultiVectorBase< Scalar > > createMembers (  const RCP< const VectorSpaceBase< Scalar > > &  vs, int  numMembers, const std::string &  label = "" )  [related]

Create a set of vector members (a MultiVectorBase) from the vector space. Calls VectorSpaceBase::createMembers() on vs but the returned MultiVectorBase object can live past vs.

template<class Scalar> RCP< MultiVectorBase< Scalar > > createMembers (  const VectorSpaceBase< Scalar > &  vs, int  numMembers, const std::string &  label = "" )  [related]

Calls createMembers(Teuchos::rcp(&vs,false),numMembers).

template<class Scalar> RCP< VectorBase< Scalar > > createMemberView (  const RCP< const VectorSpaceBase< Scalar > > &  vs, const RTOpPack::SubVectorView< Scalar > &  raw_v, const std::string &  label = "" )  [related]

Create a vector member that is a non-const view of raw data. Calls VectorSpaceBase::createMemberView() on vs but the returned VectorBase object can live past vs.

template<class Scalar> RCP< VectorBase< Scalar > > createMemberView (  const VectorSpaceBase< Scalar > &  vs, const RTOpPack::SubVectorView< Scalar > &  raw_v, const std::string &  label = "" )  [related]

Calls createMemberView(Teuchos::rcp(&vs,false),raw_v).

template<class Scalar> RCP< const VectorBase< Scalar > > createMemberView (  const RCP< const VectorSpaceBase< Scalar > > &  vs, const RTOpPack::ConstSubVectorView< Scalar > &  raw_v, const std::string &  label = "" )  [related]

Create a vector member that is a const view of raw data. Calls VectorSpaceBase::createMemberView() on vs but the returned VectorBase object can live past vs.

template<class Scalar> RCP< const VectorBase< Scalar > > createMemberView (  const VectorSpaceBase< Scalar > &  vs, const RTOpPack::ConstSubVectorView< Scalar > &  raw_v, const std::string &  label = "" )  [related]

Calls createMemberView(Teuchos::rcp(&vs,false),raw_v).

template<class Scalar> RCP< MultiVectorBase< Scalar > > createMembersView (  const RCP< const VectorSpaceBase< Scalar > > &  vs, const RTOpPack::SubMultiVectorView< Scalar > &  raw_mv, const std::string &  label = "" )  [related]

Create a multi-vector member that is a non-const view of raw data. Calls VectorSpaceBase::createMembersView() on vs but the returned MultiVectorBase object can live past vs.

template<class Scalar> RCP< MultiVectorBase< Scalar > > createMembersView (  const VectorSpaceBase< Scalar > &  vs, const RTOpPack::SubMultiVectorView< Scalar > &  raw_mv, const std::string &  label = "" )  [related]

Calls createMembersView(Teuchos::rcp(&vs,false),raw_mv).

template<class Scalar> RCP< const MultiVectorBase< Scalar > > createMembersView (  const RCP< const VectorSpaceBase< Scalar > > &  vs, const RTOpPack::ConstSubMultiVectorView< Scalar > &  raw_mv, const std::string &  label = "" )  [related]

Create a multi-vector member that is a const view of raw data. Calls VectorSpaceBase::createMembersView() on vs but the returned MultiVectorBase object can live past vs.

template<class Scalar> RCP< const MultiVectorBase< Scalar > > createMembersView (  const VectorSpaceBase< Scalar > &  vs, const RTOpPack::ConstSubMultiVectorView< Scalar > &  raw_mv, const std::string &  label = "" )  [related]

Calls createMembersView(Teuchos::rcp(&vs,false),raw_mv).

---

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

• Thyra_VectorSpaceBaseDecl.hpp
• Thyra_VectorSpaceBase.hpp

---