Stokhos Package Browser (Single Doxygen Collection)  Version of the Day
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Public Types | Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type > Class Template Reference

An operator for building pseudo-spectral coefficients using tensor-product quadrature. More...

#include <Stokhos_TensorProductPseudoSpectralOperator.hpp>

Inheritance diagram for Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >:
Inheritance graph
[legend]

Public Types

typedef ordinal_t ordinal_type
 
typedef value_t value_type
 
typedef PseudoSpectralOperator
< ordinal_type, value_type,
point_compare_type > 
base_type
 
typedef base_type::point_type point_type
 
typedef base_type::point_set_type point_set_type
 
typedef base_type::point_map_type point_map_type
 
typedef base_type::iterator iterator
 
typedef base_type::const_iterator const_iterator
 
typedef base_type::set_iterator set_iterator
 
typedef
base_type::const_set_iterator 
const_set_iterator
 
typedef MultiIndex< ordinal_typemultiindex_type
 
- Public Types inherited from Stokhos::PseudoSpectralOperator< ordinal_t, value_t, point_compare_type >
typedef ordinal_t ordinal_type
 
typedef value_t value_type
 
typedef TensorProductElement
< ordinal_type, value_type
point_type
 
typedef std::map< point_type,
std::pair< value_type,
ordinal_type >
, point_compare_type > 
point_set_type
 
typedef Teuchos::Array
< point_type
point_map_type
 
typedef point_map_type::iterator iterator
 
typedef
point_map_type::const_iterator 
const_iterator
 
typedef point_set_type::iterator set_iterator
 
typedef
point_set_type::const_iterator 
const_set_iterator
 

Public Member Functions

 TensorProductPseudoSpectralOperator (const ProductBasis< ordinal_type, value_type > &product_basis, bool use_pst_=false, multiindex_type multiindex=multiindex_type(), const point_compare_type &point_compare=point_compare_type())
 Constructor. More...
 
virtual ~TensorProductPseudoSpectralOperator ()
 Destructor. More...
 
ordinal_type point_size () const
 Number of points. More...
 
ordinal_type coeff_size () const
 Number of coefficients. More...
 
iterator begin ()
 Iterator to begining of point set. More...
 
iterator end ()
 Iterator to end of point set. More...
 
const_iterator begin () const
 Iterator to begining of point set. More...
 
const_iterator end () const
 Iterator to end of point set. More...
 
set_iterator set_begin ()
 Iterator to begining of point set. More...
 
set_iterator set_end ()
 Iterator to end of point set. More...
 
const_set_iterator set_begin () const
 Iterator to begining of point set. More...
 
const_set_iterator set_end () const
 Iterator to end of point set. More...
 
ordinal_type index (const point_type &point) const
 Get point index for given point. More...
 
const point_typepoint (ordinal_type n) const
 Get point for given index. More...
 
virtual void transformQP2PCE (const value_type &alpha, const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &input, Teuchos::SerialDenseMatrix< ordinal_type, value_type > &result, const value_type &beta, bool trans=false) const
 Transform values at quadrature points to PCE coefficients. More...
 
virtual void transformPCE2QP (const value_type &alpha, const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &input, Teuchos::SerialDenseMatrix< ordinal_type, value_type > &result, const value_type &beta, bool trans=false) const
 Transform PCE coefficients to quadrature values. More...
 
- Public Member Functions inherited from Stokhos::PseudoSpectralOperator< ordinal_t, value_t, point_compare_type >
 PseudoSpectralOperator ()
 Constructor. More...
 
virtual ~PseudoSpectralOperator ()
 Destructor. More...
 

Protected Member Functions

void apply_direct (const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &A, const value_type &alpha, const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &input, Teuchos::SerialDenseMatrix< ordinal_type, value_type > &result, const value_type &beta, bool trans) const
 Apply transformation operator using direct method. More...
 
void apply_pst (const Teuchos::Array< Teuchos::SerialDenseMatrix< ordinal_type, value_type > > &Ak, const value_type &alpha, const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &input, Teuchos::SerialDenseMatrix< ordinal_type, value_type > &result, const value_type &beta, bool trans, bool reorder_input, bool reorder_result) const
 Apply tranformation operator using PST method. More...
 

Protected Attributes

bool use_pst
 Use partial-summation-transformation. More...
 
ordinal_type coeff_sz
 Number of coefficients. More...
 
ordinal_type dim
 Dimension. More...
 
point_set_type points
 Quadrature points. More...
 
point_map_type point_map
 Map index to point term. More...
 
bool reorder
 Do we need to reorder coefficients for PST. More...
 
