Ifpack2 Templated Preconditioning Package Version 1.0
Loading...
Searching...
No Matches
Ifpack2_DenseContainer_decl.hpp
Go to the documentation of this file.
1/*@HEADER
2// ***********************************************************************
3//
4// Ifpack2: Templated Object-Oriented Algebraic Preconditioner Package
5// Copyright (2009) Sandia Corporation
6//
7// Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
8// license for use of this work by or on behalf of the U.S. Government.
9//
10// Redistribution and use in source and binary forms, with or without
11// modification, are permitted provided that the following conditions are
12// met:
13//
14// 1. Redistributions of source code must retain the above copyright
15// notice, this list of conditions and the following disclaimer.
16//
17// 2. Redistributions in binary form must reproduce the above copyright
18// notice, this list of conditions and the following disclaimer in the
19// documentation and/or other materials provided with the distribution.
20//
21// 3. Neither the name of the Corporation nor the names of the
22// contributors may be used to endorse or promote products derived from
23// this software without specific prior written permission.
24//
25// THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
26// EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
27// IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
28// PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
29// CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
30// EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
31// PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
32// PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
33// LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
34// NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
35// SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
36//
37// Questions? Contact Michael A. Heroux (maherou@sandia.gov)
38//
39// ***********************************************************************
40//@HEADER
41*/
42
43#ifndef IFPACK2_DENSECONTAINER_DECL_HPP
44#define IFPACK2_DENSECONTAINER_DECL_HPP
45
48
49#include "Ifpack2_Container.hpp"
50#include "Tpetra_MultiVector.hpp"
51#include "Tpetra_Map.hpp"
52#include "Tpetra_RowMatrix.hpp"
53#include "Teuchos_SerialDenseMatrix.hpp"
54#include "Tpetra_BlockCrsMatrix_decl.hpp"
55
56namespace Ifpack2 {
57
101
102template<class MatrixType, class LocalScalarType>
104: public ContainerImpl<MatrixType, LocalScalarType>
105{
107
108private:
115 using matrix_type = MatrixType;
117 using LSC = LocalScalarType;
118
120 using typename Container<MatrixType>::SC;
122 using typename Container<MatrixType>::LO;
124 using typename Container<MatrixType>::GO;
126 using typename Container<MatrixType>::NO;
127
128 using typename Container<MatrixType>::mv_type;
129 using local_mv_type = Tpetra::MultiVector<LSC, LO, GO, NO>;
130 using typename Container<MatrixType>::map_type;
131 using typename Container<MatrixType>::vector_type;
132 using typename Container<MatrixType>::import_type;
133
134 using typename Container<MatrixType>::ISC;
135 using typename ContainerImpl<MatrixType, LocalScalarType>::LISC;
136 using typename Container<MatrixType>::HostView;
137 using typename ContainerImpl<MatrixType, LocalScalarType>::HostViewLocal;
138 using typename ContainerImpl<MatrixType, LocalScalarType>::HostSubviewLocal;
139 using typename ContainerImpl<MatrixType, LocalScalarType>::ConstHostSubviewLocal;
140
141 static_assert(std::is_same<MatrixType, Tpetra::RowMatrix<SC, LO, GO, NO>>::value,
142 "Ifpack2::DenseContainer: Please use MatrixType = Tpetra::RowMatrix.");
143
152 using typename Container<MatrixType>::row_matrix_type;
153
154 using block_crs_matrix_type = Tpetra::BlockCrsMatrix<SC, LO, GO, NO>;
156
157public:
159
160
171 DenseContainer (const Teuchos::RCP<const row_matrix_type>& matrix,
172 const Teuchos::Array<Teuchos::Array<LO> >& partitions,
173 const Teuchos::RCP<const import_type>& importer,
174 bool pointIndexed);
175
177
179 virtual ~DenseContainer ();
180
182
184
186
188
190 virtual void initialize ();
191
193 virtual void compute ();
194
195 void clearBlocks();
196
198
200
204 virtual std::ostream& print (std::ostream& os) const;
205
207
209
211 virtual std::string description () const;
212
214 virtual void
215 describe (Teuchos::FancyOStream &out,
216 const Teuchos::EVerbosityLevel verbLevel =
217 Teuchos::Describable::verbLevel_default) const;
218
220
222 static std::string getName();
223
224private:
226 void extract();
227
231 void factor ();
232
241 void
242 solveBlock(ConstHostSubviewLocal X,
243 HostSubviewLocal Y,
244 int blockIndex,
245 Teuchos::ETransp mode,
246 const LSC alpha,
247 const LSC beta) const;
248
250 std::vector<Teuchos::SerialDenseMatrix<int, LSC>> diagBlocks_;
251
253 mutable Teuchos::Array<int> ipiv_;
254
256 bool hasBlockCrsMatrix_;
257
259 Teuchos::Array<LSC> scalars_;
260
262 Teuchos::Array<GO> scalarOffsets_;
263};
264
265} // namespace Ifpack2
266
267#endif // IFPACK2_DENSECONTAINER_DECL_HPP
The implementation of the numerical features of Container (Jacobi, Gauss-Seidel, SGS)....
Definition: Ifpack2_Container_decl.hpp:344
typename Kokkos::Details::ArithTraits< SC >::val_type ISC
Internal representation of Scalar in Kokkos::View.
Definition: Ifpack2_Container_decl.hpp:135
typename mv_type::dual_view_type::t_host HostView
Definition: Ifpack2_Container_decl.hpp:139
Store and solve a local dense linear problem.
Definition: Ifpack2_DenseContainer_decl.hpp:105
virtual std::string description() const
A one-line description of this object.
Definition: Ifpack2_DenseContainer_def.hpp:385
virtual void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel=Teuchos::Describable::verbLevel_default) const
Print the object with some verbosity level to the given FancyOStream.
Definition: Ifpack2_DenseContainer_def.hpp:408
static std::string getName()
Get the name of this container type for Details::constructContainer()
Definition: Ifpack2_DenseContainer_def.hpp:434
virtual std::ostream & print(std::ostream &os) const
Print information about this object to the given output stream.
Definition: Ifpack2_DenseContainer_def.hpp:374
virtual void compute()
Extract the local diagonal block and prepare the solver.
Definition: Ifpack2_DenseContainer_def.hpp:114
virtual ~DenseContainer()
Destructor (declared virtual for memory safety of derived classes).
Definition: Ifpack2_DenseContainer_def.hpp:93
virtual void initialize()
Do all set-up operations that only require matrix structure.
Definition: Ifpack2_DenseContainer_def.hpp:98
Preconditioners and smoothers for Tpetra sparse matrices.
Definition: Ifpack2_AdditiveSchwarz_decl.hpp:74