Thyra Version of the Day
Loading...
Searching...
No Matches
Thyra_DiagonalScalarProd_decl.hpp
1/*
2// @HEADER
3// ***********************************************************************
4//
5// Thyra: Interfaces and Support for Abstract Numerical Algorithms
6// Copyright (2004) Sandia Corporation
7//
8// Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
9// license for use of this work by or on behalf of the U.S. Government.
10//
11// Redistribution and use in source and binary forms, with or without
12// modification, are permitted provided that the following conditions are
13// met:
14//
15// 1. Redistributions of source code must retain the above copyright
16// notice, this list of conditions and the following disclaimer.
17//
18// 2. Redistributions in binary form must reproduce the above copyright
19// notice, this list of conditions and the following disclaimer in the
20// documentation and/or other materials provided with the distribution.
21//
22// 3. Neither the name of the Corporation nor the names of the
23// contributors may be used to endorse or promote products derived from
24// this software without specific prior written permission.
25//
26// THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
27// EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
28// IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
29// PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
30// CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
31// EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
32// PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
33// PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
34// LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
35// NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
36// SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
37//
38// Questions? Contact Roscoe A. Bartlett (bartlettra@ornl.gov)
39//
40// ***********************************************************************
41// @HEADER
42*/
43
44
45#ifndef THYRA_DIAGONAL_SCALAR_PROD_DECL_HPP
46#define THYRA_DIAGONAL_SCALAR_PROD_DECL_HPP
47
48
49#include "Thyra_ScalarProdBase.hpp"
50
51
52namespace Thyra {
53
54
61template<class Scalar>
62class DiagonalScalarProd : public ScalarProdBase<Scalar> {
63public:
64
67
70
72 void initialize( const RCP<const VectorBase<Scalar> > &s_diag );
73
75
76protected:
77
80
82 virtual bool isEuclideanImpl() const;
83
85 virtual void scalarProdsImpl(
87 const ArrayView<Scalar> &scalarProds_out ) const;
88
91
93
94private:
95
97
98};
99
100
105template<class Scalar>
108{
109 const RCP<DiagonalScalarProd<Scalar> > scalarProd =
111 scalarProd->initialize(s_diag);
112 return scalarProd;
113}
114
115
116
117} // end namespace Thyra
118
119
120#endif // THYRA_DIAGONAL_SCALAR_PROD_DECL_HPP
Concrete implementation of a scalar product using a diagonal vector.
virtual void scalarProdsImpl(const MultiVectorBase< Scalar > &X, const MultiVectorBase< Scalar > &Y, const ArrayView< Scalar > &scalarProds_out) const
RCP< DiagonalScalarProd< Scalar > > diagonalScalarProd(const RCP< const VectorBase< Scalar > > &s_diag)
Nonmember constructor.
virtual bool isEuclideanImpl() const
Returns false.
void initialize(const RCP< const VectorBase< Scalar > > &s_diag)
RCP< const LinearOpBase< Scalar > > getLinearOpImpl() const
Interface for a collection of column vectors called a multi-vector.
Abstract interface for scalar products.
Abstract interface for finite-dimensional dense vectors.
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)