Stokhos Package Browser (Single Doxygen Collection) Version of the Day
Loading...
Searching...
No Matches
Stokhos_OrthogPolyExpansionBase.hpp
Go to the documentation of this file.
1// @HEADER
2// ***********************************************************************
3//
4// Stokhos Package
5// Copyright (2009) Sandia Corporation
6//
7// Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
8// license for use of this work by or on behalf of the U.S. Government.
9//
10// Redistribution and use in source and binary forms, with or without
11// modification, are permitted provided that the following conditions are
12// met:
13//
14// 1. Redistributions of source code must retain the above copyright
15// notice, this list of conditions and the following disclaimer.
16//
17// 2. Redistributions in binary form must reproduce the above copyright
18// notice, this list of conditions and the following disclaimer in the
19// documentation and/or other materials provided with the distribution.
20//
21// 3. Neither the name of the Corporation nor the names of the
22// contributors may be used to endorse or promote products derived from
23// this software without specific prior written permission.
24//
25// THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
26// EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
27// IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
28// PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
29// CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
30// EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
31// PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
32// PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
33// LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
34// NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
35// SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
36//
37// Questions? Contact Eric T. Phipps (etphipp@sandia.gov).
38//
39// ***********************************************************************
40// @HEADER
41
42#ifndef STOKHOS_ORTHOGPOLYEXPANSIONBASE_HPP
43#define STOKHOS_ORTHOGPOLYEXPANSIONBASE_HPP
44
47
48#include "Teuchos_RCP.hpp"
49#include "Teuchos_ParameterList.hpp"
50
51namespace Stokhos {
52
54
58 template <typename ordinal_type, typename value_type, typename node_type>
60 public OrthogPolyExpansion<ordinal_type, value_type, node_type> {
61 public:
62
65 const Teuchos::RCP<const OrthogPolyBasis<ordinal_type, value_type> >& basis,
67 const Teuchos::RCP<Teuchos::ParameterList>& params = Teuchos::null);
68
71
73 ordinal_type size() const { return sz; }
74
76 Teuchos::RCP< const OrthogPolyBasis<ordinal_type, value_type> >
77 getBasis() const {return basis; }
78
80 virtual Teuchos::RCP<const Sparse3Tensor<ordinal_type, value_type> >
81 getTripleProduct() const { return Cijk; }
82
83 // Operations
84 void unaryMinus(
85 OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
86 const OrthogPolyApprox<ordinal_type, value_type, node_type>& a);
87
88 void plusEqual(
89 OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
90 const value_type& x);
91 void minusEqual(
92 OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
93 const value_type& x);
94 void timesEqual(
95 OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
96 const value_type& x);
97 void divideEqual(
98 OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
99 const value_type& x);
100
101 void plusEqual(
102 OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
103 const OrthogPolyApprox<ordinal_type, value_type, node_type>& x);
104 void minusEqual(
105 OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
106 const OrthogPolyApprox<ordinal_type, value_type, node_type>& x);
107 void timesEqual(
108 OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
109 const OrthogPolyApprox<ordinal_type, value_type, node_type>& x);
110 void divideEqual(
111 OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
112 const OrthogPolyApprox<ordinal_type, value_type, node_type>& x);
113
114 void plus(OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
115 const OrthogPolyApprox<ordinal_type, value_type, node_type>& a,
116 const OrthogPolyApprox<ordinal_type, value_type, node_type>& b);
117 void plus(OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
118 const value_type& a,
119 const OrthogPolyApprox<ordinal_type, value_type, node_type>& b);
120 void plus(OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
121 const OrthogPolyApprox<ordinal_type, value_type, node_type>& a,
122 const value_type& b);
123 void minus(OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
124 const OrthogPolyApprox<ordinal_type, value_type, node_type>& a,
125 const OrthogPolyApprox<ordinal_type, value_type, node_type>& b);
126 void minus(OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
127 const value_type& a,
128 const OrthogPolyApprox<ordinal_type, value_type, node_type>& b);
129 void minus(OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
130 const OrthogPolyApprox<ordinal_type, value_type, node_type>& a,
131 const value_type& b);
132 void times(OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
133 const OrthogPolyApprox<ordinal_type, value_type, node_type>& a,
134 const OrthogPolyApprox<ordinal_type, value_type, node_type>& b);
135 void times(OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
136 const value_type& a,
137 const OrthogPolyApprox<ordinal_type, value_type, node_type>& b);
138 void times(OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
139 const OrthogPolyApprox<ordinal_type, value_type, node_type>& a,
140 const value_type& b);
141 void divide(OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
142 const OrthogPolyApprox<ordinal_type, value_type, node_type>& a,
143 const OrthogPolyApprox<ordinal_type, value_type, node_type>& b);
144 void divide(OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
145 const value_type& a,
146 const OrthogPolyApprox<ordinal_type, value_type, node_type>& b);
147 void divide(OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
148 const OrthogPolyApprox<ordinal_type, value_type, node_type>& a,
149 const value_type& b);
150
151
152 void abs(OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
153 const OrthogPolyApprox<ordinal_type, value_type, node_type>& a);
154 void fabs(OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
155 const OrthogPolyApprox<ordinal_type, value_type, node_type>& a);
156 void max(OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
157 const OrthogPolyApprox<ordinal_type, value_type, node_type>& a,
158 const OrthogPolyApprox<ordinal_type, value_type, node_type>& b);
159 void max(OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
160 const value_type& a,
161 const OrthogPolyApprox<ordinal_type, value_type, node_type>& b);
162 void max(OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
163 const OrthogPolyApprox<ordinal_type, value_type, node_type>& a,
164 const value_type& b);
165 void min(OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
166 const OrthogPolyApprox<ordinal_type, value_type, node_type>& a,
167 const OrthogPolyApprox<ordinal_type, value_type, node_type>& b);
168 void min(OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
169 const value_type& a,
170 const OrthogPolyApprox<ordinal_type, value_type, node_type>& b);
171 void min(OrthogPolyApprox<ordinal_type, value_type, node_type>& c,
172 const OrthogPolyApprox<ordinal_type, value_type, node_type>& a,
173 const value_type& b);
174
175 private:
176
177 // Prohibit copying
179
180 // Prohibit Assignment
182
183 protected:
184
186 Teuchos::RCP<const OrthogPolyBasis<ordinal_type, value_type> > basis;
187
190
192 Teuchos::RCP<const Stokhos::Sparse3Tensor<ordinal_type, value_type> > Cijk;
193
195 Teuchos::RCP<Teuchos::ParameterList> params;
196
198 Teuchos::RCP<Stokhos::DivisionExpansionStrategy<ordinal_type,value_type,node_type> > division_strategy;
199
201 ordinal_type sz;
202
203 }; // class OrthogPolyExpansionBase
204
205} // namespace Stokhos
206
208
209#endif // STOKHOS_ORTHOGPOLYEXPANSIONBASE_HPP
Base class for consolidating common expansion implementations.
Teuchos::RCP< const Stokhos::Sparse3Tensor< ordinal_type, value_type > > Cijk
Triple-product tensor.
void minusEqual(OrthogPolyApprox< ordinal_type, value_type, node_type > &c, const value_type &x)
void max(OrthogPolyApprox< ordinal_type, value_type, node_type > &c, const OrthogPolyApprox< ordinal_type, value_type, node_type > &a, const OrthogPolyApprox< ordinal_type, value_type, node_type > &b)
virtual Teuchos::RCP< const Sparse3Tensor< ordinal_type, value_type > > getTripleProduct() const
Get triple product.
OrthogPolyExpansionBase(const Teuchos::RCP< const OrthogPolyBasis< ordinal_type, value_type > > &basis, const Teuchos::RCP< const Stokhos::Sparse3Tensor< ordinal_type, value_type > > &Cijk, const Teuchos::RCP< Teuchos::ParameterList > &params=Teuchos::null)
Constructor.
void fabs(OrthogPolyApprox< ordinal_type, value_type, node_type > &c, const OrthogPolyApprox< ordinal_type, value_type, node_type > &a)
void plus(OrthogPolyApprox< ordinal_type, value_type, node_type > &c, const OrthogPolyApprox< ordinal_type, value_type, node_type > &a, const OrthogPolyApprox< ordinal_type, value_type, node_type > &b)
void timesEqual(OrthogPolyApprox< ordinal_type, value_type, node_type > &c, const value_type &x)
Teuchos::RCP< Teuchos::ParameterList > params
Parameter list.
void unaryMinus(OrthogPolyApprox< ordinal_type, value_type, node_type > &c, const OrthogPolyApprox< ordinal_type, value_type, node_type > &a)
void divideEqual(OrthogPolyApprox< ordinal_type, value_type, node_type > &c, const value_type &x)
OrthogPolyExpansionBase(const OrthogPolyExpansionBase &)
void divide(OrthogPolyApprox< ordinal_type, value_type, node_type > &c, const OrthogPolyApprox< ordinal_type, value_type, node_type > &a, const OrthogPolyApprox< ordinal_type, value_type, node_type > &b)
OrthogPolyExpansionBase & operator=(const OrthogPolyExpansionBase &b)
void minus(OrthogPolyApprox< ordinal_type, value_type, node_type > &c, const OrthogPolyApprox< ordinal_type, value_type, node_type > &a, const OrthogPolyApprox< ordinal_type, value_type, node_type > &b)
Teuchos::RCP< const OrthogPolyBasis< ordinal_type, value_type > > basis
Basis.
Teuchos::RCP< Stokhos::DivisionExpansionStrategy< ordinal_type, value_type, node_type > > division_strategy
Division expansion strategy.
void min(OrthogPolyApprox< ordinal_type, value_type, node_type > &c, const OrthogPolyApprox< ordinal_type, value_type, node_type > &a, const OrthogPolyApprox< ordinal_type, value_type, node_type > &b)
ordinal_type size() const
Get expansion size.
Stokhos::Sparse3Tensor< ordinal_type, value_type > Cijk_type
Short-hand for Cijk.
void times(OrthogPolyApprox< ordinal_type, value_type, node_type > &c, const OrthogPolyApprox< ordinal_type, value_type, node_type > &a, const OrthogPolyApprox< ordinal_type, value_type, node_type > &b)
void plusEqual(OrthogPolyApprox< ordinal_type, value_type, node_type > &c, const value_type &x)
Teuchos::RCP< const OrthogPolyBasis< ordinal_type, value_type > > getBasis() const
Get basis.
void abs(OrthogPolyApprox< ordinal_type, value_type, node_type > &c, const OrthogPolyApprox< ordinal_type, value_type, node_type > &a)
Abstract base class for orthogonal polynomial-based expansions.
Data structure storing a sparse 3-tensor C(i,j,k) in a a compressed format.
Top-level namespace for Stokhos classes and functions.