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++)
73 UnitTestSetup<int,double>
setup;
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);
int main(int argc, char *argv[])
Multivariate orthogonal polynomial basis generated from a total-order complete-polynomial tensor prod...
Legendre polynomial basis.
Defines quadrature for a tensor product basis by Smolyak sparse grids.
virtual const Teuchos::Array< value_type > & getQuadWeights() const
Get quadrature weights.
An isotropic total order index set.
UnitTestSetup< int, double > setup
TEUCHOS_UNIT_TEST(Stokhos_SmolyakSparseGridQuadrature, NumPoints)
Teuchos::RCP< const Stokhos::CompletePolynomialBasis< OrdinalType, ValueType > > basis