Add a Monte Carlo simulator.

[wloh] / ziggurat.hpp

diff --git a/ziggurat.hpp b/ziggurat.hpp
new file mode 100644 (file)
index 0000000..af8272c
--- /dev/null
+++ b/ziggurat.hpp
@@ -0,0 +1,9 @@
+float r4_exp ( unsigned long int *jsr, int ke[256], float fe[256], 
+  float we[256] );
+void r4_exp_setup ( int ke[256], float fe[256], float we[256] );
+float r4_nor ( unsigned long int *jsr, int kn[128], float fn[128], 
+  float wn[128] );
+void r4_nor_setup ( int kn[128], float fn[128], float wn[128] );
+float r4_uni ( unsigned long int *jsr );
+unsigned long int shr3 ( unsigned long int *jsr );
+void timestamp ( );