#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*/ dradf2(int, int, double *, double *, double *);
int /*<<< orphan*/ dradf4(int, int, double *, double *, double *, double *,
                          double *);
int /*<<< orphan*/ dradfg(int, int, int, int, double *, double *, double *,
                          double *, double *, double *);
void drftf1(int n, double *c, double *ch, double *wa, int *ifac) {
  int i, k1, l1, l2;
  int na, kh, nf;
  int ip, iw, ido, idl1, ix2, ix3;
  nf = ifac[1];
  na = 1;
  l2 = n;
  iw = n;
  for (k1 = 0; k1 < nf; k1++) {
    kh = nf - k1;
    ip = ifac[kh + 1];
    l1 = l2 / ip;
    ido = n / l2;
    idl1 = ido * l1;
    iw -= (ip - 1) * ido;
    na = 1 - na;
    if (ip != 4)
      goto L102;
    ix2 = iw + ido;
    ix3 = ix2 + ido;
    if (na != 0)
      dradf4(ido, l1, ch, c, wa + iw - 1, wa + ix2 - 1, wa + ix3 - 1);
    else
      dradf4(ido, l1, c, ch, wa + iw - 1, wa + ix2 - 1, wa + ix3 - 1);
    goto L110;
  L102:
    if (ip != 2)
      goto L104;
    if (na != 0)
      goto L103;
    dradf2(ido, l1, c, ch, wa + iw - 1);
    goto L110;
  L103:
    dradf2(ido, l1, ch, c, wa + iw - 1);
    goto L110;
  L104:
    if (ido == 1)
      na = 1 - na;
    if (na != 0)
      goto L109;
    dradfg(ido, ip, l1, idl1, c, c, c, ch, ch, wa + iw - 1);
    na = 1;
    goto L110;
  L109:
    dradfg(ido, ip, l1, idl1, ch, ch, ch, c, c, wa + iw - 1);
    na = 0;
  L110:
    l2 = l1;
  }
  if (na == 1)
    return;
  for (i = 0; i < n; i++)
    c[i] = ch[i];
}


