コード例 #1
0
double
rangam(double a)
{
  /**
   generate gamma deviate mean a
  */
  if (a < 1.0) {
    return( randev0(a));
  }
  if (a == 1.0) {
    return( ranexp());
  }
  return( randev1(a));
}
コード例 #2
0
ファイル: gds.c プロジェクト: yiqinyang2/SOAPpop
double
rangam(double a)
{
  /**
   generate gamma deviate mean a
  */

  if (a<=0.0) {  
   fatalx("rangam called with bad param.  a:  %9.3f\n", a) ;
  }

  if (a < 1.0) {
    return( randev0(a));
  }
  if (a == 1.0) {
    return( ranexp());
  }
  return( randev1(a));
}