Teuchos::Array< ordinal_typeperm
 Permutation array when reordering for PST. More...
 
Teuchos::SerialDenseMatrix
< ordinal_type, value_type
qp2pce
 Matrix mapping points to coefficients. More...
 
Teuchos::SerialDenseMatrix
< ordinal_type, value_type
pce2qp
 Matrix mapping coefficients to points. More...
 
Teuchos::Array
< Teuchos::SerialDenseMatrix
< ordinal_type, value_type > > 
qp2pce_k
 Matrix mapping points to coefficients for each dimension for PST. More...
 
Teuchos::Array
< Teuchos::SerialDenseMatrix
< ordinal_type, value_type > > 
pce2qp_k
 Matrix mapping coefficients to points for each dimension for PST. More...
 
Teuchos::BLAS< ordinal_type,
value_type
blas
 BLAS wrappers. More...
 

Detailed Description

template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
class Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >

An operator for building pseudo-spectral coefficients using tensor-product quadrature.

Definition at line 61 of file Stokhos_TensorProductPseudoSpectralOperator.hpp.

Member Typedef Documentation

template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
typedef ordinal_t Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::ordinal_type
template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
typedef value_t Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::value_type
template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
typedef PseudoSpectralOperator<ordinal_type,value_type,point_compare_type> Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::base_type
template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
typedef base_type::point_type Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::point_type
template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
typedef base_type::point_set_type Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::point_set_type
template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
typedef base_type::point_map_type Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::point_map_type
template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
typedef base_type::iterator Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::iterator
template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
typedef base_type::const_iterator Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::const_iterator
template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
typedef base_type::set_iterator Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::set_iterator
template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
typedef base_type::const_set_iterator Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::const_set_iterator
template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
typedef MultiIndex<ordinal_type> Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::multiindex_type

Constructor & Destructor Documentation

template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::TensorProductPseudoSpectralOperator ( const ProductBasis< ordinal_type, value_type > &  product_basis,
bool  use_pst_ = false,
multiindex_type  multiindex = multiindex_type(),
const point_compare_type &  point_compare = point_compare_type() 
)
inline

Constructor.

Definition at line 79 of file Stokhos_TensorProductPseudoSpectralOperator.hpp.

template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
virtual Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::~TensorProductPseudoSpectralOperator ( )
inlinevirtual

Destructor.

Definition at line 217 of file Stokhos_TensorProductPseudoSpectralOperator.hpp.

Member Function Documentation

template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
ordinal_type Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::point_size ( ) const
inlinevirtual
template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
ordinal_type Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::coeff_size ( ) const
inlinevirtual
template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
iterator Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::begin ( )
inlinevirtual
template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
iterator Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::end ( )
inlinevirtual
template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
const_iterator Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::begin ( ) const
inlinevirtual
template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
const_iterator Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::end ( ) const
inlinevirtual
template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
set_iterator Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::set_begin ( )
inlinevirtual
template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
set_iterator Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::set_end ( )
inlinevirtual
template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
const_set_iterator Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::set_begin ( ) const
inlinevirtual
template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
const_set_iterator Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::set_end ( ) const
inlinevirtual
template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
ordinal_type Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::index ( const point_type point) const
inlinevirtual
template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
const point_type& Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::point ( ordinal_type  n) const
inlinevirtual
template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
virtual void Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::transformQP2PCE ( const value_type alpha,
const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &  input,
Teuchos::SerialDenseMatrix< ordinal_type, value_type > &  result,
const value_type beta,
bool  trans = false 
) const
inlinevirtual

Transform values at quadrature points to PCE coefficients.

input is a vector storing values of a function at the quadrature points, and result will contain the resulting polynomial chaos coefficients. input and result can have multiple columns for vector-valued functions and set trans to true if these (multi-) vectors are layed out in a transposed fashion.

Implements Stokhos::PseudoSpectralOperator< ordinal_t, value_t, point_compare_type >.

Definition at line 268 of file Stokhos_TensorProductPseudoSpectralOperator.hpp.

template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
virtual void Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::transformPCE2QP ( const value_type alpha,
const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &  input,
Teuchos::SerialDenseMatrix< ordinal_type, value_type > &  result,
const value_type beta,
bool  trans = false 
) const
inlinevirtual

Transform PCE coefficients to quadrature values.

input is a vector storing polynomial chaos coefficients and result will contain the resulting values at the quadrature points. input and result can have multiple columns for vector-valued functions and set trans to true if these (multi-) vectors are layed out in a transposed fashion.

Implements Stokhos::PseudoSpectralOperator< ordinal_t, value_t, point_compare_type >.

