Xpetra  Version of the Day
Xpetra_CrsGraphFactory.hpp
Go to the documentation of this file.
1 // @HEADER
2 //
3 // ***********************************************************************
4 //
5 // Xpetra: A linear algebra interface package
6 // Copyright 2012 Sandia Corporation
7 //
8 // Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
9 // the U.S. Government retains certain rights in this software.
10 //
11 // Redistribution and use in source and binary forms, with or without
12 // modification, are permitted provided that the following conditions are
13 // met:
14 //
15 // 1. Redistributions of source code must retain the above copyright
16 // notice, this list of conditions and the following disclaimer.
17 //
18 // 2. Redistributions in binary form must reproduce the above copyright
19 // notice, this list of conditions and the following disclaimer in the
20 // documentation and/or other materials provided with the distribution.
21 //
22 // 3. Neither the name of the Corporation nor the names of the
23 // contributors may be used to endorse or promote products derived from
24 // this software without specific prior written permission.
25 //
26 // THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
27 // EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
28 // IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
29 // PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
30 // CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
31 // EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
32 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
33 // PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
34 // LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
35 // NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
36 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
37 //
38 // Questions? Contact
39 // Jonathan Hu (jhu@sandia.gov)
40 // Andrey Prokopenko (aprokop@sandia.gov)
41 // Ray Tuminaro (rstumin@sandia.gov)
42 //
43 // ***********************************************************************
44 //
45 // @HEADER
46 #ifndef XPETRA_CRSGRAPHFACTORY_HPP
47 #define XPETRA_CRSGRAPHFACTORY_HPP
48 
49 #include "Xpetra_ConfigDefs.hpp"
50 
51 #include "Xpetra_CrsGraph.hpp"
52 
53 #ifdef HAVE_XPETRA_TPETRA
54 #include "Xpetra_TpetraCrsGraph.hpp"
55 #endif
56 
57 #ifdef HAVE_XPETRA_EPETRA
58 #include "Xpetra_EpetraCrsGraph.hpp"
59 #endif
60 
61 #include "Xpetra_Exceptions.hpp"
62 
63 namespace Xpetra {
64 
65  template <class LocalOrdinal,
66  class GlobalOrdinal,
67  class Node = KokkosClassic::DefaultNode::DefaultNodeType>
68  class CrsGraphFactory {
69  private:
71  CrsGraphFactory() {}
72 
73  public:
75  static Teuchos::RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
76  Build (const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rowMap)
77  {
78  TEUCHOS_TEST_FOR_EXCEPTION(rowMap->lib() == UseEpetra, std::logic_error,
79  "Can't create Xpetra::EpetraCrsMatrix with these scalar/LO/GO types");
80 #ifdef HAVE_XPETRA_TPETRA
81  if (rowMap->lib() == UseTpetra)
82  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node>(rowMap, 0) );
83 #endif
84 
85  XPETRA_FACTORY_END;
86  }
87 
89  static Teuchos::RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
90  Build(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &map, size_t maxNumEntriesPerRow) {
91  XPETRA_MONITOR("CrsGraphFactory::Build");
92 
93 #ifdef HAVE_XPETRA_TPETRA
94  if (map->lib() == UseTpetra)
95  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node> (map, maxNumEntriesPerRow) );
96 #endif
97 
98  XPETRA_FACTORY_ERROR_IF_EPETRA(map->lib());
99  XPETRA_FACTORY_END;
100  TEUCHOS_UNREACHABLE_RETURN(null);
101  }
102 
104  static Teuchos::RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
105  Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap,
106  const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap,
107  size_t maxNumEntriesPerRow,
108  const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null) {
109  XPETRA_MONITOR("CrsGraphFactory::Build");
110 
111 #ifdef HAVE_XPETRA_TPETRA
112  if (rowMap->lib() == UseTpetra)
113  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node>(rowMap, colMap, maxNumEntriesPerRow, plist) );
114 #endif
115 
116  XPETRA_FACTORY_ERROR_IF_EPETRA(rowMap->lib());
117  XPETRA_FACTORY_END;
118  TEUCHOS_UNREACHABLE_RETURN(null);
119  }
120 
121 
123  static Teuchos::RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
124  Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap,
125  const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap,
126  const ArrayRCP< const size_t > &NumEntriesPerRowToAlloc,
127  const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null) {
128  XPETRA_MONITOR("CrsGraphFactory::Build");
129 
130 #ifdef HAVE_XPETRA_TPETRA
131  if (rowMap->lib() == UseTpetra)
132  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node>(rowMap, colMap, NumEntriesPerRowToAlloc, plist) );
133 #endif
134 
135  XPETRA_FACTORY_ERROR_IF_EPETRA(rowMap->lib());
136  XPETRA_FACTORY_END;
137  TEUCHOS_UNREACHABLE_RETURN(null);
138  }
139 
140 
142  static Teuchos::RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
143  Build(const RCP<const CrsGraph< LocalOrdinal, GlobalOrdinal, Node > >& sourceGraph,
144  const Import< LocalOrdinal, GlobalOrdinal, Node > & importer,
145  const RCP<const Map< LocalOrdinal, GlobalOrdinal, Node >>& domainMap = Teuchos::null,
146  const RCP<const Map< LocalOrdinal, GlobalOrdinal, Node > >& rangeMap = Teuchos::null,
147  const RCP<Teuchos::ParameterList>& params = Teuchos::null) {
148 #ifdef HAVE_XPETRA_TPETRA
149  if (sourceGraph->getRowMap()->lib() == UseTpetra)
150  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node>(sourceGraph, importer, domainMap, rangeMap, params) );
151 #endif
152 
153  XPETRA_FACTORY_ERROR_IF_EPETRA(sourceGraph()->getRowMap()->lib());
154  XPETRA_FACTORY_END;
155  TEUCHOS_UNREACHABLE_RETURN(null);
156  }
157 
158 
159 
160 #ifdef HAVE_XPETRA_KOKKOS_REFACTOR
161 #ifdef HAVE_XPETRA_TPETRA
181  static Teuchos::RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
182  Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap,
183  const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap,
184  const typename Xpetra::CrsGraph<LocalOrdinal, GlobalOrdinal, Node>::local_graph_type::row_map_type& rowPointers,
185  const typename Xpetra::CrsGraph<LocalOrdinal, GlobalOrdinal, Node>::local_graph_type::entries_type::non_const_type& columnIndices,
186  const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null) {
187  XPETRA_MONITOR("CrsMatrixFactory::Build");
188 
189  if (rowMap->lib() == UseTpetra)
190  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node>(rowMap,
191  colMap,
192  rowPointers,
193  columnIndices,
194  plist));
195 
196  XPETRA_FACTORY_END;
197  }
198 
217  static Teuchos::RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
218  Build(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >& rowMap,
219  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >& colMap,
220  const typename Xpetra::CrsGraph<LocalOrdinal, GlobalOrdinal, Node>::local_graph_type& lclGraph,
221  const Teuchos::RCP<Teuchos::ParameterList>& params) {
222  XPETRA_MONITOR("CrsMatrixFactory::Build");
223 
224  if (rowMap->lib() == UseTpetra)
225  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node>(rowMap,
226  colMap,
227  lclGraph,
228  params));
229 
230  XPETRA_FACTORY_END;
231  }
232 
257  static Teuchos::RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
258  Build(const typename Xpetra::CrsGraph<LocalOrdinal, GlobalOrdinal, Node>::local_graph_type& lclGraph,
259  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >& rowMap,
260  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >& colMap,
261  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >& domainMap = Teuchos::null,
262  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >& rangeMap = Teuchos::null,
263  const Teuchos::RCP<Teuchos::ParameterList>& params = Teuchos::null) {
264  XPETRA_MONITOR("CrsMatrixFactory::Build");
265 
266  if (rowMap->lib() == UseTpetra)
267  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node>(lclGraph,
268  rowMap,
269  colMap,
270  domainMap,
271  rangeMap,
272  params));
273 
274  XPETRA_FACTORY_END;
275  }
276 #endif
277 #endif
278 
298  static Teuchos::RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
299  Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap,
300  const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap,
301  const Teuchos::ArrayRCP<size_t>& rowPointers,
302  const Teuchos::ArrayRCP<LocalOrdinal>& columnIndices,
303  const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null) {
304  XPETRA_MONITOR("CrsMatrixFactory::Build");
305 
306  if (rowMap->lib() == UseTpetra)
307  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node>(rowMap,
308  colMap,
309  rowPointers,
310  columnIndices,
311  plist));
312 
313  XPETRA_FACTORY_END;
314  }
315 
316  };
317 
318 // we need the Epetra specialization only if Epetra is enabled
319 #if (defined(HAVE_XPETRA_EPETRA) && !defined(XPETRA_EPETRA_NO_32BIT_GLOBAL_INDICES))
320 
321  template <>
322  class CrsGraphFactory<int, int, EpetraNode> {
323 
324  typedef int LocalOrdinal;
325  typedef int GlobalOrdinal;
326  typedef EpetraNode Node;
327 
328  private:
330  CrsGraphFactory() {}
331 
332  public:
334  static RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
335  Build (const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rowMap)
336  {
337  XPETRA_MONITOR("CrsMatrixFactory::Build");
338 #ifdef HAVE_XPETRA_TPETRA
339  if (rowMap->lib() == UseTpetra)
340  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node>(rowMap, 0) );
341 #endif
342 #ifdef HAVE_XPETRA_EPETRA
343  if(rowMap->lib() == UseEpetra)
344  return rcp( new EpetraCrsGraphT<int,Node>(rowMap));
345 #endif
346  XPETRA_FACTORY_END;
347  }
348 
349  static RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
350  Build(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &map, size_t maxNumEntriesPerRow) {
351  XPETRA_MONITOR("CrsGraphFactory::Build");
352 
353 #ifdef HAVE_XPETRA_TPETRA
354  if (map->lib() == UseTpetra)
355  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node> (map, maxNumEntriesPerRow) );
356 #endif
357 
358  if (map->lib() == UseEpetra)
359  return rcp( new EpetraCrsGraphT<int, Node>(map, maxNumEntriesPerRow) );
360 
361  XPETRA_FACTORY_END;
362  TEUCHOS_UNREACHABLE_RETURN(null);
363  }
364 
365  static Teuchos::RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
366  Build(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 > &plist=Teuchos::null) {
367  XPETRA_MONITOR("CrsGraphFactory::Build");
368 
369 #ifdef HAVE_XPETRA_TPETRA
370  if (rowMap->lib() == UseTpetra)
371  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node>(rowMap, colMap, NumEntriesPerRowToAlloc, plist) );
372 #endif
373 
374  if (rowMap->lib() == UseEpetra)
375  return rcp( new EpetraCrsGraphT<int, Node>(rowMap, colMap, NumEntriesPerRowToAlloc, plist) );
376 
377  XPETRA_FACTORY_END;
378  TEUCHOS_UNREACHABLE_RETURN(null);
379  }
380 
382  static Teuchos::RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
383  Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap,
384  const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap,
385  size_t maxNumEntriesPerRow,
386  const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null) {
387  XPETRA_MONITOR("CrsGraphFactory::Build");
388 
389 #ifdef HAVE_XPETRA_TPETRA
390  if (rowMap->lib() == UseTpetra)
391  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node>(rowMap, colMap, maxNumEntriesPerRow, plist) );
392 #endif
393  if (rowMap->lib() == UseEpetra)
394  return rcp( new EpetraCrsGraphT<int,Node>(rowMap, colMap, maxNumEntriesPerRow, plist) );
395 
396  XPETRA_FACTORY_END;
397  TEUCHOS_UNREACHABLE_RETURN(null);
398  }
399 
400 
402  static Teuchos::RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
403  Build(const RCP<const CrsGraph< LocalOrdinal, GlobalOrdinal, Node > >& sourceGraph,
404  const Import< LocalOrdinal, GlobalOrdinal, Node > & importer,
405  const RCP<const Map< LocalOrdinal, GlobalOrdinal, Node >>& domainMap = Teuchos::null,
406  const RCP<const Map< LocalOrdinal, GlobalOrdinal, Node > >& rangeMap = Teuchos::null,
407  const RCP<Teuchos::ParameterList>& params = Teuchos::null) {
408 #ifdef HAVE_XPETRA_TPETRA
409  if (sourceGraph->getRowMap()->lib() == UseTpetra)
410  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node>(sourceGraph, importer, domainMap, rangeMap, params) );
411 #endif
412  if (sourceGraph->getRowMap()->lib() == UseEpetra)
413  return rcp( new EpetraCrsGraphT<int, Node>(sourceGraph, importer, domainMap, rangeMap, params) );
414 
415  XPETRA_FACTORY_END;
416  TEUCHOS_UNREACHABLE_RETURN(null);
417  }
418 
419 #ifdef HAVE_XPETRA_KOKKOS_REFACTOR
420 #ifdef HAVE_XPETRA_TPETRA
440  static Teuchos::RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
441  Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap,
442  const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap,
443  const typename Xpetra::CrsGraph<LocalOrdinal, GlobalOrdinal, Node>::local_graph_type::row_map_type& rowPointers,
444  const typename Xpetra::CrsGraph<LocalOrdinal, GlobalOrdinal, Node>::local_graph_type::entries_type::non_const_type& columnIndices,
445  const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null) {
446  XPETRA_MONITOR("CrsMatrixFactory::Build");
447 
448  if (rowMap->lib() == UseTpetra)
449  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node>(rowMap,
450  colMap,
451  rowPointers,
452  columnIndices,
453  plist));
454 
455  XPETRA_FACTORY_END;
456  }
457 
476  static Teuchos::RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
477  Build(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >& rowMap,
478  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >& colMap,
479  const typename Xpetra::CrsGraph<LocalOrdinal, GlobalOrdinal, Node>::local_graph_type& lclGraph,
480  const Teuchos::RCP<Teuchos::ParameterList>& params) {
481  XPETRA_MONITOR("CrsMatrixFactory::Build");
482 
483  if (rowMap->lib() == UseTpetra)
484  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node>(rowMap,
485  colMap,
486  lclGraph,
487  params));
488 
489  XPETRA_FACTORY_END;
490  }
491 
516  static Teuchos::RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
517  Build(const typename Xpetra::CrsGraph<LocalOrdinal, GlobalOrdinal, Node>::local_graph_type& lclGraph,
518  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >& rowMap,
519  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >& colMap,
520  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >& domainMap = Teuchos::null,
521  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >& rangeMap = Teuchos::null,
522  const Teuchos::RCP<Teuchos::ParameterList>& params = Teuchos::null) {
523  XPETRA_MONITOR("CrsMatrixFactory::Build");
524 
525  if (rowMap->lib() == UseTpetra)
526  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node>(lclGraph,
527  rowMap,
528  colMap,
529  domainMap,
530  rangeMap,
531  params));
532 
533  XPETRA_FACTORY_END;
534  }
535 #endif
536 #endif
537 
557  static Teuchos::RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
558  Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap,
559  const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap,
560  const Teuchos::ArrayRCP<size_t>& rowPointers,
561  const Teuchos::ArrayRCP<LocalOrdinal>& columnIndices,
562  const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null) {
563  XPETRA_MONITOR("CrsMatrixFactory::Build");
564 
565  if (rowMap->lib() == UseTpetra)
566  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node>(rowMap,
567  colMap,
568  rowPointers,
569  columnIndices,
570  plist));
571 
572  XPETRA_FACTORY_END;
573  }
574 
575  };
576 #endif
577 
578 // we need the Epetra specialization only if Epetra is enabled
579 #if (defined(HAVE_XPETRA_EPETRA) && !defined(XPETRA_EPETRA_NO_64BIT_GLOBAL_INDICES))
580 
581  template <>
582  class CrsGraphFactory<int, long long, EpetraNode> {
583 
584  typedef int LocalOrdinal;
585  typedef long long GlobalOrdinal;
586  typedef EpetraNode Node;
587 
588  private:
590  CrsGraphFactory() {}
591 
592  public:
594  static RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
595  Build (const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rowMap)
596  {
597  XPETRA_MONITOR("CrsMatrixFactory::Build");
598 #ifdef HAVE_XPETRA_TPETRA
599  if (rowMap->lib() == UseTpetra)
600  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node>(rowMap, 0) );
601 #endif
602 #ifdef HAVE_XPETRA_EPETRA
603  if(rowMap->lib() == UseEpetra)
604  return rcp( new EpetraCrsGraphT<long long,Node>(rowMap));
605 #endif
606  XPETRA_FACTORY_END;
607  }
608 
609  static RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
610  Build(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &map, size_t maxNumEntriesPerRow) {
611  XPETRA_MONITOR("CrsGraphFactory::Build");
612 
613 #ifdef HAVE_XPETRA_TPETRA
614  if (map->lib() == UseTpetra)
615  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node> (map, maxNumEntriesPerRow) );
616 #endif
617 
618  if (map->lib() == UseEpetra)
619  return rcp( new EpetraCrsGraphT<long long, Node>(map, maxNumEntriesPerRow) );
620 
621  XPETRA_FACTORY_END;
622  TEUCHOS_UNREACHABLE_RETURN(null);
623  }
624 
625  static Teuchos::RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
626  Build(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 > &plist=Teuchos::null) {
627  XPETRA_MONITOR("CrsGraphFactory::Build");
628 
629 #ifdef HAVE_XPETRA_TPETRA
630  if (rowMap->lib() == UseTpetra)
631  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node>(rowMap, colMap, NumEntriesPerRowToAlloc, plist) );
632 #endif
633 
634  if (rowMap->lib() == UseEpetra)
635  return rcp( new EpetraCrsGraphT<long long, Node>(rowMap, colMap, NumEntriesPerRowToAlloc, plist) );
636 
637  XPETRA_FACTORY_END;
638  TEUCHOS_UNREACHABLE_RETURN(null);
639  }
640 
642  static Teuchos::RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
643  Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap,
644  const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap,
645  size_t maxNumEntriesPerRow,
646  const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null) {
647  XPETRA_MONITOR("CrsGraphFactory::Build");
648 
649 #ifdef HAVE_XPETRA_TPETRA
650  if (rowMap->lib() == UseTpetra)
651  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node>(rowMap, colMap, maxNumEntriesPerRow, plist) );
652 #endif
653  if (rowMap->lib() == UseEpetra)
654  return rcp( new EpetraCrsGraphT<long long,Node>(rowMap, colMap, maxNumEntriesPerRow, plist) );
655 
656  XPETRA_FACTORY_END;
657  TEUCHOS_UNREACHABLE_RETURN(null);
658  }
659 
661  static Teuchos::RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
662  Build(const RCP<const CrsGraph< LocalOrdinal, GlobalOrdinal, Node > >& sourceGraph,
663  const Import< LocalOrdinal, GlobalOrdinal, Node > & importer,
664  const RCP<const Map< LocalOrdinal, GlobalOrdinal, Node >>& domainMap = Teuchos::null,
665  const RCP<const Map< LocalOrdinal, GlobalOrdinal, Node > >& rangeMap = Teuchos::null,
666  const RCP<Teuchos::ParameterList>& params = Teuchos::null) {
667 #ifdef HAVE_XPETRA_TPETRA
668  if (sourceGraph->getRowMap()->lib() == UseTpetra)
669  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node>(sourceGraph, importer, domainMap, rangeMap, params) );
670 #endif
671  if (sourceGraph->getRowMap()->lib() == UseTpetra)
672  return rcp( new EpetraCrsGraphT<long long,Node><LocalOrdinal, GlobalOrdinal, Node>(sourceGraph, importer, domainMap, rangeMap, params) );
673 
674  XPETRA_FACTORY_END;
675  TEUCHOS_UNREACHABLE_RETURN(null);
676  }
677 
678 
679 #ifdef HAVE_XPETRA_KOKKOS_REFACTOR
680 #ifdef HAVE_XPETRA_TPETRA
700  static Teuchos::RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
701  Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap,
702  const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap,
703  const typename Xpetra::CrsGraph<LocalOrdinal, GlobalOrdinal, Node>::local_graph_type::row_map_type& rowPointers,
704  const typename Xpetra::CrsGraph<LocalOrdinal, GlobalOrdinal, Node>::local_graph_type::entries_type::non_const_type& columnIndices,
705  const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null) {
706  XPETRA_MONITOR("CrsMatrixFactory::Build");
707 
708  if (rowMap->lib() == UseTpetra)
709  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node>(rowMap,
710  colMap,
711  rowPointers,
712  columnIndices,
713  plist));
714 
715  XPETRA_FACTORY_END;
716  }
717 
736  static Teuchos::RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
737  Build(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >& rowMap,
738  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >& colMap,
739  const typename Xpetra::CrsGraph<LocalOrdinal, GlobalOrdinal, Node>::local_graph_type& lclGraph,
740  const Teuchos::RCP<Teuchos::ParameterList>& params) {
741  XPETRA_MONITOR("CrsMatrixFactory::Build");
742 
743  if (rowMap->lib() == UseTpetra)
744  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node>(rowMap,
745  colMap,
746  lclGraph,
747  params));
748 
749  XPETRA_FACTORY_END;
750  }
751 
776  static Teuchos::RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
777  Build(const typename Xpetra::CrsGraph<LocalOrdinal, GlobalOrdinal, Node>::local_graph_type& lclGraph,
778  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >& rowMap,
779  const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >& colMap,
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) {
783  XPETRA_MONITOR("CrsMatrixFactory::Build");
784 
785  if (rowMap->lib() == UseTpetra)
786  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node>(lclGraph,
787  rowMap,
788  colMap,
789  domainMap,
790  rangeMap,
791  params));
792 
793  XPETRA_FACTORY_END;
794  }
795 #endif
796 #endif
797 
817  static Teuchos::RCP<CrsGraph<LocalOrdinal, GlobalOrdinal, Node> >
818  Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap,
819  const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap,
820  const Teuchos::ArrayRCP<size_t>& rowPointers,
821  const Teuchos::ArrayRCP<LocalOrdinal>& columnIndices,
822  const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null) {
823  XPETRA_MONITOR("CrsMatrixFactory::Build");
824 
825  if (rowMap->lib() == UseTpetra)
826  return rcp( new TpetraCrsGraph<LocalOrdinal, GlobalOrdinal, Node>(rowMap,
827  colMap,
828  rowPointers,
829  columnIndices,
830  plist));
831 
832  XPETRA_FACTORY_END;
833  }
834 
835  };
836 #endif
837 }
838 
839 #define XPETRA_CRSGRAPHFACTORY_SHORT
840 #endif
Xpetra_ConfigDefs.hpp
XPETRA_MONITOR
#define XPETRA_MONITOR(funcName)
Definition: Xpetra_ConfigDefs.hpp:132
Xpetra_CrsGraph.hpp
Xpetra_EpetraCrsGraph.hpp
Xpetra_Exceptions.hpp
XPETRA_FACTORY_ERROR_IF_EPETRA
#define XPETRA_FACTORY_ERROR_IF_EPETRA(lib)
Definition: Xpetra_Exceptions.hpp:70
XPETRA_FACTORY_END
#define XPETRA_FACTORY_END
Definition: Xpetra_Exceptions.hpp:77
Xpetra::CrsGraphFactory< int, int, EpetraNode >::Build
static Teuchos::RCP< CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap, const Teuchos::ArrayRCP< size_t > &rowPointers, const Teuchos::ArrayRCP< LocalOrdinal > &columnIndices, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null)
Constructor specifying column Map and arrays containing the graph in sorted, local ids.
Definition: Xpetra_CrsGraphFactory.hpp:558
Xpetra::CrsGraphFactory< int, int, EpetraNode >::Build
static Teuchos::RCP< CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > Build(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 > &plist=Teuchos::null)
Constructor specifying column Map and number of entries per row.
Definition: Xpetra_CrsGraphFactory.hpp:383
Xpetra::CrsGraphFactory< int, int, EpetraNode >::CrsGraphFactory
CrsGraphFactory()
Private constructor. This is a static class.
Definition: Xpetra_CrsGraphFactory.hpp:330
Xpetra::CrsGraphFactory< int, int, EpetraNode >::Build
static RCP< CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, size_t maxNumEntriesPerRow)
Definition: Xpetra_CrsGraphFactory.hpp:350
Xpetra::CrsGraphFactory< int, int, EpetraNode >::Build
static RCP< CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > Build(const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap)
Constructor for empty graph (intended use is an import/export target - can't insert entries directly)
Definition: Xpetra_CrsGraphFactory.hpp:335
Xpetra::CrsGraphFactory< int, int, EpetraNode >::Build
static Teuchos::RCP< CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > Build(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 > &plist=Teuchos::null)
Definition: Xpetra_CrsGraphFactory.hpp:366
Xpetra::CrsGraphFactory< int, int, EpetraNode >::Build
static Teuchos::RCP< CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > Build(const RCP< const CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > &sourceGraph, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node >> &domainMap=Teuchos::null, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap=Teuchos::null, const RCP< Teuchos::ParameterList > &params=Teuchos::null)
Constructor using fused import.
Definition: Xpetra_CrsGraphFactory.hpp:403
Xpetra::CrsGraphFactory< int, long long, EpetraNode >::Build
static RCP< CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > Build(const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap)
Constructor for empty graph (intended use is an import/export target - can't insert entries directly)
Definition: Xpetra_CrsGraphFactory.hpp:595
Xpetra::CrsGraphFactory< int, long long, EpetraNode >::Build
static Teuchos::RCP< CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > Build(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 > &plist=Teuchos::null)
Constructor specifying column Map and number of entries per row.
Definition: Xpetra_CrsGraphFactory.hpp:643
Xpetra::CrsGraphFactory< int, long long, EpetraNode >::Build
static Teuchos::RCP< CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > Build(const RCP< const CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > &sourceGraph, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node >> &domainMap=Teuchos::null, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap=Teuchos::null, const RCP< Teuchos::ParameterList > &params=Teuchos::null)
Constructor using fused import.
Definition: Xpetra_CrsGraphFactory.hpp:662
Xpetra::CrsGraphFactory< int, long long, EpetraNode >::Build
static Teuchos::RCP< CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap, const Teuchos::ArrayRCP< size_t > &rowPointers, const Teuchos::ArrayRCP< LocalOrdinal > &columnIndices, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null)
Constructor specifying column Map and arrays containing the graph in sorted, local ids.
Definition: Xpetra_CrsGraphFactory.hpp:818
Xpetra::CrsGraphFactory< int, long long, EpetraNode >::Build
static Teuchos::RCP< CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > Build(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 > &plist=Teuchos::null)
Definition: Xpetra_CrsGraphFactory.hpp:626
Xpetra::CrsGraphFactory< int, long long, EpetraNode >::Build
static RCP< CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, size_t maxNumEntriesPerRow)
Definition: Xpetra_CrsGraphFactory.hpp:610
Xpetra::CrsGraphFactory< int, long long, EpetraNode >::CrsGraphFactory
CrsGraphFactory()
Private constructor. This is a static class.
Definition: Xpetra_CrsGraphFactory.hpp:590
Xpetra::CrsGraphFactory::CrsGraphFactory
CrsGraphFactory()
Private constructor. This is a static class.
Definition: Xpetra_CrsGraphFactory.hpp:71
Xpetra::CrsGraphFactory::Build
static Teuchos::RCP< CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > Build(const RCP< const CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > &sourceGraph, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node >> &domainMap=Teuchos::null, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap=Teuchos::null, const RCP< Teuchos::ParameterList > &params=Teuchos::null)
Constructor using fused import.
Definition: Xpetra_CrsGraphFactory.hpp:143
Xpetra::CrsGraphFactory::Build
static Teuchos::RCP< CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, size_t maxNumEntriesPerRow)
Constructor specifying the number of non-zeros for all rows.
Definition: Xpetra_CrsGraphFactory.hpp:90
Xpetra::CrsGraphFactory::Build
static Teuchos::RCP< CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > Build(const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap)
Constructor for empty graph (intended use is an import/export target - can't insert entries directly)
Definition: Xpetra_CrsGraphFactory.hpp:76
Xpetra::CrsGraphFactory::Build
static Teuchos::RCP< CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > Build(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 > &plist=Teuchos::null)
Constructor specifying column Map and number of entries per row.
Definition: Xpetra_CrsGraphFactory.hpp:105
Xpetra::CrsGraphFactory::Build
static Teuchos::RCP< CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap, const Teuchos::ArrayRCP< size_t > &rowPointers, const Teuchos::ArrayRCP< LocalOrdinal > &columnIndices, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null)
Constructor specifying column Map and arrays containing the graph in sorted, local ids.
Definition: Xpetra_CrsGraphFactory.hpp:299
Xpetra::CrsGraphFactory::Build
static Teuchos::RCP< CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > Build(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 > &plist=Teuchos::null)
Constructor specifying column Map and number of entries in each row.
Definition: Xpetra_CrsGraphFactory.hpp:124
TEUCHOS_TEST_FOR_EXCEPTION
#define TEUCHOS_TEST_FOR_EXCEPTION(throw_exception_test, Exception, msg)
TEUCHOS_UNREACHABLE_RETURN
#define TEUCHOS_UNREACHABLE_RETURN(dummyReturnVal)
rcp
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
Xpetra
Xpetra namespace
Definition: Xpetra_BlockedCrsMatrix.hpp:89