Ifpack2 Templated Preconditioning Package Version 1.0
Loading...
Searching...
No Matches
Ifpack2_LinearPartitioner_def.hpp
1/*@HEADER
2// ***********************************************************************
3//
4// Ifpack2: Templated Object-Oriented Algebraic Preconditioner Package
5// Copyright (2009) Sandia Corporation
6//
7// Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
8// license for use of this work by or on behalf of the U.S. Government.
9//
10// Redistribution and use in source and binary forms, with or without
11// modification, are permitted provided that the following conditions are
12// met:
13//
14// 1. Redistributions of source code must retain the above copyright
15// notice, this list of conditions and the following disclaimer.
16//
17// 2. Redistributions in binary form must reproduce the above copyright
18// notice, this list of conditions and the following disclaimer in the
19// documentation and/or other materials provided with the distribution.
20//
21// 3. Neither the name of the Corporation nor the names of the
22// contributors may be used to endorse or promote products derived from
23// this software without specific prior written permission.
24//
25// THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
26// EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
27// IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
28// PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
29// CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
30// EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
31// PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
32// PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
33// LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
34// NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
35// SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
36//
37// Questions? Contact Michael A. Heroux (maherou@sandia.gov)
38//
39// ***********************************************************************
40//@HEADER
41*/
42
43#ifndef IFPACK2_LINEAR_PARTITIONER_DEF_HPP
44#define IFPACK2_LINEAR_PARTITIONER_DEF_HPP
45#include "Ifpack2_ConfigDefs.hpp"
46#include "Ifpack2_LinearPartitioner_decl.hpp"
47
48namespace Ifpack2 {
49
50//==============================================================================
51// Constructor
52template<class GraphType>
54LinearPartitioner (const Teuchos::RCP<const row_graph_type>& graph) :
55 OverlappingPartitioner<GraphType> (graph)
56{}
57
58
59template<class GraphType>
61
62
63template<class GraphType>
64void
66setPartitionParameters (Teuchos::ParameterList& /* List */) {}
67
68
69template<class GraphType>
71{
72 using Teuchos::as;
73 // Partition_ is an array of local_ordinal_type. local_ordinal_type
74 // may be signed or unsigned. NumLocalParts_ is int, and needs to
75 // be signed, since negative values are significant. Comparisons
76 // between signed and unsigned integers often result in compiler
77 // warnings, which is why we use as() for explicit conversions
78 // below. We also use as() because in a debug build, it checks for
79 // overflow.
80 const int mod = as<int> (this->Graph_->getLocalNumRows () /
81 this->NumLocalParts_);
82 for (size_t i = 0; i < this->Graph_->getLocalNumRows (); ++i) {
83 this->Partition_[i] = as<local_ordinal_type> (i / mod);
84 if (this->Partition_[i] >= as<local_ordinal_type> (this->NumLocalParts_)) {
85 this->Partition_[i] = this->NumLocalParts_ - 1;
86 }
87 }
88}
89
90
91}// namespace Ifpack2
92
93#define IFPACK2_LINEARPARTITIONER_INSTANT(LO,GO,N) \
94 template class Ifpack2::LinearPartitioner<Tpetra::CrsGraph< LO, GO, N > >; \
95 template class Ifpack2::LinearPartitioner<Tpetra::RowGraph< LO, GO, N > >;
96
97#endif // IFPACK2_LINEARPARTITIONER_DEF_HPP
virtual ~LinearPartitioner()
Destructor.
Definition: Ifpack2_LinearPartitioner_def.hpp:60
LinearPartitioner(const Teuchos::RCP< const row_graph_type > &graph)
Constructor.
Definition: Ifpack2_LinearPartitioner_def.hpp:54
void setPartitionParameters(Teuchos::ParameterList &List)
Set the partitioner's parameters (none for linear partitioning).
Definition: Ifpack2_LinearPartitioner_def.hpp:66
void computePartitions()
Compute the partitions.
Definition: Ifpack2_LinearPartitioner_def.hpp:70
Create overlapping partitions of a local graph.
Definition: Ifpack2_OverlappingPartitioner_decl.hpp:78
Preconditioners and smoothers for Tpetra sparse matrices.
Definition: Ifpack2_AdditiveSchwarz_decl.hpp:74