void RK4(int N, double t, double h, const double X[],
         void (*F)(double t, const double X[], double R[]), double R[]) {
  double K1[N];
  F(t, X, K1);
  double T[N], K2[N];
  for (int i = 0; i < N; i++)
    T[i] = X[i] + h / 2 * K1[i];
  F(t + h / 2, T, K2);
  double K3[N];
  for (int i = 0; i < N; i++)
    T[i] = X[i] + h / 2 * K2[i];
  F(t + h / 2, T, K3);
  double K4[N];
  for (int i = 0; i < N; i++)
    T[i] = X[i] + h * K3[i];
  F(t + h, T, K4);
  for (int i = 0; i < N; i++)
    R[i] = X[i] + h / 6 * (K1[i] + 2 * (K2[i] + K3[i]) + K4[i]);
}


x86 assembly (clang, -O0):
<RK4>:
push   %rbp
mov    %rsp,%rbp
sub    $0xa0,%rsp
mov    %edi,-0x4(%rbp)
movsd  %xmm0,-0x10(%rbp)
movsd  %xmm1,-0x18(%rbp)
mov    %rsi,-0x20(%rbp)
mov    %rdx,-0x28(%rbp)
mov    %rcx,-0x30(%rbp)
mov    -0x4(%rbp),%eax
mov    %rsp,%rcx
mov    %rcx,-0x38(%rbp)
lea    0xf(,%rax,8),%rcx
and    $0xfffffffffffffff0,%rcx
mov    %rsp,%rsi
sub    %rcx,%rsi
mov    %rsi,-0x90(%rbp)
mov    %rsi,%rsp
mov    %rax,-0x40(%rbp)
mov    -0x28(%rbp),%rax
movsd  -0x10(%rbp),%xmm0
mov    -0x20(%rbp),%rdi
call   *%rax
mov    -0x4(%rbp),%eax
lea    0xf(,%rax,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0x88(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x48(%rbp)
mov    -0x4(%rbp),%eax
lea    0xf(,%rax,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0x80(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x50(%rbp)
movl   $0x0,-0x54(%rbp)
mov    -0x54(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1207 <RK4+0x107>
mov    -0x88(%rbp),%rax
mov    -0x90(%rbp),%rcx
mov    -0x20(%rbp),%rdx
movslq -0x54(%rbp),%rsi
movsd  (%rdx,%rsi,8),%xmm1
movsd  -0x18(%rbp),%xmm0
movsd  0xe25(%rip),%xmm2
divsd  %xmm2,%xmm0
movslq -0x54(%rbp),%rdx
movsd  (%rcx,%rdx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movslq -0x54(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    11a7 <RK4+0xa7>
mov    -0x80(%rbp),%rsi
mov    -0x88(%rbp),%rdi
mov    -0x28(%rbp),%rax
movsd  -0x10(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm1
movsd  0xdd8(%rip),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
call   *%rax
mov    -0x4(%rbp),%eax
lea    0xf(,%rax,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0x98(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x60(%rbp)
movl   $0x0,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    12b9 <RK4+0x1b9>
mov    -0x88(%rbp),%rax
mov    -0x80(%rbp),%rcx
mov    -0x20(%rbp),%rdx
movslq -0x64(%rbp),%rsi
movsd  (%rdx,%rsi,8),%xmm1
movsd  -0x18(%rbp),%xmm0
movsd  0xd73(%rip),%xmm2
divsd  %xmm2,%xmm0
movslq -0x64(%rbp),%rdx
movsd  (%rcx,%rdx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movslq -0x64(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    125c <RK4+0x15c>
mov    -0x98(%rbp),%rsi
mov    -0x88(%rbp),%rdi
mov    -0x28(%rbp),%rax
movsd  -0x10(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm1
movsd  0xd23(%rip),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
call   *%rax
mov    -0x4(%rbp),%eax
lea    0xf(,%rax,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0xa0(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x70(%rbp)
movl   $0x0,-0x74(%rbp)
mov    -0x74(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1365 <RK4+0x265>
mov    -0x88(%rbp),%rax
mov    -0x98(%rbp),%rcx
mov    -0x20(%rbp),%rdx
movslq -0x74(%rbp),%rsi
movsd  (%rdx,%rsi,8),%xmm1
movsd  -0x18(%rbp),%xmm0
movslq -0x74(%rbp),%rdx
movsd  (%rcx,%rdx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movslq -0x74(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x74(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x74(%rbp)
jmp    1311 <RK4+0x211>
mov    -0xa0(%rbp),%rsi
mov    -0x88(%rbp),%rdi
mov    -0x28(%rbp),%rax
movsd  -0x10(%rbp),%xmm0
addsd  -0x18(%rbp),%xmm0
call   *%rax
movl   $0x0,-0x78(%rbp)
mov    -0x78(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    141d <RK4+0x31d>
mov    -0xa0(%rbp),%rax
mov    -0x98(%rbp),%rcx
mov    -0x80(%rbp),%rdx
mov    -0x90(%rbp),%rsi
mov    -0x20(%rbp),%rdi
movslq -0x78(%rbp),%r8
movsd  (%rdi,%r8,8),%xmm1
movsd  -0x18(%rbp),%xmm0
movsd  0xc3e(%rip),%xmm2
divsd  %xmm2,%xmm0
movslq -0x78(%rbp),%rdi
movsd  (%rsi,%rdi,8),%xmm3
movslq -0x78(%rbp),%rsi
movsd  (%rdx,%rsi,8),%xmm2
movslq -0x78(%rbp),%rdx
addsd  (%rcx,%rdx,8),%xmm2
addsd  %xmm2,%xmm2
addsd  %xmm3,%xmm2
movslq -0x78(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movslq -0x78(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x78(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x78(%rbp)
jmp    138a <RK4+0x28a>
mov    -0x38(%rbp),%rax
mov    %rax,%rsp
mov    %rbp,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<RK4>:
push   %rbp
mov    %rsp,%rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
mov    %rcx,-0x78(%rbp)
mov    %rsi,%rbx
movaps %xmm1,-0x70(%rbp)
movsd  %xmm0,-0x40(%rbp)
mov    %edi,%r14d
mov    %edi,%r13d
mov    %rsp,%r15
lea    0xf(,%r13,8),%r12
and    $0xfffffffffffffff0,%r12
sub    %r12,%r15
mov    %r15,%rsp
mov    %rsi,%rdi
mov    %r15,%rsi
mov    %rdx,-0x38(%rbp)
call   *%rdx
movsd  -0x40(%rbp),%xmm1
mov    %r14d,%eax
movapd -0x70(%rbp),%xmm2
mov    %rsp,%r14
sub    %r12,%r14
mov    %r14,%rsp
mov    %rsp,%rsi
sub    %r12,%rsi
mov    %rsi,%rsp
movsd  0xe91(%rip),%xmm0
mulsd  %xmm2,%xmm0
test   %eax,%eax
jle    118f <RK4+0x8f>
mov    %r12,-0x58(%rbp)
cmp    $0x4,%eax
mov    %eax,-0x2c(%rbp)
jae    11b4 <RK4+0xb4>
xor    %eax,%eax
movsd  -0x40(%rbp),%xmm1
jmp    12b0 <RK4+0x1b0>
addsd  %xmm1,%xmm0
movapd %xmm0,-0x50(%rbp)
mov    %r14,%rdi
mov    -0x38(%rbp),%rbx
call   *%rbx
mov    %rsp,%rsi
sub    %r12,%rsi
mov    %rsi,%rsp
movapd -0x50(%rbp),%xmm0
jmp    131d <RK4+0x21d>
mov    %r13d,%eax
and    $0xfffffffc,%eax
movapd %xmm0,%xmm5
unpcklpd %xmm0,%xmm5
lea    -0x4(%rax),%rcx
mov    %rcx,%rdx
shr    $0x2,%rdx
add    $0x1,%rdx
test   %rcx,%rcx
je     1655 <RK4+0x555>
mov    %rdx,%rdi
and    $0xfffffffffffffffe,%rdi
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movupd (%rbx,%rcx,8),%xmm1
movupd 0x10(%rbx,%rcx,8),%xmm2
movapd (%r15,%rcx,8),%xmm3
mulpd  %xmm5,%xmm3
addpd  %xmm1,%xmm3
movapd 0x10(%r15,%rcx,8),%xmm1
mulpd  %xmm5,%xmm1
addpd  %xmm2,%xmm1
movapd %xmm3,(%r14,%rcx,8)
movapd %xmm1,0x10(%r14,%rcx,8)
movupd 0x20(%rbx,%rcx,8),%xmm1
movupd 0x30(%rbx,%rcx,8),%xmm2
movapd 0x20(%r15,%rcx,8),%xmm3
mulpd  %xmm5,%xmm3
addpd  %xmm1,%xmm3
movapd 0x30(%r15,%rcx,8),%xmm1
mulpd  %xmm5,%xmm1
addpd  %xmm2,%xmm1
movapd %xmm3,0x20(%r14,%rcx,8)
movapd %xmm1,0x30(%r14,%rcx,8)
add    $0x8,%rcx
add    $0xfffffffffffffffe,%rdi
jne    11f0 <RK4+0xf0>
test   $0x1,%dl
je     129d <RK4+0x19d>
movupd (%rbx,%rcx,8),%xmm1
movupd 0x10(%rbx,%rcx,8),%xmm2
movapd (%r15,%rcx,8),%xmm3
mulpd  %xmm5,%xmm3
mulpd  0x10(%r15,%rcx,8),%xmm5
addpd  %xmm1,%xmm3
addpd  %xmm2,%xmm5
movapd %xmm3,(%r14,%rcx,8)
movapd %xmm5,0x10(%r14,%rcx,8)
cmp    %r13,%rax
movsd  -0x40(%rbp),%xmm1
je     12ce <RK4+0x1ce>
nopw   0x0(%rax,%rax,1)
movsd  (%r15,%rax,8),%xmm2
mulsd  %xmm0,%xmm2
addsd  (%rbx,%rax,8),%xmm2
movsd  %xmm2,(%r14,%rax,8)
add    $0x1,%rax
cmp    %rax,%r13
jne    12b0 <RK4+0x1b0>
movapd %xmm0,-0x50(%rbp)
addsd  %xmm1,%xmm0
movapd %xmm0,-0x90(%rbp)
ff
mov    %r14,%rdi
mov    %rsi,%r12
call   *-0x38(%rbp)
mov    %rsp,%rsi
sub    -0x58(%rbp),%rsi
mov    %rsi,%rsp
mov    -0x2c(%rbp),%eax
test   %eax,%eax
jle    130d <RK4+0x20d>
cmp    $0x4,%eax
jae    134d <RK4+0x24d>
xor    %eax,%eax
mov    %r12,%rdx
movapd -0x50(%rbp),%xmm4
jmp    1440 <RK4+0x340>
mov    -0x38(%rbp),%rbx
mov    -0x58(%rbp),%r12
movapd -0x90(%rbp),%xmm0
ff
mov    %r14,%rdi
call   *%rbx
mov    %rsp,%rsi
sub    %r12,%rsi
mov    %rsi,%rsp
movapd -0x70(%rbp),%xmm1
movsd  -0x40(%rbp),%xmm0
addsd  %xmm1,%xmm0
mov    %r14,%rdi
call   *%rbx
lea    -0x28(%rbp),%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    %r13d,%eax
and    $0xfffffffc,%eax
movapd -0x50(%rbp),%xmm4
movapd %xmm4,%xmm0
unpcklpd %xmm4,%xmm0
lea    -0x4(%rax),%rcx
mov    %rcx,%r8
shr    $0x2,%r8
add    $0x1,%r8
test   %rcx,%rcx
mov    %r12,%rdx
je     1665 <RK4+0x565>
mov    %r8,%rdi
and    $0xfffffffffffffffe,%rdi
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movupd (%rbx,%rcx,8),%xmm1
movupd 0x10(%rbx,%rcx,8),%xmm2
movapd (%rdx,%rcx,8),%xmm3
mulpd  %xmm0,%xmm3
addpd  %xmm1,%xmm3
movapd 0x10(%rdx,%rcx,8),%xmm1
mulpd  %xmm0,%xmm1
addpd  %xmm2,%xmm1
movapd %xmm3,(%r14,%rcx,8)
movapd %xmm1,0x10(%r14,%rcx,8)
movupd 0x20(%rbx,%rcx,8),%xmm1
movupd 0x30(%rbx,%rcx,8),%xmm2
movapd 0x20(%rdx,%rcx,8),%xmm3
mulpd  %xmm0,%xmm3
addpd  %xmm1,%xmm3
movapd 0x30(%rdx,%rcx,8),%xmm1
mulpd  %xmm0,%xmm1
addpd  %xmm2,%xmm1
movapd %xmm3,0x20(%r14,%rcx,8)
movapd %xmm1,0x30(%r14,%rcx,8)
add    $0x8,%rcx
add    $0xfffffffffffffffe,%rdi
jne    1390 <RK4+0x290>
test   $0x1,%r8b
je     1438 <RK4+0x338>
movupd (%rbx,%rcx,8),%xmm1
movupd 0x10(%rbx,%rcx,8),%xmm2
movapd (%rdx,%rcx,8),%xmm3
mulpd  %xmm0,%xmm3
mulpd  0x10(%rdx,%rcx,8),%xmm0
addpd  %xmm1,%xmm3
addpd  %xmm2,%xmm0
movapd %xmm3,(%r14,%rcx,8)
movapd %xmm0,0x10(%r14,%rcx,8)
cmp    %r13,%rax
je     145d <RK4+0x35d>
nopl   (%rax)
movsd  (%rdx,%rax,8),%xmm0
mulsd  %xmm4,%xmm0
addsd  (%rbx,%rax,8),%xmm0
movsd  %xmm0,(%r14,%rax,8)
add    $0x1,%rax
cmp    %rax,%r13
jne    1440 <RK4+0x340>
movapd -0x90(%rbp),%xmm0
ff
mov    %r14,%rdi
mov    %rsi,-0x50(%rbp)
call   *-0x38(%rbp)
mov    %rsp,%rsi
sub    -0x58(%rbp),%rsi
mov    %rsi,%rsp
mov    -0x2c(%rbp),%eax
test   %eax,%eax
jle    1495 <RK4+0x395>
cmp    $0x4,%eax
jae    14a3 <RK4+0x3a3>
xor    %eax,%eax
movapd -0x70(%rbp),%xmm4
mov    -0x50(%rbp),%rdi
jmp    1590 <RK4+0x490>
movapd -0x70(%rbp),%xmm1
mov    -0x38(%rbp),%rbx
jmp    1330 <RK4+0x230>
mov    %r13d,%eax
and    $0xfffffffc,%eax
movapd -0x70(%rbp),%xmm4
movapd %xmm4,%xmm0
unpcklpd %xmm4,%xmm0
lea    -0x4(%rax),%rcx
mov    %rcx,%r8
shr    $0x2,%r8
add    $0x1,%r8
test   %rcx,%rcx
mov    -0x50(%rbp),%rdi
je     1676 <RK4+0x576>
mov    %r8,%rdx
and    $0xfffffffffffffffe,%rdx
xor    %ecx,%ecx
nopl   0x0(%rax,%rax,1)
movupd (%rbx,%rcx,8),%xmm1
movupd 0x10(%rbx,%rcx,8),%xmm2
movapd (%rdi,%rcx,8),%xmm3
mulpd  %xmm0,%xmm3
addpd  %xmm1,%xmm3
movapd 0x10(%rdi,%rcx,8),%xmm1
mulpd  %xmm0,%xmm1
addpd  %xmm2,%xmm1
movapd %xmm3,(%r14,%rcx,8)
movapd %xmm1,0x10(%r14,%rcx,8)
movupd 0x20(%rbx,%rcx,8),%xmm1
movupd 0x30(%rbx,%rcx,8),%xmm2
movapd 0x20(%rdi,%rcx,8),%xmm3
mulpd  %xmm0,%xmm3
addpd  %xmm1,%xmm3
movapd 0x30(%rdi,%rcx,8),%xmm1
mulpd  %xmm0,%xmm1
addpd  %xmm2,%xmm1
movapd %xmm3,0x20(%r14,%rcx,8)
movapd %xmm1,0x30(%r14,%rcx,8)
add    $0x8,%rcx
add    $0xfffffffffffffffe,%rdx
jne    14e0 <RK4+0x3e0>
test   $0x1,%r8b
je     1588 <RK4+0x488>
movupd (%rbx,%rcx,8),%xmm1
movupd 0x10(%rbx,%rcx,8),%xmm2
movapd (%rdi,%rcx,8),%xmm3
mulpd  %xmm0,%xmm3
mulpd  0x10(%rdi,%rcx,8),%xmm0
addpd  %xmm1,%xmm3
addpd  %xmm2,%xmm0
movapd %xmm3,(%r14,%rcx,8)
movapd %xmm0,0x10(%r14,%rcx,8)
cmp    %r13,%rax
je     15ad <RK4+0x4ad>
nopl   (%rax)
movsd  (%rdi,%rax,8),%xmm0
mulsd  %xmm4,%xmm0
addsd  (%rbx,%rax,8),%xmm0
movsd  %xmm0,(%r14,%rax,8)
add    $0x1,%rax
cmp    %rax,%r13
jne    1590 <RK4+0x490>
movsd  -0x40(%rbp),%xmm0
addsd  %xmm4,%xmm0
mov    %r14,%rdi
mov    %rsi,%r14
call   *-0x38(%rbp)
movapd -0x70(%rbp),%xmm3
mov    -0x2c(%rbp),%eax
test   %eax,%eax
mov    -0x78(%rbp),%r10
jle    133e <RK4+0x23e>
divsd  0xa2d(%rip),%xmm3
cmp    $0x1,%eax
je     15fa <RK4+0x4fa>
lea    (%rbx,%r13,8),%rax
cmp    %r10,%rax
jbe    1687 <RK4+0x587>
lea    (%r10,%r13,8),%rax
cmp    %rbx,%rax
jbe    1687 <RK4+0x587>
xor    %r8d,%r8d
lea    (%r10,%r8,8),%r9
lea    (%r14,%r8,8),%rdx
mov    -0x50(%rbp),%rax
lea    (%rax,%r8,8),%rsi
lea    (%r12,%r8,8),%rdi
lea    (%r15,%r8,8),%rax
lea    (%rbx,%r8,8),%rbx
sub    %r8,%r13
xor    %ecx,%ecx
xchg   %ax,%ax
movsd  (%rdi,%rcx,8),%xmm0
addsd  (%rsi,%rcx,8),%xmm0
addsd  %xmm0,%xmm0
addsd  (%rax,%rcx,8),%xmm0
addsd  (%rdx,%rcx,8),%xmm0
mulsd  %xmm3,%xmm0
addsd  (%rbx,%rcx,8),%xmm0
movsd  %xmm0,(%r9,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%r13
jne    1620 <RK4+0x520>
jmp    133e <RK4+0x23e>
xor    %ecx,%ecx
test   $0x1,%dl
jne    126c <RK4+0x16c>
jmp    129d <RK4+0x19d>
xor    %ecx,%ecx
test   $0x1,%r8b
jne    1409 <RK4+0x309>
jmp    1438 <RK4+0x338>
xor    %ecx,%ecx
test   $0x1,%r8b
jne    1559 <RK4+0x459>
jmp    1588 <RK4+0x488>
mov    %r13d,%r8d
and    $0xfffffffe,%r8d
movapd %xmm3,%xmm0
unpcklpd %xmm3,%xmm0
lea    -0x2(%r8),%rax
mov    %rax,%r9
shr    %r9
add    $0x1,%r9
test   %rax,%rax
mov    -0x50(%rbp),%rax
mov    %r12,%rdx
mov    %r14,%rdi
je     172a <RK4+0x62a>
mov    %r9,%rsi
and    $0xfffffffffffffffe,%rsi
xor    %ecx,%ecx
nopl   0x0(%rax)
movapd (%rdx,%rcx,8),%xmm1
addpd  (%rax,%rcx,8),%xmm1
movupd (%rbx,%rcx,8),%xmm2
addpd  %xmm1,%xmm1
addpd  (%r15,%rcx,8),%xmm1
addpd  (%rdi,%rcx,8),%xmm1
mulpd  %xmm0,%xmm1
addpd  %xmm2,%xmm1
movupd %xmm1,(%r10,%rcx,8)
movapd 0x10(%rdx,%rcx,8),%xmm1
addpd  0x10(%rax,%rcx,8),%xmm1
addpd  %xmm1,%xmm1
addpd  0x10(%r15,%rcx,8),%xmm1
movupd 0x10(%rbx,%rcx,8),%xmm2
addpd  0x10(%rdi,%rcx,8),%xmm1
mulpd  %xmm0,%xmm1
addpd  %xmm2,%xmm1
movupd %xmm1,0x10(%r10,%rcx,8)
add    $0x4,%rcx
add    $0xfffffffffffffffe,%rsi
jne    16c0 <RK4+0x5c0>
jmp    172c <RK4+0x62c>
xor    %ecx,%ecx
test   $0x1,%r9b
je     175e <RK4+0x65e>
movapd (%rdx,%rcx,8),%xmm1
addpd  (%rax,%rcx,8),%xmm1
addpd  %xmm1,%xmm1
addpd  (%r15,%rcx,8),%xmm1
addpd  (%rdi,%rcx,8),%xmm1
movupd (%rbx,%rcx,8),%xmm2
mulpd  %xmm1,%xmm0
addpd  %xmm2,%xmm0
movupd %xmm0,(%r10,%rcx,8)
cmp    %r13,%r8
je     133e <RK4+0x23e>
jmp    15fd <RK4+0x4fd>

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

LLVM IR (-O2):
define void @RK4(i32 %0, double %1, double %2, double* %3, void (double, double*, double*)* nocapture %4, double* nocapture writeonly %5) {
%7 = zext i32 %0 to i64
%8 = alloca double, i64 %7
%9 = alloca double, i64 %7
%10 = alloca double, i64 %7
%11 = icmp sgt i32 %0, 0
%12 = fmul double %2, 5.000000e-01
br i1 %11, label %13, label %97
13: ; preds = %6
%14 = icmp ult i32 %0, 4
br i1 %14, label %95, label %15
15: ; preds = %13
%16 = and i64 %7, 4294967292
%17 = insertelement <2 x double> poison, double %12, i64 0
%18 = shufflevector <2 x double> %17, <2 x double> poison, <2 x i32> zeroinitializer
%19 = insertelement <2 x double> poison, double %12, i64 0
%20 = shufflevector <2 x double> %19, <2 x double> poison, <2 x i32> zeroinitializer
%21 = add nsw i64 %16, -4
%22 = lshr exact i64 %21, 2
%23 = add nuw nsw i64 %22, 1
%24 = and i64 %23, 1
%25 = icmp eq i64 %21, 0
br i1 %25, label %71, label %26
26: ; preds = %15
%27 = and i64 %23, 9223372036854775806
br label %28
28: ; preds = %28, %26
%29 = phi i64 [ 0, %26 ], [ %68, %28 ]
%30 = phi i64 [ 0, %26 ], [ %69, %28 ]
%31 = getelementptr inbounds double, double* %3, i64 %29
%32 = bitcast double* %31 to <2 x double>*
%33 = load <2 x double>, <2 x double>* %32
%34 = getelementptr inbounds double, double* %31, i64 2
%35 = bitcast double* %34 to <2 x double>*
%36 = load <2 x double>, <2 x double>* %35
%37 = getelementptr inbounds double, double* %8, i64 %29
%38 = bitcast double* %37 to <2 x double>*
%39 = load <2 x double>, <2 x double>* %38
%40 = getelementptr inbounds double, double* %37, i64 2
%41 = bitcast double* %40 to <2 x double>*
%42 = load <2 x double>, <2 x double>* %41
%43 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %18, <2 x double> %39, <2 x double> %33)
%44 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %20, <2 x double> %42, <2 x double> %36)
%45 = getelementptr inbounds double, double* %9, i64 %29
%46 = bitcast double* %45 to <2 x double>*
store <2 x double> %43, <2 x double>* %46
%47 = getelementptr inbounds double, double* %45, i64 2
%48 = bitcast double* %47 to <2 x double>*
store <2 x double> %44, <2 x double>* %48
%49 = or i64 %29, 4
%50 = getelementptr inbounds double, double* %3, i64 %49
%51 = bitcast double* %50 to <2 x double>*
%52 = load <2 x double>, <2 x double>* %51
%53 = getelementptr inbounds double, double* %50, i64 2
%54 = bitcast double* %53 to <2 x double>*
%55 = load <2 x double>, <2 x double>* %54
%56 = getelementptr inbounds double, double* %8, i64 %49
%57 = bitcast double* %56 to <2 x double>*
%58 = load <2 x double>, <2 x double>* %57
%59 = getelementptr inbounds double, double* %56, i64 2
%60 = bitcast double* %59 to <2 x double>*
%61 = load <2 x double>, <2 x double>* %60
%62 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %18, <2 x double> %58, <2 x double> %52)
%63 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %20, <2 x double> %61, <2 x double> %55)
%64 = getelementptr inbounds double, double* %9, i64 %49
%65 = bitcast double* %64 to <2 x double>*
store <2 x double> %62, <2 x double>* %65
%66 = getelementptr inbounds double, double* %64, i64 2
%67 = bitcast double* %66 to <2 x double>*
store <2 x double> %63, <2 x double>* %67
%68 = add nuw i64 %29, 8
%69 = add i64 %30, 2
%70 = icmp eq i64 %69, %27
br i1 %70, label %71, label %28
71: ; preds = %28, %15
%72 = phi i64 [ 0, %15 ], [ %68, %28 ]
%73 = icmp eq i64 %24, 0
br i1 %73, label %93, label %74
74: ; preds = %71
%75 = getelementptr inbounds double, double* %3, i64 %72
%76 = bitcast double* %75 to <2 x double>*
%77 = load <2 x double>, <2 x double>* %76
%78 = getelementptr inbounds double, double* %75, i64 2
%79 = bitcast double* %78 to <2 x double>*
%80 = load <2 x double>, <2 x double>* %79
%81 = getelementptr inbounds double, double* %8, i64 %72
%82 = bitcast double* %81 to <2 x double>*
%83 = load <2 x double>, <2 x double>* %82
%84 = getelementptr inbounds double, double* %81, i64 2
%85 = bitcast double* %84 to <2 x double>*
%86 = load <2 x double>, <2 x double>* %85
%87 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %18, <2 x double> %83, <2 x double> %77)
%88 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %20, <2 x double> %86, <2 x double> %80)
%89 = getelementptr inbounds double, double* %9, i64 %72
%90 = bitcast double* %89 to <2 x double>*
store <2 x double> %87, <2 x double>* %90
%91 = getelementptr inbounds double, double* %89, i64 2
%92 = bitcast double* %91 to <2 x double>*
store <2 x double> %88, <2 x double>* %92
br label %93
93: ; preds = %71, %74
%94 = icmp eq i64 %16, %7
br i1 %94, label %100, label %95
95: ; preds = %13, %93
%96 = phi i64 [ 0, %13 ], [ %16, %93 ]
br label %187
97: ; preds = %6
%98 = fadd double %12, %1
%99 = alloca double, i64 %7
br label %197
100: ; preds = %187, %93
%101 = fadd double %12, %1
%102 = alloca double, i64 %7
br i1 %11, label %103, label %197
103: ; preds = %100
%104 = icmp ult i32 %0, 4
br i1 %104, label %185, label %105
105: ; preds = %103
%106 = and i64 %7, 4294967292
%107 = insertelement <2 x double> poison, double %12, i64 0
%108 = shufflevector <2 x double> %107, <2 x double> poison, <2 x i32> zeroinitializer
%109 = insertelement <2 x double> poison, double %12, i64 0
%110 = shufflevector <2 x double> %109, <2 x double> poison, <2 x i32> zeroinitializer
%111 = add nsw i64 %106, -4
%112 = lshr exact i64 %111, 2
%113 = add nuw nsw i64 %112, 1
%114 = and i64 %113, 1
%115 = icmp eq i64 %111, 0
br i1 %115, label %161, label %116
116: ; preds = %105
%117 = and i64 %113, 9223372036854775806
br label %118
118: ; preds = %118, %116
%119 = phi i64 [ 0, %116 ], [ %158, %118 ]
%120 = phi i64 [ 0, %116 ], [ %159, %118 ]
%121 = getelementptr inbounds double, double* %3, i64 %119
%122 = bitcast double* %121 to <2 x double>*
%123 = load <2 x double>, <2 x double>* %122
%124 = getelementptr inbounds double, double* %121, i64 2
%125 = bitcast double* %124 to <2 x double>*
%126 = load <2 x double>, <2 x double>* %125
%127 = getelementptr inbounds double, double* %10, i64 %119
%128 = bitcast double* %127 to <2 x double>*
%129 = load <2 x double>, <2 x double>* %128
%130 = getelementptr inbounds double, double* %127, i64 2
%131 = bitcast double* %130 to <2 x double>*
%132 = load <2 x double>, <2 x double>* %131
%133 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %108, <2 x double> %129, <2 x double> %123)
%134 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %110, <2 x double> %132, <2 x double> %126)
%135 = getelementptr inbounds double, double* %9, i64 %119
%136 = bitcast double* %135 to <2 x double>*
store <2 x double> %133, <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 %119, 4
%140 = getelementptr inbounds double, double* %3, 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* %10, i64 %139
%147 = bitcast double* %146 to <2 x double>*
%148 = load <2 x double>, <2 x double>* %147
%149 = getelementptr inbounds double, double* %146, i64 2
%150 = bitcast double* %149 to <2 x double>*
%151 = load <2 x double>, <2 x double>* %150
%152 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %108, <2 x double> %148, <2 x double> %142)
%153 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %110, <2 x double> %151, <2 x double> %145)
%154 = getelementptr inbounds double, double* %9, i64 %139
%155 = bitcast double* %154 to <2 x double>*
store <2 x double> %152, <2 x double>* %155
%156 = getelementptr inbounds double, double* %154, i64 2
%157 = bitcast double* %156 to <2 x double>*
store <2 x double> %153, <2 x double>* %157
%158 = add nuw i64 %119, 8
%159 = add i64 %120, 2
%160 = icmp eq i64 %159, %117
br i1 %160, label %161, label %118
161: ; preds = %118, %105
%162 = phi i64 [ 0, %105 ], [ %158, %118 ]
%163 = icmp eq i64 %114, 0
br i1 %163, label %183, label %164
164: ; preds = %161
%165 = getelementptr inbounds double, double* %3, i64 %162
%166 = bitcast double* %165 to <2 x double>*
%167 = load <2 x double>, <2 x double>* %166
%168 = getelementptr inbounds double, double* %165, i64 2
%169 = bitcast double* %168 to <2 x double>*
%170 = load <2 x double>, <2 x double>* %169
%171 = getelementptr inbounds double, double* %10, i64 %162
%172 = bitcast double* %171 to <2 x double>*
%173 = load <2 x double>, <2 x double>* %172
%174 = getelementptr inbounds double, double* %171, i64 2
%175 = bitcast double* %174 to <2 x double>*
%176 = load <2 x double>, <2 x double>* %175
%177 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %108, <2 x double> %173, <2 x double> %167)
%178 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %110, <2 x double> %176, <2 x double> %170)
%179 = getelementptr inbounds double, double* %9, i64 %162
%180 = bitcast double* %179 to <2 x double>*
store <2 x double> %177, <2 x double>* %180
%181 = getelementptr inbounds double, double* %179, i64 2
%182 = bitcast double* %181 to <2 x double>*
store <2 x double> %178, <2 x double>* %182
br label %183
183: ; preds = %161, %164
%184 = icmp eq i64 %106, %7
br i1 %184, label %201, label %185
185: ; preds = %103, %183
%186 = phi i64 [ 0, %103 ], [ %106, %183 ]
br label %287
187: ; preds = %95, %187
%188 = phi i64 [ %195, %187 ], [ %96, %95 ]
%189 = getelementptr inbounds double, double* %3, i64 %188
%190 = load double, double* %189
%191 = getelementptr inbounds double, double* %8, i64 %188
%192 = load double, double* %191
%193 = call double @llvm.fmuladd.f64(double %12, double %192, double %190)
%194 = getelementptr inbounds double, double* %9, i64 %188
store double %193, double* %194
%195 = add nuw nsw i64 %188, 1
%196 = icmp eq i64 %195, %7
br i1 %196, label %100, label %187
197: ; preds = %97, %100
%198 = phi double* [ %102, %100 ], [ %99, %97 ]
%199 = phi double [ %101, %100 ], [ %98, %97 ]
%200 = alloca double, i64 %7
br label %297
201: ; preds = %287, %183
%202 = alloca double, i64 %7
br i1 %11, label %203, label %297
203: ; preds = %201
%204 = icmp ult i32 %0, 4
br i1 %204, label %285, label %205
205: ; preds = %203
%206 = and i64 %7, 4294967292
%207 = insertelement <2 x double> poison, double %2, i64 0
%208 = shufflevector <2 x double> %207, <2 x double> poison, <2 x i32> zeroinitializer
%209 = insertelement <2 x double> poison, double %2, i64 0
%210 = shufflevector <2 x double> %209, <2 x double> poison, <2 x i32> zeroinitializer
%211 = add nsw i64 %206, -4
%212 = lshr exact i64 %211, 2
%213 = add nuw nsw i64 %212, 1
%214 = and i64 %213, 1
%215 = icmp eq i64 %211, 0
br i1 %215, label %261, label %216
216: ; preds = %205
%217 = and i64 %213, 9223372036854775806
br label %218
218: ; preds = %218, %216
%219 = phi i64 [ 0, %216 ], [ %258, %218 ]
%220 = phi i64 [ 0, %216 ], [ %259, %218 ]
%221 = getelementptr inbounds double, double* %3, i64 %219
%222 = bitcast double* %221 to <2 x double>*
%223 = load <2 x double>, <2 x double>* %222
%224 = getelementptr inbounds double, double* %221, i64 2
%225 = bitcast double* %224 to <2 x double>*
%226 = load <2 x double>, <2 x double>* %225
%227 = getelementptr inbounds double, double* %102, i64 %219
%228 = bitcast double* %227 to <2 x double>*
%229 = load <2 x double>, <2 x double>* %228
%230 = getelementptr inbounds double, double* %227, i64 2
%231 = bitcast double* %230 to <2 x double>*
%232 = load <2 x double>, <2 x double>* %231
%233 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %208, <2 x double> %229, <2 x double> %223)
%234 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %210, <2 x double> %232, <2 x double> %226)
%235 = getelementptr inbounds double, double* %9, i64 %219
%236 = bitcast double* %235 to <2 x double>*
store <2 x double> %233, <2 x double>* %236
%237 = getelementptr inbounds double, double* %235, i64 2
%238 = bitcast double* %237 to <2 x double>*
store <2 x double> %234, <2 x double>* %238
%239 = or i64 %219, 4
%240 = getelementptr inbounds double, double* %3, i64 %239
%241 = bitcast double* %240 to <2 x double>*
%242 = load <2 x double>, <2 x double>* %241
%243 = getelementptr inbounds double, double* %240, i64 2
%244 = bitcast double* %243 to <2 x double>*
%245 = load <2 x double>, <2 x double>* %244
%246 = getelementptr inbounds double, double* %102, i64 %239
%247 = bitcast double* %246 to <2 x double>*
%248 = load <2 x double>, <2 x double>* %247
%249 = getelementptr inbounds double, double* %246, i64 2
%250 = bitcast double* %249 to <2 x double>*
%251 = load <2 x double>, <2 x double>* %250
%252 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %208, <2 x double> %248, <2 x double> %242)
%253 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %210, <2 x double> %251, <2 x double> %245)
%254 = getelementptr inbounds double, double* %9, i64 %239
%255 = bitcast double* %254 to <2 x double>*
store <2 x double> %252, <2 x double>* %255
%256 = getelementptr inbounds double, double* %254, i64 2
%257 = bitcast double* %256 to <2 x double>*
store <2 x double> %253, <2 x double>* %257
%258 = add nuw i64 %219, 8
%259 = add i64 %220, 2
%260 = icmp eq i64 %259, %217
br i1 %260, label %261, label %218
261: ; preds = %218, %205
%262 = phi i64 [ 0, %205 ], [ %258, %218 ]
%263 = icmp eq i64 %214, 0
br i1 %263, label %283, label %264
264: ; preds = %261
%265 = getelementptr inbounds double, double* %3, i64 %262
%266 = bitcast double* %265 to <2 x double>*
%267 = load <2 x double>, <2 x double>* %266
%268 = getelementptr inbounds double, double* %265, i64 2
%269 = bitcast double* %268 to <2 x double>*
%270 = load <2 x double>, <2 x double>* %269
%271 = getelementptr inbounds double, double* %102, i64 %262
%272 = bitcast double* %271 to <2 x double>*
%273 = load <2 x double>, <2 x double>* %272
%274 = getelementptr inbounds double, double* %271, i64 2
%275 = bitcast double* %274 to <2 x double>*
%276 = load <2 x double>, <2 x double>* %275
%277 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %208, <2 x double> %273, <2 x double> %267)
%278 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %210, <2 x double> %276, <2 x double> %270)
%279 = getelementptr inbounds double, double* %9, i64 %262
%280 = bitcast double* %279 to <2 x double>*
store <2 x double> %277, <2 x double>* %280
%281 = getelementptr inbounds double, double* %279, i64 2
%282 = bitcast double* %281 to <2 x double>*
store <2 x double> %278, <2 x double>* %282
br label %283
283: ; preds = %261, %264
%284 = icmp eq i64 %206, %7
br i1 %284, label %300, label %285
285: ; preds = %203, %283
%286 = phi i64 [ 0, %203 ], [ %206, %283 ]
br label %400
287: ; preds = %185, %287
%288 = phi i64 [ %295, %287 ], [ %186, %185 ]
%289 = getelementptr inbounds double, double* %3, i64 %288
%290 = load double, double* %289
%291 = getelementptr inbounds double, double* %10, i64 %288
%292 = load double, double* %291
%293 = call double @llvm.fmuladd.f64(double %12, double %292, double %290)
%294 = getelementptr inbounds double, double* %9, i64 %288
store double %293, double* %294
%295 = add nuw nsw i64 %288, 1
%296 = icmp eq i64 %295, %7
br i1 %296, label %201, label %287
297: ; preds = %197, %201
%298 = phi double* [ %200, %197 ], [ %202, %201 ]
%299 = fadd double %1, %2
br label %410
300: ; preds = %400, %283
%301 = fadd double %1, %2
br i1 %11, label %302, label %410
302: ; preds = %300
%303 = fdiv double %2, 6.000000e+00
%304 = icmp eq i32 %0, 1
br i1 %304, label %398, label %305
305: ; preds = %302
%306 = getelementptr double, double* %5, i64 %7
%307 = getelementptr double, double* %3, i64 %7
%308 = icmp ugt double* %307, %5
%309 = icmp ugt double* %306, %3
%310 = and i1 %308, %309
br i1 %310, label %398, label %311
311: ; preds = %305
%312 = and i64 %7, 4294967294
%313 = insertelement <2 x double> poison, double %303, i64 0
%314 = shufflevector <2 x double> %313, <2 x double> poison, <2 x i32> zeroinitializer
%315 = add nsw i64 %312, -2
%316 = lshr exact i64 %315, 1
%317 = add nuw i64 %316, 1
%318 = and i64 %317, 1
%319 = icmp eq i64 %315, 0
br i1 %319, label %371, label %320
320: ; preds = %311
%321 = and i64 %317, -2
br label %322
322: ; preds = %322, %320
%323 = phi i64 [ 0, %320 ], [ %368, %322 ]
%324 = phi i64 [ 0, %320 ], [ %369, %322 ]
%325 = getelementptr inbounds double, double* %3, i64 %323
%326 = bitcast double* %325 to <2 x double>*
%327 = load <2 x double>, <2 x double>* %326
%328 = getelementptr inbounds double, double* %8, i64 %323
%329 = bitcast double* %328 to <2 x double>*
%330 = load <2 x double>, <2 x double>* %329
%331 = getelementptr inbounds double, double* %10, i64 %323
%332 = bitcast double* %331 to <2 x double>*
%333 = load <2 x double>, <2 x double>* %332
%334 = getelementptr inbounds double, double* %102, i64 %323
%335 = bitcast double* %334 to <2 x double>*
%336 = load <2 x double>, <2 x double>* %335
%337 = fadd <2 x double> %333, %336
%338 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %337, <2 x double> <double 2.000000e+00, double 2.000000e+00>, <2 x double> %330)
%339 = getelementptr inbounds double, double* %202, i64 %323
%340 = bitcast double* %339 to <2 x double>*
%341 = load <2 x double>, <2 x double>* %340
%342 = fadd <2 x double> %341, %338
%343 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %314, <2 x double> %342, <2 x double> %327)
%344 = getelementptr inbounds double, double* %5, i64 %323
%345 = bitcast double* %344 to <2 x double>*
store <2 x double> %343, <2 x double>* %345
%346 = or i64 %323, 2
%347 = getelementptr inbounds double, double* %3, i64 %346
%348 = bitcast double* %347 to <2 x double>*
%349 = load <2 x double>, <2 x double>* %348
%350 = getelementptr inbounds double, double* %8, i64 %346
%351 = bitcast double* %350 to <2 x double>*
%352 = load <2 x double>, <2 x double>* %351
%353 = getelementptr inbounds double, double* %10, i64 %346
%354 = bitcast double* %353 to <2 x double>*
%355 = load <2 x double>, <2 x double>* %354
%356 = getelementptr inbounds double, double* %102, i64 %346
%357 = bitcast double* %356 to <2 x double>*
%358 = load <2 x double>, <2 x double>* %357
%359 = fadd <2 x double> %355, %358
%360 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %359, <2 x double> <double 2.000000e+00, double 2.000000e+00>, <2 x double> %352)
%361 = getelementptr inbounds double, double* %202, i64 %346
%362 = bitcast double* %361 to <2 x double>*
%363 = load <2 x double>, <2 x double>* %362
%364 = fadd <2 x double> %363, %360
%365 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %314, <2 x double> %364, <2 x double> %349)
%366 = getelementptr inbounds double, double* %5, i64 %346
%367 = bitcast double* %366 to <2 x double>*
store <2 x double> %365, <2 x double>* %367
%368 = add nuw i64 %323, 4
%369 = add nuw i64 %324, 2
%370 = icmp eq i64 %369, %321
br i1 %370, label %371, label %322
371: ; preds = %322, %311
%372 = phi i64 [ 0, %311 ], [ %368, %322 ]
%373 = icmp eq i64 %318, 0
br i1 %373, label %396, label %374
374: ; preds = %371
%375 = getelementptr inbounds double, double* %3, i64 %372
%376 = bitcast double* %375 to <2 x double>*
%377 = load <2 x double>, <2 x double>* %376
%378 = getelementptr inbounds double, double* %8, i64 %372
%379 = bitcast double* %378 to <2 x double>*
%380 = load <2 x double>, <2 x double>* %379
%381 = getelementptr inbounds double, double* %10, i64 %372
%382 = bitcast double* %381 to <2 x double>*
%383 = load <2 x double>, <2 x double>* %382
%384 = getelementptr inbounds double, double* %102, i64 %372
%385 = bitcast double* %384 to <2 x double>*
%386 = load <2 x double>, <2 x double>* %385
%387 = fadd <2 x double> %383, %386
%388 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %387, <2 x double> <double 2.000000e+00, double 2.000000e+00>, <2 x double> %380)
%389 = getelementptr inbounds double, double* %202, i64 %372
%390 = bitcast double* %389 to <2 x double>*
%391 = load <2 x double>, <2 x double>* %390
%392 = fadd <2 x double> %391, %388
%393 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %314, <2 x double> %392, <2 x double> %377)
%394 = getelementptr inbounds double, double* %5, i64 %372
%395 = bitcast double* %394 to <2 x double>*
store <2 x double> %393, <2 x double>* %395
br label %396
396: ; preds = %371, %374
%397 = icmp eq i64 %312, %7
br i1 %397, label %410, label %398
398: ; preds = %305, %302, %396
%399 = phi i64 [ 0, %305 ], [ 0, %302 ], [ %312, %396 ]
br label %411
400: ; preds = %285, %400
%401 = phi i64 [ %408, %400 ], [ %286, %285 ]
%402 = getelementptr inbounds double, double* %3, i64 %401
%403 = load double, double* %402
%404 = getelementptr inbounds double, double* %102, i64 %401
%405 = load double, double* %404
%406 = call double @llvm.fmuladd.f64(double %2, double %405, double %403)
%407 = getelementptr inbounds double, double* %9, i64 %401
store double %406, double* %407
%408 = add nuw nsw i64 %401, 1
%409 = icmp eq i64 %408, %7
br i1 %409, label %300, label %400
410: ; preds = %411, %396, %297, %300
ret void
411: ; preds = %398, %411
%412 = phi i64 [ %428, %411 ], [ %399, %398 ]
%413 = getelementptr inbounds double, double* %3, i64 %412
%414 = load double, double* %413
%415 = getelementptr inbounds double, double* %8, i64 %412
%416 = load double, double* %415
%417 = getelementptr inbounds double, double* %10, i64 %412
%418 = load double, double* %417
%419 = getelementptr inbounds double, double* %102, i64 %412
%420 = load double, double* %419
%421 = fadd double %418, %420
%422 = call double @llvm.fmuladd.f64(double %421, double 2.000000e+00, double %416)
%423 = getelementptr inbounds double, double* %202, i64 %412
%424 = load double, double* %423
%425 = fadd double %424, %422
%426 = call double @llvm.fmuladd.f64(double %303, double %425, double %414)
%427 = getelementptr inbounds double, double* %5, i64 %412
store double %426, double* %427
%428 = add nuw nsw i64 %412, 1
%429 = icmp eq i64 %428, %7
br i1 %429, label %410, label %411
}