aboutsummaryrefslogtreecommitdiff
path: root/include/extent-sampling.h
diff options
context:
space:
mode:
Diffstat (limited to 'include/extent-sampling.h')
-rw-r--r--include/extent-sampling.h39
1 files changed, 39 insertions, 0 deletions
diff --git a/include/extent-sampling.h b/include/extent-sampling.h
new file mode 100644
index 0000000..8ed3fb2
--- /dev/null
+++ b/include/extent-sampling.h
@@ -0,0 +1,39 @@
+#ifndef EXTENT_SAMPLING_H
+#define EXTENT_SAMPLING_H
+
+#include <gsl/gsl_rng.h>
+#include <gsl/gsl_rstat.h>
+#include <gsl/gsl_vector.h>
+
+typedef double (*extent_oracle_t) (const gsl_vector*);
+typedef double (*integrand_t) (double, const gsl_vector*);
+
+void init_random (void);
+
+double volume
+(extent_oracle_t extent_oracle, double true_volume,
+ const gsl_rng* r, unsigned int dimension, double rtol,
+ gsl_rstat_workspace* stats);
+
+double volume_window
+(extent_oracle_t extent_oracle, double true_volume,
+ const gsl_rng* r, unsigned int dimension, double rtol,
+ unsigned int* number_of_samples);
+
+double integral
+(integrand_t integrand, extent_oracle_t extent_oracle, double true_integral,
+ const gsl_rng* r, unsigned int dimension, double rtol,
+ gsl_rstat_workspace* stats);
+
+double volume_cone
+(extent_oracle_t extent_oracle, const gsl_rng* r,
+ const gsl_vector* mean, double omega_min, double omega_max,
+ unsigned int number_of_samples, double* variance);
+
+double volume_experiment
+(extent_oracle_t extent_oracle, const gsl_rng* r,
+ const gsl_vector* mean, unsigned int samples_per_cone,
+ double solid_angle_factor, double solid_angle_threshold_exponent_factor,
+ unsigned int* number_of_samples);
+
+#endif