PhotosRandom.h

#ifndef _PhotosRandom_included_
#define _PhotosRandom_included_

/**
 * @class PhotosRandom
 *
 * @brief Photos random number generator rewritten from PHOTOS FORTRAN
 *
 * Generates uniformly distributed random numbers between 0 and 1.
 * Must be initialized by call to PhotosRandom::initialize().
 * Original authors: B. van Eijk, F. James, G. Marsaglia and A. Zaman
 *
 * @author Tomasz Przedzinski
 * @date 18th October 2010
 */

namespace Photospp
{

class PhotosRandom
{
public:
        /* Change the seed. Default is s1=1802 and s2=9373
           These values must be in range [0,31327] and [0,30080] respectively */
        static void setSeed(int s1,int s2);

        /* Initialization routine. Must be called at least once before
           the generator can be used. */
        static void initialize();

        /* Uniform distribution between 0 and 1 */
        static double randomReal();

protected:
        static bool init;
        static int iseed[2];
        static int i97;
        static int j97;
        static double uran[97];
        static double cran;
        static const double cdran;
        static const double cmran;
};

} // namespace Photospp
#endif