Subversion Repositories gng

/*
 * Generate Combined Tausworthe Generator seed
 *
 * The calling syntax is:
 *     z = ctg_seed(s)
 *
 * Input:
 *     s: seed, unsigned 32-bit integer
 *
 * Output:
 *     z: initial internal state, 3 x 1 vector of unsigned 64-bit integer
 */
 
#include "mex.h"
#include "taus176.h"
 
 
void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])
{
    /* Local variables */
    taus_state_t state;
    double *s;
    unsigned seed;
    unsigned long long *z;
 
    /* Check for proper number of arguments */
    if (nrhs != 1)
        mexErrMsgTxt("One input arguments required.");
    else if (nlhs > 1)
        mexErrMsgTxt("Too many output arguments.");
 
    if (!mxIsDouble(prhs[0]) || mxIsComplex(prhs[0]) ||
            (mxGetM(prhs[0]) != 1) || (mxGetN(prhs[0]) != 1))
        mexErrMsgTxt("ctg_seed requires that s be a unsigned 32-bit integer.");
 
    /* Create a vector for the return argument */
    plhs[0] = mxCreateNumericMatrix(3, 1, mxUINT64_CLASS, 0);
 
    /* Assign pointers to the various parameters */
    s = (unsigned *)mxGetPr(prhs[0]);
    seed = (unsigned)(*s);
    z = (unsigned long long*)mxGetData(plhs[0]);
 
    /* Call function taus_set */
    taus_set(&state, seed);
 
    /* Assign return value */
    z[0] = state.z1;
    z[1] = state.z2;
    z[2] = state.z3;
}