88 Array< RCP<const OneDOrthogPolyBasis<int,double> > > bases(d);
89 for (
int i=0; i<d; i++) {
90 bases[i] = rcp(
new HermiteBasis<int,double>(p));
92 RCP<const CompletePolynomialBasis<int,double> > basis =
93 rcp(
new CompletePolynomialBasis<int,double>(bases));
96 RCP<const Quadrature<int,double> > quad =
97 rcp(
new TensorProductQuadrature<int,double>(basis));
100 RCP<Sparse3Tensor<int,double> > Cijk =
101 basis->computeTripleProductTensor();
104 RCP<QuadOrthogPolyExpansion<int,double> > expn =
105 rcp(
new QuadOrthogPolyExpansion<int,double>(basis, Cijk, quad));
118 std::cout <<
"\tu = ";
120 std::cout <<
"\tv = ";
124 double mean = v.mean();
125 double std_dev = v.standard_deviation();
128 Teuchos::Array<double> pt(d);
129 for (
int i=0; i<d; i++)
131 double up = u.evaluate(pt);
133 double vp2 = v.evaluate(pt);
136 std::cout <<
"\tv mean = " << mean << std::endl;
137 std::cout <<
"\tv std. dev. = " << std_dev << std::endl;
138 std::cout <<
"\tv(0.25) (true) = " << vp << std::endl;
139 std::cout <<
"\tv(0.25) (pce) = " << vp2 << std::endl;
141 catch (std::exception& e) {
142 std::cout << e.what() << std::endl;