102 const OrdinalType d = 3;
103 const OrdinalType p = 5;
106 Teuchos::Array< Teuchos::RCP<const Stokhos::OneDOrthogPolyBasis<OrdinalType,ValueType> > > bases(d);
107 for (OrdinalType i=0; i<d; i++)
116 for (OrdinalType i=0; i<d; i++)
120 Teuchos::RCP<const Stokhos::Quadrature<OrdinalType,ValueType> > quad =
124 Teuchos::RCP<Stokhos::Sparse3Tensor<int,double> > Cijk =
125 basis->computeTripleProductTensor();
128 Teuchos::RCP<Teuchos::ParameterList> exp_params =
129 Teuchos::rcp(
new Teuchos::ParameterList);
130 exp_params->set(
"Use Quadrature for Times",
true);
142 Teuchos::Array< Teuchos::RCP<const Stokhos::OneDOrthogPolyBasis<OrdinalType,ValueType> > > st_bases(2);
145 Teuchos::RCP<const Stokhos::OrthogPolyBasis<OrdinalType,ValueType> > st_basis =
148 u_st.term(0, 0) =
u.
mean();
149 u_st.term(0, 1) = 1.0;
151 v_st.
term(1, 1) = 1.0;
154 Teuchos::Array<ValueType> st_points_0;
155 Teuchos::Array<ValueType> st_weights_0;
156 Teuchos::Array< Teuchos::Array<ValueType> > st_values_0;
157 st_bases[0]->getQuadPoints(p+1, st_points_0, st_weights_0, st_values_0);
158 Teuchos::Array<ValueType> st_points_1;
159 Teuchos::Array<ValueType> st_weights_1;
160 Teuchos::Array< Teuchos::Array<ValueType> > st_values_1;
161 st_bases[1]->getQuadPoints(p+1, st_points_1, st_weights_1, st_values_1);
162 Teuchos::RCP< Teuchos::Array< Teuchos::Array<ValueType> > > st_points =
163 Teuchos::rcp(
new Teuchos::Array< Teuchos::Array<ValueType> >(st_points_0.size()));
164 for (
int i=0; i<st_points_0.size(); i++) {
165 (*st_points)[i].resize(2);
166 (*st_points)[i][0] = st_points_0[i];
167 (*st_points)[i][1] = st_points_1[i];
169 Teuchos::RCP< Teuchos::Array<ValueType> > st_weights =
170 Teuchos::rcp(
new Teuchos::Array<ValueType>);
171 *st_weights = st_weights_0;
177 Teuchos::RCP<Stokhos::Sparse3Tensor<int,double> > gs_Cijk =
178 gs_basis->computeTripleProductTensor();
182 Teuchos::RCP< const Teuchos::Array< Teuchos::Array<ValueType> > > points =
184 Teuchos::RCP< const Teuchos::Array<ValueType> > weights = st_weights;