#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 */
typedef size_t k;
typedef int ic;
typedef size_t const i;

/* Variables and functions */
void radb2(size_t ido, size_t l1, const double *cc, double *ch,
           const double *wa) {
  const size_t cdim = 2;
  size_t i, k, ic;
  double ti2, tr2;
  for (k = 0; k < l1; k++) {
    ch[(0) + ido * ((k) + l1 * (0))] = cc[(0) + ido * ((0) + cdim * (k))] +
                                       cc[(ido - 1) + ido * ((1) + cdim * (k))];
    ch[(0) + ido * ((k) + l1 * (1))] = cc[(0) + ido * ((0) + cdim * (k))] -
                                       cc[(ido - 1) + ido * ((1) + cdim * (k))];
  }
  if ((ido & 1) == 0)
    for (k = 0; k < l1; k++) {
      ch[(ido - 1) + ido * ((k) + l1 * (0))] =
          2 * cc[(ido - 1) + ido * ((0) + cdim * (k))];
      ch[(ido - 1) + ido * ((k) + l1 * (1))] =
          -2 * cc[(0) + ido * ((1) + cdim * (k))];
    }
  if (ido <= 2)
    return;
  for (k = 0; k < l1; ++k)
    for (i = 2; i < ido; i += 2) {
      ic = ido - i;
      {
        ch[(i - 1) + ido * ((k) + l1 * (0))] =
            cc[(i - 1) + ido * ((0) + cdim * (k))] +
            cc[(ic - 1) + ido * ((1) + cdim * (k))];
        tr2 = cc[(i - 1) + ido * ((0) + cdim * (k))] -
              cc[(ic - 1) + ido * ((1) + cdim * (k))];
      }
      {
        ti2 = cc[(i) + ido * ((0) + cdim * (k))] +
              cc[(ic) + ido * ((1) + cdim * (k))];
        ch[(i) + ido * ((k) + l1 * (0))] = cc[(i) + ido * ((0) + cdim * (k))] -
                                           cc[(ic) + ido * ((1) + cdim * (k))];
      }
      {
        ch[(i) + ido * ((k) + l1 * (1))] =
            wa[(i - 2) + (0) * ido] * ti2 + wa[(i - 1) + (0) * ido] * tr2;
        ch[(i - 1) + ido * ((k) + l1 * (1))] =
            wa[(i - 2) + (0) * ido] * tr2 - wa[(i - 1) + (0) * ido] * ti2;
      }
    }
}


