Subversion Repositories gelsvn

/* ----------------------------------------------------------------------- *

 * This file is part of GEL, http://www.imm.dtu.dk/GEL

 * Copyright (C) the authors and DTU Informatics

 * For license and list of authors, see ../../doc/intro.pdf

 * ----------------------------------------------------------------------- */

#include "CGLA.h"

namespace

{

  unsigned int seed = 1;

  unsigned int current_rand = 1;

}

namespace CGLA

{

  void gel_srand(unsigned int s)

  {

    seed = current_rand = s;

  }

  unsigned int gel_rand(unsigned int k)

  {

    unsigned int b = 3125;

    unsigned int c = 49;

    unsigned int result = seed;

    for (;k > 0;k>>=1)

    {

      if (k & 1) result = result * b + c;

      c += b * c;

      b *= b;

    }

    return result;

  }

  unsigned int gel_rand()

  {

    const unsigned int b = 3125;

    const unsigned int c = 49;

    current_rand = current_rand*b + c;

    return current_rand;

  }

}