hpc_random.h

// -*-c++-*-
#ifndef STORMM_HPC_RANDOM_H
#define STORMM_HPC_RANDOM_H
 
#include "copyright.h"
#include "Accelerator/gpu_details.h"
#include "Accelerator/hybrid.h"
#include "DataTypes/stormm_vector_types.h"
#include "Random/random.h"
 
namespace stormm {
namespace random {
 
using card::GpuDetails;
using card::Hybrid;
  
void initXoroshiro128pArray(Hybrid<ullint2> *state_vector, int igseed, int scrub_cycles,
                            const GpuDetails &gpu);

void initXoshiro256ppArray(Hybrid<ullint2> *state_xy, Hybrid<ullint2> *state_zw, int igseed,
                           int scrub_cycles, const GpuDetails &gpu);

void fillRandomCache(ullint2* state_xy, ullint2* state_zw, double* cache, size_t length,
                     size_t depth, RandomAlgorithm method, RandomNumberKind product,
                     size_t index_start, size_t index_end, const GpuDetails &gpu);
 
void fillRandomCache(Hybrid<ullint2> *state_xy, Hybrid<ullint2> *state_zw, Hybrid<double> *cache,
                     size_t length, size_t depth, RandomAlgorithm method, RandomNumberKind product,
                     size_t index_start, size_t index_end, const GpuDetails &gpu);
 
void fillRandomCache(ullint2* state_xy, ullint2* state_zw, float* cache, size_t length,
                     size_t depth, RandomAlgorithm method, RandomNumberKind product,
                     size_t index_start, size_t index_end, const GpuDetails &gpu);
 
void fillRandomCache(Hybrid<ullint2> *state_xy, Hybrid<ullint2> *state_zw, Hybrid<float> *cache,
                     size_t length, size_t depth, RandomAlgorithm method, RandomNumberKind product,
                     size_t index_start, size_t index_end, const GpuDetails &gpu);
 
} // namespace random
} // namespace stormm
 
#endif