random.hpp

Go to the documentation of this file.

//    Copyright 2023 Jij Inc.
 
//    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.
 
#pragma once
 
#include <climits>
#include <random>
 
#ifdef USE_CUDA
 
#include <utility>
#include <cassert>
 
#include <cuda_runtime.h>
#include <curand.h>
 
#include "openjij/utility/gpu/handle_error.hpp"
#include "openjij/utility/gpu/memory.hpp"
 
#endif
 
namespace openjij {
namespace utility {
 
class Xorshift {
public:
  using result_type = uint_fast32_t;
 
  inline static constexpr unsigned min() { return 0u; }
 
  inline static constexpr unsigned max() { return UINT_MAX; }
 
  inline unsigned operator()() {
    unsigned t = x ^ (x << 11);
    x = y;
    y = z;
    z = w;
    return w = (w ^ (w >> 19)) ^ (t ^ (t >> 8));
  }
 
  Xorshift() {
    std::random_device rd;
    w = rd();
  }
 
  Xorshift(unsigned s) { w = s; }
 
private:
  unsigned x = 123456789u, y = 362436069u, z = 521288629u, w;
};
 
#ifdef USE_CUDA
namespace cuda {
template <typename FloatType>
inline curandStatus_t curand_generate_uniform_impl(curandGenerator_t generator,
                                                   FloatType *outputPtr,
                                                   size_t num) {
  static_assert(std::is_floating_point<FloatType>::value,
                "FloatType must be float or double.");
  static_assert(!std::is_same<FloatType, long double>::value,
                "long double is not supported");
  return CURAND_STATUS_SUCCESS;
}
 
template <>
inline curandStatus_t
curand_generate_uniform_impl<float>(curandGenerator_t generator,
                                    float *outputPtr, size_t num) {
  return curandGenerateUniform(generator, outputPtr, num);
}
 
template <>
inline curandStatus_t
curand_generate_uniform_impl<double>(curandGenerator_t generator,
                                     double *outputPtr, size_t num) {
  return curandGenerateUniformDouble(generator, outputPtr, num);
}
 
template <typename FloatType, curandRngType_t rng_type> class CurandWrapper {
public:
  CurandWrapper(std::uint64_t seed) {
    // generate psudo_random_number generator
    HANDLE_ERROR_CURAND(curandCreateGenerator(&_rng, rng_type));
    // set seed
    HANDLE_ERROR_CURAND(curandSetPseudoRandomGeneratorSeed(_rng, seed));
  }
 
  CurandWrapper() : CurandWrapper(std::random_device{}()) {}
 
  CurandWrapper(CurandWrapper &&obj) noexcept {
    // move curand handler
    this->_rng = obj._rng;
    obj._rng = NULL;
  }
 
  ~CurandWrapper() {
    // destroy generator
    if (_rng != NULL)
      HANDLE_ERROR_CURAND(curandDestroyGenerator(_rng));
  }
 
  inline void generate_uniform(std::size_t n,
                               cuda::unique_dev_ptr<FloatType[]> &dev_random) {
    HANDLE_ERROR_CURAND(
        curand_generate_uniform_impl(_rng, dev_random.get(), n));
  }
 
private:
  curandGenerator_t _rng;
};
} // namespace cuda
#endif
} // namespace utility
} // namespace openjij