int poisso(u)
double u;
{
  double cump, ru, ran1(), p, gasdev();
  int i = 1;
  if (u > 30.)
    return ((int)(0.5 + gasdev(u, u)));
  ru = ran1();
  p = exp(-u);
  if (ru < p)
    return (0);
  cump = p;
  while (ru > (cump += (p *= u / i)))
    i++;
  return (i);
}


x86 assembly (clang, -O0):
<poisso>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
movsd  %xmm0,-0x10(%rbp)
movl   $0x1,-0x2c(%rbp)
movsd  -0x10(%rbp),%xmm0
movsd  0xeaf(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1187 <poisso+0x57>
movsd  -0x10(%rbp),%xmm0
movsd  -0x10(%rbp),%xmm1
mov    $0x2,%al
call   1050 <gasdev@plt>
movaps %xmm0,%xmm1
movsd  0xe91(%rip),%xmm0
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x4(%rbp)
jmp    1228 <poisso+0xf8>
mov    $0x0,%al
call   1040 <ran1@plt>
movsd  %xmm0,-0x20(%rbp)
movsd  -0x10(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
call   1030 <exp@plt>
movsd  %xmm0,-0x28(%rbp)
movsd  -0x20(%rbp),%xmm1
movsd  -0x28(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    11d9 <poisso+0xa9>
movl   $0x0,-0x4(%rbp)
jmp    1228 <poisso+0xf8>
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x20(%rbp),%xmm0
movsd  -0x10(%rbp),%xmm1
cvtsi2sdl -0x2c(%rbp),%xmm2
divsd  %xmm2,%xmm1
mulsd  -0x28(%rbp),%xmm1
movsd  %xmm1,-0x28(%rbp)
addsd  -0x18(%rbp),%xmm1
movsd  %xmm1,-0x18(%rbp)
ucomisd %xmm1,%xmm0
jbe    1222 <poisso+0xf2>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11e3 <poisso+0xb3>
mov    -0x2c(%rbp),%eax
mov    %eax,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<poisso>:
push   %rbx
sub    $0x20,%rsp
ucomisd 0xec3(%rip),%xmm0
jbe    1158 <poisso+0x28>
movapd %xmm0,%xmm1
mov    $0x2,%al
call   1050 <gasdev@plt>
addsd  0xeb6(%rip),%xmm0
cvttsd2si %xmm0,%ebx
jmp    11c0 <poisso+0x90>
xor    %ebx,%ebx
xor    %eax,%eax
movapd %xmm0,0x10(%rsp)
call   1040 <ran1@plt>
movsd  %xmm0,0x8(%rsp)
movapd 0xe9b(%rip),%xmm0
xorpd  0x10(%rsp),%xmm0
call   1030 <exp@plt>
movsd  0x8(%rsp),%xmm5
movapd 0x10(%rsp),%xmm4
ucomisd %xmm5,%xmm0
ja     11c0 <poisso+0x90>
xor    %ebx,%ebx
movapd %xmm0,%xmm1
nopl   0x0(%rax,%rax,1)
add    $0x1,%ebx
xorps  %xmm2,%xmm2
cvtsi2sd %ebx,%xmm2
movapd %xmm4,%xmm3
divsd  %xmm2,%xmm3
mulsd  %xmm3,%xmm1
addsd  %xmm1,%xmm0
ucomisd %xmm0,%xmm5
ja     11a0 <poisso+0x70>
mov    %ebx,%eax
add    $0x20,%rsp
pop    %rbx
ret

LLVM IR (-O0):
define i32 @poisso(double %0) {
%2 = alloca i32
%3 = alloca double
%4 = alloca double
%5 = alloca double
%6 = alloca double
%7 = alloca i32
store double %0, double* %3
store i32 1, i32* %7
%8 = load double, double* %3
%9 = fcmp ogt double %8, 3.000000e+01
br i1 %9, label %10, label %16
10: ; preds = %1
%11 = load double, double* %3
%12 = load double, double* %3
%13 = call double (double, double, ...) bitcast (double (...)* @gasdev to double (double, double, ...)*)(double %11, double %12)
%14 = fadd double 5.000000e-01, %13
%15 = fptosi double %14 to i32
store i32 %15, i32* %2
br label %43
16: ; preds = %1
%17 = call double (...) @ran1()
store double %17, double* %5
%18 = load double, double* %3
%19 = fneg double %18
%20 = call double @exp(double %19)
store double %20, double* %6
%21 = load double, double* %5
%22 = load double, double* %6
%23 = fcmp olt double %21, %22
br i1 %23, label %24, label %25
24: ; preds = %16
store i32 0, i32* %2
br label %43
25: ; preds = %16
%26 = load double, double* %6
store double %26, double* %4
br label %27
27: ; preds = %38, %25
%28 = load double, double* %5
%29 = load double, double* %3
%30 = load i32, i32* %7
%31 = sitofp i32 %30 to double
%32 = fdiv double %29, %31
%33 = load double, double* %6
%34 = fmul double %33, %32
store double %34, double* %6
%35 = load double, double* %4
%36 = fadd double %35, %34
store double %36, double* %4
%37 = fcmp ogt double %28, %36
br i1 %37, label %38, label %41
38: ; preds = %27
%39 = load i32, i32* %7
%40 = add nsw i32 %39, 1
store i32 %40, i32* %7
br label %27
41: ; preds = %27
%42 = load i32, i32* %7
store i32 %42, i32* %2
br label %43
43: ; preds = %41, %24, %10
%44 = load i32, i32* %2
ret i32 %44
}

LLVM IR (-O2):
define i32 @poisso(double %0) {
%2 = fcmp ogt double %0, 3.000000e+01
br i1 %2, label %3, label %7
3: ; preds = %1
%4 = tail call double (double, double, ...) bitcast (double (...)* @gasdev to double (double, double, ...)*)(double %0, double %0)
%5 = fadd double %4, 5.000000e-01
%6 = fptosi double %5 to i32
br label %22
7: ; preds = %1
%8 = tail call double (...) @ran1()
%9 = fneg double %0
%10 = tail call double @exp(double %9)
%11 = fcmp olt double %8, %10
br i1 %11, label %22, label %12
12: ; preds = %7, %12
%13 = phi double [ %19, %12 ], [ %10, %7 ]
%14 = phi double [ %18, %12 ], [ %10, %7 ]
%15 = phi i32 [ %21, %12 ], [ 1, %7 ]
%16 = sitofp i32 %15 to double
%17 = fdiv double %0, %16
%18 = fmul double %14, %17
%19 = fadd double %13, %18
%20 = fcmp ogt double %8, %19
%21 = add nuw nsw i32 %15, 1
br i1 %20, label %12, label %22
22: ; preds = %12, %7, %3
%23 = phi i32 [ %6, %3 ], [ 0, %7 ], [ %15, %12 ]
ret i32 %23
}