reduce.h 1.34 KB
Newer Older
1
2
3
4
5
6
7
8
9
// Copyright (c) Microsoft Corporation.
// Licensed under the MIT License.
#pragma once

#include "common.h"

namespace tl {

struct SumOp {
10
  template <typename T> TL_DEVICE T operator()(T const &x, T const &y) {
11
12
13
14
15
    return x + y;
  }
};

struct MaxOp {
16
  template <typename T> TL_DEVICE T operator()(T const &x, T const &y) {
17
18
19
20
21
    return ck_tile::max(x, y);
  }
};

struct MinOp {
22
  template <typename T> TL_DEVICE T operator()(T const &x, T const &y) {
23
24
25
26
    return ck_tile::min(x, y);
  }
};

27
28
29
30
template <class Reducer, int threads, int scale> struct AllReduce {
  static_assert(threads == 1024 || threads == 512 || threads == 256 ||
                threads == 128 || threads == 64 || threads == 32 ||
                threads == 16 || threads == 8 || threads == 4 || threads == 2);
31
32
  static_assert(threads % scale == 0);

33
  template <typename T> static __device__ T run(T x, T *red_buf = nullptr) {
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
    constexpr int offset = threads / 2;
    constexpr int warpSize = 64;

    if constexpr (offset >= warpSize) {
      __syncthreads();
      red_buf[threadIdx.x] = x;
      __syncthreads();
      x = Reducer()(x, red_buf[threadIdx.x ^ offset]);
    } else {
      x = Reducer()(x, __shfl_xor(x, offset));
    }
    if constexpr (offset == scale) {
      return x;
    } else {
      return AllReduce<Reducer, offset, scale>::run(x, red_buf);
    }
  }
};

53
} // namespace tl