47const Teuchos::RCP<
const OrthogPolyBasis<ordinal_type,value_type> >& basis,
48const Teuchos::RCP<
const Teuchos::Array< Teuchos::Array<value_type> > >& points,
49 const Teuchos::RCP<
const Teuchos::Array<value_type> >& weights)
50 : quad_points(points),
53 ordinal_type nqp = points->size();
54 Teuchos::RCP<Teuchos::Array< Teuchos::Array<value_type> > > qv =
55 Teuchos::rcp(
new Teuchos::Array< Teuchos::Array<value_type> >(nqp));
56 for (ordinal_type i=0; i<nqp; i++) {
57 (*qv)[i].resize(basis->size());
58 basis->evaluateBases((*points)[i], (*qv)[i]);
102print(std::ostream& os)
const
104 ordinal_type nqp = quad_weights->size();
105 os <<
"Sparse Grid Quadrature with " << nqp <<
" points:"
106 << std::endl <<
"Weight : Points" << std::endl;
107 for (ordinal_type i=0; i<nqp; i++) {
108 os << i <<
": " << (*quad_weights)[i] <<
" : ";
109 for (ordinal_type
j=0;
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;
116 for (ordinal_type i=0; i<nqp; i++) {
117 os << i <<
" " <<
": ";
118 for (ordinal_type
j=0;
119 j<static_cast<ordinal_type>((*quad_values)[i].size());
121 os << (*quad_values)[i][
j] <<
" ";