44 template <
typename ordinal_type,
typename value_type>
48 const Teuchos::RCP<
const Teuchos::Array< Teuchos::Array<value_type> > >& points,
49 const Teuchos::RCP<
const Teuchos::Array<value_type> >& weights)
50 : quad_points(points),
54 Teuchos::RCP<Teuchos::Array< Teuchos::Array<value_type> > > qv =
55 Teuchos::rcp(
new Teuchos::Array< Teuchos::Array<value_type> >(nqp));
57 (*qv)[i].resize(basis->size());
58 basis->evaluateBases((*points)[i], (*qv)[i]);
63 template <
typename ordinal_type,
typename value_type>
66 const Teuchos::RCP<
const Teuchos::Array< Teuchos::Array<value_type> > >& points,
67 const Teuchos::RCP<
const Teuchos::Array<value_type> >& weights,
68 const Teuchos::RCP<
const Teuchos::Array< Teuchos::Array<value_type> > >& values)
69 : quad_points(points),
70 quad_weights(weights),
75 template <
typename ordinal_type,
typename value_type>
76 const Teuchos::Array< Teuchos::Array<value_type> >&
83 template <
typename ordinal_type,
typename value_type>
84 const Teuchos::Array<value_type>&
91 template <
typename ordinal_type,
typename value_type>
92 const Teuchos::Array< Teuchos::Array<value_type> >&
99 template <
typename ordinal_type,
typename value_type>
105 os <<
"Sparse Grid Quadrature with " << nqp <<
" points:" 106 << std::endl <<
"Weight : Points" << std::endl;
108 os << i <<
": " << (*quad_weights)[i] <<
" : ";
110 j<static_cast<ordinal_type>((*quad_points)[i].size());
112 os << (*quad_points)[i][
j] <<
" ";
115 os <<
"Basis values at quadrature points:" << std::endl;
117 os << i <<
" " <<
": ";
119 j<static_cast<ordinal_type>((*quad_values)[i].size());
121 os << (*quad_values)[i][
j] <<
" ";
Teuchos::RCP< const Teuchos::Array< Teuchos::Array< value_type > > > quad_values
Quadrature values.
virtual const Teuchos::Array< value_type > & getQuadWeights() const
Get quadrature weights.
Abstract base class for multivariate orthogonal polynomials.
virtual const Teuchos::Array< Teuchos::Array< value_type > > & getBasisAtQuadPoints() const
Get values of basis at quadrature points.
virtual const Teuchos::Array< Teuchos::Array< value_type > > & getQuadPoints() const
Get quadrature points.
virtual std::ostream & print(std::ostream &os) const
Print quadrature data.
UserDefinedQuadrature(const Teuchos::RCP< const OrthogPolyBasis< ordinal_type, value_type > > &basis, const Teuchos::RCP< const Teuchos::Array< Teuchos::Array< value_type > > > &points, const Teuchos::RCP< const Teuchos::Array< value_type > > &weights)
Constructor.