42 #ifndef THYRA_DEFAULT_MULTI_VECTOR_PRODUCT_VECTOR_HPP 43 #define THYRA_DEFAULT_MULTI_VECTOR_PRODUCT_VECTOR_HPP 46 #include "Thyra_DefaultMultiVectorProductVector_decl.hpp" 47 #include "Thyra_DefaultMultiVectorProductVectorSpace.hpp" 48 #include "Thyra_AssertOp.hpp" 49 #include "Teuchos_Assert.hpp" 58 template <
class Scalar>
65 template <
class Scalar>
72 TEUCHOS_TEST_FOR_EXCEPT(is_null(productSpace_in));
73 TEUCHOS_TEST_FOR_EXCEPT(is_null(multiVec));
75 "DefaultMultiVectorProductVector<Scalar>::initialize(productSpace,multiVec)",
76 *multiVec->range(), *productSpace_in->getBlock(0)
78 TEUCHOS_ASSERT_EQUALITY( multiVec->domain()->dim(), productSpace_in->numBlocks());
81 numBlocks_ = productSpace_in->numBlocks();
83 productSpace_ = productSpace_in;
90 template <
class Scalar>
97 TEUCHOS_TEST_FOR_EXCEPT(is_null(productSpace_in));
98 TEUCHOS_TEST_FOR_EXCEPT(is_null(multiVec));
100 "DefaultMultiVectorProductVector<Scalar>::initialize(productSpace_in,multiVec)",
101 *multiVec->range(), *productSpace_in->getBlock(0)
103 TEUCHOS_ASSERT_EQUALITY( multiVec->domain()->dim(), productSpace_in->numBlocks() );
106 numBlocks_ = productSpace_in->numBlocks();
108 productSpace_ = productSpace_in;
110 multiVec_ = multiVec;
115 template <
class Scalar>
116 RCP<MultiVectorBase<Scalar> >
119 return multiVec_.getNonconstObj();
123 template <
class Scalar>
124 RCP<const MultiVectorBase<Scalar> >
127 return multiVec_.getConstObj();
131 template <
class Scalar>
135 productSpace_ = Teuchos::null;
136 multiVec_.uninitialize();
143 template<
class Scalar>
146 std::ostringstream oss;
148 << Teuchos::Describable::description()
150 <<
"dim="<<this->space()->dim()
151 <<
",numColumns = "<<numBlocks_
156 template<
class Scalar>
158 Teuchos::FancyOStream &out_arg,
159 const Teuchos::EVerbosityLevel verbLevel
162 using Teuchos::OSTab;
163 using Teuchos::describe;
164 RCP<FancyOStream> out = rcp(&out_arg,
false);
167 case Teuchos::VERB_DEFAULT:
168 case Teuchos::VERB_LOW:
169 *out << this->description() << std::endl;
171 case Teuchos::VERB_MEDIUM:
172 case Teuchos::VERB_HIGH:
173 case Teuchos::VERB_EXTREME:
176 << Teuchos::Describable::description() <<
"{" 177 <<
"dim=" << this->space()->dim()
180 *out <<
"multiVec = " << Teuchos::describe(*multiVec_.getConstObj(),verbLevel);
184 TEUCHOS_TEST_FOR_EXCEPT(
true);
192 template <
class Scalar>
193 RCP<VectorBase<Scalar> >
197 TEUCHOS_ASSERT_IN_RANGE_UPPER_EXCLUSIVE( k, 0, numBlocks_ );
199 return multiVec_.getNonconstObj()->col(k);
203 template <
class Scalar>
204 RCP<const VectorBase<Scalar> >
208 TEUCHOS_ASSERT_IN_RANGE_UPPER_EXCLUSIVE( k, 0, numBlocks_ );
210 return multiVec_.getConstObj()->col(k);
217 template <
class Scalar>
218 RCP<const ProductVectorSpaceBase<Scalar> >
221 return productSpace_;
225 template <
class Scalar>
229 TEUCHOS_ASSERT_IN_RANGE_UPPER_EXCLUSIVE( k, 0, numBlocks_ );
231 return multiVec_.isConst();
235 template <
class Scalar>
236 RCP<MultiVectorBase<Scalar> >
239 return getNonconstVectorBlock(k);
243 template <
class Scalar>
244 RCP<const MultiVectorBase<Scalar> >
247 return getVectorBlock(k);
254 template <
class Scalar>
255 RCP< const VectorSpaceBase<Scalar> >
258 return productSpace_;
268 template <
class Scalar>
273 const Ptr<RTOpPack::ReductTarget> &reduct_obj,
277 this->getDefaultProductVector()->applyOp(
278 op, vecs, targ_vecs, reduct_obj, global_offset );
282 template <
class Scalar>
287 this->getDefaultProductVector()->acquireDetachedView(rng_in,sub_vec);
291 template <
class Scalar>
296 this->getDefaultProductVector()->releaseDetachedView(sub_vec);
300 template <
class Scalar>
305 TEUCHOS_TEST_FOR_EXCEPT(
"ToDo: Implement DefaultMultiVectorProductVector<Scalar>::acquireNonconstDetachedVectorViewImpl(...)!");
309 template <
class Scalar>
314 TEUCHOS_TEST_FOR_EXCEPT(
"ToDo: Implement DefaultMultiVectorProductVector<Scalar>::commitNonconstDetachedVectorViewImpl(...)!");
318 template <
class Scalar>
323 TEUCHOS_TEST_FOR_EXCEPT(
"ToDo: Implement DefaultMultiVectorProductVector<Scalar>::setSubVector(...)!");
330 template <
class Scalar>
333 multiVec_.getNonconstObj()->assign(alpha);
340 template <
class Scalar>
341 RCP<const DefaultProductVector<Scalar> >
351 Array<RCP<const VectorBase<Scalar> > > vecArray;
352 for (
int k = 0; k < numBlocks_; ++k) {
353 vecArray.push_back(multiVec_.getConstObj()->col(k));
356 return Thyra::defaultProductVector<Scalar>(
357 productSpace_->getDefaultProductVectorSpace(),
367 #endif // THYRA_DEFAULT_MULTI_VECTOR_PRODUCT_VECTOR_HPP void releaseDetachedVectorViewImpl(RTOpPack::ConstSubVectorView< Scalar > *sub_vec) const
#define THYRA_ASSERT_VEC_SPACES(FUNC_NAME, VS1, VS2)
This is a very useful macro that should be used to validate that two vector spaces are compatible...
void assignImpl(Scalar alpha)
bool blockIsConst(const int k) const
void acquireDetachedVectorViewImpl(const Range1D &rng, RTOpPack::ConstSubVectorView< Scalar > *sub_vec) const
Concrete implementation of a product vector which is really composed out of the columns of a multi-ve...
RCP< const VectorSpaceBase< Scalar > > space() const
void commitNonconstDetachedVectorViewImpl(RTOpPack::SubVectorView< Scalar > *sub_vec)
RCP< const MultiVectorBase< Scalar > > getMultiVector() const
Teuchos::Ordinal Ordinal
Type for the dimension of a vector space. `*.
Interface for a collection of column vectors called a multi-vector.
void acquireNonconstDetachedVectorViewImpl(const Range1D &rng, RTOpPack::SubVectorView< Scalar > *sub_vec)
Abstract interface for finite-dimensional dense vectors.
void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel) const
Standard concrete implementation of a product vector space that creates product vectors fromed implic...
void initialize(const RCP< const DefaultMultiVectorProductVectorSpace< Scalar > > &productSpace, const RCP< MultiVectorBase< Scalar > > &multiVec)
Initialize with a non-const multi-vector.
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
DefaultMultiVectorProductVector()
Construct to uninitialized.
std::string description() const
RCP< MultiVectorBase< Scalar > > getNonconstMultiVector()
RCP< const ProductVectorSpaceBase< Scalar > > productSpace() const
RCP< MultiVectorBase< Scalar > > getNonconstMultiVectorBlock(const int k)
RCP< const MultiVectorBase< Scalar > > getMultiVectorBlock(const int k) const
RCP< VectorBase< Scalar > > getNonconstVectorBlock(const int k)
RCP< const VectorBase< Scalar > > getVectorBlock(const int k) const
void setSubVectorImpl(const RTOpPack::SparseSubVectorT< Scalar > &sub_vec)