uniform_kernel.cc 2.4 KB
Newer Older
1 2 3 4 5 6 7 8 9 10 11 12 13 14
// Copyright (c) 2022 PaddlePaddle Authors. All Rights Reserved.
//
// Licensed under the Apache License, Version 2.0 (the "License");
// you may not use this file except in compliance with the License.
// You may obtain a copy of the License at
//
//     http://www.apache.org/licenses/LICENSE-2.0
//
// Unless required by applicable law or agreed to in writing, software
// distributed under the License is distributed on an "AS IS" BASIS,
// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
// See the License for the specific language governing permissions and
// limitations under the License.

15
#include "paddle/phi/kernels/uniform_kernel.h"
16

17
#include "paddle/phi/core/kernel_registry.h"
18
#include "paddle/phi/kernels/funcs/uniform_real_distribution.h"
19 20 21 22

namespace phi {

template <typename T, typename Context>
23 24 25 26 27 28 29 30 31 32
void UniformRawKernel(const Context &dev_ctx,
                      const IntArray &shape,
                      DataType dtype,
                      const Scalar &min,
                      const Scalar &max,
                      int seed,
                      int diag_num,
                      int diag_step,
                      float diag_val,
                      DenseTensor *out) {
33 34 35 36 37 38 39 40 41 42
  out->Resize(phi::make_ddim(shape.GetData()));
  T *data = dev_ctx.template Alloc<T>(out);
  auto size = out->numel();
  std::shared_ptr<std::mt19937_64> engine;
  if (seed) {
    engine = std::make_shared<std::mt19937_64>();
    engine->seed(seed);
  } else {
    engine = dev_ctx.GetGenerator()->GetCPUEngine();
  }
43 44
  UniformRealDistribution<T>(
      data, size, min.to<float>(), max.to<float>(), engine);
45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65
  if (diag_num > 0) {
    PADDLE_ENFORCE_GT(
        size,
        (diag_num - 1) * (diag_step + 1),
        phi::errors::InvalidArgument(
            "ShapeInvalid: the diagonal's elements is equal (num-1) "
            "* (step-1) with num %d, step %d,"
            "It should be smaller than %d, but received %d",
            diag_num,
            diag_step,
            (diag_num - 1) * (diag_step + 1),
            size));
    for (int64_t i = 0; i < diag_num; ++i) {
      int64_t pos = i * diag_step + i;
      data[pos] = diag_val;
    }
  }
}

}  // namespace phi

66
PD_REGISTER_KERNEL(uniform_raw,
67 68
                   CPU,
                   ALL_LAYOUT,
69
                   phi::UniformRawKernel,
70 71 72
                   float,
                   double,
                   phi::dtype::bfloat16) {}