81 const RCP<const typename Adapter::base_adapter_t> &adapter__,
82 const RCP<Teuchos::ParameterList> &pl__,
83 const RCP<
const Teuchos::Comm<int> > &comm__,
84 RCP<const Environment> &env__,
86 ) : adapter(adapter__), pl(pl__), comm(comm__), env(env__), graphFlags(graphFlags__)
96 int localOrder(
const RCP<LocalOrderingSolution<lno_t> > &solution)
102 lno_t *perm = solution->getPermutationView();
105 std::cerr <<
"perm is NULL" << std::endl;
110 const auto model = rcp(
new GraphModel<Adapter>(adapter, env, comm, graphFlags));
111 const size_t nVtx = model->getLocalNumVertices();
112 ArrayView<const gno_t> edgeIds;
113 ArrayView<const offset_t> offsets;
114 ArrayView<StridedData<lno_t, scalar_t> > wgts;
115 model->getEdgeList(edgeIds, offsets, wgts);
119 Teuchos::Array<std::pair<lno_t, offset_t> > degrees(nVtx);
121 degrees[i].first = i;
122 degrees[i].second = offsets[i+1] - offsets[i];
126 SortPairs<lno_t,offset_t> zort;
128 zort.sort(degrees, inc);
132 perm[i] = degrees[i].first;
135 solution->setHavePerm(
true);