void fft(data, ct, st, nmodes, length) float *data, *ct, *st;
int nmodes, length;
{
  double *im, *re;
  int dim[2], i;
  dim[0] = length;
  re = (double *)malloc(length * sizeof(double));
  im = (double *)malloc(length * sizeof(double));
  for (i = 0; i < length; i++) {
    im[i] = 0.0;
    re[i] = data[i];
  }
  fftn(1, dim, re, im, 1, -1);
  ct[0] = re[0];
  st[0] = 0.0;
  for (i = 1; i < nmodes; i++) {
    ct[i] = re[i] * 2.0;
    st[i] = im[i] * 2.0;
  }
  free(im);
  free(re);
}


x86 assembly (clang, -O0):
<fft>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    %eax,-0x38(%rbp)
movslq -0x20(%rbp),%rdi
shl    $0x3,%rdi
call   1050 <malloc@plt>
mov    %rax,-0x30(%rbp)
movslq -0x20(%rbp),%rdi
shl    $0x3,%rdi
call   1050 <malloc@plt>
mov    %rax,-0x28(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    11c2 <fft+0x92>
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    -0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    117a <fft+0x4a>
lea    -0x38(%rbp),%rsi
mov    -0x30(%rbp),%rdx
mov    -0x28(%rbp),%rcx
mov    $0x1,%r8d
mov    $0xffffffff,%r9d
mov    %r8d,%edi
mov    $0x0,%al
call   1040 <fftn@plt>
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x10(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x18(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax)
movl   $0x1,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1270 <fft+0x140>
mov    -0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  0xdda(%rip),%xmm0
mulsd  (%rax,%rcx,8),%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  0xdb4(%rip),%xmm0
mulsd  (%rax,%rcx,8),%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    120a <fft+0xda>
mov    -0x28(%rbp),%rdi
call   1030 <free@plt>
mov    -0x30(%rbp),%rdi
call   1030 <free@plt>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<fft>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %r8d,%r13d
mov    %ecx,0xc(%rsp)
mov    %rdx,0x18(%rsp)
mov    %rsi,0x10(%rsp)
mov    %rdi,%rbp
mov    %r8d,0x20(%rsp)
movslq %r8d,%r14
lea    0x0(,%r14,8),%r12
mov    %r12,%rdi
call   1060 <malloc@plt>
mov    %rax,%rbx
mov    %r12,%rdi
call   1060 <malloc@plt>
mov    %rax,%r12
test   %r14d,%r14d
jle    124a <fft+0x10a>
mov    %r13d,%r15d
lea    0x0(,%r15,8),%rdx
xor    %r14d,%r14d
mov    %r12,%rdi
xor    %esi,%esi
call   1050 <memset@plt>
cmp    $0x4,%r13d
jb     1230 <fft+0xf0>
mov    %r15d,%r14d
and    $0xfffffffc,%r14d
lea    -0x4(%r14),%rcx
mov    %rcx,%rax
shr    $0x2,%rax
add    $0x1,%rax
test   %rcx,%rcx
je     1447 <fft+0x307>
mov    %rax,%rdx
and    $0xfffffffffffffffe,%rdx
xor    %ecx,%ecx
nopl   0x0(%rax,%rax,1)
cvtps2pd 0x0(%rbp,%rcx,4),%xmm0
cvtps2pd 0x8(%rbp,%rcx,4),%xmm1
movups %xmm0,(%rbx,%rcx,8)
movups %xmm1,0x10(%rbx,%rcx,8)
cvtps2pd 0x10(%rbp,%rcx,4),%xmm0
cvtps2pd 0x18(%rbp,%rcx,4),%xmm1
movups %xmm0,0x20(%rbx,%rcx,8)
movups %xmm1,0x30(%rbx,%rcx,8)
add    $0x8,%rcx
add    $0xfffffffffffffffe,%rdx
jne    11e0 <fft+0xa0>
test   $0x1,%al
je     1228 <fft+0xe8>
cvtps2pd 0x0(%rbp,%rcx,4),%xmm0
cvtps2pd 0x8(%rbp,%rcx,4),%xmm1
movups %xmm0,(%rbx,%rcx,8)
movups %xmm1,0x10(%rbx,%rcx,8)
cmp    %r15,%r14
je     124a <fft+0x10a>
nopl   (%rax)
movss  0x0(%rbp,%r14,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,(%rbx,%r14,8)
add    $0x1,%r14
cmp    %r14,%r15
jne    1230 <fft+0xf0>
lea    0x20(%rsp),%rsi
mov    $0x1,%edi
mov    %rbx,%rdx
mov    %r12,%rcx
mov    $0x1,%r8d
mov    $0xffffffff,%r9d
xor    %eax,%eax
call   1040 <fftn@plt>
movsd  (%rbx),%xmm0
cvtsd2ss %xmm0,%xmm0
mov    0x10(%rsp),%rcx
movss  %xmm0,(%rcx)
mov    0x18(%rsp),%rbp
movl   $0x0,0x0(%rbp)
mov    0xc(%rsp),%eax
cmp    $0x2,%eax
jl     1428 <fft+0x2e8>
mov    %eax,%eax
lea    -0x1(%rax),%r9
mov    $0x1,%edx
cmp    $0x2,%r9
jb     1394 <fft+0x254>
lea    0x4(%rcx),%rsi
lea    0x0(,%rax,4),%rdi
add    %rbp,%rdi
cmp    %rdi,%rsi
jae    12d1 <fft+0x191>
lea    (%rcx,%rax,4),%rsi
lea    0x4(%rbp),%rdi
cmp    %rsi,%rdi
jb     1394 <fft+0x254>
mov    %r9,%r8
and    $0xfffffffffffffffe,%r8
lea    -0x2(%r8),%rsi
mov    %rsi,%rdx
shr    %rdx
add    $0x1,%rdx
test   %rsi,%rsi
je     1456 <fft+0x316>
mov    %rdx,%rdi
and    $0xfffffffffffffffe,%rdi
mov    $0x1,%esi
nopl   0x0(%rax,%rax,1)
movupd (%rbx,%rsi,8),%xmm0
addpd  %xmm0,%xmm0
cvtpd2ps %xmm0,%xmm0
movlpd %xmm0,(%rcx,%rsi,4)
movupd (%r12,%rsi,8),%xmm0
addpd  %xmm0,%xmm0
cvtpd2ps %xmm0,%xmm0
movlpd %xmm0,0x0(%rbp,%rsi,4)
movupd 0x10(%rbx,%rsi,8),%xmm0
addpd  %xmm0,%xmm0
cvtpd2ps %xmm0,%xmm0
movlpd %xmm0,0x8(%rcx,%rsi,4)
movupd 0x10(%r12,%rsi,8),%xmm0
addpd  %xmm0,%xmm0
cvtpd2ps %xmm0,%xmm0
movlpd %xmm0,0x8(%rbp,%rsi,4)
add    $0x4,%rsi
add    $0xfffffffffffffffe,%rdi
jne    1300 <fft+0x1c0>
test   $0x1,%dl
je     1384 <fft+0x244>
movupd (%rbx,%rsi,8),%xmm0
addpd  %xmm0,%xmm0
cvtpd2ps %xmm0,%xmm0
movlpd %xmm0,(%rcx,%rsi,4)
movupd (%r12,%rsi,8),%xmm0
addpd  %xmm0,%xmm0
cvtpd2ps %xmm0,%xmm0
movlpd %xmm0,0x0(%rbp,%rsi,4)
cmp    %r8,%r9
je     1428 <fft+0x2e8>
or     $0x1,%r9
mov    %r9,%rdx
mov    %rdx,%rsi
not    %rsi
test   $0x1,%al
jne    13c8 <fft+0x288>
movsd  (%rbx,%rdx,8),%xmm0
addsd  %xmm0,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,(%rcx,%rdx,4)
movsd  (%r12,%rdx,8),%xmm0
addsd  %xmm0,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0x0(%rbp,%rdx,4)
add    $0x1,%rdx
add    %rax,%rsi
je     1428 <fft+0x2e8>
nopl   (%rax)
movsd  (%rbx,%rdx,8),%xmm0
addsd  %xmm0,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,(%rcx,%rdx,4)
movsd  (%r12,%rdx,8),%xmm0
addsd  %xmm0,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0x0(%rbp,%rdx,4)
movsd  0x8(%rbx,%rdx,8),%xmm0
addsd  %xmm0,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0x4(%rcx,%rdx,4)
movsd  0x8(%r12,%rdx,8),%xmm0
addsd  %xmm0,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0x4(%rbp,%rdx,4)
add    $0x2,%rdx
cmp    %rdx,%rax
jne    13d0 <fft+0x290>
mov    %r12,%rdi
call   1030 <free@plt>
mov    %rbx,%rdi
call   1030 <free@plt>
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %ecx,%ecx
test   $0x1,%al
jne    1215 <fft+0xd5>
jmp    1228 <fft+0xe8>
mov    $0x1,%esi
test   $0x1,%dl
jne    135e <fft+0x21e>
jmp    1384 <fft+0x244>

LLVM IR (-O0):
define void @fft(float* %0, float* %1, float* %2, i32 %3, i32 %4) {
%6 = alloca float*
%7 = alloca float*
%8 = alloca float*
%9 = alloca i32
%10 = alloca i32
%11 = alloca double*
%12 = alloca double*
%13 = alloca [2 x i32]
%14 = alloca i32
store float* %0, float** %6
store float* %1, float** %7
store float* %2, float** %8
store i32 %3, i32* %9
store i32 %4, i32* %10
%15 = load i32, i32* %10
%16 = getelementptr inbounds [2 x i32], [2 x i32]* %13, i64 0, i64 0
store i32 %15, i32* %16
%17 = load i32, i32* %10
%18 = sext i32 %17 to i64
%19 = mul i64 %18, 8
%20 = call i8* @malloc(i64 %19)
%21 = bitcast i8* %20 to double*
store double* %21, double** %12
%22 = load i32, i32* %10
%23 = sext i32 %22 to i64
%24 = mul i64 %23, 8
%25 = call i8* @malloc(i64 %24)
%26 = bitcast i8* %25 to double*
store double* %26, double** %11
store i32 0, i32* %14
br label %27
27: ; preds = %46, %5
%28 = load i32, i32* %14
%29 = load i32, i32* %10
%30 = icmp slt i32 %28, %29
br i1 %30, label %31, label %49
31: ; preds = %27
%32 = load double*, double** %11
%33 = load i32, i32* %14
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds double, double* %32, i64 %34
store double 0.000000e+00, double* %35
%36 = load float*, float** %6
%37 = load i32, i32* %14
%38 = sext i32 %37 to i64
%39 = getelementptr inbounds float, float* %36, i64 %38
%40 = load float, float* %39
%41 = fpext float %40 to double
%42 = load double*, double** %12
%43 = load i32, i32* %14
%44 = sext i32 %43 to i64
%45 = getelementptr inbounds double, double* %42, i64 %44
store double %41, double* %45
br label %46
46: ; preds = %31
%47 = load i32, i32* %14
%48 = add nsw i32 %47, 1
store i32 %48, i32* %14
br label %27
49: ; preds = %27
%50 = getelementptr inbounds [2 x i32], [2 x i32]* %13, i64 0, i64 0
%51 = load double*, double** %12
%52 = load double*, double** %11
%53 = call i32 (i32, i32*, double*, double*, i32, i32, ...) bitcast (i32 (...)* @fftn to i32 (i32, i32*, double*, double*, i32, i32, ...)*)(i32 1, i32* %50, double* %51, double* %52, i32 1, i32 -1)
%54 = load double*, double** %12
%55 = getelementptr inbounds double, double* %54, i64 0
%56 = load double, double* %55
%57 = fptrunc double %56 to float
%58 = load float*, float** %7
%59 = getelementptr inbounds float, float* %58, i64 0
store float %57, float* %59
%60 = load float*, float** %8
%61 = getelementptr inbounds float, float* %60, i64 0
store float 0.000000e+00, float* %61
store i32 1, i32* %14
br label %62
62: ; preds = %89, %49
%63 = load i32, i32* %14
%64 = load i32, i32* %9
%65 = icmp slt i32 %63, %64
br i1 %65, label %66, label %92
66: ; preds = %62
%67 = load double*, double** %12
%68 = load i32, i32* %14
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds double, double* %67, i64 %69
%71 = load double, double* %70
%72 = fmul double %71, 2.000000e+00
%73 = fptrunc double %72 to float
%74 = load float*, float** %7
%75 = load i32, i32* %14
%76 = sext i32 %75 to i64
%77 = getelementptr inbounds float, float* %74, i64 %76
store float %73, float* %77
%78 = load double*, double** %11
%79 = load i32, i32* %14
%80 = sext i32 %79 to i64
%81 = getelementptr inbounds double, double* %78, i64 %80
%82 = load double, double* %81
%83 = fmul double %82, 2.000000e+00
%84 = fptrunc double %83 to float
%85 = load float*, float** %8
%86 = load i32, i32* %14
%87 = sext i32 %86 to i64
%88 = getelementptr inbounds float, float* %85, i64 %87
store float %84, float* %88
br label %89
89: ; preds = %66
%90 = load i32, i32* %14
%91 = add nsw i32 %90, 1
store i32 %91, i32* %14
br label %62
92: ; preds = %62
%93 = load double*, double** %11
%94 = bitcast double* %93 to i8*
%95 = load double*, double** %12
%96 = bitcast double* %95 to i8*
ret void
}

LLVM IR (-O2):
define void @fft(float* nocapture %0, float* nocapture writeonly %1, float* nocapture writeonly %2, i32 %3, i32 %4) {
%6 = alloca [2 x i32]
%7 = bitcast [2 x i32]* %6 to i8*
%8 = getelementptr inbounds [2 x i32], [2 x i32]* %6, i64 0, i64 0
store i32 %4, i32* %8
%9 = sext i32 %4 to i64
%10 = shl nsw i64 %9, 3
%11 = tail call i8* @malloc(i64 %10)
%12 = bitcast i8* %11 to double*
%13 = tail call i8* @malloc(i64 %10)
%14 = bitcast i8* %13 to double*
%15 = icmp sgt i32 %4, 0
br i1 %15, label %16, label %88
16: ; preds = %5
%17 = zext i32 %4 to i64
%18 = shl nuw nsw i64 %17, 3
%19 = icmp ult i32 %4, 4
br i1 %19, label %78, label %20
20: ; preds = %16
%21 = and i64 %17, 4294967292
%22 = add nsw i64 %21, -4
%23 = lshr exact i64 %22, 2
%24 = add nuw nsw i64 %23, 1
%25 = and i64 %24, 1
%26 = icmp eq i64 %22, 0
br i1 %26, label %60, label %27
27: ; preds = %20
%28 = and i64 %24, 9223372036854775806
br label %29
29: ; preds = %29, %27
%30 = phi i64 [ 0, %27 ], [ %57, %29 ]
%31 = phi i64 [ 0, %27 ], [ %58, %29 ]
%32 = getelementptr inbounds float, float* %0, i64 %30
%33 = bitcast float* %32 to <2 x float>*
%34 = load <2 x float>, <2 x float>* %33
%35 = getelementptr inbounds float, float* %32, i64 2
%36 = bitcast float* %35 to <2 x float>*
%37 = load <2 x float>, <2 x float>* %36
%38 = fpext <2 x float> %34 to <2 x double>
%39 = fpext <2 x float> %37 to <2 x double>
%40 = getelementptr inbounds double, double* %12, i64 %30
%41 = bitcast double* %40 to <2 x double>*
store <2 x double> %38, <2 x double>* %41
%42 = getelementptr inbounds double, double* %40, i64 2
%43 = bitcast double* %42 to <2 x double>*
store <2 x double> %39, <2 x double>* %43
%44 = or i64 %30, 4
%45 = getelementptr inbounds float, float* %0, i64 %44
%46 = bitcast float* %45 to <2 x float>*
%47 = load <2 x float>, <2 x float>* %46
%48 = getelementptr inbounds float, float* %45, i64 2
%49 = bitcast float* %48 to <2 x float>*
%50 = load <2 x float>, <2 x float>* %49
%51 = fpext <2 x float> %47 to <2 x double>
%52 = fpext <2 x float> %50 to <2 x double>
%53 = getelementptr inbounds double, double* %12, i64 %44
%54 = bitcast double* %53 to <2 x double>*
store <2 x double> %51, <2 x double>* %54
%55 = getelementptr inbounds double, double* %53, i64 2
%56 = bitcast double* %55 to <2 x double>*
store <2 x double> %52, <2 x double>* %56
%57 = add nuw i64 %30, 8
%58 = add i64 %31, 2
%59 = icmp eq i64 %58, %28
br i1 %59, label %60, label %29
60: ; preds = %29, %20
%61 = phi i64 [ 0, %20 ], [ %57, %29 ]
%62 = icmp eq i64 %25, 0
br i1 %62, label %76, label %63
63: ; preds = %60
%64 = getelementptr inbounds float, float* %0, i64 %61
%65 = bitcast float* %64 to <2 x float>*
%66 = load <2 x float>, <2 x float>* %65
%67 = getelementptr inbounds float, float* %64, i64 2
%68 = bitcast float* %67 to <2 x float>*
%69 = load <2 x float>, <2 x float>* %68
%70 = fpext <2 x float> %66 to <2 x double>
%71 = fpext <2 x float> %69 to <2 x double>
%72 = getelementptr inbounds double, double* %12, i64 %61
%73 = bitcast double* %72 to <2 x double>*
store <2 x double> %70, <2 x double>* %73
%74 = getelementptr inbounds double, double* %72, i64 2
%75 = bitcast double* %74 to <2 x double>*
store <2 x double> %71, <2 x double>* %75
br label %76
76: ; preds = %60, %63
%77 = icmp eq i64 %21, %17
br i1 %77, label %88, label %78
78: ; preds = %16, %76
%79 = phi i64 [ 0, %16 ], [ %21, %76 ]
br label %80
80: ; preds = %78, %80
%81 = phi i64 [ %86, %80 ], [ %79, %78 ]
%82 = getelementptr inbounds float, float* %0, i64 %81
%83 = load float, float* %82
%84 = fpext float %83 to double
%85 = getelementptr inbounds double, double* %12, i64 %81
store double %84, double* %85
%86 = add nuw nsw i64 %81, 1
%87 = icmp eq i64 %86, %17
br i1 %87, label %88, label %80
88: ; preds = %80, %76, %5
%89 = call i32 (i32, i32*, double*, double*, i32, i32, ...) bitcast (i32 (...)* @fftn to i32 (i32, i32*, double*, double*, i32, i32, ...)*)(i32 1, i32* %8, double* %12, double* %14, i32 1, i32 -1)
%90 = load double, double* %12
%91 = fptrunc double %90 to float
store float %91, float* %1
store float 0.000000e+00, float* %2
%92 = icmp sgt i32 %3, 1
br i1 %92, label %93, label %219
93: ; preds = %88
%94 = zext i32 %3 to i64
%95 = add nsw i64 %94, -1
%96 = icmp ult i64 %95, 2
br i1 %96, label %173, label %97
97: ; preds = %93
%98 = getelementptr float, float* %1, i64 1
%99 = getelementptr float, float* %1, i64 %94
%100 = getelementptr float, float* %2, i64 1
%101 = getelementptr float, float* %2, i64 %94
%102 = icmp ult float* %98, %101
%103 = icmp ult float* %100, %99
%104 = and i1 %102, %103
br i1 %104, label %173, label %105
105: ; preds = %97
%106 = and i64 %95, -2
%107 = or i64 %95, 1
%108 = add nsw i64 %106, -2
%109 = lshr exact i64 %108, 1
%110 = add nuw i64 %109, 1
%111 = and i64 %110, 1
%112 = icmp eq i64 %108, 0
br i1 %112, label %153, label %113
113: ; preds = %105
%114 = and i64 %110, -2
br label %115
115: ; preds = %115, %113
%116 = phi i64 [ 0, %113 ], [ %148, %115 ]
%117 = phi i64 [ 0, %113 ], [ %149, %115 ]
%118 = or i64 %116, 1
%119 = getelementptr inbounds double, double* %12, i64 %118
%120 = bitcast double* %119 to <2 x double>*
%121 = load <2 x double>, <2 x double>* %120
%122 = fmul <2 x double> %121, <double 2.000000e+00, double 2.000000e+00>
%123 = fptrunc <2 x double> %122 to <2 x float>
%124 = getelementptr inbounds float, float* %1, i64 %118
%125 = bitcast float* %124 to <2 x float>*
store <2 x float> %123, <2 x float>* %125
%126 = getelementptr inbounds double, double* %14, i64 %118
%127 = bitcast double* %126 to <2 x double>*
%128 = load <2 x double>, <2 x double>* %127
%129 = fmul <2 x double> %128, <double 2.000000e+00, double 2.000000e+00>
%130 = fptrunc <2 x double> %129 to <2 x float>
%131 = getelementptr inbounds float, float* %2, i64 %118
%132 = bitcast float* %131 to <2 x float>*
store <2 x float> %130, <2 x float>* %132
%133 = or i64 %116, 3
%134 = getelementptr inbounds double, double* %12, i64 %133
%135 = bitcast double* %134 to <2 x double>*
%136 = load <2 x double>, <2 x double>* %135
%137 = fmul <2 x double> %136, <double 2.000000e+00, double 2.000000e+00>
%138 = fptrunc <2 x double> %137 to <2 x float>
%139 = getelementptr inbounds float, float* %1, i64 %133
%140 = bitcast float* %139 to <2 x float>*
store <2 x float> %138, <2 x float>* %140
%141 = getelementptr inbounds double, double* %14, i64 %133
%142 = bitcast double* %141 to <2 x double>*
%143 = load <2 x double>, <2 x double>* %142
%144 = fmul <2 x double> %143, <double 2.000000e+00, double 2.000000e+00>
%145 = fptrunc <2 x double> %144 to <2 x float>
%146 = getelementptr inbounds float, float* %2, i64 %133
%147 = bitcast float* %146 to <2 x float>*
store <2 x float> %145, <2 x float>* %147
%148 = add nuw i64 %116, 4
%149 = add nuw i64 %117, 2
%150 = icmp eq i64 %149, %114
br i1 %150, label %151, label %115
151: ; preds = %115
%152 = or i64 %148, 1
br label %153
153: ; preds = %151, %105
%154 = phi i64 [ 1, %105 ], [ %152, %151 ]
%155 = icmp eq i64 %111, 0
br i1 %155, label %171, label %156
156: ; preds = %153
%157 = getelementptr inbounds double, double* %12, i64 %154
%158 = bitcast double* %157 to <2 x double>*
%159 = load <2 x double>, <2 x double>* %158
%160 = fmul <2 x double> %159, <double 2.000000e+00, double 2.000000e+00>
%161 = fptrunc <2 x double> %160 to <2 x float>
%162 = getelementptr inbounds float, float* %1, i64 %154
%163 = bitcast float* %162 to <2 x float>*
store <2 x float> %161, <2 x float>* %163
%164 = getelementptr inbounds double, double* %14, i64 %154
%165 = bitcast double* %164 to <2 x double>*
%166 = load <2 x double>, <2 x double>* %165
%167 = fmul <2 x double> %166, <double 2.000000e+00, double 2.000000e+00>
%168 = fptrunc <2 x double> %167 to <2 x float>
%169 = getelementptr inbounds float, float* %2, i64 %154
%170 = bitcast float* %169 to <2 x float>*
store <2 x float> %168, <2 x float>* %170
br label %171
171: ; preds = %153, %156
%172 = icmp eq i64 %95, %106
br i1 %172, label %219, label %173
173: ; preds = %97, %93, %171
%174 = phi i64 [ 1, %97 ], [ 1, %93 ], [ %107, %171 ]
%175 = xor i64 %174, -1
%176 = and i64 %94, 1
%177 = icmp eq i64 %176, 0
br i1 %177, label %178, label %190
178: ; preds = %173
%179 = getelementptr inbounds double, double* %12, i64 %174
%180 = load double, double* %179
%181 = fmul double %180, 2.000000e+00
%182 = fptrunc double %181 to float
%183 = getelementptr inbounds float, float* %1, i64 %174
store float %182, float* %183
%184 = getelementptr inbounds double, double* %14, i64 %174
%185 = load double, double* %184
%186 = fmul double %185, 2.000000e+00
%187 = fptrunc double %186 to float
%188 = getelementptr inbounds float, float* %2, i64 %174
store float %187, float* %188
%189 = add nuw nsw i64 %174, 1
br label %190
190: ; preds = %178, %173
%191 = phi i64 [ %174, %173 ], [ %189, %178 ]
%192 = sub nsw i64 0, %94
%193 = icmp eq i64 %175, %192
br i1 %193, label %219, label %194
194: ; preds = %190, %194
%195 = phi i64 [ %217, %194 ], [ %191, %190 ]
%196 = getelementptr inbounds double, double* %12, i64 %195
%197 = load double, double* %196
%198 = fmul double %197, 2.000000e+00
%199 = fptrunc double %198 to float
%200 = getelementptr inbounds float, float* %1, i64 %195
store float %199, float* %200
%201 = getelementptr inbounds double, double* %14, i64 %195
%202 = load double, double* %201
%203 = fmul double %202, 2.000000e+00
%204 = fptrunc double %203 to float
%205 = getelementptr inbounds float, float* %2, i64 %195
store float %204, float* %205
%206 = add nuw nsw i64 %195, 1
%207 = getelementptr inbounds double, double* %12, i64 %206
%208 = load double, double* %207
%209 = fmul double %208, 2.000000e+00
%210 = fptrunc double %209 to float
%211 = getelementptr inbounds float, float* %1, i64 %206
store float %210, float* %211
%212 = getelementptr inbounds double, double* %14, i64 %206
%213 = load double, double* %212
%214 = fmul double %213, 2.000000e+00
%215 = fptrunc double %214 to float
%216 = getelementptr inbounds float, float* %2, i64 %206
store float %215, float* %216
%217 = add nuw nsw i64 %195, 2
%218 = icmp eq i64 %217, %94
br i1 %218, label %219, label %194
219: ; preds = %190, %194, %171, %88
ret void
}