x86 assembly (clang, -O0):
<radb2>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
movq   $0x2,-0x30(%rbp)
movq   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%rax
cmp    -0x10(%rbp),%rax
jae    121e <radb2+0x11e>
mov    -0x18(%rbp),%rax
mov    -0x8(%rbp),%rcx
mov    -0x40(%rbp),%rdx
shl    $0x1,%rdx
add    $0x0,%rdx
imul   %rdx,%rcx
add    $0x0,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x8(%rbp),%rcx
sub    $0x1,%rcx
mov    -0x8(%rbp),%rdx
mov    -0x40(%rbp),%rsi
shl    $0x1,%rsi
add    $0x1,%rsi
imul   %rsi,%rdx
add    %rdx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x8(%rbp),%rcx
mov    -0x40(%rbp),%rdx
imul   $0x0,-0x10(%rbp),%rsi
add    %rsi,%rdx
imul   %rdx,%rcx
add    $0x0,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
mov    -0x8(%rbp),%rcx
mov    -0x40(%rbp),%rdx
shl    $0x1,%rdx
add    $0x0,%rdx
imul   %rdx,%rcx
add    $0x0,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x8(%rbp),%rcx
sub    $0x1,%rcx
mov    -0x8(%rbp),%rdx
mov    -0x40(%rbp),%rsi
shl    $0x1,%rsi
add    $0x1,%rsi
imul   %rsi,%rdx
add    %rdx,%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x8(%rbp),%rcx
mov    -0x40(%rbp),%rdx
mov    -0x10(%rbp),%rsi
shl    $0x0,%rsi
add    %rsi,%rdx
imul   %rdx,%rcx
add    $0x0,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    1128 <radb2+0x28>
mov    -0x8(%rbp),%rax
and    $0x1,%rax
cmp    $0x0,%rax
jne    1308 <radb2+0x208>
movq   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%rax
cmp    -0x10(%rbp),%rax
jae    1303 <radb2+0x203>
mov    -0x18(%rbp),%rax
mov    -0x8(%rbp),%rcx
sub    $0x1,%rcx
mov    -0x8(%rbp),%rdx
mov    -0x40(%rbp),%rsi
shl    $0x1,%rsi
add    $0x0,%rsi
imul   %rsi,%rdx
add    %rdx,%rcx
movsd  0xd97(%rip),%xmm0
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x8(%rbp),%rcx
sub    $0x1,%rcx
mov    -0x8(%rbp),%rdx
mov    -0x40(%rbp),%rsi
imul   $0x0,-0x10(%rbp),%rdi
add    %rdi,%rsi
imul   %rsi,%rdx
add    %rdx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
mov    -0x8(%rbp),%rcx
mov    -0x40(%rbp),%rdx
shl    $0x1,%rdx
add    $0x1,%rdx
imul   %rdx,%rcx
add    $0x0,%rcx
movsd  0xd3e(%rip),%xmm0
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x8(%rbp),%rcx
sub    $0x1,%rcx
mov    -0x8(%rbp),%rdx
mov    -0x40(%rbp),%rsi
mov    -0x10(%rbp),%rdi
shl    $0x0,%rdi
add    %rdi,%rsi
imul   %rsi,%rdx
add    %rdx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    1238 <radb2+0x138>
jmp    1308 <radb2+0x208>
cmpq   $0x2,-0x8(%rbp)
ja     1318 <radb2+0x218>
jmp    15ee <radb2+0x4ee>
movq   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%rax
cmp    -0x10(%rbp),%rax
jae    15ee <radb2+0x4ee>
movq   $0x2,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x8(%rbp),%rax
jae    15d8 <radb2+0x4d8>
mov    -0x8(%rbp),%rax
sub    -0x38(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%rcx
sub    $0x1,%rcx
mov    -0x8(%rbp),%rdx
mov    -0x40(%rbp),%rsi
shl    $0x1,%rsi
add    $0x0,%rsi
imul   %rsi,%rdx
add    %rdx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x48(%rbp),%rcx
sub    $0x1,%rcx
mov    -0x8(%rbp),%rdx
mov    -0x40(%rbp),%rsi
shl    $0x1,%rsi
add    $0x1,%rsi
imul   %rsi,%rdx
add    %rdx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%rcx
sub    $0x1,%rcx
mov    -0x8(%rbp),%rdx
mov    -0x40(%rbp),%rsi
imul   $0x0,-0x10(%rbp),%rdi
add    %rdi,%rsi
imul   %rsi,%rdx
add    %rdx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%rcx
sub    $0x1,%rcx
mov    -0x8(%rbp),%rdx
mov    -0x40(%rbp),%rsi
shl    $0x1,%rsi
add    $0x0,%rsi
imul   %rsi,%rdx
add    %rdx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x48(%rbp),%rcx
sub    $0x1,%rcx
mov    -0x8(%rbp),%rdx
mov    -0x40(%rbp),%rsi
shl    $0x1,%rsi
add    $0x1,%rsi
imul   %rsi,%rdx
add    %rdx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    -0x8(%rbp),%rdx
mov    -0x40(%rbp),%rsi
shl    $0x1,%rsi
add    $0x0,%rsi
imul   %rsi,%rdx
add    %rdx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x48(%rbp),%rcx
mov    -0x8(%rbp),%rdx
mov    -0x40(%rbp),%rsi
shl    $0x1,%rsi
add    $0x1,%rsi
imul   %rsi,%rdx
add    %rdx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    -0x8(%rbp),%rdx
mov    -0x40(%rbp),%rsi
shl    $0x1,%rsi
add    $0x0,%rsi
imul   %rsi,%rdx
add    %rdx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x48(%rbp),%rcx
mov    -0x8(%rbp),%rdx
mov    -0x40(%rbp),%rsi
shl    $0x1,%rsi
add    $0x1,%rsi
imul   %rsi,%rdx
add    %rdx,%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    -0x8(%rbp),%rdx
mov    -0x40(%rbp),%rsi
imul   $0x0,-0x10(%rbp),%rdi
add    %rdi,%rsi
imul   %rsi,%rdx
add    %rdx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
sub    $0x2,%rcx
imul   $0x0,-0x8(%rbp),%rdx
add    %rdx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x50(%rbp),%xmm2
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
sub    $0x1,%rcx
imul   $0x0,-0x8(%rbp),%rdx
add    %rdx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    -0x8(%rbp),%rdx
mov    -0x40(%rbp),%rsi
mov    -0x10(%rbp),%rdi
shl    $0x0,%rdi
add    %rdi,%rsi
imul   %rsi,%rdx
add    %rdx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
sub    $0x2,%rcx
imul   $0x0,-0x8(%rbp),%rdx
add    %rdx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x58(%rbp),%xmm2
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
sub    $0x1,%rcx
imul   $0x0,-0x8(%rbp),%rdx
add    %rdx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0x50(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%rcx
sub    $0x1,%rcx
mov    -0x8(%rbp),%rdx
mov    -0x40(%rbp),%rsi
mov    -0x10(%rbp),%rdi
shl    $0x0,%rdi
add    %rdi,%rsi
imul   %rsi,%rdx
add    %rdx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%rax
add    $0x2,%rax
mov    %rax,-0x38(%rbp)
jmp    1336 <radb2+0x236>
jmp    15dd <radb2+0x4dd>
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    1320 <radb2+0x220>
pop    %rbp
ret

x86 assembly (clang, -O2):
<radb2>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
test   %rsi,%rsi
je     13ac <radb2+0x2ac>
xor    %r9d,%r9d
cmp    $0x1,%rsi
je     11cf <radb2+0xcf>
cmp    $0x1,%rdi
jne    11cf <radb2+0xcf>
lea    (%rcx,%rsi,8),%rbx
mov    %rsi,%rax
shl    $0x4,%rax
lea    (%rcx,%rax,1),%r9
add    %rdx,%rax
cmp    %rax,%rbx
setb   %r10b
cmp    %rdx,%r9
seta   %r11b
cmp    %rcx,%rax
seta   %al
cmp    %rdx,%rbx
seta   %bl
xor    %r9d,%r9d
test   %r11b,%r10b
jne    11cf <radb2+0xcf>
and    %bl,%al
jne    11cf <radb2+0xcf>
mov    %rsi,%r9
and    $0xfffffffffffffffe,%r9
mov    %rdi,%r11
shl    $0x4,%r11
mov    %rsi,%rax
imul   %rdi,%rax
lea    (%rcx,%rax,8),%r10
xor    %eax,%eax
mov    %r9,%rbx
xchg   %ax,%ax
movupd (%rdx,%rax,2),%xmm0
movupd 0x10(%rdx,%rax,2),%xmm1
movapd %xmm0,%xmm2
unpcklpd %xmm1,%xmm2
unpckhpd %xmm1,%xmm0
addpd  %xmm2,%xmm0
movupd %xmm0,(%rcx,%rax,1)
movupd (%rdx,%rax,2),%xmm0
movupd 0x10(%rdx,%rax,2),%xmm1
movapd %xmm0,%xmm2
unpcklpd %xmm1,%xmm2
unpckhpd %xmm1,%xmm0
subpd  %xmm0,%xmm2
movupd %xmm2,(%r10,%rax,1)
add    %r11,%rax
add    $0xfffffffffffffffe,%rbx
jne    1180 <radb2+0x80>
cmp    %rsi,%r9
je     124d <radb2+0x14d>
mov    %rsi,%rbx
sub    %r9,%rbx
lea    (%r9,%rsi,1),%rax
imul   %rdi,%rax
lea    (%rcx,%rax,8),%r10
lea    0x0(,%rdi,8),%r11
mov    %r9,%r15
imul   %rdi,%r15
lea    (%rcx,%r15,8),%r14
lea    (%r9,%r9,1),%rax
add    $0x2,%rax
imul   %rdi,%rax
lea    (%rdx,%rax,8),%r9
add    $0xfffffffffffffff8,%r9
shl    $0x4,%r15
add    %rdx,%r15
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  (%r15,%rax,2),%xmm0
addsd  (%r9,%rax,2),%xmm0
movsd  %xmm0,(%r14,%rax,1)
movsd  (%r15,%rax,2),%xmm0
subsd  (%r9,%rax,2),%xmm0
movsd  %xmm0,(%r10,%rax,1)
add    %r11,%rax
add    $0xffffffffffffffff,%rbx
jne    1220 <radb2+0x120>
mov    $0x1,%r9b
test   $0x1,%dil
jne    13af <radb2+0x2af>
xor    %r10d,%r10d
cmp    $0x2,%rsi
jb     132f <radb2+0x22f>
cmp    $0x1,%rdi
jne    132f <radb2+0x22f>
lea    (%rcx,%rsi,8),%rbx
mov    %rsi,%rax
shl    $0x4,%rax
lea    (%rcx,%rax,1),%r10
add    %rdx,%rax
cmp    %rax,%rbx
setb   %r11b
cmp    %rdx,%r10
seta   %bpl
cmp    %rcx,%rax
seta   %al
cmp    %rdx,%rbx
seta   %bl
xor    %r10d,%r10d
test   %bpl,%r11b
jne    132f <radb2+0x22f>
and    %bl,%al
jne    132f <radb2+0x22f>
lea    -0x1(%rdi),%rax
mov    %rsi,%r10
and    $0xfffffffffffffffe,%r10
lea    (%rdx,%rax,8),%r12
mov    %rdi,%r14
shl    $0x4,%r14
lea    (%rcx,%rax,8),%r11
lea    0x1(%rsi),%rax
imul   %rdi,%rax
lea    (%rcx,%rax,8),%r15
add    $0xfffffffffffffff8,%r15
xor    %eax,%eax
movapd 0xd1b(%rip),%xmm0
mov    %r10,%rbx
nopl   0x0(%rax,%rax,1)
movupd (%r12,%rax,2),%xmm1
movupd 0x10(%r12,%rax,2),%xmm2
movapd %xmm1,%xmm3
unpcklpd %xmm2,%xmm3
unpckhpd %xmm2,%xmm1
addpd  %xmm3,%xmm3
movupd %xmm3,(%r11,%rax,1)
mulpd  %xmm0,%xmm1
movupd %xmm1,(%r15,%rax,1)
add    %r14,%rax
add    $0xfffffffffffffffe,%rbx
jne    12f0 <radb2+0x1f0>
cmp    %rsi,%r10
je     13af <radb2+0x2af>
mov    %rsi,%rbx
sub    %r10,%rbx
lea    (%r10,%rsi,1),%rax
add    $0x1,%rax
imul   %rdi,%rax
lea    (%rcx,%rax,8),%r11
add    $0xfffffffffffffff8,%r11
lea    0x0(,%rdi,8),%r14
lea    0x1(%r10),%rax
imul   %rdi,%rax
lea    (%rcx,%rax,8),%r15
add    $0xfffffffffffffff8,%r15
lea    (%r10,%r10,1),%rax
add    $0x1,%rax
imul   %rdi,%rax
lea    (%rdx,%rax,8),%r10
xor    %eax,%eax
movsd  0xc95(%rip),%xmm0
nopl   0x0(%rax,%rax,1)
movsd  -0x8(%r10,%rax,2),%xmm1
addsd  %xmm1,%xmm1
movsd  %xmm1,(%r15,%rax,1)
movsd  (%r10,%rax,2),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%r11,%rax,1)
add    %r14,%rax
add    $0xffffffffffffffff,%rbx
jne    1380 <radb2+0x280>
jmp    13af <radb2+0x2af>
xor    %r9d,%r9d
cmp    $0x3,%rdi
jb     1494 <radb2+0x394>
test   %r9b,%r9b
je     1494 <radb2+0x394>
mov    %rsi,%rax
imul   %rdi,%rax
lea    (%rcx,%rax,8),%r15
lea    0x0(,%rdi,8),%r9
mov    %rdi,%r10
shl    $0x4,%r10
lea    (%rdi,%rdi,1),%r11
add    $0xfffffffffffffffd,%r11
xor    %r14d,%r14d
nopw   0x0(%rax,%rax,1)
mov    $0x2,%ebx
mov    %r11,%rax
nopl   0x0(%rax,%rax,1)
movsd  -0x8(%rdx,%rbx,8),%xmm0
addsd  (%rdx,%rax,8),%xmm0
movsd  %xmm0,-0x8(%rcx,%rbx,8)
movsd  -0x8(%rdx,%rbx,8),%xmm0
movsd  (%rdx,%rbx,8),%xmm1
subsd  (%rdx,%rax,8),%xmm0
movsd  0x8(%rdx,%rax,8),%xmm2
movapd %xmm1,%xmm3
addsd  %xmm2,%xmm3
subsd  %xmm2,%xmm1
movsd  %xmm1,(%rcx,%rbx,8)
movsd  -0x8(%r8,%rbx,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  -0x10(%r8,%rbx,8),%xmm2
mulsd  %xmm3,%xmm2
addsd  %xmm1,%xmm2
movsd  %xmm2,(%r15,%rbx,8)
mulsd  -0x8(%r8,%rbx,8),%xmm3
mulsd  -0x10(%r8,%rbx,8),%xmm0
subsd  %xmm3,%xmm0
movsd  %xmm0,-0x8(%r15,%rbx,8)
add    $0x2,%rbx
add    $0xfffffffffffffffe,%rax
cmp    %rdi,%rbx
jb     1400 <radb2+0x300>
add    $0x1,%r14
add    %r9,%r15
add    %r9,%rcx
add    %r10,%rdx
cmp    %rsi,%r14
jne    13f0 <radb2+0x2f0>
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @radb2(i64 %0, i64 %1, double* %2, double* %3, double* %4) {
%6 = alloca i64
%7 = alloca i64
%8 = alloca double*
%9 = alloca double*
%10 = alloca double*
%11 = alloca i64
%12 = alloca i64
%13 = alloca i64
%14 = alloca i64
%15 = alloca double
%16 = alloca double
store i64 %0, i64* %6
store i64 %1, i64* %7
store double* %2, double** %8
store double* %3, double** %9
store double* %4, double** %10
store i64 2, i64* %11
store i64 0, i64* %13
br label %17
17: ; preds = %82, %5
%18 = load i64, i64* %13
%19 = load i64, i64* %7
%20 = icmp ult i64 %18, %19
br i1 %20, label %21, label %85
21: ; preds = %17
%22 = load double*, double** %8
%23 = load i64, i64* %6
%24 = load i64, i64* %13
%25 = mul i64 2, %24
%26 = add i64 0, %25
%27 = mul i64 %23, %26
%28 = add i64 0, %27
%29 = getelementptr inbounds double, double* %22, i64 %28
%30 = load double, double* %29
%31 = load double*, double** %8
%32 = load i64, i64* %6
%33 = sub i64 %32, 1
%34 = load i64, i64* %6
%35 = load i64, i64* %13
%36 = mul i64 2, %35
%37 = add i64 1, %36
%38 = mul i64 %34, %37
%39 = add i64 %33, %38
%40 = getelementptr inbounds double, double* %31, i64 %39
%41 = load double, double* %40
%42 = fadd double %30, %41
%43 = load double*, double** %9
%44 = load i64, i64* %6
%45 = load i64, i64* %13
%46 = load i64, i64* %7
%47 = mul i64 %46, 0
%48 = add i64 %45, %47
%49 = mul i64 %44, %48
%50 = add i64 0, %49
%51 = getelementptr inbounds double, double* %43, i64 %50
store double %42, double* %51
%52 = load double*, double** %8
%53 = load i64, i64* %6
%54 = load i64, i64* %13
%55 = mul i64 2, %54
%56 = add i64 0, %55
%57 = mul i64 %53, %56
%58 = add i64 0, %57
%59 = getelementptr inbounds double, double* %52, i64 %58
%60 = load double, double* %59
%61 = load double*, double** %8
%62 = load i64, i64* %6
%63 = sub i64 %62, 1
%64 = load i64, i64* %6
%65 = load i64, i64* %13
%66 = mul i64 2, %65
%67 = add i64 1, %66
%68 = mul i64 %64, %67
%69 = add i64 %63, %68
%70 = getelementptr inbounds double, double* %61, i64 %69
%71 = load double, double* %70
%72 = fsub double %60, %71
%73 = load double*, double** %9
%74 = load i64, i64* %6
%75 = load i64, i64* %13
%76 = load i64, i64* %7
%77 = mul i64 %76, 1
%78 = add i64 %75, %77
%79 = mul i64 %74, %78
%80 = add i64 0, %79
%81 = getelementptr inbounds double, double* %73, i64 %80
store double %72, double* %81
br label %82
82: ; preds = %21
%83 = load i64, i64* %13
%84 = add i64 %83, 1
store i64 %84, i64* %13
br label %17
85: ; preds = %17
%86 = load i64, i64* %6
%87 = and i64 %86, 1
%88 = icmp eq i64 %87, 0
br i1 %88, label %89, label %143
89: ; preds = %85
store i64 0, i64* %13
br label %90
90: ; preds = %139, %89
%91 = load i64, i64* %13
%92 = load i64, i64* %7
%93 = icmp ult i64 %91, %92
br i1 %93, label %94, label %142
94: ; preds = %90
%95 = load double*, double** %8
%96 = load i64, i64* %6
%97 = sub i64 %96, 1
%98 = load i64, i64* %6
%99 = load i64, i64* %13
%100 = mul i64 2, %99
%101 = add i64 0, %100
%102 = mul i64 %98, %101
%103 = add i64 %97, %102
%104 = getelementptr inbounds double, double* %95, i64 %103
%105 = load double, double* %104
%106 = fmul double 2.000000e+00, %105
%107 = load double*, double** %9
%108 = load i64, i64* %6
%109 = sub i64 %108, 1
%110 = load i64, i64* %6
%111 = load i64, i64* %13
%112 = load i64, i64* %7
%113 = mul i64 %112, 0
%114 = add i64 %111, %113
%115 = mul i64 %110, %114
%116 = add i64 %109, %115
%117 = getelementptr inbounds double, double* %107, i64 %116
store double %106, double* %117
%118 = load double*, double** %8
%119 = load i64, i64* %6
%120 = load i64, i64* %13
%121 = mul i64 2, %120
%122 = add i64 1, %121
%123 = mul i64 %119, %122
%124 = add i64 0, %123
%125 = getelementptr inbounds double, double* %118, i64 %124
%126 = load double, double* %125
%127 = fmul double -2.000000e+00, %126
%128 = load double*, double** %9
%129 = load i64, i64* %6
%130 = sub i64 %129, 1
%131 = load i64, i64* %6
%132 = load i64, i64* %13
%133 = load i64, i64* %7
%134 = mul i64 %133, 1
%135 = add i64 %132, %134
%136 = mul i64 %131, %135
%137 = add i64 %130, %136
%138 = getelementptr inbounds double, double* %128, i64 %137
store double %127, double* %138
br label %139
139: ; preds = %94
%140 = load i64, i64* %13
%141 = add i64 %140, 1
store i64 %141, i64* %13
br label %90
142: ; preds = %90
br label %143
143: ; preds = %142, %85
%144 = load i64, i64* %6
%145 = icmp ule i64 %144, 2
br i1 %145, label %146, label %147
146: ; preds = %143
br label %339
147: ; preds = %143
store i64 0, i64* %13
br label %148
148: ; preds = %336, %147
%149 = load i64, i64* %13
%150 = load i64, i64* %7
%151 = icmp ult i64 %149, %150
br i1 %151, label %152, label %339
152: ; preds = %148
store i64 2, i64* %12
br label %153
153: ; preds = %332, %152
%154 = load i64, i64* %12
%155 = load i64, i64* %6
%156 = icmp ult i64 %154, %155
br i1 %156, label %157, label %335
157: ; preds = %153
%158 = load i64, i64* %6
%159 = load i64, i64* %12
%160 = sub i64 %158, %159
store i64 %160, i64* %14
%161 = load double*, double** %8
%162 = load i64, i64* %12
%163 = sub i64 %162, 1
%164 = load i64, i64* %6
%165 = load i64, i64* %13
%166 = mul i64 2, %165
%167 = add i64 0, %166
%168 = mul i64 %164, %167
%169 = add i64 %163, %168
%170 = getelementptr inbounds double, double* %161, i64 %169
%171 = load double, double* %170
%172 = load double*, double** %8
%173 = load i64, i64* %14
%174 = sub i64 %173, 1
%175 = load i64, i64* %6
%176 = load i64, i64* %13
%177 = mul i64 2, %176
%178 = add i64 1, %177
%179 = mul i64 %175, %178
%180 = add i64 %174, %179
%181 = getelementptr inbounds double, double* %172, i64 %180
%182 = load double, double* %181
%183 = fadd double %171, %182
%184 = load double*, double** %9
%185 = load i64, i64* %12
%186 = sub i64 %185, 1
%187 = load i64, i64* %6
%188 = load i64, i64* %13
%189 = load i64, i64* %7
%190 = mul i64 %189, 0
%191 = add i64 %188, %190
%192 = mul i64 %187, %191
%193 = add i64 %186, %192
%194 = getelementptr inbounds double, double* %184, i64 %193
store double %183, double* %194
%195 = load double*, double** %8
%196 = load i64, i64* %12
%197 = sub i64 %196, 1
%198 = load i64, i64* %6
%199 = load i64, i64* %13
%200 = mul i64 2, %199
%201 = add i64 0, %200
%202 = mul i64 %198, %201
%203 = add i64 %197, %202
%204 = getelementptr inbounds double, double* %195, i64 %203
%205 = load double, double* %204
%206 = load double*, double** %8
%207 = load i64, i64* %14
%208 = sub i64 %207, 1
%209 = load i64, i64* %6
%210 = load i64, i64* %13
%211 = mul i64 2, %210
%212 = add i64 1, %211
%213 = mul i64 %209, %212
%214 = add i64 %208, %213
%215 = getelementptr inbounds double, double* %206, i64 %214
%216 = load double, double* %215
%217 = fsub double %205, %216
store double %217, double* %16
%218 = load double*, double** %8
%219 = load i64, i64* %12
%220 = load i64, i64* %6
%221 = load i64, i64* %13
%222 = mul i64 2, %221
%223 = add i64 0, %222
%224 = mul i64 %220, %223
%225 = add i64 %219, %224
%226 = getelementptr inbounds double, double* %218, i64 %225
%227 = load double, double* %226
%228 = load double*, double** %8
%229 = load i64, i64* %14
%230 = load i64, i64* %6
%231 = load i64, i64* %13
%232 = mul i64 2, %231
%233 = add i64 1, %232
%234 = mul i64 %230, %233
%235 = add i64 %229, %234
%236 = getelementptr inbounds double, double* %228, i64 %235
%237 = load double, double* %236
%238 = fadd double %227, %237
store double %238, double* %15
%239 = load double*, double** %8
%240 = load i64, i64* %12
%241 = load i64, i64* %6
%242 = load i64, i64* %13
%243 = mul i64 2, %242
%244 = add i64 0, %243
%245 = mul i64 %241, %244
%246 = add i64 %240, %245
%247 = getelementptr inbounds double, double* %239, i64 %246
%248 = load double, double* %247
%249 = load double*, double** %8
%250 = load i64, i64* %14
%251 = load i64, i64* %6
%252 = load i64, i64* %13
%253 = mul i64 2, %252
%254 = add i64 1, %253
%255 = mul i64 %251, %254
%256 = add i64 %250, %255
%257 = getelementptr inbounds double, double* %249, i64 %256
%258 = load double, double* %257
%259 = fsub double %248, %258
%260 = load double*, double** %9
%261 = load i64, i64* %12
%262 = load i64, i64* %6
%263 = load i64, i64* %13
%264 = load i64, i64* %7
%265 = mul i64 %264, 0
%266 = add i64 %263, %265
%267 = mul i64 %262, %266
%268 = add i64 %261, %267
%269 = getelementptr inbounds double, double* %260, i64 %268
store double %259, double* %269
%270 = load double*, double** %10
%271 = load i64, i64* %12
%272 = sub i64 %271, 2
%273 = load i64, i64* %6
%274 = mul i64 0, %273
%275 = add i64 %272, %274
%276 = getelementptr inbounds double, double* %270, i64 %275
%277 = load double, double* %276
%278 = load double, double* %15
%279 = load double*, double** %10
%280 = load i64, i64* %12
%281 = sub i64 %280, 1
%282 = load i64, i64* %6
%283 = mul i64 0, %282
%284 = add i64 %281, %283
%285 = getelementptr inbounds double, double* %279, i64 %284
%286 = load double, double* %285
%287 = load double, double* %16
%288 = fmul double %286, %287
%289 = call double @llvm.fmuladd.f64(double %277, double %278, double %288)
%290 = load double*, double** %9
%291 = load i64, i64* %12
%292 = load i64, i64* %6
%293 = load i64, i64* %13
%294 = load i64, i64* %7
%295 = mul i64 %294, 1
%296 = add i64 %293, %295
%297 = mul i64 %292, %296
%298 = add i64 %291, %297
%299 = getelementptr inbounds double, double* %290, i64 %298
store double %289, double* %299
%300 = load double*, double** %10
%301 = load i64, i64* %12
%302 = sub i64 %301, 2
%303 = load i64, i64* %6
%304 = mul i64 0, %303
%305 = add i64 %302, %304
%306 = getelementptr inbounds double, double* %300, i64 %305
%307 = load double, double* %306
%308 = load double, double* %16
%309 = load double*, double** %10
%310 = load i64, i64* %12
%311 = sub i64 %310, 1
%312 = load i64, i64* %6
%313 = mul i64 0, %312
%314 = add i64 %311, %313
%315 = getelementptr inbounds double, double* %309, i64 %314
%316 = load double, double* %315
%317 = load double, double* %15
%318 = fmul double %316, %317
%319 = fneg double %318
%320 = call double @llvm.fmuladd.f64(double %307, double %308, double %319)
%321 = load double*, double** %9
%322 = load i64, i64* %12
%323 = sub i64 %322, 1
%324 = load i64, i64* %6
%325 = load i64, i64* %13
%326 = load i64, i64* %7
%327 = mul i64 %326, 1
%328 = add i64 %325, %327
%329 = mul i64 %324, %328
%330 = add i64 %323, %329
%331 = getelementptr inbounds double, double* %321, i64 %330
store double %320, double* %331
br label %332
332: ; preds = %157
%333 = load i64, i64* %12
%334 = add i64 %333, 2
store i64 %334, i64* %12
br label %153
335: ; preds = %153
br label %336
336: ; preds = %335
%337 = load i64, i64* %13
%338 = add i64 %337, 1
store i64 %338, i64* %13
br label %148
339: ; preds = %146, %148
ret void
}

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