utilities.c File Reference

#include <complex.h>
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include "so3/so3_error.h"
#include "so3/so3_sampling.h"
#include "utilities.h"

Include dependency graph for utilities.c:

Functions
int	max (int a, int b)
void	gen_flmn_complex (complex double *flmn, const so3_parameters_t *parameters, int seed)
void	gen_flmn_real (complex double *flmn, const so3_parameters_t *parameters, int seed)
double	ran2_dp (int idum)

Function Documentation

gen_flmn_complex()

void gen_flmn_complex	(	complex double *	flmn,
		const so3_parameters_t *	parameters,
		int	seed
	)

Generate random Wigner coefficients of a complex signal.

Parameters

[out]	flmn	Random spherical harmonic coefficients generated. Provide enough memory for fully padded storage, i.e. (2*N-1)*L*L elements. Unused trailing elements will be set to zero.
[in]	parameters	A parameters object with (at least) the following fields: L0, L, N, storage, n_mode The reality flag is ignored. Use so3_test_gen_flmn_real instead for real signals.
[in]	seed	Integer seed required for random number generator.

Return values

none

Author

Martin Büttner

Jason McEwen

gen_flmn_real()

void gen_flmn_real	(	complex double *	flmn,
		const so3_parameters_t *	parameters,
		int	seed
	)

Generate random Wigner coefficients of a real signal. We only generate coefficients for n >= 0, and for n = 0, we need flm0* = (-1)^(m)*fl-m0, so that fl00 has to be real.

Parameters

[out]	flmn	Random spherical harmonic coefficients generated. Provide enough memory for fully padded, complex (!) storage, i.e. (2*N-1)*L*L elements. Unused trailing elements will be set to zero.
[in]	parameters	A parameters object with (at least) the following fields: L0, L, N, storage, n_mode The reality flag is ignored. Use so3_test_gen_flmn_complex instead for complex signals.
[in]	seed	Integer seed required for random number generator.

Return values

none

Author

Martin Büttner

Jason McEwen

max()

int max	(	int	a,
		int	b
	)

ran2_dp()

double ran2_dp

(

int

idum

)

Generate uniform deviate in range [0,1) given seed. (Using double precision.)

Note

Uniform deviate (Num rec 1992, chap 7.1), original routine said to be 'perfect'.

Parameters

[in]

idum

Seed.

Return values

ran_dp

Generated uniform deviate.

Author

Jason McEwen