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>;
91 const bool useNativeSpMV=
false);
94 setMatrix (
const Teuchos::RCP<const operator_type>& A);
97 compute (multivector_type& W,
105 using import_type = Tpetra::Import<LO, GO, NT>;
106 using export_type = Tpetra::Export<LO, GO, NT>;
108 Teuchos::RCP<const operator_type> A_op_;
109 Teuchos::RCP<const crs_matrix_type> A_crs_;
110 Teuchos::RCP<const import_type> imp_;
111 Teuchos::RCP<const export_type> exp_;
112 std::unique_ptr<vector_type> X_colMap_;
113 std::unique_ptr<multivector_type> V1_;
115 Teuchos::RCP<vector_type> W_vec_, B_vec_, X_vec_;
123 importVector (vector_type& X_domMap);
125 bool canFuse (
const multivector_type& B)
const;
128 unfusedCase (multivector_type& W,
132 const operator_type& A,
137 fusedCase (vector_type& W,
141 const crs_matrix_type& A,