mmc_rand_logistic.h

Go to the documentation of this file.

/***************************************************************************/
/***************************************************************************/
#ifndef _MCEXTREME_LOGISTIC_RAND_H
#define _MCEXTREME_LOGISTIC_RAND_H

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

#define __device__ static inline

#ifndef DOUBLE_PREC_LOGISTIC
    typedef float RandType;
#else
    typedef double RandType;
#endif

#define MCX_RNG_NAME       "Logistic-Lattice"

#define RAND_BUF_LEN       5        //register arrays
#define RAND_SEED_WORD_LEN      5        //32bit seed length (32*5=160bits)

__device__ void rand_need_more(RandType t[RAND_BUF_LEN], RandType tbuf[RAND_BUF_LEN]);
__device__ void logistic_init(RandType* t, RandType* tnew, uint seed[], uint idx);

__device__ RandType rand_uniform01(RandType v);
__device__ void rng_init(RandType t[RAND_BUF_LEN], RandType tnew[RAND_BUF_LEN], uint* n_seed, int idx);

// generate [0,1] random number for the next scattering length
__device__ float rand_next_scatlen(RandType t[RAND_BUF_LEN]);
// generate [0,1] random number for the next arimuthal angle
__device__ float rand_next_aangle(RandType t[RAND_BUF_LEN]);
// generate random number for the next zenith angle
__device__ float rand_next_zangle(RandType t[RAND_BUF_LEN]);
__device__ float rand_next_reflect(RandType t[RAND_BUF_LEN]);
__device__ float rand_do_roulette(RandType t[RAND_BUF_LEN]);

#ifdef MMC_USE_SSE_MATH
    __device__ void rand_next_aangle_sincos(RandType t[RAND_BUF_LEN], float* si, float* co);
    __device__ float rand_next_scatlen_ps(RandType t[RAND_BUF_LEN]);
#endif

#endif