Generando un doble gaussiano aleatorio en Objective-C/C

Question

Estoy intentando generar un doble gaussiano aleatorio en Objective-C (lo mismo que random.nextGaussian en Java). Sin embargo, rand_gauss() parece no funcionar. Alguien sabe una manera de lograr esto?

Answer 1

This link muestra cómo calcularlo utilizando la función estándar random().

deben tener en cuenta que es probable que tenga que hacer la rutina ranf() que convierte la salida de random() de [0,MAX_INT] ser de [0,1], pero eso no debería ser demasiado difícil.

Desde el artículo enlazado:

La forma polar de la transformación de Box-Muller es a la vez más rápido y numéricamente más robusto. La descripción algorítmica de este es: float x1, x2, w, y1, y2;

 do { 
      x1 = 2.0 * ranf() - 1.0; 
      x2 = 2.0 * ranf() - 1.0; 
      w = x1 * x1 + x2 * x2; 
    } while (w >= 1.0); 

    w = sqrt((-2.0 * ln(w))/w); 
    y1 = x1 * w; 
    y2 = x2 * w;