[c-prog] help on translating class

[Jos Timanta Tarigan]

Hi, 

i need help translating a class, i hope you guys dont mind i put the whole code 
in here. honestly speaking i dont understand what it did from top to bottom but 
i know the intention is to build a random value. the problem is i dont know if 
this random has some kind of constraint or its pure random.

so im trying to build a ray tracing program. one of the purpose is to create a 
random value to create the ray, some kind of random sampling. if you follow 
this link:

http://www.devmaster.net/articles/raytracing_series/part5.php

you'll find the explanation about the function used in one of the sniplet (the 
diffuse reflection subtopic). its the Twister.cpp file. all i need to know is 
if the value generated to some value or its pure random so i can replace it by 
c++ random function or there are some kind of constraints to the random value

hope i did not ask too much :) many thanks in advance



here is the code in the class



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

#include"twister.h"

namespace Raytracer {

unsigned short mtRand_xsubi[3] = { 723, 32761, 44444 };

#define M 397
#define MATRIX_A 0x9908b0df
#define UPPER_MASK 0x80000000
#define LOWER_MASK 0x7fffffff

#define TEMPERING_MASK_B 0x9d2c5680
#define TEMPERING_MASK_C 0xefc60000
#define TEMPERING_SHIFT_U(y)  (y >> 11)
#define TEMPERING_SHIFT_S(y)  (y << 7)
#define TEMPERING_SHIFT_T(y)  (y << 15)
#define TEMPERING_SHIFT_L(y)  (y >> 18)

void Twister::Seed( unsigned long seed )
{
    mt[0]= seed & 0xffffffff;
    for ( mti = 1; mti < mtRand_N; mti++ ) mt[mti] = (69069 * mt[mti - 1]) & 
0xffffffff;
    unsigned long s = 373737;
    for ( mti = 1; mti < mtRand_N; mti++ )
    {
mt[mti] ^= s;
s = s * 5531 + 81547;
s ^= (s >> 9) ^ (s << 19);
    }
}

float Twister::Rand()
{
    unsigned long y;
    static unsigned long mag01[2] = {0x0, MATRIX_A};
    if (mti >= mtRand_N) 
{
        int kk;
        for (kk=0;kk<mtRand_N-M;kk++) 
{
            y = (mt[kk]&UPPER_MASK)|(mt[kk+1]&LOWER_MASK);
            mt[kk] = mt[kk+M] ^ (y >> 1) ^ mag01[y & 0x1];
        }
        for (;kk<mtRand_N-1;kk++) 
{
            y = (mt[kk]&UPPER_MASK)|(mt[kk+1]&LOWER_MASK);
            mt[kk] = mt[kk+(M-mtRand_N)] ^ (y >> 1) ^ mag01[y & 0x1];
        }
        y = (mt[mtRand_N-1]&UPPER_MASK)|(mt[0]&LOWER_MASK);
        mt[mtRand_N-1] = mt[M-1] ^ (y >> 1) ^ mag01[y & 0x1];
        mti = 0;
    }
    y = mt[mti++];
    y ^= TEMPERING_SHIFT_U(y);
    y ^= TEMPERING_SHIFT_S(y) & TEMPERING_MASK_B;
    y ^= TEMPERING_SHIFT_T(y) & TEMPERING_MASK_C;
    y ^= TEMPERING_SHIFT_L(y);
    return ( (float)y * 2.3283064370807974e-10f );
}

unsigned long Twister::RandL()
{
    unsigned long y;
    static unsigned long mag01[2] = { 0x0, MATRIX_A };
    if (mti >= mtRand_N) 
{
        int kk;
        for (kk=0;kk<mtRand_N-M;kk++) 
{
            y = (mt[kk]&UPPER_MASK)|(mt[kk+1]&LOWER_MASK);
            mt[kk] = mt[kk+M] ^ (y >> 1) ^ mag01[y & 0x1];
        }
        for (;kk<mtRand_N-1;kk++) 
{
            y = (mt[kk]&UPPER_MASK)|(mt[kk+1]&LOWER_MASK);
            mt[kk] = mt[kk+(M-mtRand_N)] ^ (y >> 1) ^ mag01[y & 0x1];
        }
        y = (mt[mtRand_N-1]&UPPER_MASK)|(mt[0]&LOWER_MASK);
        mt[mtRand_N-1] = mt[M-1] ^ (y >> 1) ^ mag01[y & 0x1];
        mti = 0;
    }
    y = mt[mti++];
    y ^= TEMPERING_SHIFT_U(y);
    y ^= TEMPERING_SHIFT_S(y) & TEMPERING_MASK_B;
    y ^= TEMPERING_SHIFT_T(y) & TEMPERING_MASK_C;
    y ^= TEMPERING_SHIFT_L(y);
    return y;
}

}; // namespace Raytracer
 ================================= 
http://www.svnstrk.blogspot.com



      

[Non-text portions of this message have been removed]