67 Teuchos::RCP<Tpetra::CrsGraph<typename graph_type::local_ordinal_type, typename graph_type::global_ordinal_type, typename graph_type::node_type> >
70 typedef typename graph_type::local_ordinal_type LO;
71 typedef typename graph_type::global_ordinal_type GO;
72 typedef typename graph_type::node_type NO;
73 typedef Tpetra::Map<LO, GO, NO> map_type;
74 typedef Tpetra::CrsGraph<LO, GO, NO> crs_graph_type;
76 const size_t maxDiagEntPerRow = 1;
81 Teuchos::RCP<graph_type> diagonalGraph;
82 diagonalGraph = Teuchos::rcp(
new crs_graph_type(graph.getRowMap(), maxDiagEntPerRow));
83 const map_type& meshRowMap = *(graph.getRowMap());
85 Teuchos::Array<GO> diagGblColInds(maxDiagEntPerRow);
87 for (LO lclRowInd = meshRowMap.getMinLocalIndex(); lclRowInd <= meshRowMap.getMaxLocalIndex(); ++lclRowInd) {
88 const GO gblRowInd = meshRowMap.getGlobalElement(lclRowInd);
89 diagGblColInds[0] = gblRowInd;
90 diagonalGraph->insertGlobalIndices(gblRowInd, diagGblColInds());
93 diagonalGraph->fillComplete(graph.getDomainMap(), graph.getRangeMap());
Teuchos::RCP< Tpetra::CrsGraph< typename graph_type::local_ordinal_type, typename graph_type::global_ordinal_type, typename graph_type::node_type > > computeDiagonalGraph(graph_type const &graph)
Compute and return the graph of the diagonal of the input graph.
Definition Ifpack2_Utilities.hpp:68