void one_step_rk4(y, dt, yval, neq, tim) double dt, *tim, *yval[3], *y;
int neq;
{
  int i;
  double t = *tim, t1, t2;
  set_wieners(dt, y, t);
  rhs(t, y, yval[1], neq);
  for (i = 0; i < neq; i++) {
    yval[0][i] = y[i] + dt * yval[1][i] / 6.00;
    yval[2][i] = y[i] + dt * yval[1][i] * 0.5;
  }
  t1 = t + .5 * dt;
  rhs(t1, yval[2], yval[1], neq);
  for (i = 0; i < neq; i++) {
    yval[0][i] = yval[0][i] + dt * yval[1][i] / 3.00;
    yval[2][i] = y[i] + .5 * dt * yval[1][i];
  }
  rhs(t1, yval[2], yval[1], neq);
  for (i = 0; i < neq; i++) {
    yval[0][i] = yval[0][i] + dt * yval[1][i] / 3.000;
    yval[2][i] = y[i] + dt * yval[1][i];
  }
  t2 = t + dt;
  rhs(t2, yval[2], yval[1], neq);
  for (i = 0; i < neq; i++)
    y[i] = yval[0][i] + dt * yval[1][i] / 6.00;
  *tim = t2;
}


x86 assembly (clang, -O0):
<one_step_rk4>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x8(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
mov    %rcx,-0x28(%rbp)
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x10(%rbp),%xmm0
mov    -0x8(%rbp),%rdi
movsd  -0x38(%rbp),%xmm1
mov    $0x2,%al
call   1040 <set_wieners@plt>
movsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rdi
mov    -0x18(%rbp),%rax
mov    0x8(%rax),%rsi
mov    -0x1c(%rbp),%edx
mov    $0x1,%al
call   1030 <rhs@plt>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1221 <one_step_rk4+0x101>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x10(%rbp),%xmm1
mov    -0x18(%rbp),%rax
mov    0x8(%rax),%rax
movslq -0x2c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
movsd  0xe51(%rip),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    (%rax),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    0x8(%rax),%rax
movslq -0x2c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  0xe06(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    0x10(%rax),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1180 <one_step_rk4+0x60>
movsd  -0x38(%rbp),%xmm1
movsd  -0x10(%rbp),%xmm0
movsd  0xdcd(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    0x10(%rax),%rdi
mov    -0x18(%rbp),%rax
mov    0x8(%rax),%rsi
mov    -0x1c(%rbp),%edx
mov    $0x1,%al
call   1030 <rhs@plt>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    130a <one_step_rk4+0x1ea>
mov    -0x18(%rbp),%rax
mov    (%rax),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x10(%rbp),%xmm1
mov    -0x18(%rbp),%rax
mov    0x8(%rax),%rax
movslq -0x2c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
movsd  0xd70(%rip),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    (%rax),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  0xd33(%rip),%xmm0
mulsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    0x8(%rax),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    0x10(%rax),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1266 <one_step_rk4+0x146>
movsd  -0x40(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    0x10(%rax),%rdi
mov    -0x18(%rbp),%rax
mov    0x8(%rax),%rsi
mov    -0x1c(%rbp),%edx
mov    $0x1,%al
call   1030 <rhs@plt>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    13cc <one_step_rk4+0x2ac>
mov    -0x18(%rbp),%rax
mov    (%rax),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x10(%rbp),%xmm1
mov    -0x18(%rbp),%rax
mov    0x8(%rax),%rax
movslq -0x2c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
movsd  0xca6(%rip),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    (%rax),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    0x8(%rax),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    0x10(%rax),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1330 <one_step_rk4+0x210>
movsd  -0x38(%rbp),%xmm0
addsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    0x10(%rax),%rdi
mov    -0x18(%rbp),%rax
mov    0x8(%rax),%rsi
mov    -0x1c(%rbp),%edx
mov    $0x1,%al
call   1030 <rhs@plt>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    145e <one_step_rk4+0x33e>
mov    -0x18(%rbp),%rax
mov    (%rax),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x10(%rbp),%xmm1
mov    -0x18(%rbp),%rax
mov    0x8(%rax),%rax
movslq -0x2c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
movsd  0xbcd(%rip),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1401 <one_step_rk4+0x2e1>
movsd  -0x48(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<one_step_rk4>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %rcx,%r14
mov    %edx,%r15d
mov    %rsi,%r12
movaps %xmm0,0x30(%rsp)
mov    %rdi,%rbx
movsd  (%rcx),%xmm1
movsd  %xmm1,0x10(%rsp)
mov    $0x2,%al
call   1040 <set_wieners@plt>
mov    0x8(%r12),%rsi
movsd  0x10(%rsp),%xmm0
mov    %rbx,%rdi
mov    %r15d,%edx
mov    $0x1,%al
call   1030 <rhs@plt>
movapd 0x30(%rsp),%xmm1
movsd  0x10(%rsp),%xmm2
test   %r15d,%r15d
jle    12cd <one_step_rk4+0x1ad>
mov    (%r12),%rax
mov    0x8(%r12),%rsi
mov    0x10(%r12),%rdi
mov    %r15d,%ebp
cmp    $0x1,%r15d
mov    %r14,0x18(%rsp)
jne    1354 <one_step_rk4+0x234>
xor    %ecx,%ecx
movsd  0x10(%rsp),%xmm5
movapd 0x30(%rsp),%xmm6
movsd  0xe55(%rip),%xmm0
movsd  0xe45(%rip),%xmm1
nopl   0x0(%rax,%rax,1)
movsd  (%rsi,%rcx,8),%xmm2
mulsd  %xmm6,%xmm2
divsd  %xmm0,%xmm2
addsd  (%rbx,%rcx,8),%xmm2
movsd  %xmm2,(%rax,%rcx,8)
movsd  (%rsi,%rcx,8),%xmm2
mulsd  %xmm6,%xmm2
mulsd  %xmm1,%xmm2
addsd  (%rbx,%rcx,8),%xmm2
movsd  %xmm2,(%rdi,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%rbp
jne    11c0 <one_step_rk4+0xa0>
movsd  0xe01(%rip),%xmm0
mulsd  %xmm6,%xmm0
movapd %xmm0,0x40(%rsp)
addsd  %xmm5,%xmm0
movapd %xmm0,0x20(%rsp)
mov    %r15d,%edx
mov    $0x1,%al
call   1030 <rhs@plt>
test   %r15d,%r15d
jle    12f7 <one_step_rk4+0x1d7>
mov    (%r12),%rax
mov    0x8(%r12),%rsi
mov    0x10(%r12),%rdi
cmp    $0x1,%r15d
jne    146d <one_step_rk4+0x34d>
xor    %ecx,%ecx
movapd 0x30(%rsp),%xmm5
movaps 0x20(%rsp),%xmm0
movapd 0x40(%rsp),%xmm7
movsd  0xdb7(%rip),%xmm2
nopl   0x0(%rax)
movsd  (%rsi,%rcx,8),%xmm1
mulsd  %xmm5,%xmm1
divsd  %xmm2,%xmm1
addsd  (%rax,%rcx,8),%xmm1
movsd  %xmm1,(%rax,%rcx,8)
movsd  (%rsi,%rcx,8),%xmm1
mulsd  %xmm7,%xmm1
addsd  (%rbx,%rcx,8),%xmm1
movsd  %xmm1,(%rdi,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%rbp
jne    1260 <one_step_rk4+0x140>
mov    %r15d,%edx
mov    $0x1,%al
call   1030 <rhs@plt>
movapd 0x30(%rsp),%xmm5
test   %r15d,%r15d
movsd  0x10(%rsp),%xmm0
jle    131c <one_step_rk4+0x1fc>
mov    (%r12),%rax
mov    0x8(%r12),%rsi
mov    0x10(%r12),%rdi
cmp    $0x1,%r15d
jne    1579 <one_step_rk4+0x459>
xor    %ecx,%ecx
jmp    167b <one_step_rk4+0x55b>
mov    0x8(%r12),%rsi
mov    0x10(%r12),%rdi
movsd  0xd21(%rip),%xmm0
mulsd  %xmm1,%xmm0
addsd  %xmm2,%xmm0
movapd %xmm0,0x20(%rsp)
mov    %r15d,%edx
mov    $0x1,%al
call   1030 <rhs@plt>
movaps 0x20(%rsp),%xmm0
mov    0x8(%r12),%rsi
mov    0x10(%r12),%rdi
mov    %r15d,%edx
mov    $0x1,%al
call   1030 <rhs@plt>
movsd  0x10(%rsp),%xmm0
movapd 0x30(%rsp),%xmm5
mov    0x8(%r12),%rsi
mov    0x10(%r12),%rdi
addsd  %xmm5,%xmm0
movsd  %xmm0,0x10(%rsp)
mov    %r15d,%edx
mov    $0x1,%al
call   1030 <rhs@plt>
movsd  0x10(%rsp),%xmm4
movsd  %xmm4,(%r14)
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
lea    (%rax,%rbp,8),%rcx
lea    (%rdi,%rbp,8),%r14
lea    (%rbx,%rbp,8),%r9
cmp    %r14,%rax
setb   0x20(%rsp)
cmp    %rcx,%rdi
setb   %r11b
cmp    %r9,%rax
setb   %r13b
cmp    %rbx,%rcx
seta   0x40(%rsp)
lea    (%rsi,%rbp,8),%rdx
cmp    %rdx,%rax
setb   %r8b
cmp    %rcx,%rsi
setb   0xd(%rsp)
cmp    %r9,%rdi
setb   %r9b
cmp    %rbx,%r14
seta   0xe(%rsp)
cmp    %rdx,%rdi
setb   %r10b
cmp    %r14,%rsi
setb   %dl
xor    %ecx,%ecx
test   %r11b,0x20(%rsp)
movsd  0x10(%rsp),%xmm5
movapd 0x30(%rsp),%xmm6
jne    1662 <one_step_rk4+0x542>
and    0x40(%rsp),%r13b
jne    1662 <one_step_rk4+0x542>
and    0xd(%rsp),%r8b
mov    0x18(%rsp),%r14
jne    11ab <one_step_rk4+0x8b>
and    0xe(%rsp),%r9b
jne    11ab <one_step_rk4+0x8b>
and    %dl,%r10b
jne    11ab <one_step_rk4+0x8b>
mov    %ebp,%ecx
and    $0xfffffffe,%ecx
movapd %xmm6,%xmm0
unpcklpd %xmm6,%xmm0
xor    %edx,%edx
movapd 0xc14(%rip),%xmm1
movapd 0xc1c(%rip),%xmm2
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movupd (%rbx,%rdx,8),%xmm3
movupd (%rsi,%rdx,8),%xmm4
mulpd  %xmm0,%xmm4
divpd  %xmm1,%xmm4
addpd  %xmm3,%xmm4
movupd %xmm4,(%rax,%rdx,8)
movupd (%rbx,%rdx,8),%xmm3
movupd (%rsi,%rdx,8),%xmm4
mulpd  %xmm0,%xmm4
mulpd  %xmm2,%xmm4
addpd  %xmm3,%xmm4
movupd %xmm4,(%rdi,%rdx,8)
add    $0x2,%rdx
cmp    %rdx,%rcx
jne    1420 <one_step_rk4+0x300>
cmp    %rbp,%rcx
jne    11ab <one_step_rk4+0x8b>
jmp    11f7 <one_step_rk4+0xd7>
lea    (%rax,%rbp,8),%rcx
lea    (%rdi,%rbp,8),%r13
lea    (%rsi,%rbp,8),%rdx
lea    (%rbx,%rbp,8),%r10
cmp    %r13,%rax
setb   0xd(%rsp)
cmp    %rcx,%rdi
setb   %r11b
cmp    %rdx,%rax
setb   %r14b
cmp    %rcx,%rsi
setb   0xe(%rsp)
cmp    %r10,%rax
setb   %r8b
cmp    %rbx,%rcx
seta   0xf(%rsp)
cmp    %rdx,%rdi
setb   %r9b
cmp    %r13,%rsi
setb   %dl
cmp    %r10,%rdi
setb   %r10b
cmp    %rbx,%r13
seta   %r13b
xor    %ecx,%ecx
test   %r11b,0xd(%rsp)
jne    166c <one_step_rk4+0x54c>
and    0xe(%rsp),%r14b
jne    166c <one_step_rk4+0x54c>
and    0xf(%rsp),%r8b
mov    0x18(%rsp),%r14
movapd 0x30(%rsp),%xmm5
movaps 0x20(%rsp),%xmm0
movapd 0x40(%rsp),%xmm7
jne    1251 <one_step_rk4+0x131>
and    %dl,%r9b
jne    1251 <one_step_rk4+0x131>
and    %r13b,%r10b
jne    1251 <one_step_rk4+0x131>
mov    %ebp,%ecx
and    $0xfffffffe,%ecx
movapd %xmm5,%xmm6
unpcklpd %xmm5,%xmm6
movapd %xmm7,%xmm1
unpcklpd %xmm7,%xmm1
xor    %edx,%edx
movapd 0xb11(%rip),%xmm2
nop
movupd (%rax,%rdx,8),%xmm3
movupd (%rsi,%rdx,8),%xmm4
mulpd  %xmm6,%xmm4
divpd  %xmm2,%xmm4
addpd  %xmm3,%xmm4
movupd %xmm4,(%rax,%rdx,8)
movupd (%rbx,%rdx,8),%xmm3
movupd (%rsi,%rdx,8),%xmm4
mulpd  %xmm1,%xmm4
addpd  %xmm3,%xmm4
movupd %xmm4,(%rdi,%rdx,8)
add    $0x2,%rdx
cmp    %rdx,%rcx
jne    1530 <one_step_rk4+0x410>
cmp    %rbp,%rcx
jne    1251 <one_step_rk4+0x131>
jmp    1293 <one_step_rk4+0x173>
lea    (%rax,%rbp,8),%rcx
lea    (%rdi,%rbp,8),%r13
lea    (%rsi,%rbp,8),%rdx
lea    (%rbx,%rbp,8),%r10
cmp    %r13,%rax
setb   0x20(%rsp)
cmp    %rcx,%rdi
setb   %r11b
cmp    %rdx,%rax
setb   %r14b
cmp    %rcx,%rsi
setb   0x40(%rsp)
cmp    %r10,%rax
setb   %r8b
cmp    %rbx,%rcx
seta   0xd(%rsp)
cmp    %rdx,%rdi
setb   %r9b
cmp    %r13,%rsi
setb   %dl
cmp    %r10,%rdi
setb   %r10b
cmp    %rbx,%r13
seta   %r13b
xor    %ecx,%ecx
test   %r11b,0x20(%rsp)
jne    1676 <one_step_rk4+0x556>
and    0x40(%rsp),%r14b
jne    1676 <one_step_rk4+0x556>
and    0xd(%rsp),%r8b
mov    0x18(%rsp),%r14
jne    167b <one_step_rk4+0x55b>
and    %dl,%r9b
jne    167b <one_step_rk4+0x55b>
and    %r13b,%r10b
jne    167b <one_step_rk4+0x55b>
mov    %ebp,%ecx
and    $0xfffffffe,%ecx
movapd %xmm5,%xmm4
unpcklpd %xmm5,%xmm4
xor    %edx,%edx
movapd 0xa22(%rip),%xmm1
xchg   %ax,%ax
movupd (%rax,%rdx,8),%xmm2
movupd (%rsi,%rdx,8),%xmm3
mulpd  %xmm4,%xmm3
divpd  %xmm1,%xmm3
addpd  %xmm2,%xmm3
movupd %xmm3,(%rax,%rdx,8)
movupd (%rbx,%rdx,8),%xmm2
movupd (%rsi,%rdx,8),%xmm3
mulpd  %xmm4,%xmm3
addpd  %xmm2,%xmm3
movupd %xmm3,(%rdi,%rdx,8)
add    $0x2,%rdx
cmp    %rdx,%rcx
jne    1620 <one_step_rk4+0x500>
cmp    %rbp,%rcx
je     16c3 <one_step_rk4+0x5a3>
jmp    167b <one_step_rk4+0x55b>
mov    0x18(%rsp),%r14
jmp    11ab <one_step_rk4+0x8b>
mov    0x18(%rsp),%r14
jmp    1240 <one_step_rk4+0x120>
mov    0x18(%rsp),%r14
movsd  0x98d(%rip),%xmm2
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%rsi,%rcx,8),%xmm1
mulsd  %xmm5,%xmm1
divsd  %xmm2,%xmm1
addsd  (%rax,%rcx,8),%xmm1
movsd  %xmm1,(%rax,%rcx,8)
movsd  (%rsi,%rcx,8),%xmm1
mulsd  %xmm5,%xmm1
addsd  (%rbx,%rcx,8),%xmm1
movsd  %xmm1,(%rdi,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%rbp
jne    1690 <one_step_rk4+0x570>
addsd  %xmm5,%xmm0
movsd  %xmm0,0x10(%rsp)
mov    %r15d,%edx
mov    $0x1,%al
call   1030 <rhs@plt>
test   %r15d,%r15d
jle    133a <one_step_rk4+0x21a>
mov    (%r12),%rax
mov    0x8(%r12),%rcx
cmp    $0x1,%r15d
movsd  0x10(%rsp),%xmm4
jne    177f <one_step_rk4+0x65f>
xor    %edx,%edx
movapd 0x30(%rsp),%xmm5
mov    %rdx,%rsi
not    %rsi
test   $0x1,%bpl
je     172c <one_step_rk4+0x60c>
movsd  (%rcx,%rdx,8),%xmm0
mulsd  %xmm5,%xmm0
divsd  0x8ea(%rip),%xmm0
addsd  (%rax,%rdx,8),%xmm0
movsd  %xmm0,(%rbx,%rdx,8)
or     $0x1,%rdx
add    %rbp,%rsi
je     1340 <one_step_rk4+0x220>
movsd  0x8cb(%rip),%xmm0
nopl   (%rax)
movsd  (%rcx,%rdx,8),%xmm1
mulsd  %xmm5,%xmm1
divsd  %xmm0,%xmm1
addsd  (%rax,%rdx,8),%xmm1
movsd  %xmm1,(%rbx,%rdx,8)
movsd  0x8(%rcx,%rdx,8),%xmm1
mulsd  %xmm5,%xmm1
divsd  %xmm0,%xmm1
addsd  0x8(%rax,%rdx,8),%xmm1
movsd  %xmm1,0x8(%rbx,%rdx,8)
add    $0x2,%rdx
cmp    %rdx,%rbp
jne    1740 <one_step_rk4+0x620>
jmp    1340 <one_step_rk4+0x220>
lea    (%rbx,%rbp,8),%rdx
lea    (%rax,%rbp,8),%rsi
lea    (%rcx,%rbp,8),%rdi
cmp    %rbx,%rsi
seta   %r8b
cmp    %rdx,%rax
setb   %r9b
cmp    %rbx,%rdi
seta   %sil
cmp    %rdx,%rcx
setb   %dil
xor    %edx,%edx
test   %r9b,%r8b
movapd 0x30(%rsp),%xmm5
jne    1701 <one_step_rk4+0x5e1>
and    %dil,%sil
jne    1701 <one_step_rk4+0x5e1>
mov    %ebp,%edx
and    $0xfffffffe,%edx
movapd %xmm5,%xmm0
unpcklpd %xmm5,%xmm0
lea    -0x2(%rdx),%rsi
mov    %rsi,%r8
shr    %r8
add    $0x1,%r8
test   %rsi,%rsi
je     1845 <one_step_rk4+0x725>
mov    %r8,%rdi
and    $0xfffffffffffffffe,%rdi
xor    %esi,%esi
movapd 0x82e(%rip),%xmm1
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movupd (%rax,%rsi,8),%xmm2
movupd (%rcx,%rsi,8),%xmm3
mulpd  %xmm0,%xmm3
divpd  %xmm1,%xmm3
addpd  %xmm2,%xmm3
movupd %xmm3,(%rbx,%rsi,8)
movupd 0x10(%rax,%rsi,8),%xmm2
movupd 0x10(%rcx,%rsi,8),%xmm3
mulpd  %xmm0,%xmm3
divpd  %xmm1,%xmm3
addpd  %xmm2,%xmm3
movupd %xmm3,0x10(%rbx,%rsi,8)
add    $0x4,%rsi
add    $0xfffffffffffffffe,%rdi
jne    1800 <one_step_rk4+0x6e0>
jmp    1847 <one_step_rk4+0x727>
xor    %esi,%esi
test   $0x1,%r8b
je     186c <one_step_rk4+0x74c>
movupd (%rax,%rsi,8),%xmm1
movupd (%rcx,%rsi,8),%xmm2
mulpd  %xmm0,%xmm2
divpd  0x7bd(%rip),%xmm2
addpd  %xmm1,%xmm2
movupd %xmm2,(%rbx,%rsi,8)
cmp    %rbp,%rdx
je     1340 <one_step_rk4+0x220>
jmp    1701 <one_step_rk4+0x5e1>

LLVM IR (-O0):
define void @one_step_rk4(double* %0, double %1, double** %2, i32 %3, double* %4) {
%6 = alloca double*
%7 = alloca double
%8 = alloca double**
%9 = alloca i32
%10 = alloca double*
%11 = alloca i32
%12 = alloca double
%13 = alloca double
%14 = alloca double
store double* %0, double** %6
store double %1, double* %7
store double** %2, double*** %8
store i32 %3, i32* %9
store double* %4, double** %10
%15 = load double*, double** %10
%16 = load double, double* %15
store double %16, double* %12
%17 = load double, double* %7
%18 = load double*, double** %6
%19 = load double, double* %12
%20 = call i32 (double, double*, double, ...) bitcast (i32 (...)* @set_wieners to i32 (double, double*, double, ...)*)(double %17, double* %18, double %19)
%21 = load double, double* %12
%22 = load double*, double** %6
%23 = load double**, double*** %8
%24 = getelementptr inbounds double*, double** %23, i64 1
%25 = load double*, double** %24
%26 = load i32, i32* %9
%27 = call i32 (double, double*, double*, i32, ...) bitcast (i32 (...)* @rhs to i32 (double, double*, double*, i32, ...)*)(double %21, double* %22, double* %25, i32 %26)
store i32 0, i32* %11
br label %28
28: ; preds = %76, %5
%29 = load i32, i32* %11
%30 = load i32, i32* %9
%31 = icmp slt i32 %29, %30
br i1 %31, label %32, label %79
32: ; preds = %28
%33 = load double*, double** %6
%34 = load i32, i32* %11
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds double, double* %33, i64 %35
%37 = load double, double* %36
%38 = load double, double* %7
%39 = load double**, double*** %8
%40 = getelementptr inbounds double*, double** %39, i64 1
%41 = load double*, double** %40
%42 = load i32, i32* %11
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds double, double* %41, i64 %43
%45 = load double, double* %44
%46 = fmul double %38, %45
%47 = fdiv double %46, 6.000000e+00
%48 = fadd double %37, %47
%49 = load double**, double*** %8
%50 = getelementptr inbounds double*, double** %49, i64 0
%51 = load double*, double** %50
%52 = load i32, i32* %11
%53 = sext i32 %52 to i64
%54 = getelementptr inbounds double, double* %51, i64 %53
store double %48, double* %54
%55 = load double*, double** %6
%56 = load i32, i32* %11
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds double, double* %55, i64 %57
%59 = load double, double* %58
%60 = load double, double* %7
%61 = load double**, double*** %8
%62 = getelementptr inbounds double*, double** %61, i64 1
%63 = load double*, double** %62
%64 = load i32, i32* %11
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds double, double* %63, i64 %65
%67 = load double, double* %66
%68 = fmul double %60, %67
%69 = call double @llvm.fmuladd.f64(double %68, double 5.000000e-01, double %59)
%70 = load double**, double*** %8
%71 = getelementptr inbounds double*, double** %70, i64 2
%72 = load double*, double** %71
%73 = load i32, i32* %11
%74 = sext i32 %73 to i64
%75 = getelementptr inbounds double, double* %72, i64 %74
store double %69, double* %75
br label %76
76: ; preds = %32
%77 = load i32, i32* %11
%78 = add nsw i32 %77, 1
store i32 %78, i32* %11
br label %28
79: ; preds = %28
%80 = load double, double* %12
%81 = load double, double* %7
%82 = call double @llvm.fmuladd.f64(double 5.000000e-01, double %81, double %80)
store double %82, double* %13
%83 = load double, double* %13
%84 = load double**, double*** %8
%85 = getelementptr inbounds double*, double** %84, i64 2
%86 = load double*, double** %85
%87 = load double**, double*** %8
%88 = getelementptr inbounds double*, double** %87, i64 1
%89 = load double*, double** %88
%90 = load i32, i32* %9
%91 = call i32 (double, double*, double*, i32, ...) bitcast (i32 (...)* @rhs to i32 (double, double*, double*, i32, ...)*)(double %83, double* %86, double* %89, i32 %90)
store i32 0, i32* %11
br label %92
92: ; preds = %142, %79
%93 = load i32, i32* %11
%94 = load i32, i32* %9
%95 = icmp slt i32 %93, %94
br i1 %95, label %96, label %145
96: ; preds = %92
%97 = load double**, double*** %8
%98 = getelementptr inbounds double*, double** %97, i64 0
%99 = load double*, double** %98
%100 = load i32, i32* %11
%101 = sext i32 %100 to i64
%102 = getelementptr inbounds double, double* %99, i64 %101
%103 = load double, double* %102
%104 = load double, double* %7
%105 = load double**, double*** %8
%106 = getelementptr inbounds double*, double** %105, i64 1
%107 = load double*, double** %106
%108 = load i32, i32* %11
%109 = sext i32 %108 to i64
%110 = getelementptr inbounds double, double* %107, i64 %109
%111 = load double, double* %110
%112 = fmul double %104, %111
%113 = fdiv double %112, 3.000000e+00
%114 = fadd double %103, %113
%115 = load double**, double*** %8
%116 = getelementptr inbounds double*, double** %115, i64 0
%117 = load double*, double** %116
%118 = load i32, i32* %11
%119 = sext i32 %118 to i64
%120 = getelementptr inbounds double, double* %117, i64 %119
store double %114, double* %120
%121 = load double*, double** %6
%122 = load i32, i32* %11
%123 = sext i32 %122 to i64
%124 = getelementptr inbounds double, double* %121, i64 %123
%125 = load double, double* %124
%126 = load double, double* %7
%127 = fmul double 5.000000e-01, %126
%128 = load double**, double*** %8
%129 = getelementptr inbounds double*, double** %128, i64 1
%130 = load double*, double** %129
%131 = load i32, i32* %11
%132 = sext i32 %131 to i64
%133 = getelementptr inbounds double, double* %130, i64 %132
%134 = load double, double* %133
%135 = call double @llvm.fmuladd.f64(double %127, double %134, double %125)
%136 = load double**, double*** %8
%137 = getelementptr inbounds double*, double** %136, i64 2
%138 = load double*, double** %137
%139 = load i32, i32* %11
%140 = sext i32 %139 to i64
%141 = getelementptr inbounds double, double* %138, i64 %140
store double %135, double* %141
br label %142
142: ; preds = %96
%143 = load i32, i32* %11
%144 = add nsw i32 %143, 1
store i32 %144, i32* %11
br label %92
145: ; preds = %92
%146 = load double, double* %13
%147 = load double**, double*** %8
%148 = getelementptr inbounds double*, double** %147, i64 2
%149 = load double*, double** %148
%150 = load double**, double*** %8
%151 = getelementptr inbounds double*, double** %150, i64 1
%152 = load double*, double** %151
%153 = load i32, i32* %9
%154 = call i32 (double, double*, double*, i32, ...) bitcast (i32 (...)* @rhs to i32 (double, double*, double*, i32, ...)*)(double %146, double* %149, double* %152, i32 %153)
store i32 0, i32* %11
br label %155
155: ; preds = %204, %145
%156 = load i32, i32* %11
%157 = load i32, i32* %9
%158 = icmp slt i32 %156, %157
br i1 %158, label %159, label %207
159: ; preds = %155
%160 = load double**, double*** %8
%161 = getelementptr inbounds double*, double** %160, i64 0
%162 = load double*, double** %161
%163 = load i32, i32* %11
%164 = sext i32 %163 to i64
%165 = getelementptr inbounds double, double* %162, i64 %164
%166 = load double, double* %165
%167 = load double, double* %7
%168 = load double**, double*** %8
%169 = getelementptr inbounds double*, double** %168, i64 1
%170 = load double*, double** %169
%171 = load i32, i32* %11
%172 = sext i32 %171 to i64
%173 = getelementptr inbounds double, double* %170, i64 %172
%174 = load double, double* %173
%175 = fmul double %167, %174
%176 = fdiv double %175, 3.000000e+00
%177 = fadd double %166, %176
%178 = load double**, double*** %8
%179 = getelementptr inbounds double*, double** %178, i64 0
%180 = load double*, double** %179
%181 = load i32, i32* %11
%182 = sext i32 %181 to i64
%183 = getelementptr inbounds double, double* %180, i64 %182
store double %177, double* %183
%184 = load double*, double** %6
%185 = load i32, i32* %11
%186 = sext i32 %185 to i64
%187 = getelementptr inbounds double, double* %184, i64 %186
%188 = load double, double* %187
%189 = load double, double* %7
%190 = load double**, double*** %8
%191 = getelementptr inbounds double*, double** %190, i64 1
%192 = load double*, double** %191
%193 = load i32, i32* %11
%194 = sext i32 %193 to i64
%195 = getelementptr inbounds double, double* %192, i64 %194
%196 = load double, double* %195
%197 = call double @llvm.fmuladd.f64(double %189, double %196, double %188)
%198 = load double**, double*** %8
%199 = getelementptr inbounds double*, double** %198, i64 2
%200 = load double*, double** %199
%201 = load i32, i32* %11
%202 = sext i32 %201 to i64
%203 = getelementptr inbounds double, double* %200, i64 %202
store double %197, double* %203
br label %204
204: ; preds = %159
%205 = load i32, i32* %11
%206 = add nsw i32 %205, 1
store i32 %206, i32* %11
br label %155
207: ; preds = %155
%208 = load double, double* %12
%209 = load double, double* %7
%210 = fadd double %208, %209
store double %210, double* %14
%211 = load double, double* %14
%212 = load double**, double*** %8
%213 = getelementptr inbounds double*, double** %212, i64 2
%214 = load double*, double** %213
%215 = load double**, double*** %8
%216 = getelementptr inbounds double*, double** %215, i64 1
%217 = load double*, double** %216
%218 = load i32, i32* %9
%219 = call i32 (double, double*, double*, i32, ...) bitcast (i32 (...)* @rhs to i32 (double, double*, double*, i32, ...)*)(double %211, double* %214, double* %217, i32 %218)
store i32 0, i32* %11
br label %220
220: ; preds = %247, %207
%221 = load i32, i32* %11
%222 = load i32, i32* %9
%223 = icmp slt i32 %221, %222
br i1 %223, label %224, label %250
224: ; preds = %220
%225 = load double**, double*** %8
%226 = getelementptr inbounds double*, double** %225, i64 0
%227 = load double*, double** %226
%228 = load i32, i32* %11
%229 = sext i32 %228 to i64
%230 = getelementptr inbounds double, double* %227, i64 %229
%231 = load double, double* %230
%232 = load double, double* %7
%233 = load double**, double*** %8
%234 = getelementptr inbounds double*, double** %233, i64 1
%235 = load double*, double** %234
%236 = load i32, i32* %11
%237 = sext i32 %236 to i64
%238 = getelementptr inbounds double, double* %235, i64 %237
%239 = load double, double* %238
%240 = fmul double %232, %239
%241 = fdiv double %240, 6.000000e+00
%242 = fadd double %231, %241
%243 = load double*, double** %6
%244 = load i32, i32* %11
%245 = sext i32 %244 to i64
%246 = getelementptr inbounds double, double* %243, i64 %245
store double %242, double* %246
br label %247
247: ; preds = %224
%248 = load i32, i32* %11
%249 = add nsw i32 %248, 1
store i32 %249, i32* %11
br label %220
250: ; preds = %220
%251 = load double, double* %14
%252 = load double*, double** %10
store double %251, double* %252
ret void
}

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