void
Stokhos::SmolyakPseudoSpectralOperator<ordinal_type,value_type,point_compare_type>::
transformQP2PCE(
const value_type& alpha,
const Teuchos::SerialDenseMatrix<ordinal_type,value_type>& input,
Teuchos::SerialDenseMatrix<ordinal_type,value_type>& result,
const value_type& beta,
bool trans) const {
if (use_smolyak)
transformQP2PCE_smolyak(alpha, input, result, beta, trans);
else
apply_direct(qp2pce, alpha, input, result, beta, trans);
}
void
Stokhos::SmolyakPseudoSpectralOperator<ordinal_type,value_type,point_compare_type>::
transformPCE2QP(
const value_type& alpha,
const Teuchos::SerialDenseMatrix<ordinal_type,value_type>& input,
Teuchos::SerialDenseMatrix<ordinal_type,value_type>& result,
const value_type& beta,
bool trans) const {
// Currently we use the direct method for mapping PCE->QP because the
// current implementation doesn't work. Need to evaluate tensor bases
// on all quad points, not just the quad points associated with that
// basis.
// if (use_smolyak)
// transformPCE2QP_smolyak(alpha, input, result, beta, trans);
// else
apply_direct(pce2qp, alpha, input, result, beta, trans);
}
Image * ReplacedImages::apply(Image * image, Image * src_image)
{
Image * ret = apply_direct(image, src_image);
index = 0;
return ret;
}
