Subversion Repositories gng

/*

 * Generate inverse of the normal cumulative distribution function

 *

 * The calling syntax is:

 *     x = icdf_gen(r)

 *

 * Input:

 *     z: pseudorandom numbers

 *

 * Output:

 *     x: Gaussian random numbers, n x 1 vector of 32-bit integer (s<16,11>)

 */

/*

 * Copyright (C) 2014, Guangxi Liu <guangxi.liu@opencores.org>

 *

 * This source file may be used and distributed without restriction provided

 * that this copyright statement is not removed from the file and that any

 * derivative work contains the original copyright notice and the associated

 * disclaimer.

 *

 * This source file is free software; you can redistribute it and/or modify it

 * under the terms of the GNU Lesser General Public License as published by

 * the Free Software Foundation; either version 2.1 of the License,

 * or (at your option) any later version.

 *

 * This source is distributed in the hope that it will be useful, but

 * WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY

 * or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public

 * License for more details.

 *

 * You should have received a copy of the GNU Lesser General Public License

 * along with this source; if not, download it from

 * http://www.opencores.org/lgpl.shtml

 */

#include "mex.h"

#include "icdf.h"

void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])

{

    /* Local variables */

    unsigned long long *r;

    unsigned len;

    long *x;

    unsigned i;

    /* Check for proper number of arguments */

    if (nrhs != 1)

        mexErrMsgTxt("One input arguments required.");

    if (!mxIsUint64(prhs[0]) || (mxGetN(prhs[0]) != 1))

        mexErrMsgTxt("icdf_gen requires that r be a column vector of unsigned 64-bit integer.");

    len = mxGetM(prhs[0]);

    if (len < 1)

        mexErrMsgTxt("icdf_gen requires that vector length be a positive integer.");

    /* Create a vector for the return argument */

    plhs[0] = mxCreateNumericMatrix(len, 1, mxINT32_CLASS, 0);

    /* Assign pointers to the various parameters */

    r = (unsigned long long*)mxGetData(prhs[0]);

    x = (long*)mxGetData(plhs[0]);

    /* Call function icdf and assign return value */

    for (i = 0; i < len; i++)

        x[i] = icdf(r[i]);

}