Xpetra  Version of the Day
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Pages
Xpetra_TpetraBlockCrsMatrix_def.hpp
Go to the documentation of this file.
1 // @HEADER
2 // *****************************************************************************
3 // Xpetra: A linear algebra interface package
4 //
5 // Copyright 2012 NTESS and the Xpetra contributors.
6 // SPDX-License-Identifier: BSD-3-Clause
7 // *****************************************************************************
8 // @HEADER
9 
10 #ifndef XPETRA_TPETRABLOCKCRSMATRIX_DEF_HPP
11 #define XPETRA_TPETRABLOCKCRSMATRIX_DEF_HPP
12 
13 #include "Xpetra_TpetraBlockCrsMatrix_decl.hpp"
14 #include "Xpetra_TpetraCrsGraph.hpp"
15 
16 namespace Xpetra {
17 
19 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
20 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
21  TpetraBlockCrsMatrix(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rowMap,
22  size_t maxNumEntriesPerRow,
23  const Teuchos::RCP<Teuchos::ParameterList> &params) {
24  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in" + std::string(__FILE__) + ":" + std::to_string(__LINE__));
25 }
26 
28 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
29 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
30  TpetraBlockCrsMatrix(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rowMap,
31  const ArrayRCP<const size_t> &NumEntriesPerRowToAlloc,
32  const Teuchos::RCP<Teuchos::ParameterList> &params) {
33  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in" + std::string(__FILE__) + ":" + std::to_string(__LINE__));
34 }
35 
37 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
38 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
39  TpetraBlockCrsMatrix(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rowMap,
40  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &colMap,
41  size_t maxNumEntriesPerRow,
42  const Teuchos::RCP<Teuchos::ParameterList> &params) {
43  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in" + std::string(__FILE__) + ":" + std::to_string(__LINE__));
44 }
45 
47 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
48 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
49  TpetraBlockCrsMatrix(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rowMap,
50  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &colMap,
51  const ArrayRCP<const size_t> &NumEntriesPerRowToAlloc,
52  const Teuchos::RCP<Teuchos::ParameterList> &params) {
53  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in" + std::string(__FILE__) + ":" + std::to_string(__LINE__));
54 }
55 
57 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
58 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
59  TpetraBlockCrsMatrix(const Teuchos::RCP<const CrsGraph<LocalOrdinal, GlobalOrdinal, Node> > &graph,
60  const Teuchos::RCP<Teuchos::ParameterList> &params)
61 // : mtx_(Teuchos::rcp(new Tpetra::BlockCrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >(toTpetra(graph), params)))
62 // * there is no Tpetra::BlockCrsMatrix(graph, params) c'tor. We throw anyways here so no need to set mtx_.
63 {
64  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in" + std::string(__FILE__) + ":" + std::to_string(__LINE__));
65 }
66 
68 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
69 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
70  TpetraBlockCrsMatrix(const Teuchos::RCP<const CrsGraph<LocalOrdinal, GlobalOrdinal, Node> > &graph,
71  const LocalOrdinal blockSize)
72  : mtx_(Teuchos::rcp(new Tpetra::BlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(*toTpetra(graph), blockSize))) {}
73 
75 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
76 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
77  TpetraBlockCrsMatrix(const Teuchos::RCP<const CrsGraph<LocalOrdinal, GlobalOrdinal, Node> > &graph,
78  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &pointDomainMap,
79  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &pointRangeMap,
80  const LocalOrdinal blockSize)
81  : mtx_(Teuchos::rcp(new Tpetra::BlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(*toTpetra(graph), *toTpetra(pointDomainMap), *toTpetra(pointRangeMap), blockSize))) {}
82 
84 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
85 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
86  TpetraBlockCrsMatrix(const Teuchos::RCP<const Tpetra::BlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > &sourceMatrix,
87  const Import<LocalOrdinal, GlobalOrdinal, Node> &importer,
88  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &domainMap,
89  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rangeMap,
90  const Teuchos::RCP<Teuchos::ParameterList> &params) {
91  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in" + std::string(__FILE__) + ":" + std::to_string(__LINE__));
92 }
93 
95 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
96 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
97  TpetraBlockCrsMatrix(const Teuchos::RCP<const Tpetra::BlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > &sourceMatrix,
98  const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter,
99  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &domainMap,
100  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rangeMap,
101  const Teuchos::RCP<Teuchos::ParameterList> &params) {
102  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in" + std::string(__FILE__) + ":" + std::to_string(__LINE__));
103 }
104 
106 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
107 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
108  TpetraBlockCrsMatrix(const Teuchos::RCP<const Tpetra::BlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > &sourceMatrix,
109  const Import<LocalOrdinal, GlobalOrdinal, Node> &RowImporter,
110  const Teuchos::RCP<const Import<LocalOrdinal, GlobalOrdinal, Node> > DomainImporter,
111  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &domainMap,
112  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rangeMap,
113  const Teuchos::RCP<Teuchos::ParameterList> &params) {
114  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in" + std::string(__FILE__) + ":" + std::to_string(__LINE__));
115 }
116 
118 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
119 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
120  TpetraBlockCrsMatrix(const Teuchos::RCP<const Tpetra::BlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > &sourceMatrix,
121  const Export<LocalOrdinal, GlobalOrdinal, Node> &RowExporter,
122  const Teuchos::RCP<const Export<LocalOrdinal, GlobalOrdinal, Node> > DomainExporter,
123  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &domainMap,
124  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rangeMap,
125  const Teuchos::RCP<Teuchos::ParameterList> &params) {
126  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in" + std::string(__FILE__) + ":" + std::to_string(__LINE__));
127 }
128 
130 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
131 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
132  ~TpetraBlockCrsMatrix() {}
133 
135 
137 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
138 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
139  insertGlobalValues(GlobalOrdinal globalRow,
140  const ArrayView<const GlobalOrdinal> &cols,
141  const ArrayView<const Scalar> &vals) {
142  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in" + std::string(__FILE__) + ":" + std::to_string(__LINE__));
143 }
144 
146 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
147 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
148  insertLocalValues(LocalOrdinal localRow,
149  const ArrayView<const LocalOrdinal> &cols,
150  const ArrayView<const Scalar> &vals) {
151  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in" + std::string(__FILE__) + ":" + std::to_string(__LINE__));
152 }
153 
155 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
156 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
157  replaceGlobalValues(GlobalOrdinal globalRow,
158  const ArrayView<const GlobalOrdinal> &cols,
159  const ArrayView<const Scalar> &vals) {
160  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in" + std::string(__FILE__) + ":" + std::to_string(__LINE__));
161 }
162 
164 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
165 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
166  replaceLocalValues(LocalOrdinal localRow, const ArrayView<const LocalOrdinal> &cols, const ArrayView<const Scalar> &vals) {
167  XPETRA_MONITOR("TpetraBlockCrsMatrix::replaceLocalValues");
168  mtx_->replaceLocalValues(localRow, cols.getRawPtr(), vals.getRawPtr(), cols.size());
169 }
170 
172 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
173 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
174  setAllToScalar(const Scalar &alpha) {
175  XPETRA_MONITOR("TpetraBlockCrsMatrix::setAllToScalar");
176  mtx_->setAllToScalar(alpha);
177 }
178 
180 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
181 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
182  scale(const Scalar &alpha) {
183  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in " + std::string(__FILE__) + ":" + std::to_string(__LINE__));
184 }
185 
187 //** \warning This is an expert-only routine and should not be called from user code. (not implemented)
188 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
189 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
190  allocateAllValues(size_t numNonZeros, ArrayRCP<size_t> &rowptr, ArrayRCP<LocalOrdinal> &colind, ArrayRCP<Scalar> &values) {
191  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in " + std::string(__FILE__) + ":" + std::to_string(__LINE__));
192 }
193 
195 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
196 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
197  setAllValues(const ArrayRCP<size_t> &rowptr, const ArrayRCP<LocalOrdinal> &colind, const ArrayRCP<Scalar> &values) {
198  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in " + std::string(__FILE__) + ":" + std::to_string(__LINE__));
199 }
200 
202 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
203 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
204  getAllValues(ArrayRCP<const size_t> &rowptr,
205  ArrayRCP<const LocalOrdinal> &colind,
206  ArrayRCP<const Scalar> &values) const {
207  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in " + std::string(__FILE__) + ":" + std::to_string(__LINE__));
208 }
209 
211 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
212 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
213  getAllValues(ArrayRCP<Scalar> &values) {
214  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in " + std::string(__FILE__) + ":" + std::to_string(__LINE__));
215 }
216 
218 
219 // Transformational Methods
221 
222 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
223 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
224  resumeFill(const RCP<ParameterList> &params) {
225  /*noop*/
226 }
227 
228 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
229 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
230  fillComplete(const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &domainMap,
231  const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rangeMap,
232  const RCP<ParameterList> &params) {
233  /*noop*/
234 }
235 
236 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
237 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
238  fillComplete(const RCP<ParameterList> &params) {
239  /*noop*/
240 }
241 
242 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
243 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
244  replaceDomainMapAndImporter(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &newDomainMap,
245  Teuchos::RCP<const Import<LocalOrdinal, GlobalOrdinal, Node> > &newImporter) {
246  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in " + std::string(__FILE__) + ":" + std::to_string(__LINE__));
247 }
248 
249 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
250 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
251  expertStaticFillComplete(const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &domainMap,
252  const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rangeMap,
253  const RCP<const Import<LocalOrdinal, GlobalOrdinal, Node> > &importer,
254  const RCP<const Export<LocalOrdinal, GlobalOrdinal, Node> > &exporter,
255  const RCP<ParameterList> &params) {
256  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in " + std::string(__FILE__) + ":" + std::to_string(__LINE__));
257 }
258 
260 
262 
263 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
264 const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >
265 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
266  getRowMap() const {
267  XPETRA_MONITOR("TpetraBlockCrsMatrix::getRowMap");
268  return toXpetra(mtx_->getRowMap());
269 }
270 
271 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
272 const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >
273 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
274  getColMap() const {
275  XPETRA_MONITOR("TpetraBlockCrsMatrix::getColMap");
276  return toXpetra(mtx_->getColMap());
277 }
278 
279 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
280 RCP<const CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
281 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
282  getCrsGraph() const {
283  XPETRA_MONITOR("TpetraBlockCrsMatrix::getCrsGraph");
284  using G_t = Tpetra::CrsGraph<LocalOrdinal, GlobalOrdinal, Node>;
285  using G_x = TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node>;
286  RCP<G_t> t_graph = Teuchos::rcp_const_cast<G_t>(Teuchos::rcpFromRef(mtx_->getCrsGraph()));
287  RCP<const G_x> x_graph = rcp(new G_x(t_graph));
288  return x_graph;
289 }
290 
291 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
292 global_size_t
293 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
294  getGlobalNumRows() const {
295  XPETRA_MONITOR("TpetraBlockCrsMatrix::getGlobalNumRows");
296  return mtx_->getGlobalNumRows();
297 }
298 
299 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
300 global_size_t
301 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
302  getGlobalNumCols() const {
303  XPETRA_MONITOR("TpetraBlockCrsMatrix::getGlobalNumCols");
304  return mtx_->getGlobalNumCols();
305 }
306 
307 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
308 size_t
309 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
310  getLocalNumRows() const {
311  XPETRA_MONITOR("TpetraBlockCrsMatrix::getLocalNumRows");
312  return mtx_->getLocalNumRows();
313 }
314 
315 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
316 size_t
317 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
318  getLocalNumCols() const {
319  XPETRA_MONITOR("TpetraBlockCrsMatrix::getLocalNumCols");
320  return mtx_->getLocalNumCols();
321 }
322 
323 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
324 global_size_t
325 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
326  getGlobalNumEntries() const {
327  XPETRA_MONITOR("TpetraBlockCrsMatrix::getGlobalNumEntries");
328  return mtx_->getGlobalNumEntries();
329 }
330 
331 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
332 size_t
333 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
334  getLocalNumEntries() const {
335  XPETRA_MONITOR("TpetraBlockCrsMatrix::getLocalNumEntries");
336  return mtx_->getLocalNumEntries();
337 }
338 
339 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
340 size_t
341 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
342  getNumEntriesInLocalRow(LocalOrdinal localRow) const {
343  XPETRA_MONITOR("TpetraBlockCrsMatrix::getNumEntriesInLocalRow");
344  return mtx_->getNumEntriesInLocalRow(localRow);
345 }
346 
347 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
348 size_t
349 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
350  getNumEntriesInGlobalRow(GlobalOrdinal globalRow) const {
351  XPETRA_MONITOR("TpetraBlockCrsMatrix::getNumEntriesInGlobalRow");
352  return mtx_->getNumEntriesInGlobalRow(globalRow);
353 }
354 
355 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
356 size_t TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::getGlobalMaxNumRowEntries() const {
357  XPETRA_MONITOR("TpetraBlockCrsMatrix::getGlobalMaxNumRowEntries");
358  return mtx_->getGlobalMaxNumRowEntries();
359 }
360 
361 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
362 size_t TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::getLocalMaxNumRowEntries() const {
363  XPETRA_MONITOR("TpetraBlockCrsMatrix::getLocalMaxNumRowEntries");
364  return mtx_->getLocalMaxNumRowEntries();
365 }
366 
367 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
368 bool TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::isLocallyIndexed() const {
369  XPETRA_MONITOR("TpetraBlockCrsMatrix::isLocallyIndexed");
370  return mtx_->isLocallyIndexed();
371 }
372 
373 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
374 bool TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::isGloballyIndexed() const {
375  XPETRA_MONITOR("TpetraBlockCrsMatrix::isGloballyIndexed");
376  return mtx_->isGloballyIndexed();
377 }
378 
379 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
380 bool TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::isFillComplete() const {
381  XPETRA_MONITOR("TpetraBlockCrsMatrix::isFillComplete");
382  return mtx_->isFillComplete();
383 }
384 
385 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
386 bool TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::isFillActive() const {
387  XPETRA_MONITOR("TpetraBlockCrsMatrix::isFillActive");
388  return false;
389 }
390 
391 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
392 typename ScalarTraits<Scalar>::magnitudeType TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::getFrobeniusNorm() const {
393  XPETRA_MONITOR("TpetraBlockCrsMatrix::getFrobeniusNorm");
394  return mtx_->getFrobeniusNorm();
395 }
396 
397 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
398 bool TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::supportsRowViews() const {
399  XPETRA_MONITOR("TpetraBlockCrsMatrix::supportsRowViews");
400  return mtx_->supportsRowViews();
401 }
402 
403 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
404 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
405  getLocalRowCopy(LocalOrdinal LocalRow,
406  const ArrayView<LocalOrdinal> &Indices,
407  const ArrayView<Scalar> &Values,
408  size_t &NumEntries) const {
409  XPETRA_MONITOR("TpetraBlockCrsMatrix::getLocalRowCopy");
410  typename Tpetra::CrsGraph<LocalOrdinal, GlobalOrdinal, Node>::nonconst_local_inds_host_view_type indices("indices", Indices.size());
411  typename Tpetra::CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::nonconst_values_host_view_type values("values", Values.size());
412 
413  mtx_->getLocalRowCopy(LocalRow, indices, values, NumEntries);
414  for (size_t i = 0; i < NumEntries; ++i) {
415  Indices[i] = indices(i);
416  Values[i] = values(i);
417  }
418 }
419 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
420 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
421  getLocalRowView(LocalOrdinal LocalRow, ArrayView<const LocalOrdinal> &Indices,
422  ArrayView<const Scalar> &Values) const {
423  XPETRA_MONITOR("TpetraBlockCrsMatrix::getLocalRowView");
424  typename Tpetra::CrsGraph<LocalOrdinal, GlobalOrdinal, Node>::local_inds_host_view_type indices;
425  typename Tpetra::CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::values_host_view_type values;
426 
427  mtx_->getLocalRowView(LocalRow, indices, values);
428  Indices = ArrayView<const LocalOrdinal>(indices.data(), indices.extent(0));
429  Values = ArrayView<const Scalar>(reinterpret_cast<const Scalar *>(values.data()), values.extent(0));
430 }
431 
432 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
433 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
434  getGlobalRowView(GlobalOrdinal GlobalRow,
435  ArrayView<const GlobalOrdinal> &Indices,
436  ArrayView<const Scalar> &Values) const {
437  XPETRA_MONITOR("TpetraBlockCrsMatrix::getGlobalRowView");
438  typename Tpetra::CrsGraph<LocalOrdinal, GlobalOrdinal, Node>::global_inds_host_view_type indices;
439  typename Tpetra::CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::values_host_view_type values;
440 
441  mtx_->getGlobalRowView(GlobalRow, indices, values);
442  Indices = ArrayView<const GlobalOrdinal>(indices.data(), indices.extent(0));
443  Values = ArrayView<const Scalar>(reinterpret_cast<const Scalar *>(values.data()), values.extent(0));
444 }
445 
446 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
447 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
448  getGlobalRowCopy(GlobalOrdinal GlobalRow,
449  const ArrayView<GlobalOrdinal> &Indices,
450  const ArrayView<Scalar> &Values,
451  size_t &NumEntries) const {
452  XPETRA_MONITOR("TpetraBlockCrsMatrix::getGlobalRowCopy");
453  typename Tpetra::CrsGraph<LocalOrdinal, GlobalOrdinal, Node>::nonconst_global_inds_host_view_type indices("indices", Indices.size());
454  typename Tpetra::CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::nonconst_values_host_view_type values("values", Values.size());
455 
456  mtx_->getGlobalRowCopy(GlobalRow, indices, values, NumEntries);
457  for (size_t i = 0; i < NumEntries; ++i) {
458  Indices[i] = indices(i);
459  Values[i] = values(i);
460  }
461 }
462 
463 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
464 bool TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
465  haveGlobalConstants() const { return true; }
466 
467 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
468 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
469  apply(const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &X,
470  MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &Y,
471  Teuchos::ETransp mode,
472  Scalar alpha,
473  Scalar beta) const {
474  XPETRA_MONITOR("TpetraBlockCrsMatrix::apply");
475  mtx_->apply(toTpetra(X), toTpetra(Y), mode, alpha, beta);
476 }
477 
478 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
479 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
480  apply(const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &X,
481  MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &Y,
482  Teuchos::ETransp mode,
483  Scalar alpha,
484  Scalar beta,
485  bool sumInterfaceValues,
486  const RCP<Import<LocalOrdinal, GlobalOrdinal, Node> > &regionInterfaceImporter,
487  const Teuchos::ArrayRCP<LocalOrdinal> &regionInterfaceLIDs) const {}
488 
489 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
490 const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >
491 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
492  getDomainMap() const {
493  XPETRA_MONITOR("TpetraBlockCrsMatrix::getDomainMap");
494  return toXpetra(mtx_->getDomainMap());
495 }
496 
497 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
498 const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >
499 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
500  getRangeMap() const {
501  XPETRA_MONITOR("TpetraBlockCrsMatrix::getRangeMap");
502  return toXpetra(mtx_->getRangeMap());
503 }
504 
505 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
506 std::string
507 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
508  description() const {
509  XPETRA_MONITOR("TpetraBlockCrsMatrix::description");
510  return mtx_->description();
511 }
512 
513 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
514 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
515  describe(Teuchos::FancyOStream &out,
516  const Teuchos::EVerbosityLevel verbLevel) const {
517  XPETRA_MONITOR("TpetraBlockCrsMatrix::describe");
518  mtx_->describe(out, verbLevel);
519 }
520 
521 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
522 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
523  setObjectLabel(const std::string &objectLabel) {
524  XPETRA_MONITOR("TpetraCrsMatrix::setObjectLabel");
525  Teuchos::LabeledObject::setObjectLabel(objectLabel);
526  mtx_->setObjectLabel(objectLabel);
527 }
528 
529 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
530 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
531  getLocalDiagCopy(Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &diag) const {
532  XPETRA_MONITOR("TpetraBlockCrsMatrix::getLocalDiagCopy");
533  XPETRA_DYNAMIC_CAST(TpetraVectorClass,
534  diag,
535  tDiag,
536  "Xpetra::TpetraBlockCrsMatrix.getLocalDiagCopy() only accept Xpetra::TpetraVector as input arguments.");
537  mtx_->getLocalDiagCopy(*tDiag.getTpetra_Vector());
538 }
539 
541 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
542 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
543  getLocalDiagCopy(Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &diag,
544  const Teuchos::ArrayView<const size_t> &offsets) const {
545  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in " + std::string(__FILE__) + ":" + std::to_string(__LINE__));
546 }
547 
549 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
550 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
551  getLocalDiagCopy(Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &diag,
552  const Kokkos::View<const size_t *, typename Node::device_type, Kokkos::MemoryUnmanaged> &offsets) const {
553  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in " + std::string(__FILE__) + ":" + std::to_string(__LINE__));
554 }
555 
556 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
557 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
558  getLocalDiagOffsets(Teuchos::ArrayRCP<size_t> &offsets) const {
559  XPETRA_MONITOR("TpetraBlockCrsMatrix::getLocalDiagOffsets");
560 
561  const size_t lclNumRows = mtx_->getGraph()->getLocalNumRows();
562  if (static_cast<size_t>(offsets.size()) < lclNumRows) {
563  offsets.resize(lclNumRows);
564  }
565 
566  // The input ArrayRCP must always be a host pointer. Thus, if
567  // device_type::memory_space is Kokkos::HostSpace, it's OK for us
568  // to write to that allocation directly as a Kokkos::View.
569  typedef typename Node::device_type device_type;
570  typedef typename device_type::memory_space memory_space;
571  if (std::is_same<memory_space, Kokkos::HostSpace>::value) {
572  // It is always syntactically correct to assign a raw host
573  // pointer to a device View, so this code will compile correctly
574  // even if this branch never runs.
575  typedef Kokkos::View<size_t *, device_type, Kokkos::MemoryUnmanaged> output_type;
576  output_type offsetsOut(offsets.getRawPtr(), offsets.size());
577  mtx_->getLocalDiagOffsets(offsetsOut);
578  } else {
579  Kokkos::View<size_t *, device_type> offsetsTmp("diagOffsets", offsets.size());
580  mtx_->getLocalDiagOffsets(offsetsTmp);
581  typedef Kokkos::View<size_t *, Kokkos::HostSpace, Kokkos::MemoryUnmanaged> output_type;
582  output_type offsetsOut(offsets.getRawPtr(), offsets.size());
583  Kokkos::deep_copy(offsetsOut, offsetsTmp);
584  }
585 }
586 
587 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
588 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
589  replaceDiag(const Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &diag) {
590  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix::replaceDiag: function not implemented in " + std::string(__FILE__) + ":" + std::to_string(__LINE__));
591 }
592 
593 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
594 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
595  leftScale(const Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &x) {
596  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in " + std::string(__FILE__) + ":" + std::to_string(__LINE__));
597 }
598 
599 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
600 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
601  rightScale(const Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &x) {
602  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in " + std::string(__FILE__) + ":" + std::to_string(__LINE__));
603 }
604 
605 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
606 Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >
607 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
608  getMap() const {
609  XPETRA_MONITOR("TpetraBlockCrsMatrix::getMap");
610  return rcp(new TpetraMap<LocalOrdinal, GlobalOrdinal, Node>(mtx_->getMap()));
611 }
612 
614 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
615 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
616  doImport(const DistObject<char, LocalOrdinal, GlobalOrdinal, Node> &source,
617  const Import<LocalOrdinal, GlobalOrdinal, Node> &importer, CombineMode CM) {
618  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in " + std::string(__FILE__) + ":" + std::to_string(__LINE__));
619 }
620 
622 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
623 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
624  doExport(const DistObject<char, LocalOrdinal, GlobalOrdinal, Node> &dest,
625  const Import<LocalOrdinal, GlobalOrdinal, Node> &importer, CombineMode CM) {
626  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in " + std::string(__FILE__) + ":" + std::to_string(__LINE__));
627 }
628 
630 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
631 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
632  doImport(const DistObject<char, LocalOrdinal, GlobalOrdinal, Node> &source,
633  const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter, CombineMode CM) {
634  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in " + std::string(__FILE__) + ":" + std::to_string(__LINE__));
635 }
636 
638 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
639 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
640  doExport(const DistObject<char, LocalOrdinal, GlobalOrdinal, Node> &dest,
641  const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter, CombineMode CM) {
642  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in " + std::string(__FILE__) + ":" + std::to_string(__LINE__));
643 }
644 
645 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
646 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
647  removeEmptyProcessesInPlace(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &newMap) {
648  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix function not implemented in " + std::string(__FILE__) + ":" + std::to_string(__LINE__));
649 }
650 
651 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
652 bool TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
653  hasMatrix() const {
654  return !mtx_.is_null();
655 }
656 
657 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
658 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
659  TpetraBlockCrsMatrix(const Teuchos::RCP<Tpetra::BlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > &mtx)
660  : mtx_(mtx) {}
661 
662 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
663 RCP<const Tpetra::BlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> >
664 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
665  getTpetra_BlockCrsMatrix() const {
666  return mtx_;
667 }
668 
669 // TODO: remove
670 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
671 RCP<Tpetra::BlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> >
672 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
673  getTpetra_BlockCrsMatrixNonConst() const {
674  return mtx_;
675 }
676 
677 #ifdef HAVE_XPETRA_TPETRA
678 
679 // was: typedef typename Xpetra::CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::local_matrix_type local_matrix_type;
680 // using local_matrix_type = typename CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::local_matrix_type;
681 
682 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
683 typename CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::local_matrix_type
684 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
685  getLocalMatrixDevice() const {
686  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix does not support getLocalMatrix due to missing Kokkos::CrsMatrix in Tpetra's experimental implementation in " + std::string(__FILE__) + ":" + std::to_string(__LINE__));
687 
688 #ifndef __NVCC__
689  local_matrix_type ret;
690 #endif // __NVCC__
691 
692  TEUCHOS_UNREACHABLE_RETURN(ret);
693 }
694 
695 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
696 typename CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::local_matrix_type::HostMirror
697 TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
698  getLocalMatrixHost() const {
699  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix does not support getLocalMatrix due to missing Kokkos::CrsMatrix in Tpetra's experimental implementation in " + std::string(__FILE__) + ":" + std::to_string(__LINE__));
700 
701 #ifndef __NVCC__
702  typename local_matrix_type::HostMirror ret;
703 #endif // __NVCC__
704 
705  TEUCHOS_UNREACHABLE_RETURN(ret);
706 }
707 
708 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
709 void TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
710  setAllValues(const typename local_matrix_type::row_map_type &ptr,
711  const typename local_matrix_type::StaticCrsGraphType::entries_type::non_const_type &ind,
712  const typename local_matrix_type::values_type &val) {
713  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix does not support setAllValues due to missing Kokkos::CrsMatrix in Tpetra's experimental implementation in " + std::string(__FILE__) + ":" + std::to_string(__LINE__));
714 }
715 
716 #endif // HAVE_XPETRA_TPETRA
717 
718 #ifdef HAVE_XPETRA_EPETRA
719 
720 #if ((defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_OPENMP) || !defined(HAVE_TPETRA_INST_INT_INT))) || \
721  (!defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_SERIAL) || !defined(HAVE_TPETRA_INST_INT_INT))))
722 
723 // specialization of TpetraBlockCrsMatrix for GO=LO=int and Node=EpetraNode
724 template <class Scalar>
725 class TpetraBlockCrsMatrix<Scalar, int, int, EpetraNode>
726  : public CrsMatrix<Scalar, int, int, EpetraNode> //, public TpetraRowMatrix<Scalar,int,int,Node>
727 {
728  // The following typedef are used by the XPETRA_DYNAMIC_CAST() macro.
729  typedef int LocalOrdinal;
730  typedef int GlobalOrdinal;
731  typedef EpetraNode Node;
732  typedef TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> TpetraBlockCrsMatrixClass;
733  typedef TpetraVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> TpetraVectorClass;
734  typedef TpetraImport<LocalOrdinal, GlobalOrdinal, Node> TpetraImportClass;
735  typedef TpetraExport<LocalOrdinal, GlobalOrdinal, Node> TpetraExportClass;
736 
737  public:
739 
741  TpetraBlockCrsMatrix(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rowMap, size_t maxNumEntriesPerRow, const Teuchos::RCP<Teuchos::ParameterList> &params = Teuchos::null) {
742  XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "int", typeid(EpetraNode).name());
743  }
744 
746  TpetraBlockCrsMatrix(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rowMap, const ArrayRCP<const size_t> &NumEntriesPerRowToAlloc, const Teuchos::RCP<Teuchos::ParameterList> &params = Teuchos::null) {
747  XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "int", typeid(EpetraNode).name());
748  }
749 
751  TpetraBlockCrsMatrix(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rowMap, const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &colMap, size_t maxNumEntriesPerRow, const Teuchos::RCP<Teuchos::ParameterList> &params = Teuchos::null) {
752  XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "int", typeid(EpetraNode).name());
753  }
754 
756  TpetraBlockCrsMatrix(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rowMap, const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &colMap, const ArrayRCP<const size_t> &NumEntriesPerRowToAlloc, const Teuchos::RCP<Teuchos::ParameterList> &params = Teuchos::null) {
757  XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "int", typeid(EpetraNode).name());
758  }
759 
761  TpetraBlockCrsMatrix(const Teuchos::RCP<const CrsGraph<LocalOrdinal, GlobalOrdinal, Node> > &graph, const Teuchos::RCP<Teuchos::ParameterList> &params = Teuchos::null) {
762  XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "int", typeid(EpetraNode).name());
763  }
764 
766  TpetraBlockCrsMatrix(const Teuchos::RCP<const CrsGraph<LocalOrdinal, GlobalOrdinal, Node> > &graph, const LocalOrdinal blockSize) {
767  XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "int", typeid(EpetraNode).name());
768  }
769 
771  TpetraBlockCrsMatrix(const Teuchos::RCP<const Tpetra::BlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > &sourceMatrix,
772  const Import<LocalOrdinal, GlobalOrdinal, Node> &importer,
773  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &domainMap = Teuchos::null,
774  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rangeMap = Teuchos::null,
775  const Teuchos::RCP<Teuchos::ParameterList> &params = Teuchos::null) { XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "int", typeid(EpetraNode).name()); }
776 
778  TpetraBlockCrsMatrix(const Teuchos::RCP<const Tpetra::BlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > &sourceMatrix,
779  const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter,
780  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &domainMap = Teuchos::null,
781  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rangeMap = Teuchos::null,
782  const Teuchos::RCP<Teuchos::ParameterList> &params = Teuchos::null) { XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "int", typeid(EpetraNode).name()); }
783 
785  TpetraBlockCrsMatrix(const Teuchos::RCP<const Tpetra::BlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > &sourceMatrix,
786  const Import<LocalOrdinal, GlobalOrdinal, Node> &RowImporter,
787  const Teuchos::RCP<const Import<LocalOrdinal, GlobalOrdinal, Node> > DomainImporter,
788  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &domainMap,
789  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rangeMap,
790  const Teuchos::RCP<Teuchos::ParameterList> &params) { XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "int", typeid(EpetraNode).name()); }
791 
793  TpetraBlockCrsMatrix(const Teuchos::RCP<const Tpetra::BlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > &sourceMatrix,
794  const Export<LocalOrdinal, GlobalOrdinal, Node> &RowExporter,
795  const Teuchos::RCP<const Export<LocalOrdinal, GlobalOrdinal, Node> > DomainExporter,
796  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &domainMap,
797  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rangeMap,
798  const Teuchos::RCP<Teuchos::ParameterList> &params) {
799  XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "int", typeid(EpetraNode).name());
800  }
801 
803  ~TpetraBlockCrsMatrix() {}
804 
806 
808  void insertGlobalValues(GlobalOrdinal globalRow, const ArrayView<const GlobalOrdinal> &cols, const ArrayView<const Scalar> &vals) {}
809 
811  void insertLocalValues(LocalOrdinal localRow, const ArrayView<const LocalOrdinal> &cols, const ArrayView<const Scalar> &vals) {}
812 
814  void replaceGlobalValues(GlobalOrdinal globalRow, const ArrayView<const GlobalOrdinal> &cols, const ArrayView<const Scalar> &vals) {}
815 
817  void replaceLocalValues(LocalOrdinal localRow, const ArrayView<const LocalOrdinal> &cols, const ArrayView<const Scalar> &vals) {}
818 
820  void setAllToScalar(const Scalar &alpha) {}
821 
823  void scale(const Scalar &alpha) {}
824 
826  //** \warning This is an expert-only routine and should not be called from user code. (not implemented)
827  void allocateAllValues(size_t numNonZeros, ArrayRCP<size_t> &rowptr, ArrayRCP<LocalOrdinal> &colind, ArrayRCP<Scalar> &values) {}
828 
830  void setAllValues(const ArrayRCP<size_t> &rowptr, const ArrayRCP<LocalOrdinal> &colind, const ArrayRCP<Scalar> &values) {}
831 
833  void getAllValues(ArrayRCP<const size_t> &rowptr, ArrayRCP<const LocalOrdinal> &colind, ArrayRCP<const Scalar> &values) const {}
834 
836  void getAllValues(ArrayRCP<Scalar> &values) {}
837 
839 
841  void resumeFill(const RCP<ParameterList> &params = null) { /*noop*/
842  }
843 
845  void fillComplete(const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &domainMap, const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rangeMap, const RCP<ParameterList> &params = null) { /*noop*/
846  }
847 
849  void fillComplete(const RCP<ParameterList> &params = null) { /*noop*/
850  }
851 
853  void replaceDomainMapAndImporter(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &newDomainMap, Teuchos::RCP<const Import<LocalOrdinal, GlobalOrdinal, Node> > &newImporter) {}
854 
856  void expertStaticFillComplete(const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &domainMap,
857  const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rangeMap,
858  const RCP<const Import<LocalOrdinal, GlobalOrdinal, Node> > &importer = Teuchos::null,
859  const RCP<const Export<LocalOrdinal, GlobalOrdinal, Node> > &exporter = Teuchos::null,
860  const RCP<ParameterList> &params = Teuchos::null) {}
861 
863 
865  const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > getRowMap() const { return Teuchos::null; }
866 
868  const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > getColMap() const { return Teuchos::null; }
869 
871  RCP<const CrsGraph<LocalOrdinal, GlobalOrdinal, Node> > getCrsGraph() const { return Teuchos::null; }
872 
874  global_size_t getGlobalNumRows() const { return 0; }
875 
877  global_size_t getGlobalNumCols() const { return 0; }
878 
880  size_t getLocalNumRows() const { return 0; }
881 
883  size_t getLocalNumCols() const { return 0; }
884 
886  global_size_t getGlobalNumEntries() const { return 0; }
887 
889  size_t getLocalNumEntries() const { return 0; }
890 
892  size_t getNumEntriesInLocalRow(LocalOrdinal localRow) const { return 0; }
893 
895  size_t getNumEntriesInGlobalRow(GlobalOrdinal globalRow) const { return 0; }
896 
898  size_t getGlobalMaxNumRowEntries() const { return 0; }
899 
901  size_t getLocalMaxNumRowEntries() const { return 0; }
902 
904  bool isLocallyIndexed() const { return false; }
905 
907  bool isGloballyIndexed() const { return false; }
908 
910  bool isFillComplete() const { return false; }
911 
913  bool isFillActive() const { return false; }
914 
916  typename ScalarTraits<Scalar>::magnitudeType getFrobeniusNorm() const { return Teuchos::ScalarTraits<Scalar>::magnitude(Teuchos::ScalarTraits<Scalar>::zero()); }
917 
919  bool supportsRowViews() const { return false; }
920 
922  void getLocalRowCopy(LocalOrdinal LocalRow, const ArrayView<LocalOrdinal> &Indices, const ArrayView<Scalar> &Values, size_t &NumEntries) const {}
923 
925  void getGlobalRowView(GlobalOrdinal GlobalRow, ArrayView<const GlobalOrdinal> &indices, ArrayView<const Scalar> &values) const {}
926 
928  void getGlobalRowCopy(GlobalOrdinal GlobalRow, const ArrayView<GlobalOrdinal> &indices, const ArrayView<Scalar> &values, size_t &numEntries) const {}
929 
931  void getLocalRowView(LocalOrdinal LocalRow, ArrayView<const LocalOrdinal> &indices, ArrayView<const Scalar> &values) const {}
932 
934  bool haveGlobalConstants() const { return false; }
935 
937 
939  void apply(const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &X, MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &Y, Teuchos::ETransp mode = Teuchos::NO_TRANS, Scalar alpha = ScalarTraits<Scalar>::one(), Scalar beta = ScalarTraits<Scalar>::zero()) const {}
940 
942  const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > getDomainMap() const { return Teuchos::null; }
943 
945  const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > getRangeMap() const { return Teuchos::null; }
946 
948 
950  std::string description() const { return std::string(""); }
951 
953  void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel = Teuchos::Describable::verbLevel_default) const {}
954 
956  TpetraBlockCrsMatrix(const TpetraBlockCrsMatrix &matrix) {}
957 
959  void getLocalDiagCopy(Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &diag) const {}
960 
962  void getLocalDiagOffsets(Teuchos::ArrayRCP<size_t> &offsets) const {}
963 
965  void getLocalDiagCopy(Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &diag, const Teuchos::ArrayView<const size_t> &offsets) const {}
966 
968  void getLocalDiagCopy(Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &diag, const Kokkos::View<const size_t *, typename Node::device_type, Kokkos::MemoryUnmanaged> &offsets) const {}
969 
970  void replaceDiag(const Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &diag) {}
971 
972  void leftScale(const Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &x) {}
973  void rightScale(const Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &x) {}
974 
976 
978  Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > getMap() const { return Teuchos::null; }
979 
981  void doImport(const DistObject<char, LocalOrdinal, GlobalOrdinal, Node> &source,
982  const Import<LocalOrdinal, GlobalOrdinal, Node> &importer, CombineMode CM) {}
983 
985  void doExport(const DistObject<char, LocalOrdinal, GlobalOrdinal, Node> &dest,
986  const Import<LocalOrdinal, GlobalOrdinal, Node> &importer, CombineMode CM) {}
987 
989  void doImport(const DistObject<char, LocalOrdinal, GlobalOrdinal, Node> &source,
990  const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter, CombineMode CM) {}
991 
993  void doExport(const DistObject<char, LocalOrdinal, GlobalOrdinal, Node> &dest,
994  const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter, CombineMode CM) {}
995 
996  void removeEmptyProcessesInPlace(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &newMap) {}
997 
999 
1001  bool hasMatrix() const { return false; }
1002 
1004  TpetraBlockCrsMatrix(const Teuchos::RCP<Tpetra::BlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > &mtx) {
1005  XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "int", typeid(EpetraNode).name());
1006  }
1007 
1009  RCP<const Tpetra::BlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > getTpetra_BlockCrsMatrix() const { return Teuchos::null; }
1010 
1012  RCP<Tpetra::BlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > getTpetra_BlockCrsMatrixNonConst() const { return Teuchos::null; }
1013 
1014 #ifdef HAVE_XPETRA_TPETRA
1015  typedef typename Xpetra::CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::local_matrix_type local_matrix_type;
1016 
1017  local_matrix_type getLocalMatrix() const {
1018  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix does not support getLocalMatrix due to missing Kokkos::CrsMatrix in Tpetra's experimental implementation in " + std::string(__FILE__) + ":" + std::to_string(__LINE__));
1019  local_matrix_type ret;
1020  return ret; // make compiler happy
1021  }
1022 
1023  void setAllValues(const typename local_matrix_type::row_map_type &ptr,
1024  const typename local_matrix_type::StaticCrsGraphType::entries_type::non_const_type &ind,
1025  const typename local_matrix_type::values_type &val) {
1026  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix does not support setAllValues due to missing Kokkos::CrsMatrix in Tpetra's experimental implementation in " + std::string(__FILE__) + ":" + std::to_string(__LINE__));
1027  }
1028 #endif // HAVE_XPETRA_TPETRA
1029 
1030 }; // TpetraBlockCrsMatrix class
1031 
1032 #endif // #if ((defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_OPENMP) || !defined(HAVE_TPETRA_INST_INT_INT)))
1033 
1034 #if ((defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_OPENMP) || !defined(HAVE_TPETRA_INST_INT_LONG_LONG))) || \
1035  (!defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_SERIAL) || !defined(HAVE_TPETRA_INST_INT_LONG_LONG))))
1036 
1037 // specialization of TpetraBlockCrsMatrix for GO=long long and Node=EpetraNode
1038 template <class Scalar>
1039 class TpetraBlockCrsMatrix<Scalar, int, long long, EpetraNode>
1040  : public CrsMatrix<Scalar, int, long long, EpetraNode> //, public TpetraRowMatrix<Scalar,int,int,Node>
1041 {
1042  // The following typedef are used by the XPETRA_DYNAMIC_CAST() macro.
1043  typedef int LocalOrdinal;
1044  typedef long long GlobalOrdinal;
1045  typedef EpetraNode Node;
1046  typedef TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> TpetraBlockCrsMatrixClass;
1047  typedef TpetraVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> TpetraVectorClass;
1048  typedef TpetraImport<LocalOrdinal, GlobalOrdinal, Node> TpetraImportClass;
1049  typedef TpetraExport<LocalOrdinal, GlobalOrdinal, Node> TpetraExportClass;
1050 
1051  public:
1053 
1055  TpetraBlockCrsMatrix(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rowMap, size_t maxNumEntriesPerRow, const Teuchos::RCP<Teuchos::ParameterList> &params = Teuchos::null) { XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "long long", typeid(EpetraNode).name()); }
1056 
1058  TpetraBlockCrsMatrix(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rowMap, const ArrayRCP<const size_t> &NumEntriesPerRowToAlloc, const Teuchos::RCP<Teuchos::ParameterList> &params = Teuchos::null) { XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "long long", typeid(EpetraNode).name()); }
1059 
1061  TpetraBlockCrsMatrix(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rowMap, const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &colMap, size_t maxNumEntriesPerRow, const Teuchos::RCP<Teuchos::ParameterList> &params = Teuchos::null) { XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "long long", typeid(EpetraNode).name()); }
1062 
1064  TpetraBlockCrsMatrix(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rowMap, const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &colMap, const ArrayRCP<const size_t> &NumEntriesPerRowToAlloc, const Teuchos::RCP<Teuchos::ParameterList> &params = Teuchos::null) { XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "long long", typeid(EpetraNode).name()); }
1065 
1067  TpetraBlockCrsMatrix(const Teuchos::RCP<const CrsGraph<LocalOrdinal, GlobalOrdinal, Node> > &graph, const Teuchos::RCP<Teuchos::ParameterList> &params = Teuchos::null) { XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "long long", typeid(EpetraNode).name()); }
1068 
1070  TpetraBlockCrsMatrix(const Teuchos::RCP<const CrsGraph<LocalOrdinal, GlobalOrdinal, Node> > &graph, const LocalOrdinal blockSize) { XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "long long", typeid(EpetraNode).name()); }
1071 
1073  TpetraBlockCrsMatrix(const Teuchos::RCP<const Tpetra::BlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > &sourceMatrix,
1074  const Import<LocalOrdinal, GlobalOrdinal, Node> &importer,
1075  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &domainMap = Teuchos::null,
1076  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rangeMap = Teuchos::null,
1077  const Teuchos::RCP<Teuchos::ParameterList> &params = Teuchos::null) { XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "long long", typeid(EpetraNode).name()); }
1078 
1080  TpetraBlockCrsMatrix(const Teuchos::RCP<const Tpetra::BlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > &sourceMatrix,
1081  const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter,
1082  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &domainMap = Teuchos::null,
1083  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rangeMap = Teuchos::null,
1084  const Teuchos::RCP<Teuchos::ParameterList> &params = Teuchos::null) { XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "long long", typeid(EpetraNode).name()); }
1085 
1087  TpetraBlockCrsMatrix(const Teuchos::RCP<const Tpetra::BlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > &sourceMatrix,
1088  const Import<LocalOrdinal, GlobalOrdinal, Node> &RowImporter,
1089  const Teuchos::RCP<const Import<LocalOrdinal, GlobalOrdinal, Node> > DomainImporter,
1090  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &domainMap,
1091  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rangeMap,
1092  const Teuchos::RCP<Teuchos::ParameterList> &params) { XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "long long", typeid(EpetraNode).name()); }
1093 
1095  TpetraBlockCrsMatrix(const Teuchos::RCP<const Tpetra::BlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > &sourceMatrix,
1096  const Export<LocalOrdinal, GlobalOrdinal, Node> &RowExporter,
1097  const Teuchos::RCP<const Export<LocalOrdinal, GlobalOrdinal, Node> > DomainExporter,
1098  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &domainMap,
1099  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rangeMap,
1100  const Teuchos::RCP<Teuchos::ParameterList> &params) { XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "long long", typeid(EpetraNode).name()); }
1101 
1103  ~TpetraBlockCrsMatrix() {}
1104 
1106 
1108  void insertGlobalValues(GlobalOrdinal globalRow, const ArrayView<const GlobalOrdinal> &cols, const ArrayView<const Scalar> &vals) {}
1109 
1111  void insertLocalValues(LocalOrdinal localRow, const ArrayView<const LocalOrdinal> &cols, const ArrayView<const Scalar> &vals) {}
1112 
1114  void replaceGlobalValues(GlobalOrdinal globalRow, const ArrayView<const GlobalOrdinal> &cols, const ArrayView<const Scalar> &vals) {}
1115 
1117  void replaceLocalValues(LocalOrdinal localRow, const ArrayView<const LocalOrdinal> &cols, const ArrayView<const Scalar> &vals) {}
1118 
1120  void setAllToScalar(const Scalar &alpha) {}
1121 
1123  void scale(const Scalar &alpha) {}
1124 
1126  //** \warning This is an expert-only routine and should not be called from user code. (not implemented)
1127  void allocateAllValues(size_t numNonZeros, ArrayRCP<size_t> &rowptr, ArrayRCP<LocalOrdinal> &colind, ArrayRCP<Scalar> &values) {}
1128 
1130  void setAllValues(const ArrayRCP<size_t> &rowptr, const ArrayRCP<LocalOrdinal> &colind, const ArrayRCP<Scalar> &values) {}
1131 
1133  void getAllValues(ArrayRCP<const size_t> &rowptr, ArrayRCP<const LocalOrdinal> &colind, ArrayRCP<const Scalar> &values) const {}
1134 
1136  void getAllValues(ArrayRCP<Scalar> &values) {}
1137 
1139 
1141  void resumeFill(const RCP<ParameterList> &params = null) { /*noop*/
1142  }
1143 
1145  void fillComplete(const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &domainMap, const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rangeMap, const RCP<ParameterList> &params = null) { /*noop*/
1146  }
1147 
1149  void fillComplete(const RCP<ParameterList> &params = null) { /*noop*/
1150  }
1151 
1153  void replaceDomainMapAndImporter(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &newDomainMap, Teuchos::RCP<const Import<LocalOrdinal, GlobalOrdinal, Node> > &newImporter) {}
1154 
1156  void expertStaticFillComplete(const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &domainMap,
1157  const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rangeMap,
1158  const RCP<const Import<LocalOrdinal, GlobalOrdinal, Node> > &importer = Teuchos::null,
1159  const RCP<const Export<LocalOrdinal, GlobalOrdinal, Node> > &exporter = Teuchos::null,
1160  const RCP<ParameterList> &params = Teuchos::null) {}
1161 
1163 
1165  const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > getRowMap() const { return Teuchos::null; }
1166 
1168  const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > getColMap() const { return Teuchos::null; }
1169 
1171  RCP<const CrsGraph<LocalOrdinal, GlobalOrdinal, Node> > getCrsGraph() const { return Teuchos::null; }
1172 
1174  global_size_t getGlobalNumRows() const { return 0; }
1175 
1177  global_size_t getGlobalNumCols() const { return 0; }
1178 
1180  size_t getLocalNumRows() const { return 0; }
1181 
1183  size_t getLocalNumCols() const { return 0; }
1184 
1186  global_size_t getGlobalNumEntries() const { return 0; }
1187 
1189  size_t getLocalNumEntries() const { return 0; }
1190 
1192  size_t getNumEntriesInLocalRow(LocalOrdinal localRow) const { return 0; }
1193 
1195  size_t getNumEntriesInGlobalRow(GlobalOrdinal globalRow) const { return 0; }
1196 
1198  size_t getGlobalMaxNumRowEntries() const { return 0; }
1199 
1201  size_t getLocalMaxNumRowEntries() const { return 0; }
1202 
1204  bool isLocallyIndexed() const { return false; }
1205 
1207  bool isGloballyIndexed() const { return false; }
1208 
1210  bool isFillComplete() const { return false; }
1211 
1213  bool isFillActive() const { return false; }
1214 
1216  typename ScalarTraits<Scalar>::magnitudeType getFrobeniusNorm() const { return Teuchos::ScalarTraits<Scalar>::magnitude(Teuchos::ScalarTraits<Scalar>::zero()); }
1217 
1219  bool supportsRowViews() const { return false; }
1220 
1222  void getLocalRowCopy(LocalOrdinal LocalRow, const ArrayView<LocalOrdinal> &Indices, const ArrayView<Scalar> &Values, size_t &NumEntries) const {}
1223 
1225  void getGlobalRowView(GlobalOrdinal GlobalRow, ArrayView<const GlobalOrdinal> &indices, ArrayView<const Scalar> &values) const {}
1226 
1228  void getGlobalRowCopy(GlobalOrdinal GlobalRow, const ArrayView<GlobalOrdinal> &indices, const ArrayView<Scalar> &values, size_t &numEntries) const {}
1229 
1231  void getLocalRowView(LocalOrdinal LocalRow, ArrayView<const LocalOrdinal> &indices, ArrayView<const Scalar> &values) const {}
1232 
1234  bool haveGlobalConstants() const { return true; }
1235 
1237 
1239  void apply(const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &X, MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &Y, Teuchos::ETransp mode = Teuchos::NO_TRANS, Scalar alpha = ScalarTraits<Scalar>::one(), Scalar beta = ScalarTraits<Scalar>::zero()) const {}
1240 
1242  const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > getDomainMap() const { return Teuchos::null; }
1243 
1245  const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > getRangeMap() const { return Teuchos::null; }
1246 
1248 
1250  std::string description() const { return std::string(""); }
1251 
1253  void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel = Teuchos::Describable::verbLevel_default) const {}
1254 
1256  TpetraBlockCrsMatrix(const TpetraBlockCrsMatrix &matrix) {}
1257 
1259  void getLocalDiagCopy(Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &diag) const {}
1260 
1262  void getLocalDiagOffsets(Teuchos::ArrayRCP<size_t> &offsets) const {}
1263 
1265  void getLocalDiagCopy(Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &diag, const Teuchos::ArrayView<const size_t> &offsets) const {}
1266 
1268  void getLocalDiagCopy(Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &diag, const Kokkos::View<const size_t *, typename Node::device_type, Kokkos::MemoryUnmanaged> &offsets) const {}
1269 
1270  void replaceDiag(Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &diag) const {}
1271 
1272  void leftScale(const Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &x) {}
1273  void rightScale(const Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &x) {}
1274 
1276 
1278  Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > getMap() const { return Teuchos::null; }
1279 
1281  void doImport(const DistObject<char, LocalOrdinal, GlobalOrdinal, Node> &source,
1282  const Import<LocalOrdinal, GlobalOrdinal, Node> &importer, CombineMode CM) {}
1283 
1285  void doExport(const DistObject<char, LocalOrdinal, GlobalOrdinal, Node> &dest,
1286  const Import<LocalOrdinal, GlobalOrdinal, Node> &importer, CombineMode CM) {}
1287 
1289  void doImport(const DistObject<char, LocalOrdinal, GlobalOrdinal, Node> &source,
1290  const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter, CombineMode CM) {}
1291 
1293  void doExport(const DistObject<char, LocalOrdinal, GlobalOrdinal, Node> &dest,
1294  const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter, CombineMode CM) {}
1295 
1296  void removeEmptyProcessesInPlace(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &newMap) {}
1297 
1299 
1301  bool hasMatrix() const { return false; }
1302 
1304  TpetraBlockCrsMatrix(const Teuchos::RCP<Tpetra::BlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > &mtx) {
1305  XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraBlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "long long", typeid(EpetraNode).name());
1306  }
1307 
1309  RCP<const Tpetra::BlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > getTpetra_BlockCrsMatrix() const { return Teuchos::null; }
1310 
1312  RCP<Tpetra::BlockCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > getTpetra_BlockCrsMatrixNonConst() const { return Teuchos::null; }
1313 
1314 #ifdef HAVE_XPETRA_TPETRA
1315  typedef typename Xpetra::CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::local_matrix_type local_matrix_type;
1316 
1317  local_matrix_type getLocalMatrix() const {
1318  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix does not support getLocalMatrix due to missing Kokkos::CrsMatrix in Tpetra's experimental implementation in " + std::string(__FILE__) + ":" + std::to_string(__LINE__));
1319  local_matrix_type ret;
1320  TEUCHOS_UNREACHABLE_RETURN(ret);
1321  }
1322 
1323  void setAllValues(const typename local_matrix_type::row_map_type &ptr,
1324  const typename local_matrix_type::StaticCrsGraphType::entries_type::non_const_type &ind,
1325  const typename local_matrix_type::values_type &val) {
1326  throw std::runtime_error("Xpetra::TpetraBlockCrsMatrix does not support setAllValues due to missing Kokkos::CrsMatrix in Tpetra's experimental implementation in " + std::string(__FILE__) + ":" + std::to_string(__LINE__));
1327  }
1328 #endif // HAVE_XPETRA_TPETRA
1329 
1330 }; // TpetraBlockCrsMatrix class
1331 
1332 #endif // IF ((defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_OPENMP) || !defined(HAVE_TPETRA_INST_INT_LONG_LONG)))
1333 
1334 #endif // HAVE_XPETRA_EPETRA
1335 
1336 } // namespace Xpetra
1337 
1338 #endif // XPETRA_TPETRABLOCKCRSMATRIX_DEF_HPP
Xpetra::TpetraBlockCrsMatrix::isGloballyIndexed
bool isGloballyIndexed() const
If matrix indices are in the global range, this function returns true. Otherwise, this function retur...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:374
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::doExport
void doExport(const DistObject< char, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Export (using an Importer).
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1293
Xpetra::TpetraBlockCrsMatrix::scale
void scale(const Scalar &alpha)
Scale the current values of a matrix, this = alpha*this (not implemented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:182
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::TpetraBlockCrsMatrix
TpetraBlockCrsMatrix(const TpetraBlockCrsMatrix &matrix)
Deep copy constructor.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1256
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::fillComplete
void fillComplete(const RCP< ParameterList > &params=null)
Signal that data entry is complete.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1149
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::replaceGlobalValues
void replaceGlobalValues(GlobalOrdinal globalRow, const ArrayView< const GlobalOrdinal > &cols, const ArrayView< const Scalar > &vals)
Replace matrix entries, using global IDs (not implemented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1114
Xpetra::TpetraBlockCrsMatrix::getLocalRowCopy
void getLocalRowCopy(LocalOrdinal LocalRow, const ArrayView< LocalOrdinal > &Indices, const ArrayView< Scalar > &Values, size_t &NumEntries) const
Extract a list of entries in a specified local row of the matrix. Put into storage allocated by calli...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:405
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::removeEmptyProcessesInPlace
void removeEmptyProcessesInPlace(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &newMap)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1296
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::replaceDiag
void replaceDiag(Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &diag) const
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1270
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getMap
Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > getMap() const
Implements DistObject interface.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:978
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::setAllToScalar
void setAllToScalar(const Scalar &alpha)
Set all matrix entries equal to scalarThis.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1120
Xpetra::TpetraBlockCrsMatrix::getGlobalRowCopy
void getGlobalRowCopy(GlobalOrdinal GlobalRow, const ArrayView< GlobalOrdinal > &indices, const ArrayView< Scalar > &values, size_t &numEntries) const
Extract a list of entries in a specified global row of this matrix. Put into pre-allocated storage...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:448
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::supportsRowViews
bool supportsRowViews() const
Returns true if getLocalRowView() and getGlobalRowView() are valid for this class.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1219
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getDomainMap
const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > getDomainMap() const
Returns the Map associated with the domain of this operator. This will be null until fillComplete() i...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1242
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::insertGlobalValues
void insertGlobalValues(GlobalOrdinal globalRow, const ArrayView< const GlobalOrdinal > &cols, const ArrayView< const Scalar > &vals)
Insert matrix entries, using global IDs (not implemented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:808
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::fillComplete
void fillComplete(const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &domainMap, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap, const RCP< ParameterList > &params=null)
Signal that data entry is complete, specifying domain and range maps.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1145
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getCrsGraph
RCP< const CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > getCrsGraph() const
Returns the CrsGraph associated with this matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:871
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getDomainMap
const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > getDomainMap() const
Returns the Map associated with the domain of this operator. This will be null until fillComplete() i...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:942
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getLocalRowCopy
void getLocalRowCopy(LocalOrdinal LocalRow, const ArrayView< LocalOrdinal > &Indices, const ArrayView< Scalar > &Values, size_t &NumEntries) const
Extract a list of entries in a specified local row of the matrix. Put into storage allocated by calli...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:922
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::rightScale
void rightScale(const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &x)
Right scale matrix using the given vector entries.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1273
Xpetra::TpetraBlockCrsMatrix::setAllToScalar
void setAllToScalar(const Scalar &alpha)
Set all matrix entries equal to scalarThis.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:174
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::TpetraExportClass
TpetraExport< LocalOrdinal, GlobalOrdinal, Node > TpetraExportClass
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:735
Xpetra::TpetraBlockCrsMatrix::allocateAllValues
void allocateAllValues(size_t numNonZeros, ArrayRCP< size_t > &rowptr, ArrayRCP< LocalOrdinal > &colind, ArrayRCP< Scalar > &values)
Allocates and returns ArrayRCPs of the Crs arrays — This is an Xpetra-only routine.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:190
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getLocalDiagCopy
void getLocalDiagCopy(Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &diag) const
Get a copy of the diagonal entries owned by this node, with local row idices.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:959
Xpetra::TpetraBlockCrsMatrix::getLocalMatrixHost
local_matrix_type::HostMirror getLocalMatrixHost() const
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:698
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::insertGlobalValues
void insertGlobalValues(GlobalOrdinal globalRow, const ArrayView< const GlobalOrdinal > &cols, const ArrayView< const Scalar > &vals)
Insert matrix entries, using global IDs (not implemented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1108
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::replaceDomainMapAndImporter
void replaceDomainMapAndImporter(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &newDomainMap, Teuchos::RCP< const Import< LocalOrdinal, GlobalOrdinal, Node > > &newImporter)
Replaces the current domainMap and importer with the user-specified objects.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:853
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getRowMap
const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > getRowMap() const
Returns the Map that describes the row distribution in this matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:865
Xpetra::TpetraBlockCrsMatrix::getCrsGraph
RCP< const CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > getCrsGraph() const
Returns the CrsGraph associated with this matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:282
Xpetra::TpetraBlockCrsMatrix::insertLocalValues
void insertLocalValues(LocalOrdinal localRow, const ArrayView< const LocalOrdinal > &cols, const ArrayView< const Scalar > &vals)
Insert matrix entries, using local IDs (not implemented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:148
Xpetra::TpetraBlockCrsMatrix::insertGlobalValues
void insertGlobalValues(GlobalOrdinal globalRow, const ArrayView< const GlobalOrdinal > &cols, const ArrayView< const Scalar > &vals)
Insert matrix entries, using global IDs (not implemented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:139
Xpetra::TpetraBlockCrsMatrix::getTpetra_BlockCrsMatrixNonConst
RCP< Tpetra::BlockCrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getTpetra_BlockCrsMatrixNonConst() const
Get the underlying Tpetra matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:673
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getLocalRowCopy
void getLocalRowCopy(LocalOrdinal LocalRow, const ArrayView< LocalOrdinal > &Indices, const ArrayView< Scalar > &Values, size_t &NumEntries) const
Extract a list of entries in a specified local row of the matrix. Put into storage allocated by calli...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1222
Xpetra::TpetraBlockCrsMatrix::replaceLocalValues
void replaceLocalValues(LocalOrdinal localRow, const ArrayView< const LocalOrdinal > &cols, const ArrayView< const Scalar > &vals)
Replace matrix entries, using local IDs.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:166
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getLocalDiagCopy
void getLocalDiagCopy(Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &diag) const
Get a copy of the diagonal entries owned by this node, with local row idices.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1259
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::expertStaticFillComplete
void expertStaticFillComplete(const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &domainMap, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap, const RCP< const Import< LocalOrdinal, GlobalOrdinal, Node > > &importer=Teuchos::null, const RCP< const Export< LocalOrdinal, GlobalOrdinal, Node > > &exporter=Teuchos::null, const RCP< ParameterList > &params=Teuchos::null)
Expert static fill complete.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1156
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::describe
void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel=Teuchos::Describable::verbLevel_default) const
Print the object with some verbosity level to an FancyOStream object.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1253
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::local_matrix_type
Xpetra::CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::local_matrix_type local_matrix_type
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1315
Xpetra::TpetraBlockCrsMatrix::rightScale
void rightScale(const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &x)
Right scale operator with given vector values.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:601
Xpetra::TpetraBlockCrsMatrix::isLocallyIndexed
bool isLocallyIndexed() const
If matrix indices are in the local range, this function returns true. Otherwise, this function return...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:368
Xpetra::TpetraBlockCrsMatrix::getGlobalNumCols
global_size_t getGlobalNumCols() const
Number of global columns in the matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:302
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::isGloballyIndexed
bool isGloballyIndexed() const
If matrix indices are in the global range, this function returns true. Otherwise, this function retur...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:907
Xpetra::TpetraBlockCrsMatrix::hasMatrix
bool hasMatrix() const
Does this have an underlying matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:653
Xpetra::TpetraBlockCrsMatrix::getLocalNumEntries
size_t getLocalNumEntries() const
Returns the local number of entries in this matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:334
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::doExport
void doExport(const DistObject< char, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Export.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1285
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getTpetra_BlockCrsMatrixNonConst
RCP< Tpetra::BlockCrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getTpetra_BlockCrsMatrixNonConst() const
Get the underlying Tpetra matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1012
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::TpetraBlockCrsMatrix
TpetraBlockCrsMatrix(const Teuchos::RCP< const Tpetra::BlockCrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &domainMap=Teuchos::null, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap=Teuchos::null, const Teuchos::RCP< Teuchos::ParameterList > &params=Teuchos::null)
Constructor for a fused import ( not implemented )
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:771
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::TpetraBlockCrsMatrix
TpetraBlockCrsMatrix(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap, size_t maxNumEntriesPerRow, const Teuchos::RCP< Teuchos::ParameterList > &params=Teuchos::null)
Constructor specifying column Map and fixed number of entries for each row (not implemented) ...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1061
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::leftScale
void leftScale(const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &x)
Left scale matrix using the given vector entries.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:972
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getLocalDiagOffsets
void getLocalDiagOffsets(Teuchos::ArrayRCP< size_t > &offsets) const
Get offsets of the diagonal entries in the matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1262
Xpetra::TpetraBlockCrsMatrix::getGlobalNumRows
global_size_t getGlobalNumRows() const
Number of global elements in the row map of this matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:294
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getGlobalNumCols
global_size_t getGlobalNumCols() const
Number of global columns in the matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:877
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getLocalMaxNumRowEntries
size_t getLocalMaxNumRowEntries() const
Returns the maximum number of entries across all rows/columns on this node.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:901
Xpetra::TpetraBlockCrsMatrix::getNumEntriesInLocalRow
size_t getNumEntriesInLocalRow(LocalOrdinal localRow) const
Returns the current number of entries on this node in the specified local row.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:342
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::scale
void scale(const Scalar &alpha)
Scale the current values of a matrix, this = alpha*this (not implemented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1123
Xpetra::TpetraBlockCrsMatrix::expertStaticFillComplete
void expertStaticFillComplete(const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &domainMap, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap, const RCP< const Import< LocalOrdinal, GlobalOrdinal, Node > > &importer=Teuchos::null, const RCP< const Export< LocalOrdinal, GlobalOrdinal, Node > > &exporter=Teuchos::null, const RCP< ParameterList > &params=Teuchos::null)
Expert static fill complete.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:251
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::TpetraBlockCrsMatrix
TpetraBlockCrsMatrix(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap, size_t maxNumEntriesPerRow, const Teuchos::RCP< Teuchos::ParameterList > &params=Teuchos::null)
Constructor specifying column Map and fixed number of entries for each row (not implemented) ...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:751
Xpetra::TpetraBlockCrsMatrix::setObjectLabel
void setObjectLabel(const std::string &objectLabel)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:523
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getTpetra_BlockCrsMatrix
RCP< const Tpetra::BlockCrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getTpetra_BlockCrsMatrix() const
Get the underlying Tpetra matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1009
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getLocalNumCols
size_t getLocalNumCols() const
Returns the number of columns connected to the locally owned rows of this matrix. ...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:883
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::doImport
void doImport(const DistObject< char, LocalOrdinal, GlobalOrdinal, Node > &source, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Import (using an Exporter).
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:989
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getLocalNumCols
size_t getLocalNumCols() const
Returns the number of columns connected to the locally owned rows of this matrix. ...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1183
XPETRA_TPETRA_ETI_EXCEPTION
#define XPETRA_TPETRA_ETI_EXCEPTION(cl, obj, go, node)
Definition: Xpetra_Exceptions.hpp:43
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::doImport
void doImport(const DistObject< char, LocalOrdinal, GlobalOrdinal, Node > &source, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Import (using an Exporter).
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1289
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::allocateAllValues
void allocateAllValues(size_t numNonZeros, ArrayRCP< size_t > &rowptr, ArrayRCP< LocalOrdinal > &colind, ArrayRCP< Scalar > &values)
Allocates and returns ArrayRCPs of the Crs arrays — This is an Xpetra-only routine.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:827
Xpetra::TpetraBlockCrsMatrix::getRangeMap
const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > getRangeMap() const
Returns the Map associated with the range of this operator, which must be compatible with Y...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:500
Xpetra::TpetraBlockCrsMatrix::getLocalNumRows
size_t getLocalNumRows() const
Returns the number of matrix rows owned on the calling node.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:310
Xpetra::TpetraBlockCrsMatrix::removeEmptyProcessesInPlace
void removeEmptyProcessesInPlace(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &newMap)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:647
Xpetra::TpetraBlockCrsMatrix::describe
void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel=Teuchos::Describable::verbLevel_default) const
Print the object with some verbosity level to an FancyOStream object.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:515
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getFrobeniusNorm
ScalarTraits< Scalar >::magnitudeType getFrobeniusNorm() const
Returns the Frobenius norm of the matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1216
Xpetra::TpetraBlockCrsMatrix::getGlobalNumEntries
global_size_t getGlobalNumEntries() const
Returns the global number of entries in this matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:326
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getLocalNumRows
size_t getLocalNumRows() const
Returns the number of matrix rows owned on the calling node.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:880
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getGlobalNumEntries
global_size_t getGlobalNumEntries() const
Returns the global number of entries in this matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:886
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::isLocallyIndexed
bool isLocallyIndexed() const
If matrix indices are in the local range, this function returns true. Otherwise, this function return...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:904
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::description
std::string description() const
A simple one-line description of this object.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1250
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::setAllValues
void setAllValues(const ArrayRCP< size_t > &rowptr, const ArrayRCP< LocalOrdinal > &colind, const ArrayRCP< Scalar > &values)
Sets the 1D pointer arrays of the graph (not impelmented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:830
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::TpetraBlockCrsMatrix
TpetraBlockCrsMatrix(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, size_t maxNumEntriesPerRow, const Teuchos::RCP< Teuchos::ParameterList > &params=Teuchos::null)
Constructor specifying fixed number of entries for each row (not implemented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1055
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::fillComplete
void fillComplete(const RCP< ParameterList > &params=null)
Signal that data entry is complete.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:849
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::TpetraBlockCrsMatrix
TpetraBlockCrsMatrix(const Teuchos::RCP< const CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > &graph, const LocalOrdinal blockSize)
Constructor specifying a previously constructed graph &amp; blocksize.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1070
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getLocalNumRows
size_t getLocalNumRows() const
Returns the number of matrix rows owned on the calling node.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1180
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::TpetraBlockCrsMatrix
TpetraBlockCrsMatrix(const Teuchos::RCP< Tpetra::BlockCrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &mtx)
TpetraBlockCrsMatrix constructor to wrap a Tpetra::BlockCrsMatrix object.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1004
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getLocalDiagCopy
void getLocalDiagCopy(Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &diag, const Kokkos::View< const size_t *, typename Node::device_type, Kokkos::MemoryUnmanaged > &offsets) const
Get a copy of the diagonal entries owned by this node, with local row indices.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:968
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::doImport
void doImport(const DistObject< char, LocalOrdinal, GlobalOrdinal, Node > &source, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Import.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:981
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getGlobalRowCopy
void getGlobalRowCopy(GlobalOrdinal GlobalRow, const ArrayView< GlobalOrdinal > &indices, const ArrayView< Scalar > &values, size_t &numEntries) const
Extract a list of entries in a specified global row of this matrix. Put into pre-allocated storage...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:928
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::replaceGlobalValues
void replaceGlobalValues(GlobalOrdinal globalRow, const ArrayView< const GlobalOrdinal > &cols, const ArrayView< const Scalar > &vals)
Replace matrix entries, using global IDs (not implemented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:814
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::haveGlobalConstants
bool haveGlobalConstants() const
Returns true if globalConstants have been computed; false otherwise.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1234
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::TpetraBlockCrsMatrixClass
TpetraBlockCrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > TpetraBlockCrsMatrixClass
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1046
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::scale
void scale(const Scalar &alpha)
Scale the current values of a matrix, this = alpha*this (not implemented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:823
Xpetra::TpetraBlockCrsMatrix::haveGlobalConstants
virtual bool haveGlobalConstants() const
Returns true if globalConstants have been computed; false otherwise.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:465
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::supportsRowViews
bool supportsRowViews() const
Returns true if getLocalRowView() and getGlobalRowView() are valid for this class.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:919
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getAllValues
void getAllValues(ArrayRCP< Scalar > &values)
Gets the 1D pointer arrays of the graph (not implemented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1136
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::setAllValues
void setAllValues(const ArrayRCP< size_t > &rowptr, const ArrayRCP< LocalOrdinal > &colind, const ArrayRCP< Scalar > &values)
Sets the 1D pointer arrays of the graph (not impelmented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1130
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::TpetraBlockCrsMatrixClass
TpetraBlockCrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > TpetraBlockCrsMatrixClass
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:732
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::TpetraBlockCrsMatrix
TpetraBlockCrsMatrix(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap, const ArrayRCP< const size_t > &NumEntriesPerRowToAlloc, const Teuchos::RCP< Teuchos::ParameterList > &params=Teuchos::null)
Constructor specifying column Map and number of entries in each row (not implemented) ...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1064
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::TpetraImportClass
TpetraImport< LocalOrdinal, GlobalOrdinal, Node > TpetraImportClass
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:734
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::replaceDiag
void replaceDiag(const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &diag)
Replace the diagonal entries of the matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:970
Xpetra::TpetraBlockCrsMatrix::getTpetra_BlockCrsMatrix
RCP< const Tpetra::BlockCrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getTpetra_BlockCrsMatrix() const
Get the underlying Tpetra matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:665
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getCrsGraph
RCP< const CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > getCrsGraph() const
Returns the CrsGraph associated with this matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1171
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::TpetraBlockCrsMatrix
TpetraBlockCrsMatrix(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, const ArrayRCP< const size_t > &NumEntriesPerRowToAlloc, const Teuchos::RCP< Teuchos::ParameterList > &params=Teuchos::null)
Constructor specifying (possibly different) number of entries in each row (not implemented) ...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:746
Xpetra::TpetraBlockCrsMatrix::apply
void apply(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &X, MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &Y, Teuchos::ETransp mode=Teuchos::NO_TRANS, Scalar alpha=ScalarTraits< Scalar >::one(), Scalar beta=ScalarTraits< Scalar >::zero()) const
Computes the sparse matrix-multivector multiplication.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:469
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::description
std::string description() const
A simple one-line description of this object.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:950
Xpetra::TpetraBlockCrsMatrix::getLocalDiagCopy
void getLocalDiagCopy(Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &diag) const
Get a copy of the diagonal entries owned by this node, with local row idices.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:531
Xpetra::TpetraBlockCrsMatrix::description
std::string description() const
A simple one-line description of this object.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:508
Xpetra::TpetraBlockCrsMatrix::setAllValues
void setAllValues(const ArrayRCP< size_t > &rowptr, const ArrayRCP< LocalOrdinal > &colind, const ArrayRCP< Scalar > &values)
Sets the 1D pointer arrays of the graph (not impelmented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:197
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getTpetra_BlockCrsMatrixNonConst
RCP< Tpetra::BlockCrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getTpetra_BlockCrsMatrixNonConst() const
Get the underlying Tpetra matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1312
Xpetra::TpetraBlockCrsMatrix::getLocalNumCols
size_t getLocalNumCols() const
Returns the number of columns connected to the locally owned rows of this matrix. ...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:318
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getGlobalNumRows
global_size_t getGlobalNumRows() const
Number of global elements in the row map of this matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1174
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::TpetraBlockCrsMatrix
TpetraBlockCrsMatrix(const Teuchos::RCP< const Tpetra::BlockCrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix, const Export< LocalOrdinal, GlobalOrdinal, Node > &RowExporter, const Teuchos::RCP< const Export< LocalOrdinal, GlobalOrdinal, Node > > DomainExporter, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &domainMap, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap, const Teuchos::RCP< Teuchos::ParameterList > &params)
Constructor for a fused export (not implemented(.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:793
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::describe
void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel=Teuchos::Describable::verbLevel_default) const
Print the object with some verbosity level to an FancyOStream object.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:953
Xpetra::toTpetra
RCP< const Tpetra::CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > toTpetra(const RCP< const CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > &graph)
Definition: Xpetra_TpetraCrsGraph_decl.hpp:372
Xpetra::TpetraBlockCrsMatrix::isFillComplete
bool isFillComplete() const
Returns true if the matrix is in compute mode, i.e. if fillComplete() has been called.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:380
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getNumEntriesInGlobalRow
size_t getNumEntriesInGlobalRow(GlobalOrdinal globalRow) const
Returns the current number of entries in the (locally owned) global row.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1195
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getLocalDiagCopy
void getLocalDiagCopy(Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &diag, const Teuchos::ArrayView< const size_t > &offsets) const
Get a copy of the diagonal entries owned by this node, with local row indices.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1265
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::TpetraBlockCrsMatrix
TpetraBlockCrsMatrix(const Teuchos::RCP< const Tpetra::BlockCrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix, const Import< LocalOrdinal, GlobalOrdinal, Node > &RowImporter, const Teuchos::RCP< const Import< LocalOrdinal, GlobalOrdinal, Node > > DomainImporter, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &domainMap, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap, const Teuchos::RCP< Teuchos::ParameterList > &params)
Constructor for a fused import (not implemented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1087
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getLocalDiagOffsets
void getLocalDiagOffsets(Teuchos::ArrayRCP< size_t > &offsets) const
Get offsets of the diagonal entries in the matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:962
Xpetra::TpetraBlockCrsMatrix::fillComplete
void fillComplete(const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &domainMap, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap, const RCP< ParameterList > &params=null)
Signal that data entry is complete, specifying domain and range maps.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:230
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::setAllToScalar
void setAllToScalar(const Scalar &alpha)
Set all matrix entries equal to scalarThis.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:820
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::rightScale
void rightScale(const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &x)
Right scale matrix using the given vector entries.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:973
Xpetra::TpetraBlockCrsMatrix::getNumEntriesInGlobalRow
size_t getNumEntriesInGlobalRow(GlobalOrdinal globalRow) const
Returns the current number of entries in the (locally owned) global row.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:350
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getGlobalNumCols
global_size_t getGlobalNumCols() const
Number of global columns in the matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1177
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getGlobalRowView
void getGlobalRowView(GlobalOrdinal GlobalRow, ArrayView< const GlobalOrdinal > &indices, ArrayView< const Scalar > &values) const
Extract a const, non-persisting view of global indices in a specified row of the matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1225
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::isFillActive
bool isFillActive() const
Returns true if the matrix is in edit mode.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1213
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::replaceDomainMapAndImporter
void replaceDomainMapAndImporter(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &newDomainMap, Teuchos::RCP< const Import< LocalOrdinal, GlobalOrdinal, Node > > &newImporter)
Replaces the current domainMap and importer with the user-specified objects.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1153
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::TpetraBlockCrsMatrix
TpetraBlockCrsMatrix(const Teuchos::RCP< const Tpetra::BlockCrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &domainMap=Teuchos::null, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap=Teuchos::null, const Teuchos::RCP< Teuchos::ParameterList > &params=Teuchos::null)
Constructor for a fused export (not implemented(.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:778
XPETRA_DYNAMIC_CAST
#define XPETRA_DYNAMIC_CAST(type, obj, newObj, exceptionMsg)
Definition: Xpetra_Exceptions.hpp:22
Xpetra::global_size_t
size_t global_size_t
Global size_t object.
Definition: Xpetra_ConfigDefs.hpp:138
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getAllValues
void getAllValues(ArrayRCP< Scalar > &values)
Gets the 1D pointer arrays of the graph (not implemented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:836
Xpetra::TpetraBlockCrsMatrix::getFrobeniusNorm
ScalarTraits< Scalar >::magnitudeType getFrobeniusNorm() const
Returns the Frobenius norm of the matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:392
Xpetra::TpetraBlockCrsMatrix::getGlobalMaxNumRowEntries
size_t getGlobalMaxNumRowEntries() const
Returns the maximum number of entries across all rows/columns on all nodes.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:356
Xpetra::TpetraBlockCrsMatrix::doImport
void doImport(const DistObject< char, LocalOrdinal, GlobalOrdinal, Node > &source, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Import.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:616
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::TpetraBlockCrsMatrix
TpetraBlockCrsMatrix(const TpetraBlockCrsMatrix &matrix)
Deep copy constructor.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:956
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getColMap
const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > getColMap() const
Returns the Map that describes the column distribution in this matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1168
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::removeEmptyProcessesInPlace
void removeEmptyProcessesInPlace(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &newMap)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:996
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getLocalNumEntries
size_t getLocalNumEntries() const
Returns the local number of entries in this matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1189
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::TpetraBlockCrsMatrix
TpetraBlockCrsMatrix(const Teuchos::RCP< const Tpetra::BlockCrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &domainMap=Teuchos::null, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap=Teuchos::null, const Teuchos::RCP< Teuchos::ParameterList > &params=Teuchos::null)
Constructor for a fused export (not implemented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1080
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::TpetraBlockCrsMatrix
TpetraBlockCrsMatrix(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap, const ArrayRCP< const size_t > &NumEntriesPerRowToAlloc, const Teuchos::RCP< Teuchos::ParameterList > &params=Teuchos::null)
Constructor specifying column Map and number of entries in each row (not implemented) ...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:756
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::allocateAllValues
void allocateAllValues(size_t numNonZeros, ArrayRCP< size_t > &rowptr, ArrayRCP< LocalOrdinal > &colind, ArrayRCP< Scalar > &values)
Allocates and returns ArrayRCPs of the Crs arrays — This is an Xpetra-only routine.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1127
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::isFillComplete
bool isFillComplete() const
Returns true if the matrix is in compute mode, i.e. if fillComplete() has been called.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1210
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::isFillComplete
bool isFillComplete() const
Returns true if the matrix is in compute mode, i.e. if fillComplete() has been called.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:910
Xpetra::TpetraBlockCrsMatrix::getMap
Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > getMap() const
Implements DistObject interface.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:608
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::TpetraBlockCrsMatrix
TpetraBlockCrsMatrix(const Teuchos::RCP< const CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > &graph, const Teuchos::RCP< Teuchos::ParameterList > &params=Teuchos::null)
Constructor specifying a previously constructed graph ( not implemented )
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:761
Xpetra::TpetraBlockCrsMatrix::supportsRowViews
bool supportsRowViews() const
Returns true if getLocalRowView() and getGlobalRowView() are valid for this class.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:398
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::doExport
void doExport(const DistObject< char, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Export (using an Importer).
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:993
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getMap
Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > getMap() const
Implements DistObject interface.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1278
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getGlobalMaxNumRowEntries
size_t getGlobalMaxNumRowEntries() const
Returns the maximum number of entries across all rows/columns on all nodes.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:898
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::leftScale
void leftScale(const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &x)
Left scale matrix using the given vector entries.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1272
Xpetra::EpetraNode
Tpetra::KokkosCompat::KokkosSerialWrapperNode EpetraNode
Definition: Xpetra_Map_decl.hpp:35
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getRangeMap
const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > getRangeMap() const
Returns the Map associated with the range of this operator, which must be compatible with Y...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1245
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::TpetraBlockCrsMatrix
TpetraBlockCrsMatrix(const Teuchos::RCP< const CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > &graph, const Teuchos::RCP< Teuchos::ParameterList > &params=Teuchos::null)
Constructor specifying a previously constructed graph ( not implemented )
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1067
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::setAllValues
void setAllValues(const typename local_matrix_type::row_map_type &ptr, const typename local_matrix_type::StaticCrsGraphType::entries_type::non_const_type &ind, const typename local_matrix_type::values_type &val)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1023
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::TpetraBlockCrsMatrix
TpetraBlockCrsMatrix(const Teuchos::RCP< const Tpetra::BlockCrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &domainMap=Teuchos::null, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap=Teuchos::null, const Teuchos::RCP< Teuchos::ParameterList > &params=Teuchos::null)
Constructor for a fused import (not implemented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1073
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::doImport
void doImport(const DistObject< char, LocalOrdinal, GlobalOrdinal, Node > &source, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Import.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1281
Xpetra::TpetraBlockCrsMatrix::TpetraBlockCrsMatrix
TpetraBlockCrsMatrix(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, size_t maxNumEntriesPerRow, const Teuchos::RCP< Teuchos::ParameterList > &params=Teuchos::null)
Constructor specifying fixed number of entries for each row (not implemented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:21
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::TpetraBlockCrsMatrix
TpetraBlockCrsMatrix(const Teuchos::RCP< const CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > &graph, const LocalOrdinal blockSize)
Constructor specifying a previously constructed graph &amp; blocksize.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:766
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::TpetraVectorClass
TpetraVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > TpetraVectorClass
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:733
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getGlobalRowView
void getGlobalRowView(GlobalOrdinal GlobalRow, ArrayView< const GlobalOrdinal > &indices, ArrayView< const Scalar > &values) const
Extract a const, non-persisting view of global indices in a specified row of the matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:925
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getTpetra_BlockCrsMatrix
RCP< const Tpetra::BlockCrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getTpetra_BlockCrsMatrix() const
Get the underlying Tpetra matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1309
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getLocalRowView
void getLocalRowView(LocalOrdinal LocalRow, ArrayView< const LocalOrdinal > &indices, ArrayView< const Scalar > &values) const
Extract a const, non-persisting view of local indices in a specified row of the matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1231
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::replaceLocalValues
void replaceLocalValues(LocalOrdinal localRow, const ArrayView< const LocalOrdinal > &cols, const ArrayView< const Scalar > &vals)
Replace matrix entries, using local IDs.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1117
Xpetra_TpetraBlockCrsMatrix_decl.hpp
Xpetra::TpetraBlockCrsMatrix::getGlobalRowView
void getGlobalRowView(GlobalOrdinal GlobalRow, ArrayView< const GlobalOrdinal > &indices, ArrayView< const Scalar > &values) const
Extract a const, non-persisting view of global indices in a specified row of the matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:434
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getAllValues
void getAllValues(ArrayRCP< const size_t > &rowptr, ArrayRCP< const LocalOrdinal > &colind, ArrayRCP< const Scalar > &values) const
Gets the 1D pointer arrays of the graph (not implemented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1133
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::isLocallyIndexed
bool isLocallyIndexed() const
If matrix indices are in the local range, this function returns true. Otherwise, this function return...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1204
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getNumEntriesInLocalRow
size_t getNumEntriesInLocalRow(LocalOrdinal localRow) const
Returns the current number of entries on this node in the specified local row.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:892
Xpetra::TpetraBlockCrsMatrix::leftScale
void leftScale(const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &x)
Left scale operator with given vector values.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:595
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::doExport
void doExport(const DistObject< char, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Export.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:985
Xpetra::TpetraBlockCrsMatrix::getAllValues
void getAllValues(ArrayRCP< const size_t > &rowptr, ArrayRCP< const LocalOrdinal > &colind, ArrayRCP< const Scalar > &values) const
Gets the 1D pointer arrays of the graph (not implemented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:204
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getGlobalNumEntries
global_size_t getGlobalNumEntries() const
Returns the global number of entries in this matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1186
Xpetra::TpetraBlockCrsMatrix::getDomainMap
const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > getDomainMap() const
Returns the Map associated with the domain of this operator. This will be null until fillComplete() i...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:492
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getLocalMaxNumRowEntries
size_t getLocalMaxNumRowEntries() const
Returns the maximum number of entries across all rows/columns on this node.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1201
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getColMap
const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > getColMap() const
Returns the Map that describes the column distribution in this matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:868
Xpetra::TpetraBlockCrsMatrix::replaceDiag
void replaceDiag(const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &diag)
Replace the diagonal entries of the matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:589
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::TpetraVectorClass
TpetraVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > TpetraVectorClass
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1047
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::replaceLocalValues
void replaceLocalValues(LocalOrdinal localRow, const ArrayView< const LocalOrdinal > &cols, const ArrayView< const Scalar > &vals)
Replace matrix entries, using local IDs.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:817
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::fillComplete
void fillComplete(const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &domainMap, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap, const RCP< ParameterList > &params=null)
Signal that data entry is complete, specifying domain and range maps.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:845
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::expertStaticFillComplete
void expertStaticFillComplete(const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &domainMap, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap, const RCP< const Import< LocalOrdinal, GlobalOrdinal, Node > > &importer=Teuchos::null, const RCP< const Export< LocalOrdinal, GlobalOrdinal, Node > > &exporter=Teuchos::null, const RCP< ParameterList > &params=Teuchos::null)
Expert static fill complete.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:856
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::setAllValues
void setAllValues(const typename local_matrix_type::row_map_type &ptr, const typename local_matrix_type::StaticCrsGraphType::entries_type::non_const_type &ind, const typename local_matrix_type::values_type &val)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1323
Xpetra::CombineMode
CombineMode
Xpetra::Combine Mode enumerable type.
Definition: Xpetra_ConfigDefs.hpp:176
Xpetra::TpetraBlockCrsMatrix::getColMap
const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > getColMap() const
Returns the Map that describes the column distribution in this matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:274
Xpetra::TpetraBlockCrsMatrix::replaceGlobalValues
void replaceGlobalValues(GlobalOrdinal globalRow, const ArrayView< const GlobalOrdinal > &cols, const ArrayView< const Scalar > &vals)
Replace matrix entries, using global IDs (not implemented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:157
XPETRA_MONITOR
#define XPETRA_MONITOR(funcName)
Definition: Xpetra_ConfigDefs.hpp:96
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::apply
void apply(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &X, MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &Y, Teuchos::ETransp mode=Teuchos::NO_TRANS, Scalar alpha=ScalarTraits< Scalar >::one(), Scalar beta=ScalarTraits< Scalar >::zero()) const
Computes the sparse matrix-multivector multiplication.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:939
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getLocalRowView
void getLocalRowView(LocalOrdinal LocalRow, ArrayView< const LocalOrdinal > &indices, ArrayView< const Scalar > &values) const
Extract a const, non-persisting view of local indices in a specified row of the matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:931
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::TpetraBlockCrsMatrix
TpetraBlockCrsMatrix(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, size_t maxNumEntriesPerRow, const Teuchos::RCP< Teuchos::ParameterList > &params=Teuchos::null)
Constructor specifying fixed number of entries for each row (not implemented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:741
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::haveGlobalConstants
bool haveGlobalConstants() const
Returns true if globalConstants have been computed; false otherwise.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:934
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::TpetraBlockCrsMatrix
TpetraBlockCrsMatrix(const Teuchos::RCP< Tpetra::BlockCrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &mtx)
TpetraBlockCrsMatrix constructor to wrap a Tpetra::BlockCrsMatrix object.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1304
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getRowMap
const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > getRowMap() const
Returns the Map that describes the row distribution in this matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1165
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getLocalDiagCopy
void getLocalDiagCopy(Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &diag, const Teuchos::ArrayView< const size_t > &offsets) const
Get a copy of the diagonal entries owned by this node, with local row indices.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:965
Xpetra::TpetraBlockCrsMatrix::getLocalMaxNumRowEntries
size_t getLocalMaxNumRowEntries() const
Returns the maximum number of entries across all rows/columns on this node.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:362
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::isGloballyIndexed
bool isGloballyIndexed() const
If matrix indices are in the global range, this function returns true. Otherwise, this function retur...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1207
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getNumEntriesInLocalRow
size_t getNumEntriesInLocalRow(LocalOrdinal localRow) const
Returns the current number of entries on this node in the specified local row.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1192
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::isFillActive
bool isFillActive() const
Returns true if the matrix is in edit mode.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:913
Xpetra::TpetraBlockCrsMatrix::getLocalDiagOffsets
void getLocalDiagOffsets(Teuchos::ArrayRCP< size_t > &offsets) const
Get offsets of the diagonal entries in the matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:558
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getLocalNumEntries
size_t getLocalNumEntries() const
Returns the local number of entries in this matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:889
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::TpetraBlockCrsMatrix
TpetraBlockCrsMatrix(const Teuchos::RCP< const Tpetra::BlockCrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix, const Export< LocalOrdinal, GlobalOrdinal, Node > &RowExporter, const Teuchos::RCP< const Export< LocalOrdinal, GlobalOrdinal, Node > > DomainExporter, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &domainMap, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap, const Teuchos::RCP< Teuchos::ParameterList > &params)
Constructor for a fused export (not implemented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1095
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::apply
void apply(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &X, MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &Y, Teuchos::ETransp mode=Teuchos::NO_TRANS, Scalar alpha=ScalarTraits< Scalar >::one(), Scalar beta=ScalarTraits< Scalar >::zero()) const
Computes the sparse matrix-multivector multiplication.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1239
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getGlobalMaxNumRowEntries
size_t getGlobalMaxNumRowEntries() const
Returns the maximum number of entries across all rows/columns on all nodes.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1198
Xpetra::CrsMatrix::local_matrix_type
KokkosSparse::CrsMatrix< impl_scalar_type, LocalOrdinal, execution_space, void, typename local_graph_type::size_type > local_matrix_type
The specialization of Kokkos::CrsMatrix that represents the part of the sparse matrix on each MPI pro...
Definition: Xpetra_CrsMatrix.hpp:289
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::local_matrix_type
Xpetra::CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::local_matrix_type local_matrix_type
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1015
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getGlobalRowCopy
void getGlobalRowCopy(GlobalOrdinal GlobalRow, const ArrayView< GlobalOrdinal > &indices, const ArrayView< Scalar > &values, size_t &numEntries) const
Extract a list of entries in a specified global row of this matrix. Put into pre-allocated storage...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1228
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::insertLocalValues
void insertLocalValues(LocalOrdinal localRow, const ArrayView< const LocalOrdinal > &cols, const ArrayView< const Scalar > &vals)
Insert matrix entries, using local IDs (not implemented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1111
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getGlobalNumRows
global_size_t getGlobalNumRows() const
Number of global elements in the row map of this matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:874
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::TpetraBlockCrsMatrix
TpetraBlockCrsMatrix(const Teuchos::RCP< const Tpetra::BlockCrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix, const Import< LocalOrdinal, GlobalOrdinal, Node > &RowImporter, const Teuchos::RCP< const Import< LocalOrdinal, GlobalOrdinal, Node > > DomainImporter, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &domainMap, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap, const Teuchos::RCP< Teuchos::ParameterList > &params)
Constructor for a fused import ( not implemented )
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:785
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getAllValues
void getAllValues(ArrayRCP< const size_t > &rowptr, ArrayRCP< const LocalOrdinal > &colind, ArrayRCP< const Scalar > &values) const
Gets the 1D pointer arrays of the graph (not implemented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:833
Xpetra::TpetraBlockCrsMatrix::replaceDomainMapAndImporter
void replaceDomainMapAndImporter(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &newDomainMap, Teuchos::RCP< const Import< LocalOrdinal, GlobalOrdinal, Node > > &newImporter)
Replaces the current domainMap and importer with the user-specified objects.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:244
Xpetra::TpetraBlockCrsMatrix::isFillActive
bool isFillActive() const
Returns true if the matrix is in edit mode.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:386
Xpetra::TpetraBlockCrsMatrix::doExport
void doExport(const DistObject< char, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Export.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:624
Xpetra::TpetraBlockCrsMatrix::getRowMap
const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > getRowMap() const
Returns the Map that describes the row distribution in this matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:266
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::insertLocalValues
void insertLocalValues(LocalOrdinal localRow, const ArrayView< const LocalOrdinal > &cols, const ArrayView< const Scalar > &vals)
Insert matrix entries, using local IDs (not implemented)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:811
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::getLocalDiagCopy
void getLocalDiagCopy(Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &diag, const Kokkos::View< const size_t *, typename Node::device_type, Kokkos::MemoryUnmanaged > &offsets) const
Get a copy of the diagonal entries owned by this node, with local row indices.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1268
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getRangeMap
const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > getRangeMap() const
Returns the Map associated with the range of this operator, which must be compatible with Y...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:945
Xpetra::TpetraBlockCrsMatrix::getLocalRowView
void getLocalRowView(LocalOrdinal LocalRow, ArrayView< const LocalOrdinal > &indices, ArrayView< const Scalar > &values) const
Extract a const, non-persisting view of local indices in a specified row of the matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:421
Xpetra::TpetraBlockCrsMatrix::resumeFill
void resumeFill(const RCP< ParameterList > &params=null)
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:224
Xpetra::TpetraBlockCrsMatrix< Scalar, int, long long, EpetraNode >::TpetraBlockCrsMatrix
TpetraBlockCrsMatrix(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, const ArrayRCP< const size_t > &NumEntriesPerRowToAlloc, const Teuchos::RCP< Teuchos::ParameterList > &params=Teuchos::null)
Constructor specifying (possibly different) number of entries in each row (not implemented) ...
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:1058
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getNumEntriesInGlobalRow
size_t getNumEntriesInGlobalRow(GlobalOrdinal globalRow) const
Returns the current number of entries in the (locally owned) global row.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:895
Xpetra::TpetraBlockCrsMatrix< Scalar, int, int, EpetraNode >::getFrobeniusNorm
ScalarTraits< Scalar >::magnitudeType getFrobeniusNorm() const
Returns the Frobenius norm of the matrix.
Definition: Xpetra_TpetraBlockCrsMatrix_def.hpp:916
Xpetra::toXpetra
RCP< const CrsGraph< int, GlobalOrdinal, Node > > toXpetra(const Epetra_CrsGraph &g)
Definition: Xpetra_EpetraCrsGraph.cpp:29