Stokhos Package Browser (Single Doxygen Collection)  Version of the Day
Stokhos_RecurrenceBasisImp.hpp
Go to the documentation of this file.
1 // $Id$
2 // $Source$
3 // @HEADER
4 // ***********************************************************************
5 //
6 // Stokhos Package
7 // Copyright (2009) Sandia Corporation
8 //
9 // Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
10 // license for use of this work by or on behalf of the U.S. Government.
11 //
12 // Redistribution and use in source and binary forms, with or without
13 // modification, are permitted provided that the following conditions are
14 // met:
15 //
16 // 1. Redistributions of source code must retain the above copyright
17 // notice, this list of conditions and the following disclaimer.
18 //
19 // 2. Redistributions in binary form must reproduce the above copyright
20 // notice, this list of conditions and the following disclaimer in the
21 // documentation and/or other materials provided with the distribution.
22 //
23 // 3. Neither the name of the Corporation nor the names of the
24 // contributors may be used to endorse or promote products derived from
25 // this software without specific prior written permission.
26 //
27 // THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
28 // EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
29 // IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
30 // PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
31 // CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
32 // EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
33 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
34 // PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
35 // LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
36 // NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
37 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
38 //
39 // Questions? Contact Eric T. Phipps (etphipp@sandia.gov).
40 //
41 // ***********************************************************************
42 // @HEADER
43 
44 #include "Teuchos_LAPACK.hpp"
45 #include "Teuchos_SerialDenseMatrix.hpp"
46 
47 template <typename ordinal_type, typename value_type>
48 Stokhos::RecurrenceBasis<ordinal_type, value_type>::
49 RecurrenceBasis(const std::string& name_, ordinal_type p_, bool normalize_,
50  Stokhos::GrowthPolicy growth_) :
51  name(name_),
52  p(p_),
53  normalize(normalize_),
54  growth(growth_),
55  quad_zero_tol(1.0e-14),
56 #ifdef HAVE_STOKHOS_DAKOTA
57  sparse_grid_growth_rule(webbur::level_to_order_linear_nn),
58 #else
59  sparse_grid_growth_rule(NULL),
60 #endif
61  alpha(p+1, value_type(0.0)),
62  beta(p+1, value_type(0.0)),
63  delta(p+1, value_type(0.0)),
64  gamma(p+1, value_type(0.0)),
65  norms(p+1, value_type(0.0))
66 {
67 }
68 
69 template <typename ordinal_type, typename value_type>
70 Stokhos::RecurrenceBasis<ordinal_type, value_type>::
71 RecurrenceBasis(ordinal_type p_, const RecurrenceBasis& basis) :
72  name(basis.name),
73  p(p_),
74  normalize(basis.normalize),
75  growth(basis.growth),
76  quad_zero_tol(basis.quad_zero_tol),
77  sparse_grid_growth_rule(basis.sparse_grid_growth_rule),
78  alpha(p+1, value_type(0.0)),
79  beta(p+1, value_type(0.0)),
80  delta(p+1, value_type(0.0)),
81  gamma(p+1, value_type(0.0)),
82  norms(p+1, value_type(0.0))
83 {
84 }
85 
86 template <typename ordinal_type, typename value_type>
87 void
88 Stokhos::RecurrenceBasis<ordinal_type, value_type>::
89 setup()
90 {
91  bool is_normalized =
92  computeRecurrenceCoefficients(p+1, alpha, beta, delta, gamma);
93  if (normalize && !is_normalized) {
94  normalizeRecurrenceCoefficients(alpha, beta, delta, gamma);
95  }
96 
97 
98  // Compute norms -- when polynomials are normalized, this should work
99  // out to one (norms[0] == 1, delta[k] == 1, beta[k] == gamma[k]
100  norms[0] = beta[0]/(gamma[0]*gamma[0]);
101  for (ordinal_type k=1; k<=p; k++) {
102  norms[k] = (beta[k]/gamma[k])*(delta[k-1]/delta[k])*norms[k-1];
103  }
104 }
105 
106 template <typename ordinal_type, typename value_type>
107 Stokhos::RecurrenceBasis<ordinal_type, value_type>::
108 ~RecurrenceBasis()
109 {
110 }
111 
112 template <typename ordinal_type, typename value_type>
113 ordinal_type
114 Stokhos::RecurrenceBasis<ordinal_type, value_type>::
115 order() const
116 {
117  return p;
118 }
119 
120 template <typename ordinal_type, typename value_type>
121 ordinal_type
122 Stokhos::RecurrenceBasis<ordinal_type, value_type>::
123 size() const
124 {
125  return p+1;
126 }
127 
128 template <typename ordinal_type, typename value_type>
129 const Teuchos::Array<value_type>&
130 Stokhos::RecurrenceBasis<ordinal_type, value_type>::
131 norm_squared() const
132 {
133  return norms;
134 }
135 
136 template <typename ordinal_type, typename value_type>
137 const value_type&
138 Stokhos::RecurrenceBasis<ordinal_type, value_type>::
139 norm_squared(ordinal_type i) const
140 {
141  return norms[i];
142 }
143 
144 template <typename ordinal_type, typename value_type>
145 Teuchos::RCP< Stokhos::Dense3Tensor<ordinal_type, value_type> >
146 Stokhos::RecurrenceBasis<ordinal_type, value_type>::
147 computeTripleProductTensor() const
148 {
149  // Compute Cijk = < \Psi_i \Psi_j \Psi_k >
150  ordinal_type sz = size();
151  Teuchos::RCP< Stokhos::Dense3Tensor<ordinal_type, value_type> > Cijk =
152  Teuchos::rcp(new Dense3Tensor<ordinal_type, value_type>(sz));
153  Teuchos::Array<value_type> points, weights;
154  Teuchos::Array< Teuchos::Array<value_type> > values;
155  getQuadPoints(3*p, points, weights, values);
156 
157  for (ordinal_type i=0; i<sz; i++) {
158  for (ordinal_type j=0; j<sz; j++) {
159  for (ordinal_type k=0; k<sz; k++) {
160  value_type triple_product = 0;
161  for (ordinal_type l=0; l<static_cast<ordinal_type>(points.size());
162  l++){
163  triple_product +=
164  weights[l]*(values[l][i])*(values[l][j])*(values[l][k]);
165  }
166  (*Cijk)(i,j,k) = triple_product;
167  }
168  }
169  }
170 
171  return Cijk;
172 }
173 
174 template <typename ordinal_type, typename value_type>
175 Teuchos::RCP< Stokhos::Sparse3Tensor<ordinal_type, value_type> >
176 Stokhos::RecurrenceBasis<ordinal_type, value_type>::
177 computeSparseTripleProductTensor(ordinal_type theOrder) const
178 {
179  // Compute Cijk = < \Psi_i \Psi_j \Psi_k >
180  value_type sparse_tol = 1.0e-15;
181  ordinal_type sz = size();
182  Teuchos::RCP< Stokhos::Sparse3Tensor<ordinal_type, value_type> > Cijk =
183  Teuchos::rcp(new Sparse3Tensor<ordinal_type, value_type>());
184  Teuchos::Array<value_type> points, weights;
185  Teuchos::Array< Teuchos::Array<value_type> > values;
186  getQuadPoints(3*p, points, weights, values);
187 
188  for (ordinal_type i=0; i<sz; i++) {
189  for (ordinal_type j=0; j<sz; j++) {
190  for (ordinal_type k=0; k<theOrder; k++) {
191  value_type triple_product = 0;
192  for (ordinal_type l=0; l<static_cast<ordinal_type>(points.size());
193  l++){
194  triple_product +=
195  weights[l]*(values[l][i])*(values[l][j])*(values[l][k]);
196  }
197  if (std::abs(triple_product/norms[i]) > sparse_tol)
198  Cijk->add_term(i,j,k,triple_product);
199  }
200  }
201  }
202  Cijk->fillComplete();
203 
204  return Cijk;
205 }
206 
207 template <typename ordinal_type, typename value_type>
208 Teuchos::RCP< Teuchos::SerialDenseMatrix<ordinal_type, value_type> >
209 Stokhos::RecurrenceBasis<ordinal_type, value_type>::
210 computeDerivDoubleProductTensor() const
211 {
212  // Compute Bij = < \Psi_i' \Psi_j >
213  Teuchos::Array<value_type> points, weights;
214  Teuchos::Array< Teuchos::Array<value_type> > values, derivs;
215  getQuadPoints(2*p, points, weights, values);
216  ordinal_type nqp = weights.size();
217  derivs.resize(nqp);
218  ordinal_type sz = size();
219  for (ordinal_type i=0; i<nqp; i++) {
220  derivs[i].resize(sz);
221  evaluateBasesAndDerivatives(points[i], values[i], derivs[i]);
222  }
223  Teuchos::RCP< Teuchos::SerialDenseMatrix<ordinal_type, value_type> > Bij =
224  Teuchos::rcp(new Teuchos::SerialDenseMatrix<ordinal_type, value_type>(sz,sz));
225  for (ordinal_type i=0; i<sz; i++) {
226  for (ordinal_type j=0; j<sz; j++) {
227  value_type b = value_type(0.0);
228  for (int qp=0; qp<nqp; qp++)
229  b += weights[qp]*derivs[qp][i]*values[qp][j];
230  (*Bij)(i,j) = b;
231  }
232  }
233 
234  return Bij;
235 }
236 
237 template <typename ordinal_type, typename value_type>
238 void
239 Stokhos::RecurrenceBasis<ordinal_type,value_type>::
240 evaluateBases(const value_type& x, Teuchos::Array<value_type>& basis_pts) const
241 {
242  // Evaluate basis polynomials P(x) using 3 term recurrence
243  // P_0(x) = 1/gamma[0], P_-1 = 0
244  // P_i(x) = [(delta[i-1]*x-alpha[i-1])*P_{i-1}(x) -
245  // beta[i-1]*P_{i-2}(x)]/gamma[i],
246  // i=2,3,...
247  basis_pts[0] = value_type(1)/gamma[0];
248  if (p >= 1)
249  basis_pts[1] = (delta[0]*x-alpha[0])*basis_pts[0]/gamma[1];
250  for (ordinal_type i=2; i<=p; i++)
251  basis_pts[i] = ((delta[i-1]*x-alpha[i-1])*basis_pts[i-1] -
252  beta[i-1]*basis_pts[i-2])/gamma[i];
253 }
254 
255 template <typename ordinal_type, typename value_type>
256 void
257 Stokhos::RecurrenceBasis<ordinal_type,value_type>::
258 evaluateBasesAndDerivatives(const value_type& x,
259  Teuchos::Array<value_type>& vals,
260  Teuchos::Array<value_type>& derivs) const
261 {
262  evaluateBases(x, vals);
263  derivs[0] = 0.0;
264  if (p >= 1)
265  derivs[1] = delta[0]/(gamma[0]*gamma[1]);
266  for (ordinal_type i=2; i<=p; i++)
267  derivs[i] = (delta[i-1]*vals[i-1] + (delta[i-1]*x-alpha[i-1])*derivs[i-1] -
268  beta[i-1]*derivs[i-2])/gamma[i];
269 }
270 
271 template <typename ordinal_type, typename value_type>
272 value_type
273 Stokhos::RecurrenceBasis<ordinal_type,value_type>::
274 evaluate(const value_type& x, ordinal_type k) const
275 {
276  // Evaluate basis polynomials P(x) using 3 term recurrence
277  // P_0(x) = 1/gamma[0], P_-1 = 0
278  // P_i(x) = [(delta[i-1]*x-alpha[i-1])*P_{i-1}(x) -
279  // beta[i-1]*P_{i-2}(x)]/gamma[i],
280  // i=2,3,...
281 
282  value_type v0 = value_type(1.0)/gamma[0];
283  if (k == 0)
284  return v0;
285 
286  value_type v1 = (x*delta[0]-alpha[0])*v0/gamma[1];
287  if (k == 1)
288  return v1;
289 
290  value_type v2 = value_type(0.0);
291  for (ordinal_type i=2; i<=k; i++) {
292  v2 = ((delta[i-1]*x-alpha[i-1])*v1 - beta[i-1]*v0)/gamma[i];
293  v0 = v1;
294  v1 = v2;
295  }
296 
297  return v2;
298 }
299 
300 template <typename ordinal_type, typename value_type>
301 void
302 Stokhos::RecurrenceBasis<ordinal_type, value_type>::
303 print(std::ostream& os) const
304 {
305  os << name << " basis of order " << p << "." << std::endl;
306 
307  os << "Alpha recurrence coefficients:\n\t";
308  for (ordinal_type i=0; i<=p; i++)
309  os << alpha[i] << " ";
310  os << std::endl;
311 
312  os << "Beta recurrence coefficients:\n\t";
313  for (ordinal_type i=0; i<=p; i++)
314  os << beta[i] << " ";
315  os << std::endl;
316 
317  os << "Delta recurrence coefficients:\n\t";
318  for (ordinal_type i=0; i<=p; i++)
319  os << delta[i] << " ";
320  os << std::endl;
321 
322  os << "Gamma recurrence coefficients:\n\t";
323  for (ordinal_type i=0; i<=p; i++)
324  os << gamma[i] << " ";
325  os << std::endl;
326 
327  os << "Basis polynomial norms (squared):\n\t";
328  for (ordinal_type i=0; i<=p; i++)
329  os << norms[i] << " ";
330  os << std::endl;
331 }
332 
333 template <typename ordinal_type, typename value_type>
334 const std::string&
335 Stokhos::RecurrenceBasis<ordinal_type, value_type>::
336 getName() const
337 {
338  return name;
339 }
340 
341 template <typename ordinal_type, typename value_type>
342 void
343 Stokhos::RecurrenceBasis<ordinal_type,value_type>::
344 getQuadPoints(ordinal_type quad_order,
345  Teuchos::Array<value_type>& quad_points,
346  Teuchos::Array<value_type>& quad_weights,
347  Teuchos::Array< Teuchos::Array<value_type> >& quad_values) const
348 {
349 
350  //This is a transposition into C++ of Gautschi's code for taking the first
351  // N recurrance coefficient and generating a N point quadrature rule.
352  // The MATLAB version is available at
353  // http://www.cs.purdue.edu/archives/2002/wxg/codes/gauss.m
354  ordinal_type num_points =
355  static_cast<ordinal_type>(std::ceil((quad_order+1)/2.0));
356  Teuchos::Array<value_type> a(num_points,0);
357  Teuchos::Array<value_type> b(num_points,0);
358  Teuchos::Array<value_type> c(num_points,0);
359  Teuchos::Array<value_type> d(num_points,0);
360 
361  // If we don't have enough recurrance coefficients, get some more.
362  if(num_points > p+1){
363  bool is_normalized = computeRecurrenceCoefficients(num_points, a, b, c, d);
364  if (!is_normalized)
365  normalizeRecurrenceCoefficients(a, b, c, d);
366  }
367  else { //else just take the ones we already have.
368  for(ordinal_type n = 0; n<num_points; n++){
369  a[n] = alpha[n];
370  b[n] = beta[n];
371  c[n] = delta[n];
372  d[n] = gamma[n];
373  }
374  if (!normalize)
375  normalizeRecurrenceCoefficients(a, b, c, d);
376  }
377 
378  // With normalized coefficients, A is symmetric and tri-diagonal, with
379  // diagonal = a, and off-diagonal = b, starting at b[1]
380  Teuchos::SerialDenseMatrix<ordinal_type,value_type> eig_vectors(num_points,
381  num_points);
382  Teuchos::Array<value_type> workspace(2*num_points);
383  ordinal_type info_flag;
384  Teuchos::LAPACK<ordinal_type,value_type> my_lapack;
385 
386  // compute the eigenvalues (stored in a) and right eigenvectors.
387  if (num_points == 1)
388  my_lapack.STEQR('I', num_points, &a[0], &b[0], eig_vectors.values(),
389  num_points, &workspace[0], &info_flag);
390  else
391  my_lapack.STEQR('I', num_points, &a[0], &b[1], eig_vectors.values(),
392  num_points, &workspace[0], &info_flag);
393 
394  // eigenvalues are sorted by STEQR
395  quad_points.resize(num_points);
396  quad_weights.resize(num_points);
397  for (ordinal_type i=0; i<num_points; i++) {
398  quad_points[i] = a[i];
399  if (std::abs(quad_points[i]) < quad_zero_tol)
400  quad_points[i] = 0.0;
401  quad_weights[i] = beta[0]*eig_vectors[i][0]*eig_vectors[i][0];
402  }
403 
404  // Evalute basis at gauss points
405  quad_values.resize(num_points);
406  for (ordinal_type i=0; i<num_points; i++) {
407  quad_values[i].resize(p+1);
408  evaluateBases(quad_points[i], quad_values[i]);
409  }
410 }
411 
412 template <typename ordinal_type, typename value_type>
413 ordinal_type
414 Stokhos::RecurrenceBasis<ordinal_type,value_type>::
415 quadDegreeOfExactness(ordinal_type n) const
416 {
417  return ordinal_type(2)*n-ordinal_type(1);
418 }
419 
420 template <typename ordinal_type, typename value_type>
421 ordinal_type
422 Stokhos::RecurrenceBasis<ordinal_type,value_type>::
423 coefficientGrowth(ordinal_type n) const
424 {
425  if (growth == SLOW_GROWTH)
426  return n;
427 
428  // else moderate
429  return ordinal_type(2)*n;
430 }
431 
432 template <typename ordinal_type, typename value_type>
433 ordinal_type
434 Stokhos::RecurrenceBasis<ordinal_type,value_type>::
435 pointGrowth(ordinal_type n) const
436 {
437  if (growth == SLOW_GROWTH) {
438  if (n % ordinal_type(2) == ordinal_type(1))
439  return n + ordinal_type(1);
440  else
441  return n;
442  }
443 
444  // else moderate
445  return n;
446 }
447 
448 template <typename ordinal_type, typename value_type>
449 void
450 Stokhos::RecurrenceBasis<ordinal_type,value_type>::
451 getRecurrenceCoefficients(Teuchos::Array<value_type>& a,
452  Teuchos::Array<value_type>& b,
453  Teuchos::Array<value_type>& c,
454  Teuchos::Array<value_type>& g) const
455 {
456  a = alpha;
457  b = beta;
458  c = delta;
459  g = gamma;
460 }
461 
462 template <typename ordinal_type, typename value_type>
463 void
464 Stokhos::RecurrenceBasis<ordinal_type,value_type>::
465 normalizeRecurrenceCoefficients(
466  Teuchos::Array<value_type>& a,
467  Teuchos::Array<value_type>& b,
468  Teuchos::Array<value_type>& c,
469  Teuchos::Array<value_type>& g) const
470 {
471  ordinal_type n = a.size();
472  a[0] = a[0] / c[0];
473  b[0] = std::sqrt(b[0]);
474  g[0] = b[0];
475  for (ordinal_type k=1; k<n; k++) {
476  a[k] = a[k] / c[k];
477  b[k] = std::sqrt((b[k]*g[k])/(c[k]*c[k-1]));
478  g[k] = b[k];
479  }
480  for (ordinal_type k=0; k<n; k++)
481  c[k] = value_type(1);
482 }
g
ScalarType g(const Teuchos::Array< ScalarType > &x, const ScalarType &y)
Definition: gram_schmidt_example3.cpp:109
Stokhos::RecurrenceBasis::evaluate
virtual value_type evaluate(const value_type &point, ordinal_type order) const
Evaluate basis polynomial given by order order at given point point.
Definition: Stokhos_RecurrenceBasisImp.hpp:274
Stokhos::RecurrenceBasis::~RecurrenceBasis
virtual ~RecurrenceBasis()
Destructor.
Definition: Stokhos_RecurrenceBasisImp.hpp:108
Teuchos::Array::size
size_type size() const
Stokhos::RecurrenceBasis::computeTripleProductTensor
virtual Teuchos::RCP< Stokhos::Dense3Tensor< ordinal_type, value_type > > computeTripleProductTensor() const
Compute triple product tensor.
Definition: Stokhos_RecurrenceBasisImp.hpp:147
Sacado::UQ::abs
KOKKOS_INLINE_FUNCTION PCE< Storage > abs(const PCE< Storage > &a)
Definition: Sacado_UQ_PCE_Imp.hpp:1181
Teuchos::Array::resize
void resize(size_type new_size, const value_type &x=value_type())
Teuchos::SerialDenseMatrix::values
ScalarType * values() const
Stokhos::RecurrenceBasis::RecurrenceBasis
RecurrenceBasis(const std::string &name, ordinal_type p, bool normalize, GrowthPolicy growth=SLOW_GROWTH)
Constructor to be called by derived classes.
Definition: Stokhos_RecurrenceBasisImp.hpp:49
Teuchos::LAPACK::STEQR
void STEQR(const char &COMPZ, const OrdinalType &n, ScalarType *D, ScalarType *E, ScalarType *Z, const OrdinalType &ldz, ScalarType *WORK, OrdinalType *info) const
TotalOrderBasisUnitTest::value_type
double value_type
Definition: Stokhos_LexicographicTreeBasisUnitTest.cpp:70
Teuchos::rcp
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
Stokhos::GrowthPolicy
GrowthPolicy
Enumerated type for determining Smolyak growth policies.
Definition: Stokhos_RecurrenceBasis.hpp:50
Stokhos::RecurrenceBasis::computeDerivDoubleProductTensor
virtual Teuchos::RCP< Teuchos::SerialDenseMatrix< ordinal_type, value_type > > computeDerivDoubleProductTensor() const
Compute derivative double product tensor.
Definition: Stokhos_RecurrenceBasisImp.hpp:210
Stokhos::RecurrenceBasis::normalizeRecurrenceCoefficients
void normalizeRecurrenceCoefficients(Teuchos::Array< value_type > &alpha, Teuchos::Array< value_type > &beta, Teuchos::Array< value_type > &delta, Teuchos::Array< value_type > &gamma) const
Normalize coefficients.
Definition: Stokhos_RecurrenceBasisImp.hpp:465
Stokhos::Sparse3Tensor
Data structure storing a sparse 3-tensor C(i,j,k) in a a compressed format.
Definition: Stokhos_Sparse3Tensor.hpp:56
Teuchos_LAPACK.hpp
Teuchos::RCP
Stokhos::RecurrenceBasis::norm_squared
virtual const Teuchos::Array< value_type > & norm_squared() const
Return array storing norm-squared of each basis polynomial.
Definition: Stokhos_RecurrenceBasisImp.hpp:131
Stokhos::SLOW_GROWTH
Definition: Stokhos_RecurrenceBasis.hpp:51
Teuchos::Array< value_type >
Stokhos::RecurrenceBasis::computeSparseTripleProductTensor
virtual Teuchos::RCP< Stokhos::Sparse3Tensor< ordinal_type, value_type > > computeSparseTripleProductTensor(ordinal_type order) const
Compute triple product tensor.
Definition: Stokhos_RecurrenceBasisImp.hpp:177
Stokhos::RecurrenceBasis::getRecurrenceCoefficients
virtual void getRecurrenceCoefficients(Teuchos::Array< value_type > &alpha, Teuchos::Array< value_type > &beta, Teuchos::Array< value_type > &delta, Teuchos::Array< value_type > &gamma) const
Return recurrence coefficients defined by above formula.
Definition: Stokhos_RecurrenceBasisImp.hpp:451
TotalOrderBasisUnitTest::ordinal_type
int ordinal_type
Definition: Stokhos_LexicographicTreeBasisUnitTest.cpp:69
Stokhos::RecurrenceBasis::quadDegreeOfExactness
virtual ordinal_type quadDegreeOfExactness(ordinal_type n) const
Definition: Stokhos_RecurrenceBasisImp.hpp:415
Stokhos::RecurrenceBasis::print
virtual void print(std::ostream &os) const
Print basis to stream os.
Definition: Stokhos_RecurrenceBasisImp.hpp:303
Stokhos::Dense3Tensor
Data structure storing a dense 3-tensor C(i,j,k).
Definition: Stokhos_Dense3Tensor.hpp:57
Stokhos::RecurrenceBasis
Implementation of OneDOrthogPolyBasis based on the general three-term recurrence relationship:
Definition: Stokhos_RecurrenceBasis.hpp:85
Sacado::UQ::sqrt
KOKKOS_INLINE_FUNCTION PCE< Storage > sqrt(const PCE< Storage > &a)
Definition: Sacado_UQ_PCE_Imp.hpp:878
Stokhos::RecurrenceBasis::pointGrowth
virtual ordinal_type pointGrowth(ordinal_type n) const
Evaluate point growth rule for Smolyak-type bases.
Definition: Stokhos_RecurrenceBasisImp.hpp:435
Stokhos::RecurrenceBasis::size
virtual ordinal_type size() const
Return total size of basis (given by order() + 1).
Definition: Stokhos_RecurrenceBasisImp.hpp:123
Stokhos::RecurrenceBasis::evaluateBasesAndDerivatives
virtual void evaluateBasesAndDerivatives(const value_type &point, Teuchos::Array< value_type > &vals, Teuchos::Array< value_type > &derivs) const
Evaluate basis polynomials and their derivatives at given point point.
Definition: Stokhos_RecurrenceBasisImp.hpp:258
Teuchos_SerialDenseMatrix.hpp
j
j
Definition: Sacado_Fad_Exp_MP_Vector.hpp:527
Stokhos::RecurrenceBasis::setup
virtual void setup()
Setup basis after computing recurrence coefficients.
Definition: Stokhos_RecurrenceBasisImp.hpp:89
Stokhos::RecurrenceBasis::getName
virtual const std::string & getName() const
Return string name of basis.
Definition: Stokhos_RecurrenceBasisImp.hpp:336
Stokhos::RecurrenceBasis::order
virtual ordinal_type order() const
Return order of basis (largest monomial degree ).
Definition: Stokhos_RecurrenceBasisImp.hpp:115
cusp::detail::device::x
const IndexType const IndexType const IndexType const IndexType const ValueType const ValueType * x
Definition: csr_vector.h:260
Stokhos::RecurrenceBasis::evaluateBases
virtual void evaluateBases(const value_type &point, Teuchos::Array< value_type > &basis_pts) const
Evaluate each basis polynomial at given point point.
Definition: Stokhos_RecurrenceBasisImp.hpp:240
Stokhos::RecurrenceBasis::coefficientGrowth
virtual ordinal_type coefficientGrowth(ordinal_type n) const
Evaluate coefficient growth rule for Smolyak-type bases.
Definition: Stokhos_RecurrenceBasisImp.hpp:423
Stokhos::RecurrenceBasis::getQuadPoints
virtual void getQuadPoints(ordinal_type quad_order, Teuchos::Array< value_type > &points, Teuchos::Array< value_type > &weights, Teuchos::Array< Teuchos::Array< value_type > > &values) const
Compute quadrature points, weights, and values of basis polynomials at given set of points points.
Definition: Stokhos_RecurrenceBasisImp.hpp:344
Teuchos::SerialDenseMatrix< ordinal_type, value_type >
Teuchos::LAPACK< ordinal_type, value_type >
n
int n
Generated by   doxygen 1.8.16