34 #ifndef ANASAZI_THYRA_ADAPTER_HPP 35 #define ANASAZI_THYRA_ADAPTER_HPP 41 #include <Thyra_DetachedMultiVectorView.hpp> 42 #include <Thyra_MultiVectorBase.hpp> 43 #include <Thyra_MultiVectorStdOps.hpp> 44 #include <Thyra_VectorStdOps.hpp> 46 #include <Teuchos_Ptr.hpp> 47 #include <Teuchos_ArrayRCP.hpp> 48 #include <Teuchos_ArrayView.hpp> 66 template<
class ScalarType>
70 typedef Thyra::MultiVectorBase<ScalarType> TMVB;
96 const int numvecs = mv.domain()->dim();
100 Thyra::assign (Teuchos::outArg (*cc), mv);
111 const int numvecs = index.size();
117 Thyra::assign (Teuchos::outArg (*cc), *view);
124 const int numVecs = index.
size();
130 Thyra::assign (Teuchos::outArg (*cc), *view);
141 int numvecs = index.size();
157 for (
int i=0; i<numvecs; i++) {
158 if (lb+i != index[i]) contig =
false;
163 const Thyra::Range1D rng(lb,lb+numvecs-1);
165 cc = mv.subView(rng);
169 cc = mv.subView( Teuchos::arrayViewFromVector( index ) );
182 return mv.subView (index);
192 int numvecs = index.size();
208 for (
int i=0; i<numvecs; i++) {
209 if (lb+i != index[i]) contig =
false;
214 const Thyra::Range1D rng(lb,lb+numvecs-1);
216 cc = mv.subView(rng);
220 cc = mv.subView(Teuchos::arrayViewFromVector( index ) );
233 return mv.subView (index);
243 {
return mv.range()->dim(); }
247 {
return mv.domain()->dim(); }
258 const ScalarType beta, TMVB & mv )
264 B_thyra = Thyra::createMembersView(
266 RTOpPack::ConstSubMultiVectorView<ScalarType>(
268 Teuchos::arcpFromArrayView(Teuchos::arrayView(&B(0,0), B.
stride()*B.
numCols())), B.
stride()
272 A.apply(Thyra::NOTRANS,*B_thyra,Teuchos::outArg (mv),alpha,beta);
277 static void MvAddMv(
const ScalarType alpha,
const TMVB & A,
278 const ScalarType beta,
const TMVB & B, TMVB & mv )
280 using Teuchos::tuple;
using Teuchos::ptrInArg;
using Teuchos::inoutArg;
282 Thyra::linear_combination<ScalarType>(
283 tuple(alpha, beta)(), tuple(ptrInArg(A), ptrInArg(B))(), ST::zero(), inoutArg(mv));
288 static void MvTransMv(
const ScalarType alpha,
const TMVB & A,
const TMVB & mv,
292 int m = A.domain()->dim();
293 int n = mv.domain()->dim();
296 B_thyra = Thyra::createMembersView(
298 RTOpPack::SubMultiVectorView<ScalarType>(
300 Teuchos::arcpFromArrayView(Teuchos::arrayView(&B(0,0), B.
stride()*B.
numCols())), B.
stride()
309 static void MvDot(
const TMVB & mv,
const TMVB & A, std::vector<ScalarType> &b )
310 { Thyra::dots(mv,A,Teuchos::arrayViewFromVector( b )); }
316 const ScalarType alpha)
318 Thyra::scale (alpha, Teuchos::inOutArg (mv));
325 const std::vector<ScalarType>& alpha)
327 for (
unsigned int i=0; i<alpha.size(); i++) {
328 Thyra::scale (alpha[i], mv.col(i).ptr());
341 { Thyra::norms_2(mv,Teuchos::arrayViewFromVector( normvec )); }
350 static void SetBlock(
const TMVB & A,
const std::vector<int>& index, TMVB & mv )
353 int numvecs = index.size();
354 std::vector<int> indexA(numvecs);
355 int numAcols = A.domain()->dim();
356 for (
int i=0; i<numvecs; i++) {
361 if ( numAcols < numvecs ) {
366 else if ( numAcols > numvecs ) {
368 indexA.resize( numAcols );
375 Thyra::assign (Teuchos::outArg (*reldest), *relsource);
381 const int numColsA = A.domain()->dim();
382 const int numColsMv = mv.domain()->dim();
384 const bool validIndex = index.
lbound() >= 0 && index.
ubound() < numColsMv;
386 const bool validSource = index.
size() <= numColsA;
388 if (! validIndex || ! validSource)
390 std::ostringstream os;
391 os <<
"Anasazi::MultiVecTraits<Scalar, Thyra::MultiVectorBase<Scalar> " 392 ">::SetBlock(A, [" << index.
lbound() <<
", " << index.
ubound()
395 os.str() <<
"Range lower bound must be nonnegative.");
397 os.str() <<
"Range upper bound must be less than " 398 "the number of columns " << numColsA <<
" in the " 399 "'mv' output argument.");
401 os.str() <<
"Range must have no more elements than" 402 " the number of columns " << numColsA <<
" in the " 403 "'A' input argument.");
411 if (index.
lbound() == 0 && index.
ubound()+1 == numColsMv)
412 mv_view = Teuchos::rcpFromRef (mv);
414 mv_view = mv.subView (index);
420 if (index.
size() == numColsA)
421 A_view = Teuchos::rcpFromRef (A);
426 Thyra::assign (Teuchos::outArg (*mv_view), *A_view);
430 Assign (
const TMVB& A, TMVB& mv)
436 const int numColsA = A.domain()->dim();
437 const int numColsMv = mv.domain()->dim();
438 if (numColsA > numColsMv) {
439 const std::string prefix (
"Anasazi::MultiVecTraits<Scalar, " 440 "Thyra::MultiVectorBase<Scalar>" 441 " >::Assign(A, mv): ");
443 prefix <<
"Input multivector 'A' has " 444 << numColsA <<
" columns, but output multivector " 445 "'mv' has only " << numColsMv <<
" columns.");
449 if (numColsA == numColsMv) {
450 Thyra::assign (outArg (mv), A);
454 Thyra::assign (outArg (*mv_view), A);
466 Teuchos::outArg (mv));
475 Thyra::assign (Teuchos::outArg (mv), alpha);
485 static void MvPrint(
const TMVB & mv, std::ostream& os )
488 out->setf(std::ios_base::scientific);
513 template <
class ScalarType>
514 class OperatorTraits < ScalarType, Thyra::MultiVectorBase<ScalarType>, Thyra::LinearOpBase<ScalarType> >
521 static void Apply (
const Thyra::LinearOpBase< ScalarType >& Op,
const Thyra::MultiVectorBase< ScalarType > & x, Thyra::MultiVectorBase< ScalarType > & y )
static void MvNorm(const TMVB &mv, std::vector< typename Teuchos::ScalarTraits< ScalarType >::magnitudeType > &normvec)
Compute the 2-norm of each individual vector of mv. Upon return, normvec[i] holds the value of ...
static Teuchos::RCP< TMVB > CloneCopy(const TMVB &mv, const std::vector< int > &index)
Creates a new MultiVectorBase and copies the selected contents of mv into the new vector (deep copy)...
static void MvDot(const TMVB &mv, const TMVB &A, std::vector< ScalarType > &b)
Compute a vector b where the components are the individual dot-products of the i-th columns of A and ...
Declaration of basic traits for the multivector type.
#define TEUCHOS_TEST_FOR_EXCEPTION(throw_exception_test, Exception, msg)
Virtual base class which defines basic traits for the operator type.
static void MvTimesMatAddMv(const ScalarType alpha, const TMVB &A, const Teuchos::SerialDenseMatrix< int, ScalarType > &B, const ScalarType beta, TMVB &mv)
Update mv with .
static void Assign(const MV &A, MV &mv)
mv := A
Namespace Anasazi contains the classes, structs, enums and utilities used by the Anasazi package...
static void MvPrint(const TMVB &mv, std::ostream &os)
Print the mv multi-vector to the os output stream.
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
OrdinalType numRows() const
static void MvTransMv(const ScalarType alpha, const TMVB &A, const TMVB &mv, Teuchos::SerialDenseMatrix< int, ScalarType > &B)
Compute a dense matrix B through the matrix-matrix multiply .
Traits class which defines basic operations on multivectors.
static Teuchos::RCP< MV > CloneCopy(const MV &mv)
Creates a new MV and copies contents of mv into the new vector (deep copy).
Virtual base class which defines basic traits for the operator type.
static void MvAddMv(const ScalarType alpha, const TMVB &A, const ScalarType beta, const TMVB &B, TMVB &mv)
Replace mv with .
Anasazi header file which uses auto-configuration information to include necessary C++ headers...
static void MvRandom(TMVB &mv)
Replace the vectors in mv with random vectors.
static Teuchos::RCP< MV > CloneViewNonConst(MV &mv, const std::vector< int > &index)
Creates a new MV that shares the selected contents of mv (shallow copy).
static Teuchos::RCP< const MV > CloneView(const MV &mv, const std::vector< int > &index)
Creates a new const MV that shares the selected contents of mv (shallow copy).
static void SetBlock(const MV &A, const std::vector< int > &index, MV &mv)
Copy the vectors in A to a set of vectors in mv indicated by the indices given in index...
static int GetNumberVecs(const TMVB &mv)
Obtain the number of vectors in mv.
static void SetBlock(const TMVB &A, const std::vector< int > &index, TMVB &mv)
Copy the vectors in A to a set of vectors in mv indicated by the indices given in index...
static Teuchos::RCP< TMVB > CloneViewNonConst(TMVB &mv, const std::vector< int > &index)
Creates a new MultiVectorBase that shares the selected contents of mv (shallow copy).
static void MvInit(TMVB &mv, ScalarType alpha=Teuchos::ScalarTraits< ScalarType >::zero())
Replace each element of the vectors in mv with alpha.
static Teuchos::RCP< TMVB > CloneCopy(const TMVB &mv)
Creates a new MultiVectorBase and copies contents of mv into the new vector (deep copy)...
OrdinalType stride() const
static void MvScale(TMVB &mv, const std::vector< ScalarType > &alpha)
Scale each element of the i-th vector in *this with alpha[i].
static ptrdiff_t GetGlobalLength(const TMVB &mv)
Obtain the vector length of mv.
static void Apply(const Thyra::LinearOpBase< ScalarType > &Op, const Thyra::MultiVectorBase< ScalarType > &x, Thyra::MultiVectorBase< ScalarType > &y)
This method takes the MultiVectorBase x and applies the LinearOpBase Op to it resulting in the MultiV...
OrdinalType numCols() const
static Teuchos::RCP< const TMVB > CloneView(const TMVB &mv, const std::vector< int > &index)
Creates a new const MultiVectorBase that shares the selected contents of mv (shallow copy)...
static void MvScale(TMVB &mv, const ScalarType alpha)
Scale each element of the vectors in *this with alpha.
static Teuchos::RCP< TMVB > Clone(const TMVB &mv, const int numvecs)
Creates a new empty MultiVectorBase containing numvecs columns.