#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ Reverse(long, long, double *);
void FFT1DOnce(long direction, long M, long N, double *u, double *x) {
  long j;
  long k;
  long q;
  long L;
  long r;
  long Lstar;
  double *u1;
  double *x1;
  double *x2;
  double omega_r;
  double omega_c;
  double tau_r;
  double tau_c;
  double x_r;
  double x_c;
  Reverse(N, M, x);
  for (q = 1; q <= M; q++) {
    L = 1 << q;
    r = N / L;
    Lstar = L / 2;
    u1 = &u[2 * (Lstar - 1)];
    for (k = 0; k < r; k++) {
      x1 = &x[2 * (k * L)];
      x2 = &x[2 * (k * L + Lstar)];
      for (j = 0; j < Lstar; j++) {
        omega_r = u1[2 * j];
        omega_c = direction * u1[2 * j + 1];
        x_r = x2[2 * j];
        x_c = x2[2 * j + 1];
        tau_r = omega_r * x_r - omega_c * x_c;
        tau_c = omega_r * x_c + omega_c * x_r;
        x_r = x1[2 * j];
        x_c = x1[2 * j + 1];
        x2[2 * j] = x_r - tau_r;
        x2[2 * j + 1] = x_c - tau_c;
        x1[2 * j] = x_r + tau_r;
        x1[2 * j + 1] = x_c + tau_c;
      }
    }
  }
}


