42 #include "Teuchos_UnitTestHarness.hpp" 43 #include "Teuchos_TestingHelpers.hpp" 44 #include "Teuchos_UnitTestRepository.hpp" 45 #include "Teuchos_GlobalMPISession.hpp" 53 template <
typename OrdinalType,
typename ValueType>
55 Teuchos::RCP<const Stokhos::CompletePolynomialBasis<OrdinalType,ValueType> >
basis;
62 Teuchos::Array< Teuchos::RCP<const Stokhos::OneDOrthogPolyBasis<OrdinalType,ValueType> > > bases(
d);
63 for (OrdinalType i=0; i<
d; i++)
75 #ifdef HAVE_STOKHOS_DAKOTA 78 const Stokhos::SparseGridQuadrature<int,double> quad(
79 setup.basis,
setup.p, 1e-12, Pecos::MODERATE_RESTRICTED_GROWTH);
80 const Teuchos::Array<double>& weights = quad.getQuadWeights();
81 int nqp = weights.size();
90 <<
"Check: quad_weight.size() = " << nqp <<
" == " << nqp_gold
92 if (success) out <<
"Passed.";
103 setup.basis, index_set, 1e-12);
105 int nqp = weights.size();
114 <<
"Check: quad_weight.size() = " << nqp <<
" == " << nqp_gold
116 if (success) out <<
"Passed.";
125 Teuchos::GlobalMPISession mpiSession(&argc, &
argv);
126 return Teuchos::UnitTestRepository::runUnitTestsFromMain(argc,
argv);
Teuchos::RCP< const Stokhos::CompletePolynomialBasis< OrdinalType, ValueType > > basis
Defines quadrature for a tensor product basis by Smolyak sparse grids.
int main(int argc, char *argv[])
virtual const Teuchos::Array< value_type > & getQuadWeights() const
Get quadrature weights.
Multivariate orthogonal polynomial basis generated from a total-order complete-polynomial tensor prod...
An isotropic total order index set.
TEUCHOS_UNIT_TEST(Stokhos_SmolyakSparseGridQuadrature, NumPoints)
UnitTestSetup< int, double > setup