#ifndef GAUSSIAN_ND_RANDOM_H
#define GAUSSIAN_ND_RANDOM_H

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

double planar_angle_to_standard_deviation
(double mean, double theta_max, double truncation, unsigned int dimension);

unsigned int shifted_gaussian_random_vector
(const gsl_rng* r, const gsl_vector* mean,
 double theta_max, double truncation, gsl_vector* x);

double shifted_gaussian_pdf
(double theta, double mean, double theta_max,
 double truncation, unsigned int dimension, gsl_integration_workspace* w);

#endif