Zoltan2
Loading...
Searching...
No Matches
Zoltan2_AlgNatural.hpp
Go to the documentation of this file.
1// @HEADER
2//
3// ***********************************************************************
4//
5// Zoltan2: A package of combinatorial algorithms for scientific computing
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 Karen Devine (kddevin@sandia.gov)
39// Erik Boman (egboman@sandia.gov)
40// Siva Rajamanickam (srajama@sandia.gov)
41//
42// ***********************************************************************
43//
44// @HEADER
45#ifndef _ZOLTAN2_ALGNATURAL_HPP_
46#define _ZOLTAN2_ALGNATURAL_HPP_
47
48#include <Zoltan2_Algorithm.hpp>
51
52
53namespace Zoltan2{
54
59
60template <typename Adapter>
61class AlgNatural : public Algorithm<Adapter>
62{
63 private:
64
65 const RCP<const typename Adapter::base_adapter_t> adapter;
66 const RCP<Teuchos::ParameterList> pl;
67 const RCP<const Teuchos::Comm<int> > comm;
68 const RCP<const Environment> env;
69
70 public:
71
72 typedef typename Adapter::lno_t lno_t;
73 typedef typename Adapter::gno_t gno_t;
74
76 const RCP<const typename Adapter::base_adapter_t> &adapter__,
77 const RCP<Teuchos::ParameterList> &pl__,
78 const RCP<const Teuchos::Comm<int> > &comm__,
79 const RCP<const Environment> &env__
80 ) : adapter(adapter__), pl(pl__), comm(comm__), env(env__)
81 {
82 }
83
84 int globalOrder(const RCP<GlobalOrderingSolution<gno_t> > &/* solution */) {
85 throw std::logic_error("AlgNatural does not yet support global ordering.");
86 }
87
88 int localOrder(const RCP<LocalOrderingSolution<lno_t> > &solution)
89 {
90
91 int ierr= 0;
92
93 HELLO;
94
95 // Local permutation only for now.
96
97 // Set identity permutation.
98 modelFlag_t modelFlags;
99 const auto model = rcp(new IdentifierModel<Adapter>(adapter, env, comm, modelFlags));
100 const size_t n = model->getLocalNumIdentifiers();
101 lno_t *perm = solution->getPermutationView();
102 if (perm){
103 for (size_t i=0; i<n; i++){
104 perm[i] = i;
105 }
106 }
107 else
108 // TODO: throw exception?
109 ierr = -1;
110
111 solution->setHavePerm(true);
112 return ierr;
113
114 }
115
116};
117}
118#endif
Defines the IdentifierModel interface.
Defines the OrderingSolution class.
#define HELLO
int localOrder(const RCP< LocalOrderingSolution< lno_t > > &solution)
Ordering method.
int globalOrder(const RCP< GlobalOrderingSolution< gno_t > > &)
Ordering method.
AlgNatural(const RCP< const typename Adapter::base_adapter_t > &adapter__, const RCP< Teuchos::ParameterList > &pl__, const RCP< const Teuchos::Comm< int > > &comm__, const RCP< const Environment > &env__)
Algorithm defines the base class for all algorithms.
Created by mbenlioglu on Aug 31, 2020.
std::bitset< NUM_MODEL_FLAGS > modelFlag_t