aboutsummaryrefslogtreecommitdiff
path: root/src
diff options
context:
space:
mode:
authorArun Isaac2021-05-07 17:18:25 +0530
committerArun Isaac2021-05-07 17:18:25 +0530
commit7ec6e79a83868fb618b844f7eb46a693c3f021e1 (patch)
treee73b3cf4d087d85693dd01481bdff5d16969c2dc /src
parentd90f2982c2de41dee7dd91c20698b4c30085eeb3 (diff)
downloadnsmc-7ec6e79a83868fb618b844f7eb46a693c3f021e1.tar.gz
nsmc-7ec6e79a83868fb618b844f7eb46a693c3f021e1.tar.lz
nsmc-7ec6e79a83868fb618b844f7eb46a693c3f021e1.zip
Implement spheroid importance sampling experiments.
* src/spheroid.sc: New file. * CMakeLists.txt: Add spheroid executable.
Diffstat (limited to 'src')
-rw-r--r--src/spheroid.sc59
1 files changed, 59 insertions, 0 deletions
diff --git a/src/spheroid.sc b/src/spheroid.sc
new file mode 100644
index 0000000..830fa7c
--- /dev/null
+++ b/src/spheroid.sc
@@ -0,0 +1,59 @@
+(sc-include "macros/macros")
+
+(pre-include "stdio.h")
+(pre-include "stdlib.h")
+(pre-include "gsl/gsl_blas.h")
+(pre-include "extent-sampling.h")
+(pre-include "oracles.h")
+(pre-include "utils.h")
+
+(pre-define DIMENSION-START 10)
+(pre-define DIMENSION-STEP 10)
+(pre-define DIMENSION-STOP 100)
+
+(pre-define RTOL 0.1)
+(pre-define SAMPLES-PER-CONE 1000)
+(pre-define SOLID-ANGLE-FACTOR 1.1)
+
+(define (experiment eccentricity solid-angle-threshold-exponent-factor fp)
+ (void double double FILE*)
+ (declare params (struct-variable spheroid-params
+ (eccentricity eccentricity)))
+ (declare oracle (struct-variable extent-oracle-t
+ (oracle spheroid-extent-oracle)
+ (params (address-of params))))
+ (fprintf fp "dimension\tsamples\n")
+ (with-rng r
+ (for-i-step dim DIMENSION-START DIMENSION-STOP DIMENSION-STEP
+ (with-vector mean dim
+ (gsl-vector-set-basis mean 0)
+ (gsl-vector-scale mean (/ (gsl-blas-dnrm2 mean)))
+ (declare samples (unsigned int))
+ (define volume double (* 2 (volume-importance (address-of oracle)
+ r
+ mean
+ SAMPLES-PER-CONE
+ SOLID-ANGLE-FACTOR
+ solid-angle-threshold-exponent-factor
+ (address-of samples))))
+ (let* ((true-volume double (spheroid-true-volume dim (address-of params))))
+ (unless (rtol? volume true-volume RTOL)
+ (fprintf stderr
+ "Volume estimation failed to converge: dimension=%d, true volume=%g, estimated volume=%g, rtol=%g\n"
+ dim
+ true-volume
+ volume
+ RTOL)
+ (exit EXIT_FAILURE)))
+ (fprintf fp "%d\t%d\n" dim samples)
+ (fflush NULL)))))
+
+(sc-define-syntax* (run-experiment eccentricity solid-angle-threshold-exponent-factor)
+ `(begin
+ (with-data-file fp ,(format #f "spheroid-eccentricity-~a.dat" eccentricity) "w"
+ (experiment ,eccentricity ,solid-angle-threshold-exponent-factor fp))))
+
+(define (main) (int)
+ (run-experiment 10 2)
+ (run-experiment 100 6)
+ (return 0))