44 #include "Teuchos_UnitTestHarness.hpp" 45 #include "Teuchos_TestingHelpers.hpp" 46 #include "Teuchos_UnitTestRepository.hpp" 47 #include "Teuchos_GlobalMPISession.hpp" 57 template <
typename OrdinalType,
typename ValueType>
80 Teuchos::Array< Teuchos::RCP<const Stokhos::OneDOrthogPolyBasis<ordinal_type,value_type> > > bases(
setup.d);
91 TEUCHOS_TEST_EQUALITY(sz, basis->size(), out, success);
93 std::ostream_iterator<ordinal_type> out_iterator(out,
" ");
98 out <<
"term " << term <<
" <= " <<
setup.p <<
" : ";
101 is_less = is_less && term[
j] <=
setup.p;
103 out <<
"passed" << std::endl;
105 out <<
"failed" << std::endl;
117 Teuchos::Array< Teuchos::RCP<const Stokhos::OneDOrthogPolyBasis<ordinal_type,value_type> > > bases(
setup.d);
128 TEUCHOS_TEST_EQUALITY(sz, basis->size(), out, success);
130 std::ostream_iterator<ordinal_type> out_iterator(out,
" ");
135 out <<
"term " << term <<
" <= " <<
setup.p <<
" : ";
138 is_less = is_less && term[
j] <=
setup.p;
140 out <<
"passed" << std::endl;
142 out <<
"failed" << std::endl;
154 Teuchos::Array< Teuchos::RCP<const Stokhos::OneDOrthogPolyBasis<ordinal_type,value_type> > > bases(
setup.d);
158 Teuchos::RCP< Stokhos::Sparse3Tensor<ordinal_type, value_type> > Cijk =
159 basis->computeTripleProductTensor();
169 Teuchos::Array< Teuchos::RCP<const Stokhos::OneDOrthogPolyBasis<ordinal_type,value_type> > > bases(
setup.d);
173 Teuchos::RCP< Stokhos::Sparse3Tensor<ordinal_type, value_type> > Cijk =
174 basis->computeLinearTripleProductTensor();
184 Teuchos::GlobalMPISession mpiSession(&argc, &
argv);
185 int res = Teuchos::UnitTestRepository::runUnitTestsFromMain(argc,
argv);
TEUCHOS_UNIT_TEST(Coefficients, Isotropic)
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)
UnitTestSetup< ordinal_type, value_type > setup
int main(int argc, char *argv[])
Multivariate orthogonal polynomial basis generated from a tensor product of univariate polynomials...
Legendre polynomial basis.