Definition at line 288 of file Stokhos_TensorProductPseudoSpectralOperator.hpp.

template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
void Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::apply_direct ( const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &  A,
const value_type alpha,
const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &  input,
Teuchos::SerialDenseMatrix< ordinal_type, value_type > &  result,
const value_type beta,
bool  trans 
) const
inlineprotected

Apply transformation operator using direct method.

Definition at line 312 of file Stokhos_TensorProductPseudoSpectralOperator.hpp.

template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
void Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::apply_pst ( const Teuchos::Array< Teuchos::SerialDenseMatrix< ordinal_type, value_type > > &  Ak,
const value_type alpha,
const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &  input,
Teuchos::SerialDenseMatrix< ordinal_type, value_type > &  result,
const value_type beta,
bool  trans,
bool  reorder_input,
bool  reorder_result 
) const
inlineprotected

Apply tranformation operator using PST method.

For k=1,...,d, let A_k be the m_k-by-n_k quadrature operator for dimension k and A = A_1 ... A_d be m--n where m = m_1...m_d and n = n_1...n_d. For any two matrices B and C and vector x, (B C)x = vec( C*X*B^T ) = vec( C*(B*X^T)^T ) where x = vec(X), and the vec() operator makes a vector out of a matrix by stacking columns. Applying this formula recursively to A yields the simple algorithm for computing y = A*x (x is n-by-1 and y is m-by-1): X = x; for k=1:d n = n / n_k; X = reshape(X, n, n_k); X = A_k*X'; n = n * m_k; end y = reshape(X, m, 1);

When x has p columns, it is somehwat more complicated because the standard transpose above isn't correct. Instead the transpose needs to be applied to each p block in X.

The strategy here for dealing with transposed input and result is to transpose them when copying to/from the temporary buffers used in the algorithm.

Definition at line 366 of file Stokhos_TensorProductPseudoSpectralOperator.hpp.

Member Data Documentation

template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
bool Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::use_pst
protected

Use partial-summation-transformation.

Definition at line 498 of file Stokhos_TensorProductPseudoSpectralOperator.hpp.

template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
ordinal_type Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::coeff_sz
protected

Number of coefficients.

Definition at line 501 of file Stokhos_TensorProductPseudoSpectralOperator.hpp.

template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
ordinal_type Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::dim
protected

Dimension.

Definition at line 504 of file Stokhos_TensorProductPseudoSpectralOperator.hpp.

template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
point_set_type Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::points
protected

Quadrature points.

Definition at line 507 of file Stokhos_TensorProductPseudoSpectralOperator.hpp.

template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
point_map_type Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::point_map
protected

Map index to point term.

Definition at line 510 of file Stokhos_TensorProductPseudoSpectralOperator.hpp.

template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
bool Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::reorder
protected

Do we need to reorder coefficients for PST.

Definition at line 513 of file Stokhos_TensorProductPseudoSpectralOperator.hpp.

template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
Teuchos::Array<ordinal_type> Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::perm
protected

Permutation array when reordering for PST.

Definition at line 516 of file Stokhos_TensorProductPseudoSpectralOperator.hpp.

template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
Teuchos::SerialDenseMatrix<ordinal_type,value_type> Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::qp2pce
protected

Matrix mapping points to coefficients.

Definition at line 519 of file Stokhos_TensorProductPseudoSpectralOperator.hpp.

template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
Teuchos::SerialDenseMatrix<ordinal_type,value_type> Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::pce2qp
protected

Matrix mapping coefficients to points.

Definition at line 522 of file Stokhos_TensorProductPseudoSpectralOperator.hpp.

template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
Teuchos::Array< Teuchos::SerialDenseMatrix<ordinal_type,value_type> > Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::qp2pce_k
protected

Matrix mapping points to coefficients for each dimension for PST.

Definition at line 525 of file Stokhos_TensorProductPseudoSpectralOperator.hpp.

template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
Teuchos::Array< Teuchos::SerialDenseMatrix<ordinal_type,value_type> > Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::pce2qp_k
protected

Matrix mapping coefficients to points for each dimension for PST.

Definition at line 528 of file Stokhos_TensorProductPseudoSpectralOperator.hpp.

template<typename ordinal_t, typename value_t, typename point_compare_type = typename DefaultPointCompare<ordinal_t,value_t>::type>
Teuchos::BLAS<ordinal_type,value_type> Stokhos::TensorProductPseudoSpectralOperator< ordinal_t, value_t, point_compare_type >::blas
protected

BLAS wrappers.

Definition at line 531 of file Stokhos_TensorProductPseudoSpectralOperator.hpp.


The documentation for this class was generated from the following file: