void twiddleFFT_2(int sofarRadix, int radix, int remainRadix, int groupOffset,
                  double twiddleRe[], double twiddleIm[], double yRe[],
                  double yIm[]) {
  int yIndx[2];
  int groupNo, blockNo, adr;
  double yr, yi, twr, twi;
  double ar0, ar1;
  double ai0, ai1;
  adr = groupOffset;
  for (groupNo = 0; groupNo < remainRadix; groupNo++) {
    for (blockNo = 0; blockNo < radix; blockNo++) {
      yIndx[blockNo] = adr;
      adr = adr + sofarRadix;
    }
    ar0 = yRe[yIndx[0]];
    ai0 = yIm[yIndx[0]];
    yr = yRe[yIndx[1]];
    yi = yIm[yIndx[1]];
    twr = twiddleRe[1];
    twi = twiddleIm[1];
    ar1 = twr * yr - twi * yi;
    ai1 = twr * yi + twi * yr;
    yRe[yIndx[0]] = ar0 + ar1;
    yRe[yIndx[1]] = ar0 - ar1;
    yIm[yIndx[0]] = ai0 + ai1;
    yIm[yIndx[1]] = ai0 - ai1;
    groupOffset = groupOffset + sofarRadix * radix;
    adr = groupOffset;
  }
}


x86 assembly (clang, -O0):
<twiddleFFT_2>:
push   %rbp
mov    %rsp,%rbp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %ecx,-0x10(%rbp)
mov    %r8,-0x18(%rbp)
mov    %r9,-0x20(%rbp)
mov    -0x10(%rbp),%eax
mov    %eax,-0x34(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12aa <twiddleFFT_2+0x1aa>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    116e <twiddleFFT_2+0x6e>
mov    -0x34(%rbp),%ecx
movslq -0x30(%rbp),%rax
mov    %ecx,-0x28(%rbp,%rax,4)
mov    -0x34(%rbp),%eax
add    -0x4(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1140 <twiddleFFT_2+0x40>
mov    0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x20(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x50(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm2
movsd  -0x58(%rbp),%xmm1
mulsd  -0x48(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x50(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm2
movsd  -0x58(%rbp),%xmm1
mulsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x60(%rbp),%xmm0
addsd  -0x68(%rbp),%xmm0
mov    0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x60(%rbp),%xmm0
subsd  -0x68(%rbp),%xmm0
mov    0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm0
addsd  -0x78(%rbp),%xmm0
mov    0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm0
subsd  -0x78(%rbp),%xmm0
mov    0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%eax
mov    -0x4(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0x10(%rbp)
mov    -0x10(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    112d <twiddleFFT_2+0x2d>
pop    %rbp
ret

x86 assembly (clang, -O2):
<twiddleFFT_2>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %r9,-0x8(%rsp)
mov    %r8,-0x10(%rsp)
mov    %edx,-0x3c(%rsp)
test   %edx,%edx
jle    1374 <twiddleFFT_2+0x274>
mov    0x40(%rsp),%r10
mov    0x38(%rsp),%r14
mov    %esi,%edx
imul   %edi,%edx
mov    %esi,%r8d
mov    %r8d,%eax
and    $0xfffffff8,%eax
lea    -0x8(%rax),%rbp
mov    %rbp,-0x18(%rsp)
shr    $0x3,%rbp
add    $0x1,%rbp
mov    %rax,-0x38(%rsp)
imul   %edi,%eax
mov    %rax,-0x30(%rsp)
movd   %edi,%xmm0
pshufd $0x0,%xmm0,%xmm0
movdqa 0xe99(%rip),%xmm1
pmuludq %xmm0,%xmm1
pshufd $0xe8,%xmm1,%xmm2
movdqa 0xe98(%rip),%xmm1
pmuludq %xmm0,%xmm1
pshufd $0xe8,%xmm1,%xmm8
punpckldq %xmm2,%xmm8
pslld  $0x2,%xmm0
mov    %ebp,%r11d
and    $0x3,%r11d
and    $0xfffffffffffffffc,%rbp
mov    %rbp,-0x28(%rsp)
mov    %r11,-0x20(%rsp)
shl    $0x5,%r11
xor    %r13d,%r13d
mov    -0x10(%rsp),%rbx
mov    -0x8(%rsp),%r9
jmp    1250 <twiddleFFT_2+0x150>
nopw   0x0(%rax,%rax,1)
mov    -0x44(%rsp),%r12d
mov    -0x40(%rsp),%r15d
movslq %r12d,%rax
movsd  (%r14,%rax,8),%xmm2
movsd  (%r10,%rax,8),%xmm3
movslq %r15d,%rbp
movsd  (%r14,%rbp,8),%xmm4
movsd  (%r10,%rbp,8),%xmm5
movsd  0x8(%rbx),%xmm1
movsd  0x8(%r9),%xmm7
movapd %xmm1,%xmm6
mulsd  %xmm5,%xmm1
mulsd  %xmm7,%xmm5
mulsd  %xmm4,%xmm6
subsd  %xmm5,%xmm6
mulsd  %xmm7,%xmm4
addsd  %xmm4,%xmm1
movapd %xmm2,%xmm4
addsd  %xmm6,%xmm4
movsd  %xmm4,(%r14,%rax,8)
subsd  %xmm6,%xmm2
movsd  %xmm2,(%r14,%rbp,8)
movapd %xmm3,%xmm2
addsd  %xmm1,%xmm2
movsd  %xmm2,(%r10,%rax,8)
subsd  %xmm1,%xmm3
movsd  %xmm3,(%r10,%rbp,8)
add    %edx,%ecx
add    $0x1,%r13d
cmp    -0x3c(%rsp),%r13d
je     1374 <twiddleFFT_2+0x274>
test   %esi,%esi
jle    11ca <twiddleFFT_2+0xca>
cmp    $0x8,%esi
jae    1270 <twiddleFFT_2+0x170>
xor    %eax,%eax
mov    %ecx,%ebp
jmp    1360 <twiddleFFT_2+0x260>
cs nopw 0x0(%rax,%rax,1)
movd   %ecx,%xmm2
pshufd $0x0,%xmm2,%xmm2
paddd  %xmm8,%xmm2
cmpq   $0x18,-0x18(%rsp)
jae    128a <twiddleFFT_2+0x18a>
xor    %eax,%eax
jmp    12fa <twiddleFFT_2+0x1fa>
mov    -0x28(%rsp),%rbp
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movdqu %xmm2,-0x44(%rsp,%rax,4)
paddd  %xmm0,%xmm2
movdqu %xmm2,-0x34(%rsp,%rax,4)
paddd  %xmm0,%xmm2
movdqu %xmm2,-0x24(%rsp,%rax,4)
paddd  %xmm0,%xmm2
movdqu %xmm2,-0x14(%rsp,%rax,4)
paddd  %xmm0,%xmm2
movdqu %xmm2,-0x4(%rsp,%rax,4)
paddd  %xmm0,%xmm2
movdqu %xmm2,0xc(%rsp,%rax,4)
paddd  %xmm0,%xmm2
movdqu %xmm2,0x1c(%rsp,%rax,4)
paddd  %xmm0,%xmm2
movdqu %xmm2,0x2c(%rsp,%rax,4)
add    $0x20,%rax
paddd  %xmm0,%xmm2
add    $0xfffffffffffffffc,%rbp
jne    12a0 <twiddleFFT_2+0x1a0>
cmpq   $0x0,-0x20(%rsp)
je     133c <twiddleFFT_2+0x23c>
lea    -0x34(%rsp),%rbp
lea    0x0(,%rax,4),%rax
add    %rbp,%rax
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movdqu %xmm2,-0x10(%rax,%rbp,1)
paddd  %xmm0,%xmm2
movdqu %xmm2,(%rax,%rbp,1)
paddd  %xmm0,%xmm2
add    $0x20,%rbp
cmp    %rbp,%r11
jne    1320 <twiddleFFT_2+0x220>
cmp    %r8,-0x38(%rsp)
je     11c0 <twiddleFFT_2+0xc0>
mov    -0x30(%rsp),%rax
lea    (%rcx,%rax,1),%ebp
mov    -0x38(%rsp),%rax
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    %ebp,-0x44(%rsp,%rax,4)
add    %edi,%ebp
add    $0x1,%rax
cmp    %rax,%r8
jne    1360 <twiddleFFT_2+0x260>
jmp    11c0 <twiddleFFT_2+0xc0>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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