double iglsadfft(double x, double *c, const int m, double *d) {
  int i;
  double y;
  y = x + d[0];
  d[m] = c[m] * x;
  for (i = m - 1; i >= 1; i--)
    d[i] += c[i] * x;
  for (i = 0; i < m; i++)
    d[i] = d[i + 1];
  return (y);
}


x86 assembly (clang, -O0):
<iglsadfft>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x8(%rbp)
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    %rdx,-0x20(%rbp)
movsd  -0x8(%rbp),%xmm0
mov    -0x20(%rbp),%rax
addsd  (%rax),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x8(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x14(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x1,-0x24(%rbp)
jl     1192 <iglsadfft+0x92>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x8(%rbp),%xmm2
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
jmp    114e <iglsadfft+0x4e>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11d2 <iglsadfft+0xd2>
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1199 <iglsadfft+0x99>
movsd  -0x30(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<iglsadfft>:
push   %rbx
sub    $0x20,%rsp
movsd  (%rdx),%xmm6
movslq %esi,%rax
movsd  (%rdi,%rax,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%rax,8)
mov    %esi,%r10d
cmp    $0x2,%eax
jl     11e0 <iglsadfft+0xd0>
xor    %r9d,%r9d
mov    %r10,%rcx
sub    $0x2,%rcx
mov    $0x0,%r8d
cmovae %rcx,%r8
add    $0x1,%r8
mov    %r10,%rax
cmp    $0x4,%r8
jb     11b9 <iglsadfft+0xa9>
sub    $0x2,%rax
mov    $0x0,%ecx
cmovae %rax,%rcx
mov    %r10d,%eax
sub    $0x2,%eax
cmovb  %r9d,%eax
lea    -0x1(%rsi),%r9d
cmp    %eax,%r9d
jb     11b6 <iglsadfft+0xa6>
shr    $0x20,%rcx
jne    11b6 <iglsadfft+0xa6>
cmp    $0x2,%r10
mov    $0x2,%ecx
cmovb  %r10,%rcx
add    %r9,%rcx
sub    %r10,%rcx
lea    (%rdx,%rcx,8),%r11
lea    (%rdi,%r9,8),%rax
add    $0x8,%rax
cmp    %rax,%r11
jae    121d <iglsadfft+0x10d>
lea    (%rdx,%r9,8),%rax
add    $0x8,%rax
lea    (%rdi,%rcx,8),%rcx
cmp    %rax,%rcx
jae    121d <iglsadfft+0x10d>
mov    %r10,%rax
add    $0x1,%rax
nopl   (%rax)
lea    -0x2(%rax),%ecx
movsd  (%rdi,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rdx,%rcx,8),%xmm1
movsd  %xmm1,(%rdx,%rcx,8)
add    $0xffffffffffffffff,%rax
cmp    $0x2,%rax
ja     11c0 <iglsadfft+0xb0>
test   %esi,%esi
jle    120f <iglsadfft+0xff>
lea    0x8(%rdx),%rsi
shl    $0x3,%r10
mov    %rdx,%rdi
mov    %r10,%rdx
movapd %xmm0,0x10(%rsp)
movsd  %xmm6,0x8(%rsp)
call   1030 <memmove@plt>
movsd  0x8(%rsp),%xmm6
movapd 0x10(%rsp),%xmm0
addsd  %xmm0,%xmm6
movapd %xmm6,%xmm0
add    $0x20,%rsp
pop    %rbx
ret
mov    %r8,%r9
and    $0xfffffffffffffffc,%r9
mov    %r10,%rax
sub    %r9,%rax
movapd %xmm0,%xmm5
unpcklpd %xmm0,%xmm5
lea    -0x1(%r10),%rcx
shufpd $0x1,%xmm5,%xmm5
mov    %r9,%r11
xchg   %ax,%ax
mov    %ecx,%ebx
shl    $0x3,%rbx
movupd -0x18(%rdi,%rbx,1),%xmm1
movupd -0x8(%rdi,%rbx,1),%xmm2
movupd -0x18(%rdx,%rbx,1),%xmm3
movupd -0x8(%rdx,%rbx,1),%xmm4
mulpd  %xmm5,%xmm2
addpd  %xmm4,%xmm2
movupd %xmm2,-0x8(%rdx,%rbx,1)
mulpd  %xmm5,%xmm1
addpd  %xmm3,%xmm1
movupd %xmm1,-0x18(%rdx,%rbx,1)
add    $0xfffffffffffffffc,%rcx
add    $0xfffffffffffffffc,%r11
jne    1240 <iglsadfft+0x130>
cmp    %r9,%r8
je     11e0 <iglsadfft+0xd0>
jmp    11b9 <iglsadfft+0xa9>

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

LLVM IR (-O2):
define double @iglsadfft(double %0, double* nocapture %1, i32 %2, double* nocapture %3) {
%5 = bitcast double* %3 to i8*
%6 = load double, double* %3
%7 = sext i32 %2 to i64
%8 = getelementptr inbounds double, double* %1, i64 %7
%9 = load double, double* %8
%10 = fmul double %9, %0
%11 = getelementptr inbounds double, double* %3, i64 %7
store double %10, double* %11
%12 = icmp sgt i32 %2, 1
br i1 %12, label %13, label %85
13: ; preds = %4
%14 = zext i32 %2 to i64
%15 = add nuw nsw i64 %14, 1
%16 = call i64 @llvm.umin.i64(i64 %14, i64 2)
%17 = sub nsw i64 %15, %16
%18 = icmp ult i64 %17, 4
br i1 %18, label %83, label %19
19: ; preds = %13
%20 = call i64 @llvm.umin.i64(i64 %14, i64 2)
%21 = sub nsw i64 %14, %20
%22 = add i32 %2, -1
%23 = trunc i64 %21 to i32
%24 = icmp ult i32 %22, %23
%25 = icmp ugt i64 %21, 4294967295
%26 = or i1 %24, %25
br i1 %26, label %83, label %27
27: ; preds = %19
%28 = call i64 @llvm.umin.i64(i64 %14, i64 2)
%29 = add i32 %2, -1
%30 = zext i32 %29 to i64
%31 = add nuw nsw i64 %28, %30
%32 = sub nsw i64 %31, %14
%33 = getelementptr double, double* %3, i64 %32
%34 = add nuw nsw i64 %30, 1
%35 = getelementptr double, double* %3, i64 %34
%36 = getelementptr double, double* %1, i64 %32
%37 = getelementptr double, double* %1, i64 %34
%38 = icmp ult double* %33, %37
%39 = icmp ult double* %36, %35
%40 = and i1 %38, %39
br i1 %40, label %83, label %41
41: ; preds = %27
%42 = and i64 %17, -4
%43 = sub nsw i64 %14, %42
%44 = insertelement <2 x double> poison, double %0, i64 0
%45 = shufflevector <2 x double> %44, <2 x double> poison, <2 x i32> zeroinitializer
%46 = insertelement <2 x double> poison, double %0, i64 0
%47 = shufflevector <2 x double> %46, <2 x double> poison, <2 x i32> zeroinitializer
br label %48
48: ; preds = %48, %41
%49 = phi i64 [ 0, %41 ], [ %79, %48 ]
%50 = xor i64 %49, -1
%51 = add i64 %50, %14
%52 = and i64 %51, 4294967295
%53 = getelementptr inbounds double, double* %1, i64 %52
%54 = getelementptr inbounds double, double* %53, i64 -1
%55 = bitcast double* %54 to <2 x double>*
%56 = load <2 x double>, <2 x double>* %55
%57 = shufflevector <2 x double> %56, <2 x double> poison, <2 x i32> <i32 1, i32 0>
%58 = getelementptr inbounds double, double* %53, i64 -2
%59 = getelementptr inbounds double, double* %58, i64 -1
%60 = bitcast double* %59 to <2 x double>*
%61 = load <2 x double>, <2 x double>* %60
%62 = shufflevector <2 x double> %61, <2 x double> poison, <2 x i32> <i32 1, i32 0>
%63 = getelementptr inbounds double, double* %3, i64 %52
%64 = getelementptr inbounds double, double* %63, i64 -1
%65 = bitcast double* %64 to <2 x double>*
%66 = load <2 x double>, <2 x double>* %65
%67 = shufflevector <2 x double> %66, <2 x double> poison, <2 x i32> <i32 1, i32 0>
%68 = getelementptr inbounds double, double* %63, i64 -2
%69 = getelementptr inbounds double, double* %68, i64 -1
%70 = bitcast double* %69 to <2 x double>*
%71 = load <2 x double>, <2 x double>* %70
%72 = shufflevector <2 x double> %71, <2 x double> poison, <2 x i32> <i32 1, i32 0>
%73 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %57, <2 x double> %45, <2 x double> %67)
%74 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %62, <2 x double> %47, <2 x double> %72)
%75 = shufflevector <2 x double> %73, <2 x double> poison, <2 x i32> <i32 1, i32 0>
%76 = bitcast double* %64 to <2 x double>*
store <2 x double> %75, <2 x double>* %76
%77 = shufflevector <2 x double> %74, <2 x double> poison, <2 x i32> <i32 1, i32 0>
%78 = bitcast double* %69 to <2 x double>*
store <2 x double> %77, <2 x double>* %78
%79 = add nuw i64 %49, 4
%80 = icmp eq i64 %79, %42
br i1 %80, label %81, label %48
81: ; preds = %48
%82 = icmp eq i64 %17, %42
br i1 %82, label %85, label %83
83: ; preds = %27, %19, %13, %81
%84 = phi i64 [ %14, %27 ], [ %14, %19 ], [ %14, %13 ], [ %43, %81 ]
br label %92
85: ; preds = %92, %81, %4
%86 = icmp sgt i32 %2, 0
br i1 %86, label %87, label %102
87: ; preds = %85
%88 = getelementptr double, double* %3, i64 1
%89 = bitcast double* %88 to i8*
%90 = zext i32 %2 to i64
%91 = shl nuw nsw i64 %90, 3
br label %102
92: ; preds = %83, %92
%93 = phi i64 [ %94, %92 ], [ %84, %83 ]
%94 = add nsw i64 %93, -1
%95 = and i64 %94, 4294967295
%96 = getelementptr inbounds double, double* %1, i64 %95
%97 = load double, double* %96
%98 = getelementptr inbounds double, double* %3, i64 %95
%99 = load double, double* %98
%100 = tail call double @llvm.fmuladd.f64(double %97, double %0, double %99)
store double %100, double* %98
%101 = icmp ugt i64 %93, 2
br i1 %101, label %92, label %85
102: ; preds = %87, %85
%103 = fadd double %6, %0
ret double %103
}