#ifndef ND_RANDOM_H
#define ND_RANDOM_H

#include <gsl/gsl_integration.h>
#include <gsl/gsl_rng.h>
#include <gsl/gsl_vector.h>

void random_direction_vector (const gsl_rng* r, gsl_vector* x);
void hollow_cone_random_vector (const gsl_rng* r, const gsl_vector* mean, double theta_min, double theta_max, gsl_vector* x);
void subsampling_random_vector (const gsl_rng* r, const gsl_vector* mean, double theta_max, gsl_vector* x);

double planar_angle_to_solid_angle (double planar_angle, unsigned int dimension);
double solid_angle_to_planar_angle (double solid_angle, unsigned int dimension);

#endif