Stokhos Package Browser (Single Doxygen Collection) Version of the Day
Loading...
Searching...
No Matches
Stokhos_UserDefinedQuadratureImp.hpp
Go to the documentation of this file.
1// $Id$
2// $Source$
3// @HEADER
4// ***********************************************************************
5//
6// Stokhos Package
7// Copyright (2009) Sandia Corporation
8//
9// Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
10// license for use of this work by or on behalf of the U.S. Government.
11//
12// Redistribution and use in source and binary forms, with or without
13// modification, are permitted provided that the following conditions are
14// met:
15//
16// 1. Redistributions of source code must retain the above copyright
17// notice, this list of conditions and the following disclaimer.
18//
19// 2. Redistributions in binary form must reproduce the above copyright
20// notice, this list of conditions and the following disclaimer in the
21// documentation and/or other materials provided with the distribution.
22//
23// 3. Neither the name of the Corporation nor the names of the
24// contributors may be used to endorse or promote products derived from
25// this software without specific prior written permission.
26//
27// THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
28// EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
29// IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
30// PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
31// CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
32// EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
33// PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
34// PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
35// LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
36// NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
37// SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
38//
39// Questions? Contact Eric T. Phipps (etphipp@sandia.gov).
40//
41// ***********************************************************************
42// @HEADER
43
44template <typename ordinal_type, typename value_type>
47const Teuchos::RCP<const OrthogPolyBasis<ordinal_type,value_type> >& basis,
48const Teuchos::RCP<const Teuchos::Array< Teuchos::Array<value_type> > >& points,
49 const Teuchos::RCP<const Teuchos::Array<value_type> >& weights)
50 : quad_points(points),
51 quad_weights(weights)
52{
53 ordinal_type nqp = points->size();
54 Teuchos::RCP<Teuchos::Array< Teuchos::Array<value_type> > > qv =
55 Teuchos::rcp(new Teuchos::Array< Teuchos::Array<value_type> >(nqp));
56 for (ordinal_type i=0; i<nqp; i++) {
57 (*qv)[i].resize(basis->size());
58 basis->evaluateBases((*points)[i], (*qv)[i]);
59 }
60 quad_values = qv;
61}
62
63template <typename ordinal_type, typename value_type>
66const Teuchos::RCP<const Teuchos::Array< Teuchos::Array<value_type> > >& points,
67const Teuchos::RCP<const Teuchos::Array<value_type> >& weights,
68const Teuchos::RCP<const Teuchos::Array< Teuchos::Array<value_type> > >& values)
69 : quad_points(points),
70 quad_weights(weights),
71 quad_values(values)
72{
73}
74
75template <typename ordinal_type, typename value_type>
76const Teuchos::Array< Teuchos::Array<value_type> >&
82
83template <typename ordinal_type, typename value_type>
84const Teuchos::Array<value_type>&
90
91template <typename ordinal_type, typename value_type>
92const Teuchos::Array< Teuchos::Array<value_type> >&
98
99template <typename ordinal_type, typename value_type>
100std::ostream&
102print(std::ostream& os) const
103{
104 ordinal_type nqp = quad_weights->size();
105 os << "Sparse Grid Quadrature with " << nqp << " points:"
106 << std::endl << "Weight : Points" << std::endl;
107 for (ordinal_type i=0; i<nqp; i++) {
108 os << i << ": " << (*quad_weights)[i] << " : ";
109 for (ordinal_type j=0;
110 j<static_cast<ordinal_type>((*quad_points)[i].size());
111 j++)
112 os << (*quad_points)[i][j] << " ";
113 os << std::endl;
114 }
115 os << "Basis values at quadrature points:" << std::endl;
116 for (ordinal_type i=0; i<nqp; i++) {
117 os << i << " " << ": ";
118 for (ordinal_type j=0;
119 j<static_cast<ordinal_type>((*quad_values)[i].size());
120 j++)
121 os << (*quad_values)[i][j] << " ";
122 os << std::endl;
123 }
124
125 return os;
126}
virtual const Teuchos::Array< Teuchos::Array< value_type > > & getQuadPoints() const
Get quadrature points.
virtual const Teuchos::Array< Teuchos::Array< value_type > > & getBasisAtQuadPoints() const
Get values of basis at quadrature points.
Teuchos::RCP< const Teuchos::Array< Teuchos::Array< value_type > > > quad_values
Quadrature values.
UserDefinedQuadrature(const Teuchos::RCP< const OrthogPolyBasis< ordinal_type, value_type > > &basis, const Teuchos::RCP< const Teuchos::Array< Teuchos::Array< value_type > > > &points, const Teuchos::RCP< const Teuchos::Array< value_type > > &weights)
Constructor.
virtual std::ostream & print(std::ostream &os) const
Print quadrature data.
virtual const Teuchos::Array< value_type > & getQuadWeights() const
Get quadrature weights.