double glaguerre(long m, double b, double x) {
  long i;
  double gln, glnm1, glnp1;
  if (m == 0)
    return 1.0;
  else {
    if (m == 1)
      return 1.0 + b - x;
    else {
      gln = 1.0 + b - x;
      glnm1 = 1.0;
      for (i = 2; i <= m; i++) {
        glnp1 = ((2 * (i - 1) + b + 1.0 - x) * gln - (i - 1 + b) * glnm1) / i;
        glnm1 = gln;
        gln = glnp1;
      }
      return gln;
    }
  }
}


x86 assembly (clang, -O0):
<glaguerre>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
movsd  %xmm1,-0x20(%rbp)
cmpq   $0x0,-0x10(%rbp)
jne    112f <glaguerre+0x2f>
movsd  0xedb(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1230 <glaguerre+0x130>
cmpq   $0x1,-0x10(%rbp)
jne    1156 <glaguerre+0x56>
movsd  0xebe(%rip),%xmm0
addsd  -0x18(%rbp),%xmm0
subsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1230 <glaguerre+0x130>
movsd  0xea2(%rip),%xmm0
addsd  -0x18(%rbp),%xmm0
subsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  0xe8b(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
movq   $0x2,-0x28(%rbp)
mov    -0x28(%rbp),%rax
cmp    -0x10(%rbp),%rax
jg     1226 <glaguerre+0x126>
mov    -0x28(%rbp),%rax
sub    $0x1,%rax
shl    $0x1,%rax
cvtsi2sd %rax,%xmm0
addsd  -0x18(%rbp),%xmm0
movsd  0xe52(%rip),%xmm1
addsd  %xmm1,%xmm0
subsd  -0x20(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
mov    -0x28(%rbp),%rax
sub    $0x1,%rax
cvtsi2sd %rax,%xmm1
addsd  -0x18(%rbp),%xmm1
mulsd  -0x38(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsi2sdq -0x28(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    1182 <glaguerre+0x82>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<glaguerre>:
test   %rdi,%rdi
je     1118 <glaguerre+0x18>
cmp    $0x1,%rdi
jne    1121 <glaguerre+0x21>
addsd  0xeed(%rip),%xmm0
subsd  %xmm1,%xmm0
ret
movsd  0xee0(%rip),%xmm0
ret
movsd  0xed7(%rip),%xmm3
addsd  %xmm0,%xmm3
subsd  %xmm1,%xmm3
jg     1138 <glaguerre+0x38>
movapd %xmm3,%xmm0
ret
movsd  0xec0(%rip),%xmm4
mov    $0x1,%eax
mov    $0x2,%ecx
movapd %xmm4,%xmm5
xchg   %ax,%ax
xorps  %xmm2,%xmm2
cvtsi2sd %rcx,%xmm2
addsd  %xmm0,%xmm2
addsd  %xmm4,%xmm2
subsd  %xmm1,%xmm2
xorps  %xmm6,%xmm6
cvtsi2sd %rax,%xmm6
add    $0x1,%rax
addsd  %xmm0,%xmm6
mulsd  %xmm5,%xmm6
mulsd  %xmm3,%xmm2
xorps  %xmm5,%xmm5
cvtsi2sd %rax,%xmm5
subsd  %xmm6,%xmm2
divsd  %xmm5,%xmm2
add    $0x2,%rcx
movapd %xmm3,%xmm5
movapd %xmm2,%xmm3
cmp    %rax,%rdi
jne    1150 <glaguerre+0x50>
movapd %xmm2,%xmm0
ret

LLVM IR (-O0):
define double @glaguerre(i64 %0, double %1, double %2) {
%4 = alloca double
%5 = alloca i64
%6 = alloca double
%7 = alloca double
%8 = alloca i64
%9 = alloca double
%10 = alloca double
%11 = alloca double
store i64 %0, i64* %5
store double %1, double* %6
store double %2, double* %7
%12 = load i64, i64* %5
%13 = icmp eq i64 %12, 0
br i1 %13, label %14, label %15
14: ; preds = %3
store double 1.000000e+00, double* %4
br label %62
15: ; preds = %3
%16 = load i64, i64* %5
%17 = icmp eq i64 %16, 1
br i1 %17, label %18, label %23
18: ; preds = %15
%19 = load double, double* %6
%20 = fadd double 1.000000e+00, %19
%21 = load double, double* %7
%22 = fsub double %20, %21
store double %22, double* %4
br label %62
23: ; preds = %15
%24 = load double, double* %6
%25 = fadd double 1.000000e+00, %24
%26 = load double, double* %7
%27 = fsub double %25, %26
store double %27, double* %9
store double 1.000000e+00, double* %10
store i64 2, i64* %8
br label %28
28: ; preds = %57, %23
%29 = load i64, i64* %8
%30 = load i64, i64* %5
%31 = icmp sle i64 %29, %30
br i1 %31, label %32, label %60
32: ; preds = %28
%33 = load i64, i64* %8
%34 = sub nsw i64 %33, 1
%35 = mul nsw i64 2, %34
%36 = sitofp i64 %35 to double
%37 = load double, double* %6
%38 = fadd double %36, %37
%39 = fadd double %38, 1.000000e+00
%40 = load double, double* %7
%41 = fsub double %39, %40
%42 = load double, double* %9
%43 = load i64, i64* %8
%44 = sub nsw i64 %43, 1
%45 = sitofp i64 %44 to double
%46 = load double, double* %6
%47 = fadd double %45, %46
%48 = load double, double* %10
%49 = fmul double %47, %48
%50 = fneg double %49
%51 = call double @llvm.fmuladd.f64(double %41, double %42, double %50)
%52 = load i64, i64* %8
%53 = sitofp i64 %52 to double
%54 = fdiv double %51, %53
store double %54, double* %11
%55 = load double, double* %9
store double %55, double* %10
%56 = load double, double* %11
store double %56, double* %9
br label %57
57: ; preds = %32
%58 = load i64, i64* %8
%59 = add nsw i64 %58, 1
store i64 %59, i64* %8
br label %28
60: ; preds = %28
%61 = load double, double* %9
store double %61, double* %4
br label %62
62: ; preds = %60, %18, %14
%63 = load double, double* %4
ret double %63
}

LLVM IR (-O2):
define double @glaguerre(i64 %0, double %1, double %2) {
switch i64 %0, label %7 [
i64 0, label %30
i64 1, label %4
]
4: ; preds = %3
%5 = fadd double %1, 1.000000e+00
%6 = fsub double %5, %2
br label %30
7: ; preds = %3
%8 = fadd double %1, 1.000000e+00
%9 = fsub double %8, %2
%10 = icmp slt i64 %0, 2
br i1 %10, label %30, label %11
11: ; preds = %7, %11
%12 = phi double [ %13, %11 ], [ 1.000000e+00, %7 ]
%13 = phi double [ %27, %11 ], [ %9, %7 ]
%14 = phi i64 [ %28, %11 ], [ 2, %7 ]
%15 = add nsw i64 %14, -1
%16 = shl nsw i64 %15, 1
%17 = sitofp i64 %16 to double
%18 = fadd double %17, %1
%19 = fadd double %18, 1.000000e+00
%20 = fsub double %19, %2
%21 = sitofp i64 %15 to double
%22 = fadd double %21, %1
%23 = fneg double %22
%24 = fmul double %12, %23
%25 = tail call double @llvm.fmuladd.f64(double %20, double %13, double %24)
%26 = sitofp i64 %14 to double
%27 = fdiv double %25, %26
%28 = add nuw i64 %14, 1
%29 = icmp eq i64 %14, %0
br i1 %29, label %30, label %11
30: ; preds = %11, %7, %3, %4
%31 = phi double [ %6, %4 ], [ 1.000000e+00, %3 ], [ %9, %7 ], [ %27, %11 ]
ret double %31
}