Stokhos Package Browser (Single Doxygen Collection) Version of the Day
Loading...
Searching...
No Matches
Stokhos_FullyAssembledOperator.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
44#ifndef STOKHOS_FULLY_ASSEMBLED_OPERATOR_HPP
45#define STOKHOS_FULLY_ASSEMBLED_OPERATOR_HPP
46
48#include "EpetraExt_MultiComm.h"
51#include "Teuchos_ParameterList.hpp"
52#include "EpetraExt_BlockCrsMatrix.h"
53#include <vector>
54
55namespace Stokhos {
56
62 public Stokhos::SGOperator,
63 public EpetraExt::BlockCrsMatrix {
64
65 public:
66
69 const Teuchos::RCP<const EpetraExt::MultiComm>& sg_comm,
70 const Teuchos::RCP<const Stokhos::OrthogPolyBasis<int,double> >& sg_basis,
71 const Teuchos::RCP<const Stokhos::EpetraSparse3Tensor>& epetraCijk,
72 const Teuchos::RCP<const Epetra_CrsGraph>& base_graph,
73 const Teuchos::RCP<const Epetra_Map>& domain_sg_map,
74 const Teuchos::RCP<const Epetra_Map>& range_sg_map,
75 const Teuchos::RCP<Teuchos::ParameterList>& params);
76
79
82
84 virtual void setupOperator(
85 const Teuchos::RCP<Stokhos::EpetraOperatorOrthogPoly >& poly);
86
88 virtual Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly >
90
92 virtual Teuchos::RCP<const Stokhos::EpetraOperatorOrthogPoly >
93 getSGPolynomial() const;
94
96
99
103 virtual int Apply(const Epetra_MultiVector& Input,
104 Epetra_MultiVector& Result) const;
105
107
108 private:
109
112
115
116 protected:
117
119 Teuchos::RCP<const EpetraExt::MultiComm> sg_comm;
120
122 Teuchos::RCP<const Stokhos::OrthogPolyBasis<int,double> > sg_basis;
123
125 Teuchos::RCP<const Stokhos::EpetraSparse3Tensor> epetraCijk;
126
128 Teuchos::RCP<const Epetra_Map> domain_sg_map;
129
131 Teuchos::RCP<const Epetra_Map> range_sg_map;
132
135
137 Teuchos::RCP<const Cijk_type> Cijk;
138
140 Teuchos::RCP<Stokhos::EpetraOperatorOrthogPoly > block_ops;
141
144
147
150
151 }; // class FullyAssembledOperator
152
153} // namespace Stokhos
154
155#endif // STOKHOS_FULLY_ASSEMBLED_OPERATOR_HPP
An Epetra operator representing the block stochastic Galerkin operator generated by fully assembling ...
Teuchos::RCP< const Epetra_Map > domain_sg_map
Stores domain SG map.
Stokhos::Sparse3Tensor< int, double > Cijk_type
Short-hand for Cijk.
bool scale_op
Flag indicating whether operator be scaled with <\psi_i^2>
Teuchos::RCP< const Stokhos::EpetraSparse3Tensor > epetraCijk
Stores Epetra Cijk tensor.
Teuchos::RCP< const EpetraExt::MultiComm > sg_comm
Stores SG parallel communicator.
virtual int Apply(const Epetra_MultiVector &Input, Epetra_MultiVector &Result) const
Wrap Apply() to add a timer.
FullyAssembledOperator(const FullyAssembledOperator &)
Private to prohibit copying.
Teuchos::RCP< const Epetra_Map > range_sg_map
Stores range SG map.
Teuchos::RCP< const Cijk_type > Cijk
Stores triple product tensor.
bool include_mean
Flag indicating whether to include mean term.
virtual Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly > getSGPolynomial()
Get SG polynomial.
FullyAssembledOperator & operator=(const FullyAssembledOperator &)
Private to prohibit copying.
Teuchos::RCP< const Stokhos::OrthogPolyBasis< int, double > > sg_basis
Stochastic Galerking basis.
bool only_use_linear
Flag indicating whether to only use linear terms.
Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly > block_ops
Stores operators.
virtual void setupOperator(const Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly > &poly)
Setup operator.
Abstract base class for multivariate orthogonal polynomials.
An abstract class to represent a generic stochastic Galerkin operator as an Epetra_Operator.
Data structure storing a sparse 3-tensor C(i,j,k) in a a compressed format.
Top-level namespace for Stokhos classes and functions.