DocsTutorial/example1_kernel.cl

Device code for Example 1: Monte Carlo

See also

DocsTutorial/example1.c

/*
 * Authors:
 *
 *   David Munger <mungerd@iro.umontreal.ca>        (2015)
 *   Pierre L'Ecuyer <lecuyer@iro.umontreal.ca>     (2015)
 *
 */

#define CLRNG_ENABLE_SUBSTREAMS
#include <clRNG/mrg31k3p.clh>

#define StreamType     clrngMrg31k3pStream
#define nextCoordinate clrngMrg31k3pRandomU01
#include <clQMC/DocsTutorial/common.clh>


__kernel void simulateWithMC(
        __global const clrngMrg31k3pHostStream* streams,
        uint points_per_work_item,
        __global clqmc_fptype* out)
{
  uint gsize = get_global_size(0);
  uint gid   = get_global_id(0);

  clrngMrg31k3pStream stream;
  clrngMrg31k3pCopyOverStreamsFromGlobal(1, &stream, &streams[gid]);

  clqmc_fptype sum = 0.0;

  for (uint i = 0; i < points_per_work_item; i++) {
    sum += simulateOneRun(&stream);
    clrngMrg31k3pForwardToNextSubstreams(1, &stream);
  }

  out[gid] = sum / points_per_work_item;
}

/*
vim: ft=c
*/