x86 assembly (clang, -O0):
<drftf1>:
push   %rbp
mov    %rsp,%rbp
push   %r14
push   %rbx
sub    $0x80,%rsp
mov    %edi,-0x14(%rbp)
mov    %rsi,-0x20(%rbp)
mov    %rdx,-0x28(%rbp)
mov    %rcx,-0x30(%rbp)
mov    %r8,-0x38(%rbp)
mov    -0x38(%rbp),%rax
mov    0x4(%rax),%eax
mov    %eax,-0x54(%rbp)
movl   $0x1,-0x4c(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x5c(%rbp)
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x54(%rbp),%eax
jge    1402 <drftf1+0x2d2>
mov    -0x54(%rbp),%eax
sub    -0x40(%rbp),%eax
mov    %eax,-0x50(%rbp)
mov    -0x38(%rbp),%rax
mov    -0x50(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x58(%rbp)
mov    -0x48(%rbp),%eax
cltd
idivl  -0x58(%rbp)
mov    %eax,-0x44(%rbp)
mov    -0x14(%rbp),%eax
cltd
idivl  -0x48(%rbp)
mov    %eax,-0x60(%rbp)
mov    -0x60(%rbp),%eax
imul   -0x44(%rbp),%eax
mov    %eax,-0x64(%rbp)
mov    -0x58(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x60(%rbp),%ecx
mov    -0x5c(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x5c(%rbp)
mov    $0x1,%eax
sub    -0x4c(%rbp),%eax
mov    %eax,-0x4c(%rbp)
cmpl   $0x4,-0x58(%rbp)
je     11e7 <drftf1+0xb7>
jmp    12ad <drftf1+0x17d>
mov    -0x5c(%rbp),%eax
add    -0x60(%rbp),%eax
mov    %eax,-0x68(%rbp)
mov    -0x68(%rbp),%eax
add    -0x60(%rbp),%eax
mov    %eax,-0x6c(%rbp)
cmpl   $0x0,-0x4c(%rbp)
je     1258 <drftf1+0x128>
mov    -0x60(%rbp),%edi
mov    -0x44(%rbp),%esi
mov    -0x28(%rbp),%rdx
mov    -0x20(%rbp),%rcx
mov    -0x30(%rbp),%r8
movslq -0x5c(%rbp),%rax
shl    $0x3,%rax
add    %rax,%r8
add    $0xfffffffffffffff8,%r8
mov    -0x30(%rbp),%r9
movslq -0x68(%rbp),%rax
shl    $0x3,%rax
add    %rax,%r9
add    $0xfffffffffffffff8,%r9
mov    -0x30(%rbp),%rax
movslq -0x6c(%rbp),%r10
shl    $0x3,%r10
add    %r10,%rax
add    $0xfffffffffffffff8,%rax
mov    %rax,(%rsp)
call   1050 <dradf4@plt>
jmp    12a8 <drftf1+0x178>
mov    -0x60(%rbp),%edi
mov    -0x44(%rbp),%esi
mov    -0x20(%rbp),%rdx
mov    -0x28(%rbp),%rcx
mov    -0x30(%rbp),%r8
movslq -0x5c(%rbp),%rax
shl    $0x3,%rax
add    %rax,%r8
add    $0xfffffffffffffff8,%r8
mov    -0x30(%rbp),%r9
movslq -0x68(%rbp),%rax
shl    $0x3,%rax
add    %rax,%r9
add    $0xfffffffffffffff8,%r9
mov    -0x30(%rbp),%rax
movslq -0x6c(%rbp),%r10
shl    $0x3,%r10
add    %r10,%rax
add    $0xfffffffffffffff8,%rax
mov    %rax,(%rsp)
call   1050 <dradf4@plt>
jmp    13ee <drftf1+0x2be>
cmpl   $0x2,-0x58(%rbp)
je     12bc <drftf1+0x18c>
jmp    1321 <drftf1+0x1f1>
cmpl   $0x0,-0x4c(%rbp)
je     12cb <drftf1+0x19b>
jmp    12f6 <drftf1+0x1c6>
mov    -0x60(%rbp),%edi
mov    -0x44(%rbp),%esi
mov    -0x20(%rbp),%rdx
mov    -0x28(%rbp),%rcx
mov    -0x30(%rbp),%r8
movslq -0x5c(%rbp),%rax
shl    $0x3,%rax
add    %rax,%r8
add    $0xfffffffffffffff8,%r8
call   1030 <dradf2@plt>
jmp    13ee <drftf1+0x2be>
mov    -0x60(%rbp),%edi
mov    -0x44(%rbp),%esi
mov    -0x28(%rbp),%rdx
mov    -0x20(%rbp),%rcx
mov    -0x30(%rbp),%r8
movslq -0x5c(%rbp),%rax
shl    $0x3,%rax
add    %rax,%r8
add    $0xfffffffffffffff8,%r8
call   1030 <dradf2@plt>
jmp    13ee <drftf1+0x2be>
cmpl   $0x1,-0x60(%rbp)
jne    1336 <drftf1+0x206>
mov    $0x1,%eax
sub    -0x4c(%rbp),%eax
mov    %eax,-0x4c(%rbp)
cmpl   $0x0,-0x4c(%rbp)
je     1345 <drftf1+0x215>
jmp    139c <drftf1+0x26c>
mov    -0x60(%rbp),%edi
mov    -0x58(%rbp),%esi
mov    -0x44(%rbp),%edx
mov    -0x64(%rbp),%ecx
mov    -0x20(%rbp),%r8
mov    -0x20(%rbp),%r9
mov    -0x20(%rbp),%rbx
mov    -0x28(%rbp),%r11
mov    -0x28(%rbp),%r10
mov    -0x30(%rbp),%rax
movslq -0x5c(%rbp),%r14
shl    $0x3,%r14
add    %r14,%rax
add    $0xfffffffffffffff8,%rax
mov    %rbx,(%rsp)
mov    %r11,0x8(%rsp)
mov    %r10,0x10(%rsp)
mov    %rax,0x18(%rsp)
call   1040 <dradfg@plt>
movl   $0x1,-0x4c(%rbp)
jmp    13ee <drftf1+0x2be>
mov    -0x60(%rbp),%edi
mov    -0x58(%rbp),%esi
mov    -0x44(%rbp),%edx
mov    -0x64(%rbp),%ecx
mov    -0x28(%rbp),%r8
mov    -0x28(%rbp),%r9
mov    -0x28(%rbp),%rbx
mov    -0x20(%rbp),%r11
mov    -0x20(%rbp),%r10
mov    -0x30(%rbp),%rax
movslq -0x5c(%rbp),%r14
shl    $0x3,%r14
add    %r14,%rax
add    $0xfffffffffffffff8,%rax
mov    %rbx,(%rsp)
mov    %r11,0x8(%rsp)
mov    %r10,0x10(%rsp)
mov    %rax,0x18(%rsp)
call   1040 <dradfg@plt>
movl   $0x0,-0x4c(%rbp)
mov    -0x44(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1175 <drftf1+0x45>
cmpl   $0x1,-0x4c(%rbp)
jne    1411 <drftf1+0x2e1>
jmp    144c <drftf1+0x31c>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    144c <drftf1+0x31c>
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1418 <drftf1+0x2e8>
add    $0x80,%rsp
pop    %rbx
pop    %r14
pop    %rbp
ret

x86 assembly (clang, -O2):
<drftf1>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %r8,0x20(%rsp)
mov    0x4(%r8),%r15d
test   %r15d,%r15d
jle    14d4 <drftf1+0x3a4>
mov    %rcx,%rbp
mov    %rdx,%r13
mov    %rsi,%r10
mov    %edi,%ebx
add    $0xfffffffffffffff8,%rbp
mov    $0x1,%edi
mov    %ebx,%r14d
mov    %ebx,%r12d
mov    %ebx,0xc(%rsp)
mov    %rbp,0x18(%rsp)
mov    %rsi,0x10(%rsp)
jmp    11a2 <drftf1+0x72>
mov    %eax,%edi
mov    %r14d,%esi
mov    %r10,%rdx
mov    %r13,%rcx
call   1030 <dradf2@plt>
mov    0x10(%rsp),%r10
mov    0x18(%rsp),%rbp
xor    %edi,%edi
mov    $0x1,%al
add    $0xffffffffffffffff,%r15
je     132c <drftf1+0x1fc>
mov    %r14d,%ecx
lea    0x1(%r15),%eax
cltq
mov    0x20(%rsp),%rdx
mov    (%rdx,%rax,4),%esi
mov    %r14d,%eax
cltd
idiv   %esi
mov    %eax,%r14d
mov    %ebx,%eax
cltd
idiv   %ecx
lea    -0x1(%rsi),%ecx
imul   %eax,%ecx
sub    %ecx,%r12d
movslq %r12d,%rdx
cmp    $0x2,%esi
je     1240 <drftf1+0x110>
cmp    $0x4,%esi
jne    1270 <drftf1+0x140>
lea    (%r12,%rax,1),%esi
mov    %eax,%ecx
add    %esi,%ecx
lea    0x0(,%rdx,8),%r8
add    %rbp,%r8
movslq %esi,%rdx
lea    0x0(,%rdx,8),%r9
add    %rbp,%r9
movslq %ecx,%rcx
lea    0x0(,%rcx,8),%rbx
add    %rbp,%rbx
sub    $0x8,%rsp
cmp    $0x1,%edi
jne    12bd <drftf1+0x18d>
mov    %eax,%edi
mov    %r14d,%esi
mov    %r10,%rdx
mov    %r13,%rcx
push   %rbx
mov    0x1c(%rsp),%ebx
call   1050 <dradf4@plt>
mov    0x20(%rsp),%r10
mov    0x28(%rsp),%rbp
add    $0x10,%rsp
jmp    1194 <drftf1+0x64>
nop
lea    0x0(,%rdx,8),%r8
add    %rbp,%r8
cmp    $0x1,%edi
je     117a <drftf1+0x4a>
mov    %eax,%edi
mov    %r14d,%esi
mov    %r13,%rdx
mov    %r10,%rcx
call   1030 <dradf2@plt>
mov    0x10(%rsp),%r10
mov    0x18(%rsp),%rbp
jmp    12e0 <drftf1+0x1b0>
mov    %r14d,%ecx
imul   %eax,%ecx
mov    $0x1,%ebx
sub    %edi,%ebx
cmp    $0x1,%eax
cmove  %edi,%ebx
test   %ebx,%ebx
je     12ec <drftf1+0x1bc>
lea    0x0(,%rdx,8),%rbx
add    %rbp,%rbx
mov    %eax,%edi
mov    %r14d,%edx
mov    %r13,%r8
mov    %r13,%r9
push   %rbx
push   %r10
push   %r10
push   %r13
call   1040 <dradfg@plt>
mov    0x30(%rsp),%r10
mov    0x38(%rsp),%rbp
add    $0x20,%rsp
xor    %edi,%edi
mov    $0x1,%al
jmp    1323 <drftf1+0x1f3>
mov    %eax,%edi
mov    %r14d,%esi
mov    %r13,%rdx
mov    %r10,%rcx
push   %rbx
mov    0x1c(%rsp),%ebx
call   1050 <dradf4@plt>
mov    0x20(%rsp),%r10
mov    0x28(%rsp),%rbp
add    $0x10,%rsp
mov    $0x1,%edi
xor    %eax,%eax
jmp    1198 <drftf1+0x68>
lea    0x0(,%rdx,8),%rbx
add    %rbp,%rbx
mov    %eax,%edi
mov    %r14d,%edx
mov    %r10,%r8
mov    %r10,%r9
push   %rbx
push   %r13
push   %r13
push   %r10
call   1040 <dradfg@plt>
mov    0x30(%rsp),%r10
mov    0x38(%rsp),%rbp
add    $0x20,%rsp
mov    $0x1,%edi
xor    %eax,%eax
mov    0xc(%rsp),%ebx
jmp    1198 <drftf1+0x68>
test   %ebx,%ebx
setle  %cl
xor    $0x1,%al
or     %cl,%al
jne    14d4 <drftf1+0x3a4>
mov    %ebx,%eax
cmp    $0x4,%ebx
jb     1363 <drftf1+0x233>
lea    0x0(,%rax,8),%rcx
add    %r13,%rcx
cmp    %r10,%rcx
jbe    13f5 <drftf1+0x2c5>
lea    (%r10,%rax,8),%rcx
cmp    %r13,%rcx
jbe    13f5 <drftf1+0x2c5>
xor    %ecx,%ecx
mov    %rcx,%rdx
not    %rdx
add    %rax,%rdx
mov    %rax,%rsi
and    $0x3,%rsi
je     1397 <drftf1+0x267>
nopw   0x0(%rax,%rax,1)
movsd  0x0(%r13,%rcx,8),%xmm0
movsd  %xmm0,(%r10,%rcx,8)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rsi
jne    1380 <drftf1+0x250>
cmp    $0x3,%rdx
jb     14d4 <drftf1+0x3a4>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  0x0(%r13,%rcx,8),%xmm0
movsd  %xmm0,(%r10,%rcx,8)
movsd  0x8(%r13,%rcx,8),%xmm0
movsd  %xmm0,0x8(%r10,%rcx,8)
movsd  0x10(%r13,%rcx,8),%xmm0
movsd  %xmm0,0x10(%r10,%rcx,8)
movsd  0x18(%r13,%rcx,8),%xmm0
movsd  %xmm0,0x18(%r10,%rcx,8)
add    $0x4,%rcx
cmp    %rcx,%rax
jne    13b0 <drftf1+0x280>
jmp    14d4 <drftf1+0x3a4>
mov    %eax,%ecx
and    $0xfffffffc,%ecx
lea    -0x4(%rcx),%rdi
mov    %rdi,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
mov    %esi,%edx
and    $0x3,%edx
cmp    $0xc,%rdi
jae    1418 <drftf1+0x2e8>
xor    %edi,%edi
jmp    1489 <drftf1+0x359>
and    $0xfffffffffffffffc,%rsi
xor    %edi,%edi
xchg   %ax,%ax
movups 0x0(%r13,%rdi,8),%xmm0
movups 0x10(%r13,%rdi,8),%xmm1
movups %xmm0,(%r10,%rdi,8)
movups %xmm1,0x10(%r10,%rdi,8)
movups 0x20(%r13,%rdi,8),%xmm0
movups 0x30(%r13,%rdi,8),%xmm1
movups %xmm0,0x20(%r10,%rdi,8)
movups %xmm1,0x30(%r10,%rdi,8)
movups 0x40(%r13,%rdi,8),%xmm0
movups 0x50(%r13,%rdi,8),%xmm1
movups %xmm0,0x40(%r10,%rdi,8)
movups %xmm1,0x50(%r10,%rdi,8)
movups 0x60(%r13,%rdi,8),%xmm0
movups 0x70(%r13,%rdi,8),%xmm1
movups %xmm0,0x60(%r10,%rdi,8)
movups %xmm1,0x70(%r10,%rdi,8)
add    $0x10,%rdi
add    $0xfffffffffffffffc,%rsi
jne    1420 <drftf1+0x2f0>
test   %rdx,%rdx
je     14cb <drftf1+0x39b>
lea    (%r10,%rdi,8),%rsi
add    $0x10,%rsi
lea    0x10(,%rdi,8),%rdi
add    %r13,%rdi
shl    $0x5,%rdx
xor    %ebp,%ebp
nopw   0x0(%rax,%rax,1)
movups -0x10(%rdi,%rbp,1),%xmm0
movups (%rdi,%rbp,1),%xmm1
movups %xmm0,-0x10(%rsi,%rbp,1)
movups %xmm1,(%rsi,%rbp,1)
add    $0x20,%rbp
cmp    %rbp,%rdx
jne    14b0 <drftf1+0x380>
cmp    %rax,%rcx
jne    1365 <drftf1+0x235>
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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