diff options
Diffstat (limited to 'include/extent-sampling.h')
-rw-r--r-- | include/extent-sampling.h | 39 |
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 |