expdev.c

来自「IBM(原)数据生成器和源代码」· C语言 代码 · 共 15 行

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

float ExpDist::expdev(void)
// Returns an exponentially distributed, positive, random deviate of
// unit mean, using ran1(idum) as the source of uniform deviates.
{
  float dum;
  
  do
    dum = (*ran1)();
  while (dum == 0.0);
  return -log(dum);
}