x86 assembly (clang, -O0):
<FFT1DOnce>:
push   %rbp
mov    %rsp,%rbp
sub    $0xa0,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0x18(%rbp),%rdi
mov    -0x10(%rbp),%rsi
mov    -0x28(%rbp),%rdx
call   1030 <Reverse@plt>
movq   $0x1,-0x40(%rbp)
mov    -0x40(%rbp),%rax
cmp    -0x10(%rbp),%rax
jg     13d3 <FFT1DOnce+0x2c3>
mov    -0x40(%rbp),%rax
mov    %eax,%ecx
mov    $0x1,%eax
shl    %cl,%eax
cltq
mov    %rax,-0x48(%rbp)
mov    -0x18(%rbp),%rax
cqto
idivq  -0x48(%rbp)
mov    %rax,-0x50(%rbp)
mov    -0x48(%rbp),%rax
mov    $0x2,%ecx
cqto
idiv   %rcx
mov    %rax,-0x58(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x58(%rbp),%rcx
sub    $0x1,%rcx
shl    $0x1,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x60(%rbp)
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x50(%rbp),%rax
jge    13bd <FFT1DOnce+0x2ad>
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
imul   -0x48(%rbp),%rcx
shl    $0x1,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x68(%rbp)
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
imul   -0x48(%rbp),%rcx
add    -0x58(%rbp),%rcx
shl    $0x1,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x70(%rbp)
movq   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x58(%rbp),%rax
jge    13a7 <FFT1DOnce+0x297>
mov    -0x60(%rbp),%rax
mov    -0x30(%rbp),%rcx
shl    $0x1,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
cvtsi2sdq -0x8(%rbp),%xmm0
mov    -0x60(%rbp),%rax
mov    -0x30(%rbp),%rcx
shl    $0x1,%rcx
mulsd  0x8(%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x70(%rbp),%rax
mov    -0x30(%rbp),%rcx
shl    $0x1,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
mov    -0x70(%rbp),%rax
mov    -0x30(%rbp),%rcx
shl    $0x1,%rcx
movsd  0x8(%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x78(%rbp),%xmm0
movsd  -0x98(%rbp),%xmm2
ff
movsd  -0x80(%rbp),%xmm1
mulsd  -0xa0(%rbp),%xmm1
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x78(%rbp),%xmm0
movsd  -0xa0(%rbp),%xmm2
ff
movsd  -0x80(%rbp),%xmm1
mulsd  -0x98(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0x68(%rbp),%rax
mov    -0x30(%rbp),%rcx
shl    $0x1,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
mov    -0x68(%rbp),%rax
mov    -0x30(%rbp),%rcx
shl    $0x1,%rcx
movsd  0x8(%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x98(%rbp),%xmm0
ff
subsd  -0x88(%rbp),%xmm0
ff
mov    -0x70(%rbp),%rax
mov    -0x30(%rbp),%rcx
shl    $0x1,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xa0(%rbp),%xmm0
ff
subsd  -0x90(%rbp),%xmm0
ff
mov    -0x70(%rbp),%rax
mov    -0x30(%rbp),%rcx
shl    $0x1,%rcx
movsd  %xmm0,0x8(%rax,%rcx,8)
movsd  -0x98(%rbp),%xmm0
ff
addsd  -0x88(%rbp),%xmm0
ff
mov    -0x68(%rbp),%rax
mov    -0x30(%rbp),%rcx
shl    $0x1,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xa0(%rbp),%xmm0
ff
addsd  -0x90(%rbp),%xmm0
ff
mov    -0x68(%rbp),%rax
mov    -0x30(%rbp),%rcx
shl    $0x1,%rcx
movsd  %xmm0,0x8(%rax,%rcx,8)
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    11fe <FFT1DOnce+0xee>
jmp    13ac <FFT1DOnce+0x29c>
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    11ac <FFT1DOnce+0x9c>
jmp    13c2 <FFT1DOnce+0x2b2>
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    1148 <FFT1DOnce+0x38>
add    $0xa0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<FFT1DOnce>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
mov    %rcx,(%rsp)
mov    %rdx,%rbp
mov    %rsi,%rbx
mov    %rdi,%r14
mov    %rdx,%rdi
mov    %r8,0x18(%rsp)
mov    %r8,%rdx
call   1030 <Reverse@plt>
mov    %rbx,0x30(%rsp)
test   %rbx,%rbx
jle    1468 <FFT1DOnce+0x358>
cvtsi2sd %r14,%xmm0
mov    (%rsp),%rax
lea    -0x10(%rax),%rcx
mov    %rcx,0x8(%rsp)
movapd %xmm0,%xmm8
unpcklpd %xmm0,%xmm8
add    $0xfffffffffffffff8,%rax
mov    %rax,0x20(%rsp)
mov    $0x1,%ecx
movapd 0xe84(%rip),%xmm9
mov    %rbp,0x28(%rsp)
jmp    11ac <FFT1DOnce+0x9c>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    0x38(%rsp),%rcx
lea    0x1(%rcx),%rax
cmp    0x30(%rsp),%rcx
mov    %rax,%rcx
mov    0x28(%rsp),%rbp
je     1468 <FFT1DOnce+0x358>
mov    $0x1,%r11d
shl    %cl,%r11d
movslq %r11d,%r8
mov    %rbp,%rax
or     %r8,%rax
shr    $0x20,%rax
mov    %rcx,0x38(%rsp)
je     11e0 <FFT1DOnce+0xd0>
mov    %rbp,%rax
cqto
idiv   %r8
test   %rax,%rax
jg     11ec <FFT1DOnce+0xdc>
jmp    1190 <FFT1DOnce+0x80>
nopl   0x0(%rax,%rax,1)
mov    %ebp,%eax
xor    %edx,%edx
div    %r11d
test   %rax,%rax
jle    1190 <FFT1DOnce+0x80>
mov    %r11d,%ecx
shr    $0x1f,%ecx
add    %r11d,%ecx
sar    %ecx
movslq %ecx,%r14
lea    (%r14,%r14,1),%rsi
mov    %r14,%rdx
shl    $0x4,%rdx
mov    (%rsp),%rcx
add    %rdx,%rcx
add    $0xfffffffffffffff0,%rcx
mov    %rcx,0x58(%rsp)
lea    (%r8,%r8,1),%rcx
mov    %rcx,0x50(%rsp)
lea    0x0(,%r14,4),%rcx
mov    %rcx,0x48(%rsp)
mov    %r14,%rcx
shl    $0x5,%rcx
add    0x8(%rsp),%rcx
mov    %rcx,0x40(%rsp)
mov    %r14,%rcx
and    $0xfffffffffffffffe,%rcx
mov    %rcx,0x10(%rsp)
shl    $0x4,%r8
mov    0x20(%rsp),%rcx
lea    (%rcx,%rsi,8),%r15
mov    0x18(%rsp),%r13
mov    %rsi,0x60(%rsp)
lea    0x0(,%rsi,8),%r9
add    %r13,%r9
xor    %ecx,%ecx
jmp    1283 <FFT1DOnce+0x173>
xchg   %ax,%ax
add    $0x1,%rcx
add    %r8,%r13
add    %r8,%r9
cmp    %rax,%rcx
je     1190 <FFT1DOnce+0x80>
cmp    $0x2,%r11d
jl     1270 <FFT1DOnce+0x160>
cmp    $0x2,%r14d
jb     13f0 <FFT1DOnce+0x2e0>
mov    0x50(%rsp),%rbp
imul   %rcx,%rbp
mov    0x60(%rsp),%rsi
lea    (%rsi,%rbp,1),%rbx
mov    0x48(%rsp),%rsi
lea    (%rsi,%rbp,1),%rdi
mov    0x18(%rsp),%rsi
lea    (%rsi,%rbp,8),%rbp
lea    (%rsi,%rbx,8),%rbx
lea    (%rsi,%rdi,8),%r10
mov    0x40(%rsp),%rdi
cmp    %rdi,%rbp
setb   %bpl
mov    0x58(%rsp),%rsi
cmp    %rbx,%rsi
setb   %r12b
cmp    %rdi,%rbx
setb   %bl
cmp    %r10,%rsi
setb   %r10b
test   %r12b,%bpl
jne    13f0 <FFT1DOnce+0x2e0>
mov    $0x0,%ebp
and    %r10b,%bl
jne    13f2 <FFT1DOnce+0x2e2>
mov    0x10(%rsp),%r12
mov    0x8(%rsp),%r10
mov    %r13,%rbp
nopl   0x0(%rax,%rax,1)
movupd (%r10,%rdx,1),%xmm4
movupd 0x10(%r10,%rdx,1),%xmm3
movapd %xmm4,%xmm5
unpcklpd %xmm3,%xmm5
unpckhpd %xmm3,%xmm4
mulpd  %xmm8,%xmm4
movupd 0x0(%rbp,%rdx,1),%xmm3
movupd 0x10(%rbp,%rdx,1),%xmm6
movapd %xmm3,%xmm7
unpcklpd %xmm6,%xmm7
unpckhpd %xmm6,%xmm3
movapd %xmm3,%xmm6
mulpd  %xmm4,%xmm6
mulpd  %xmm5,%xmm3
mulpd  %xmm7,%xmm5
subpd  %xmm6,%xmm5
mulpd  %xmm4,%xmm7
addpd  %xmm7,%xmm3
movupd 0x0(%rbp),%xmm4
movupd 0x10(%rbp),%xmm6
movapd %xmm4,%xmm7
unpcklpd %xmm6,%xmm7
unpckhpd %xmm6,%xmm4
movapd %xmm7,%xmm6
subpd  %xmm5,%xmm6
movapd %xmm4,%xmm1
subpd  %xmm3,%xmm1
movapd %xmm6,%xmm2
unpckhpd %xmm1,%xmm2
unpcklpd %xmm1,%xmm6
movupd %xmm6,0x0(%rbp,%rdx,1)
movupd %xmm2,0x10(%rbp,%rdx,1)
addpd  %xmm5,%xmm7
addpd  %xmm3,%xmm4
movapd %xmm7,%xmm1
unpcklpd %xmm4,%xmm1
movhlps %xmm7,%xmm4
movups %xmm4,0x10(%rbp)
movupd %xmm1,0x0(%rbp)
add    $0x20,%rbp
add    $0x20,%r10
add    $0xfffffffffffffffe,%r12
jne    1310 <FFT1DOnce+0x200>
mov    0x10(%rsp),%rsi
mov    %rsi,%rbp
cmp    %r14,%rsi
je     1270 <FFT1DOnce+0x160>
jmp    13f2 <FFT1DOnce+0x2e2>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
xor    %ebp,%ebp
mov    %r14,%rbx
sub    %rbp,%rbx
shl    $0x4,%rbp
nopl   0x0(%rax)
movsd  (%r15,%rbp,1),%xmm1
mulsd  %xmm0,%xmm1
movupd (%r9,%rbp,1),%xmm2
movapd %xmm1,%xmm3
xorpd  %xmm9,%xmm3
unpcklpd %xmm1,%xmm3
movsd  -0x8(%r15,%rbp,1),%xmm1
unpcklpd %xmm1,%xmm1
mulpd  %xmm2,%xmm1
shufpd $0x1,%xmm2,%xmm2
mulpd  %xmm3,%xmm2
addpd  %xmm2,%xmm1
movupd 0x0(%r13,%rbp,1),%xmm2
movapd %xmm2,%xmm3
subpd  %xmm1,%xmm3
movupd %xmm3,(%r9,%rbp,1)
addpd  %xmm2,%xmm1
movupd %xmm1,0x0(%r13,%rbp,1)
add    $0x10,%rbp
add    $0xffffffffffffffff,%rbx
jne    1400 <FFT1DOnce+0x2f0>
jmp    1270 <FFT1DOnce+0x160>
add    $0x68,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @FFT1DOnce(i64 %0, i64 %1, i64 %2, double* %3, double* %4) {
%6 = alloca i64
%7 = alloca i64
%8 = alloca i64
%9 = alloca double*
%10 = alloca double*
%11 = alloca i64
%12 = alloca i64
%13 = alloca i64
%14 = alloca i64
%15 = alloca i64
%16 = alloca i64
%17 = alloca double*
%18 = alloca double*
%19 = alloca double*
%20 = alloca double
%21 = alloca double
%22 = alloca double
%23 = alloca double
%24 = alloca double
%25 = alloca double
store i64 %0, i64* %6
store i64 %1, i64* %7
store i64 %2, i64* %8
store double* %3, double** %9
store double* %4, double** %10
%26 = load i64, i64* %8
%27 = load i64, i64* %7
%28 = load double*, double** %10
%29 = call i32 @Reverse(i64 %26, i64 %27, double* %28)
store i64 1, i64* %13
br label %30
30: ; preds = %160, %5
%31 = load i64, i64* %13
%32 = load i64, i64* %7
%33 = icmp sle i64 %31, %32
br i1 %33, label %34, label %163
34: ; preds = %30
%35 = load i64, i64* %13
%36 = trunc i64 %35 to i32
%37 = shl i32 1, %36
%38 = sext i32 %37 to i64
store i64 %38, i64* %14
%39 = load i64, i64* %8
%40 = load i64, i64* %14
%41 = sdiv i64 %39, %40
store i64 %41, i64* %15
%42 = load i64, i64* %14
%43 = sdiv i64 %42, 2
store i64 %43, i64* %16
%44 = load double*, double** %9
%45 = load i64, i64* %16
%46 = sub nsw i64 %45, 1
%47 = mul nsw i64 2, %46
%48 = getelementptr inbounds double, double* %44, i64 %47
store double* %48, double** %17
store i64 0, i64* %12
br label %49
49: ; preds = %156, %34
%50 = load i64, i64* %12
%51 = load i64, i64* %15
%52 = icmp slt i64 %50, %51
br i1 %52, label %53, label %159
53: ; preds = %49
%54 = load double*, double** %10
%55 = load i64, i64* %12
%56 = load i64, i64* %14
%57 = mul nsw i64 %55, %56
%58 = mul nsw i64 2, %57
%59 = getelementptr inbounds double, double* %54, i64 %58
store double* %59, double** %18
%60 = load double*, double** %10
%61 = load i64, i64* %12
%62 = load i64, i64* %14
%63 = mul nsw i64 %61, %62
%64 = load i64, i64* %16
%65 = add nsw i64 %63, %64
%66 = mul nsw i64 2, %65
%67 = getelementptr inbounds double, double* %60, i64 %66
store double* %67, double** %19
store i64 0, i64* %11
br label %68
68: ; preds = %152, %53
%69 = load i64, i64* %11
%70 = load i64, i64* %16
%71 = icmp slt i64 %69, %70
br i1 %71, label %72, label %155
72: ; preds = %68
%73 = load double*, double** %17
%74 = load i64, i64* %11
%75 = mul nsw i64 2, %74
%76 = getelementptr inbounds double, double* %73, i64 %75
%77 = load double, double* %76
store double %77, double* %20
%78 = load i64, i64* %6
%79 = sitofp i64 %78 to double
%80 = load double*, double** %17
%81 = load i64, i64* %11
%82 = mul nsw i64 2, %81
%83 = add nsw i64 %82, 1
%84 = getelementptr inbounds double, double* %80, i64 %83
%85 = load double, double* %84
%86 = fmul double %79, %85
store double %86, double* %21
%87 = load double*, double** %19
%88 = load i64, i64* %11
%89 = mul nsw i64 2, %88
%90 = getelementptr inbounds double, double* %87, i64 %89
%91 = load double, double* %90
store double %91, double* %24
%92 = load double*, double** %19
%93 = load i64, i64* %11
%94 = mul nsw i64 2, %93
%95 = add nsw i64 %94, 1
%96 = getelementptr inbounds double, double* %92, i64 %95
%97 = load double, double* %96
store double %97, double* %25
%98 = load double, double* %20
%99 = load double, double* %24
%100 = load double, double* %21
%101 = load double, double* %25
%102 = fmul double %100, %101
%103 = fneg double %102
%104 = call double @llvm.fmuladd.f64(double %98, double %99, double %103)
store double %104, double* %22
%105 = load double, double* %20
%106 = load double, double* %25
%107 = load double, double* %21
%108 = load double, double* %24
%109 = fmul double %107, %108
%110 = call double @llvm.fmuladd.f64(double %105, double %106, double %109)
store double %110, double* %23
%111 = load double*, double** %18
%112 = load i64, i64* %11
%113 = mul nsw i64 2, %112
%114 = getelementptr inbounds double, double* %111, i64 %113
%115 = load double, double* %114
store double %115, double* %24
%116 = load double*, double** %18
%117 = load i64, i64* %11
%118 = mul nsw i64 2, %117
%119 = add nsw i64 %118, 1
%120 = getelementptr inbounds double, double* %116, i64 %119
%121 = load double, double* %120
store double %121, double* %25
%122 = load double, double* %24
%123 = load double, double* %22
%124 = fsub double %122, %123
%125 = load double*, double** %19
%126 = load i64, i64* %11
%127 = mul nsw i64 2, %126
%128 = getelementptr inbounds double, double* %125, i64 %127
store double %124, double* %128
%129 = load double, double* %25
%130 = load double, double* %23
%131 = fsub double %129, %130
%132 = load double*, double** %19
%133 = load i64, i64* %11
%134 = mul nsw i64 2, %133
%135 = add nsw i64 %134, 1
%136 = getelementptr inbounds double, double* %132, i64 %135
store double %131, double* %136
%137 = load double, double* %24
%138 = load double, double* %22
%139 = fadd double %137, %138
%140 = load double*, double** %18
%141 = load i64, i64* %11
%142 = mul nsw i64 2, %141
%143 = getelementptr inbounds double, double* %140, i64 %142
store double %139, double* %143
%144 = load double, double* %25
%145 = load double, double* %23
%146 = fadd double %144, %145
%147 = load double*, double** %18
%148 = load i64, i64* %11
%149 = mul nsw i64 2, %148
%150 = add nsw i64 %149, 1
%151 = getelementptr inbounds double, double* %147, i64 %150
store double %146, double* %151
br label %152
152: ; preds = %72
%153 = load i64, i64* %11
%154 = add nsw i64 %153, 1
store i64 %154, i64* %11
br label %68
155: ; preds = %68
br label %156
156: ; preds = %155
%157 = load i64, i64* %12
%158 = add nsw i64 %157, 1
store i64 %158, i64* %12
br label %49
159: ; preds = %49
br label %160
160: ; preds = %159
%161 = load i64, i64* %13
%162 = add nsw i64 %161, 1
store i64 %162, i64* %13
br label %30
163: ; preds = %30
ret void
}

LLVM IR (-O2):
define void @FFT1DOnce(i64 %0, i64 %1, i64 %2, double* nocapture %3, double* %4) {
%6 = tail call i32 @Reverse(i64 %2, i64 %1, double* %4)
%7 = icmp slt i64 %1, 1
br i1 %7, label %134, label %8
8: ; preds = %5
%9 = sitofp i64 %0 to double
%10 = getelementptr double, double* %3, i64 -2
%11 = insertelement <2 x double> poison, double %9, i64 0
%12 = shufflevector <2 x double> %11, <2 x double> poison, <2 x i32> zeroinitializer
br label %13
13: ; preds = %8, %131
%14 = phi i64 [ 1, %8 ], [ %132, %131 ]
%15 = trunc i64 %14 to i32
%16 = shl nuw i32 1, %15
%17 = sext i32 %16 to i64
%18 = sdiv i64 %2, %17
%19 = sdiv i32 %16, 2
%20 = sext i32 %19 to i64
%21 = shl nsw i64 %20, 1
%22 = add nsw i64 %21, -2
%23 = getelementptr double, double* %3, i64 %22
%24 = icmp sgt i64 %18, 0
br i1 %24, label %25, label %131
25: ; preds = %13
%26 = icmp sgt i32 %16, 1
%27 = shl nsw i64 %17, 1
%28 = shl nsw i64 %20, 2
%29 = getelementptr double, double* %10, i64 %28
%30 = icmp ult i32 %19, 2
%31 = and i64 %20, -2
%32 = icmp eq i64 %31, %20
br label %33
33: ; preds = %25, %128
%34 = phi i64 [ 0, %25 ], [ %129, %128 ]
%35 = mul i64 %27, %34
%36 = getelementptr double, double* %4, i64 %35
%37 = add i64 %21, %35
%38 = getelementptr double, double* %4, i64 %37
%39 = add i64 %28, %35
%40 = getelementptr double, double* %4, i64 %39
%41 = mul nsw i64 %34, %17
%42 = shl nsw i64 %41, 1
%43 = getelementptr inbounds double, double* %4, i64 %42
%44 = add nsw i64 %41, %20
%45 = shl nsw i64 %44, 1
%46 = getelementptr inbounds double, double* %4, i64 %45
br i1 %26, label %47, label %128
47: ; preds = %33
br i1 %30, label %97, label %48
48: ; preds = %47
%49 = icmp ult double* %36, %29
%50 = icmp ult double* %23, %38
%51 = and i1 %49, %50
%52 = icmp ult double* %38, %29
%53 = icmp ult double* %23, %40
%54 = and i1 %52, %53
%55 = or i1 %51, %54
br i1 %55, label %97, label %56
56: ; preds = %48
%57 = getelementptr inbounds double, double* %46, i64 -1
%58 = getelementptr inbounds double, double* %43, i64 -1
br label %59
59: ; preds = %59, %56
%60 = phi i64 [ 0, %56 ], [ %94, %59 ]
%61 = shl nuw nsw i64 %60, 1
%62 = getelementptr inbounds double, double* %23, i64 %61
%63 = bitcast double* %62 to <4 x double>*
%64 = load <4 x double>, <4 x double>* %63
%65 = shufflevector <4 x double> %64, <4 x double> poison, <2 x i32> <i32 0, i32 2>
%66 = shufflevector <4 x double> %64, <4 x double> poison, <2 x i32> <i32 1, i32 3>
%67 = or i64 %61, 1
%68 = fmul <2 x double> %66, %12
%69 = getelementptr inbounds double, double* %46, i64 %61
%70 = bitcast double* %69 to <4 x double>*
%71 = load <4 x double>, <4 x double>* %70
%72 = shufflevector <4 x double> %71, <4 x double> poison, <2 x i32> <i32 0, i32 2>
%73 = shufflevector <4 x double> %71, <4 x double> poison, <2 x i32> <i32 1, i32 3>
%74 = fneg <2 x double> %68
%75 = fmul <2 x double> %73, %74
%76 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %65, <2 x double> %72, <2 x double> %75)
%77 = fmul <2 x double> %68, %72
%78 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %65, <2 x double> %73, <2 x double> %77)
%79 = getelementptr inbounds double, double* %43, i64 %61
%80 = bitcast double* %79 to <4 x double>*
%81 = load <4 x double>, <4 x double>* %80
%82 = shufflevector <4 x double> %81, <4 x double> poison, <2 x i32> <i32 0, i32 2>
%83 = shufflevector <4 x double> %81, <4 x double> poison, <2 x i32> <i32 1, i32 3>
%84 = fsub <2 x double> %82, %76
%85 = fsub <2 x double> %83, %78
%86 = getelementptr inbounds double, double* %57, i64 %67
%87 = bitcast double* %86 to <4 x double>*
%88 = shufflevector <2 x double> %84, <2 x double> %85, <4 x i32> <i32 0, i32 2, i32 1, i32 3>
store <4 x double> %88, <4 x double>* %87
%89 = fadd <2 x double> %82, %76
%90 = fadd <2 x double> %78, %83
%91 = getelementptr inbounds double, double* %58, i64 %67
%92 = bitcast double* %91 to <4 x double>*
%93 = shufflevector <2 x double> %89, <2 x double> %90, <4 x i32> <i32 0, i32 2, i32 1, i32 3>
store <4 x double> %93, <4 x double>* %92
%94 = add nuw i64 %60, 2
%95 = icmp eq i64 %94, %31
br i1 %95, label %96, label %59
96: ; preds = %59
br i1 %32, label %128, label %97
97: ; preds = %48, %47, %96
%98 = phi i64 [ 0, %48 ], [ 0, %47 ], [ %31, %96 ]
br label %99
99: ; preds = %97, %99
%100 = phi i64 [ %126, %99 ], [ %98, %97 ]
%101 = shl nuw nsw i64 %100, 1
%102 = getelementptr inbounds double, double* %23, i64 %101
%103 = load double, double* %102
%104 = or i64 %101, 1
%105 = getelementptr inbounds double, double* %23, i64 %104
%106 = load double, double* %105
%107 = fmul double %106, %9
%108 = getelementptr inbounds double, double* %46, i64 %101
%109 = bitcast double* %108 to <2 x double>*
%110 = load <2 x double>, <2 x double>* %109
%111 = fneg double %107
%112 = insertelement <2 x double> poison, double %107, i64 0
%113 = insertelement <2 x double> %112, double %111, i64 1
%114 = fmul <2 x double> %110, %113
%115 = shufflevector <2 x double> %114, <2 x double> poison, <2 x i32> <i32 1, i32 0>
%116 = insertelement <2 x double> poison, double %103, i64 0
%117 = shufflevector <2 x double> %116, <2 x double> poison, <2 x i32> zeroinitializer
%118 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %117, <2 x double> %110, <2 x double> %115)
%119 = getelementptr inbounds double, double* %43, i64 %101
%120 = bitcast double* %119 to <2 x double>*
%121 = load <2 x double>, <2 x double>* %120
%122 = fsub <2 x double> %121, %118
%123 = bitcast double* %108 to <2 x double>*
store <2 x double> %122, <2 x double>* %123
%124 = fadd <2 x double> %121, %118
%125 = bitcast double* %119 to <2 x double>*
store <2 x double> %124, <2 x double>* %125
%126 = add nuw nsw i64 %100, 1
%127 = icmp eq i64 %126, %20
br i1 %127, label %128, label %99
128: ; preds = %99, %96, %33
%129 = add nuw nsw i64 %34, 1
%130 = icmp eq i64 %129, %18
br i1 %130, label %131, label %33
131: ; preds = %128, %13
%132 = add nuw i64 %14, 1
%133 = icmp eq i64 %14, %1
br i1 %133, label %134, label %13
134: ; preds = %131, %5
ret void
}