48#ifndef IFPACK2_DETAILS_SCALEDDAMPEDRESIDUAL_DECL_HPP
49#define IFPACK2_DETAILS_SCALEDDAMPEDRESIDUAL_DECL_HPP
51#include "Ifpack2_config.h"
52#include "Tpetra_CrsMatrix_fwd.hpp"
53#include "Tpetra_MultiVector_fwd.hpp"
54#include "Tpetra_Operator_fwd.hpp"
55#include "Tpetra_Vector_fwd.hpp"
56#include "Tpetra_Export_fwd.hpp"
57#include "Tpetra_Import_fwd.hpp"
58#include "Teuchos_RCP.hpp"
76template<
class TpetraOperatorType>
79 using SC =
typename TpetraOperatorType::scalar_type;
80 using LO =
typename TpetraOperatorType::local_ordinal_type;
81 using GO =
typename TpetraOperatorType::global_ordinal_type;
82 using NT =
typename TpetraOperatorType::node_type;
84 using crs_matrix_type = Tpetra::CrsMatrix<SC, LO, GO, NT>;
85 using multivector_type = Tpetra::MultiVector<SC, LO, GO, NT>;
86 using operator_type = Tpetra::Operator<SC, LO, GO, NT>;
87 using vector_type = Tpetra::Vector<SC, LO, GO, NT>;
93 setMatrix (
const Teuchos::RCP<const operator_type>& A);
96 compute (multivector_type& W,
104 using import_type = Tpetra::Import<LO, GO, NT>;
105 using export_type = Tpetra::Export<LO, GO, NT>;
107 Teuchos::RCP<const operator_type> A_op_;
108 Teuchos::RCP<const crs_matrix_type> A_crs_;
109 Teuchos::RCP<const import_type> imp_;
110 Teuchos::RCP<const export_type> exp_;
111 std::unique_ptr<vector_type> X_colMap_;
112 std::unique_ptr<multivector_type> V1_;
114 Teuchos::RCP<vector_type> W_vec_, B_vec_, X_vec_;
118 importVector (vector_type& X_domMap);
120 bool canFuse (
const multivector_type& B)
const;
123 unfusedCase (multivector_type& W,
127 const operator_type& A,
132 fusedCase (vector_type& W,
136 const crs_matrix_type& A,
Compute scaled damped residual for Chebyshev.
Definition: Ifpack2_Details_ScaledDampedResidual_decl.hpp:77
Ifpack2 implementation details.
Preconditioners and smoothers for Tpetra sparse matrices.
Definition: Ifpack2_AdditiveSchwarz_decl.hpp:74