//  Cauchy random numbers a la RPP
//  G. Cowan

#include "random.h"
#include <math.h>

float CauchyRan(int* SeedPtr){
  float rsq, v1, v2;
  do {
    v1 =  2.0 * random(*SeedPtr) - 1.0;
    v2 =  2.0 * random(*SeedPtr) - 1.0;
    rsq = v1*v1 + v2*v2;
  } while (rsq > 1.0 || v2 == 0 );
  return v1/v2;
}