42 #include "Teuchos_UnitTestHarness.hpp" 43 #include "Teuchos_TestingHelpers.hpp" 44 #include "Teuchos_UnitTestRepository.hpp" 45 #include "Teuchos_GlobalMPISession.hpp" 55 template <
typename OrdinalType,
typename ValueType>
78 Teuchos::Array< Teuchos::RCP<const Stokhos::OneDOrthogPolyBasis<ordinal_type,value_type> > > bases(
setup.d);
87 TEUCHOS_TEST_EQUALITY(sz, basis->size(), out, success);
89 std::ostream_iterator<ordinal_type> out_iterator(out,
" ");
94 out <<
"term " << term <<
" <= " <<
setup.p <<
" : ";
97 is_less = is_less && term[
j] <=
setup.p;
99 out <<
"passed" << std::endl;
101 out <<
"failed" << std::endl;
113 Teuchos::Array< Teuchos::RCP<const Stokhos::OneDOrthogPolyBasis<ordinal_type,value_type> > > bases(
setup.d);
123 TEUCHOS_TEST_EQUALITY(sz, basis->size(), out, success);
125 std::ostream_iterator<ordinal_type> out_iterator(out,
" ");
130 out <<
"term " << term <<
" <= " <<
setup.p <<
" : ";
133 is_less = is_less && term[
j] <=
j+1;
135 out <<
"passed" << std::endl;
137 out <<
"failed" << std::endl;
149 Teuchos::Array< Teuchos::RCP<const Stokhos::OneDOrthogPolyBasis<ordinal_type,value_type> > > bases(
setup.d);
153 Teuchos::RCP< Stokhos::Sparse3Tensor<ordinal_type, value_type> > Cijk =
154 basis->computeTripleProductTensor();
165 Teuchos::Array< Teuchos::RCP<const Stokhos::OneDOrthogPolyBasis<ordinal_type,value_type> > > bases(
setup.d);
169 Teuchos::RCP< Stokhos::Sparse3Tensor<ordinal_type, value_type> > Cijk =
170 basis->computeLinearTripleProductTensor();
180 Teuchos::GlobalMPISession mpiSession(&argc, &
argv);
181 int res = Teuchos::UnitTestRepository::runUnitTestsFromMain(argc,
argv);
ordinal_type n_choose_k(const ordinal_type &n, const ordinal_type &k)
Compute bionomial coefficient (n ; k) = n!/( k! (n-k)! )
Multivariate orthogonal polynomial basis generated from a total order tensor product of univariate po...
int main(int argc, char *argv[])
bool testSparse3Tensor(const Stokhos::Sparse3Tensor< ordinal_type, scalar_type > &Cijk, const Stokhos::ProductBasis< ordinal_type, scalar_type > &basis, const scalar_type &sparse_tol, const scalar_type &rel_tol, const scalar_type &abs_tol, Teuchos::FancyOStream &out, bool linear=false)
TEUCHOS_UNIT_TEST(LexicographicTreeCoefficients, Isotropic)
Legendre polynomial basis.
UnitTestSetup< ordinal_type, value_type > setup