void UncoupledAggregationAlgorithm<LocalOrdinal, GlobalOrdinal, Node, LocalMatOps>::RandomReorder(ArrayRCP<LO> list) const {
   //TODO: replace int
   int n = list.size();
   for(int i = 0; i < n-1; i++) {
     std::swap(list[i], list[RandomOrdinal(i,n-1)]);
   }
 }
 void AggregationPhase1Algorithm_kokkos<LocalOrdinal, GlobalOrdinal, Node>::RandomReorder(ArrayRCP<LO> list) const {
   //TODO: replace int
   int n = list.size();
   for(int i = 0; i < n-1; i++)
     std::swap(list[i], list[RandomOrdinal(i,n-1)]);
 }