metis_cpu.cpp 1.15 KB
Newer Older
rusty1s's avatar
rusty1s committed
1
#include "metis_cpu.h"
rusty1s's avatar
rusty1s committed
2

rusty1s's avatar
rusty1s committed
3
#ifdef WITH_METIS
rusty1s's avatar
rusty1s committed
4
#include <metis.h>
rusty1s's avatar
rusty1s committed
5
#endif
rusty1s's avatar
rusty1s committed
6
7
8

#include "utils.h"

9
torch::Tensor partition_cpu(torch::Tensor rowptr, torch::Tensor col,int64_t num_parts,
10
                            torch::optional<torch::Tensor> edge_wgt,  bool recursive) {
rusty1s's avatar
rusty1s committed
11
#ifdef WITH_METIS
rusty1s's avatar
rusty1s committed
12
13
14
15
  CHECK_CPU(rowptr);
  CHECK_CPU(col);

  int64_t nvtxs = rowptr.numel() - 1;
rusty1s's avatar
rusty1s committed
16

17
  auto part = torch::empty(nvtxs, rowptr.options());
rusty1s's avatar
rusty1s committed
18
19
  auto *xadj = rowptr.data_ptr<int64_t>();
  auto *adjncy = col.data_ptr<int64_t>();
20
21
22
23
24
25
26
27

  int64_t * adjwgt;
  if (edge_wgt.has_value()){
      adjwgt=edge_wgt.value().data_ptr<int64_t>();
      adjwgt=(idx_t*) adjwgt;
  }else{
       adjwgt=nullptr;
  }
bowendeng's avatar
bowendeng committed
28

rusty1s's avatar
rusty1s committed
29
30
31
32
  int64_t ncon = 1;
  int64_t objval = -1;
  auto part_data = part.data_ptr<int64_t>();

rusty1s's avatar
rusty1s committed
33
  if (recursive) {
bwdeng20's avatar
bwdeng20 committed
34
    METIS_PartGraphRecursive(&nvtxs, &ncon, xadj, adjncy, NULL, NULL, adjwgt,
rusty1s's avatar
rusty1s committed
35
36
                             &num_parts, NULL, NULL, NULL, &objval, part_data);
  } else {
bwdeng20's avatar
bwdeng20 committed
37
    METIS_PartGraphKway(&nvtxs, &ncon, xadj, adjncy, NULL, NULL, adjwgt,
rusty1s's avatar
rusty1s committed
38
39
                        &num_parts, NULL, NULL, NULL, &objval, part_data);
  }
rusty1s's avatar
rusty1s committed
40
41

  return part;
rusty1s's avatar
rusty1s committed
42
43
44
#else
  AT_ERROR("Not compiled with METIS support");
#endif
bowendeng's avatar
bowendeng committed
45
}