Changeset 2433 for trunk/psLib/src/math/psRandom.c

Timestamp:

Nov 24, 2004, 11:59:43 AM (22 years ago)

Author:

gusciora

Message:

Changed the order of the arguments for the polynomial eval functions.

File:

• trunk/psLib/src/math/psRandom.c (modified) (4 diffs)

trunk/psLib/src/math/psRandom.c

@@ -10 +10 @@
 *  @author GLG, MHPCC
 *
-*  @version $Revision: 1.1 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2004-11-24 20:55:01 $
+*  @version $Revision: 1.2 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2004-11-24 21:59:43 $
 *
 *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -20 +20 @@
 #include <float.h>
 #include <math.h>
+#include <gsl/gsl_rng.h>
+#include <gsl/gsl_randist.h>
 
+#include "psMemory.h"
 #include "psRandom.h"
 #include "psScalar.h"
+#include "psError.h"
+#include "psTrace.h"
+#include "psLogMsg.h"
+#include "psConstants.h"
+#include "psDataManipErrors.h"
 
 /*
@@ -39 +47 @@
                         psU64 seed)
 {
+    gsl_rng *r = NULL;
+    psRandom *myRNG = (psRandom *) psAlloc(sizeof(psRandom));
+
+    switch (type) {
+    case PS_RANDOM_TAUS:
+        r = gsl_rng_alloc(gsl_rng_taus);
+        myRNG->gsl = r;
+        return(myRNG);
+
+    default:
+        psError(PS_ERR_UNEXPECTED_NULL,
+                true,
+                PS_ERRORTEXT_psRandom_UNKNOWN_RANDFOM_NUMBER_GENERATOR_TYPE);
+    }
     return(NULL);
 }
@@ -44 +66 @@
 void psRandomReset(psRandom *rand,
                    psU64 seed)
-{}
+{
+    gsl_rng_set(rand->gsl, seed);
+}
 
 psF64 psRandomUniform(const psRandom *r)
 {
-    return(0.0);
+    return(gsl_rng_uniform(r->gsl));
 }
 
 psF64 psRandomGaussian(const psRandom *r)
 {
-    return(0.0);
+    // XXX: What should sigma be?
+    return(gsl_ran_gaussian(r->gsl, 1.0));
 }
 
 psF64 psRandomPoisson(const psRandom *r, psF64 mean)
 {
-    return(0.0);
+    return((psF64) gsl_ran_poisson(r->gsl, mean));
 }