42#ifndef THYRA_TPETRA_EUCLIDEAN_SCALAR_PROD_DEF_HPP
43#define THYRA_TPETRA_EUCLIDEAN_SCALAR_PROD_DEF_HPP
45#include "Thyra_TpetraEuclideanScalarProd_decl.hpp"
46#include "Thyra_TpetraMultiVector.hpp"
47#include "Thyra_TpetraVector.hpp"
53template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
60 typedef Tpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> TMV;
64 if (nonnull(X_tpetra) && nonnull(Y_tpetra)) {
71 X_tpetra->dot(*Y_tpetra, scalarProds_out);
78template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
83 using Teuchos::rcp_dynamic_cast;
89 return tmv->getConstTpetraMultiVector();
92 RCP<const TV> tv = rcp_dynamic_cast<const TV>(mv);
94 return tv->getConstTpetraVector();
virtual void scalarProdsImpl(const MultiVectorBase< Scalar > &X, const MultiVectorBase< Scalar > &Y, const ArrayView< Scalar > &scalarProds) const
Simply calls dots(X,Y,scalar_prods).
Interface for a collection of column vectors called a multi-vector.
Extends concrete implementation of a Euclidean scalar product for specifically Tpetra vectors/multive...
virtual void scalarProdsImpl(const MultiVectorBase< Scalar > &X, const MultiVectorBase< Scalar > &Y, const ArrayView< Scalar > &scalarProds) const
If X and Y are both Tpetra wrappers, computes the pair-wise scalar products directly with Tpetra call...
Concrete implementation of Thyra::MultiVector in terms of Tpetra::MultiVector.
Concrete Thyra::SpmdVectorBase using Tpetra::Vector.
T_To & dyn_cast(T_From &from)