Thyra_DiagonalScalarProd_decl.hpp

// @HEADER
// *****************************************************************************
//    Thyra: Interfaces and Support for Abstract Numerical Algorithms
//
// Copyright 2004 NTESS and the Thyra contributors.
// SPDX-License-Identifier: BSD-3-Clause
// *****************************************************************************
// @HEADER
 
#ifndef THYRA_DIAGONAL_SCALAR_PROD_DECL_HPP
#define THYRA_DIAGONAL_SCALAR_PROD_DECL_HPP
 
 
#include "Thyra_ScalarProdBase.hpp"
 
 
namespace Thyra {
 

template<class Scalar>
class DiagonalScalarProd : public ScalarProdBase<Scalar> {
public:
  

  DiagonalScalarProd();

  void initialize( const RCP<const VectorBase<Scalar> > &s_diag );

 
protected:
  

  virtual bool isEuclideanImpl() const;
  
  virtual void scalarProdsImpl(
    const MultiVectorBase<Scalar>& X, const MultiVectorBase<Scalar>& Y,
    const ArrayView<Scalar> &scalarProds_out ) const;

  RCP<const LinearOpBase<Scalar> > getLinearOpImpl() const;

 
private:
 
  RCP<const VectorBase<Scalar> > s_diag_;
 
};
 

template<class Scalar>
RCP<DiagonalScalarProd<Scalar> >
diagonalScalarProd(const RCP<const VectorBase<Scalar> > &s_diag)
{
  const RCP<DiagonalScalarProd<Scalar> > scalarProd =
    Teuchos::rcp(new DiagonalScalarProd<Scalar>());
  scalarProd->initialize(s_diag);
  return scalarProd;
}
 
 
 
} // end namespace Thyra
 
 
#endif  // THYRA_DIAGONAL_SCALAR_PROD_DECL_HPP