void jacobn(x, y, dfdx, dermat, eps, work, n) double x, *y, *dermat, *dfdx, eps,
    *work;
int n;
{
  int i, j;
  double r;
  double *yval, *ynew, ytemp;
  yval = work;
  ynew = work + n;
  rhs(x, y, yval, n);
  r = eps * ((eps) > (fabs(x)) ? (eps) : (fabs(x)));
  rhs(x + r, y, ynew, n);
  for (i = 0; i < n; i++) {
    dfdx[i] = (ynew[i] - yval[i]) / r;
  }
  for (i = 0; i < n; i++) {
    ytemp = y[i];
    r = eps * ((eps) > (fabs(ytemp)) ? (eps) : (fabs(ytemp)));
    y[i] = ytemp + r;
    rhs(x, y, ynew, n);
    for (j = 0; j < n; j++) {
      dermat[j * n + i] = (ynew[j] - yval[j]) / r;
    }
    y[i] = ytemp;
  }
}


x86 assembly (clang, -O0):
<jacobn>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
movsd  %xmm0,-0x8(%rbp)
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
movsd  %xmm1,-0x28(%rbp)
mov    %rcx,-0x30(%rbp)
mov    %r8d,-0x34(%rbp)
mov    -0x30(%rbp),%rax
mov    %rax,-0x50(%rbp)
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x58(%rbp)
movsd  -0x8(%rbp),%xmm0
mov    -0x10(%rbp),%rdi
mov    -0x50(%rbp),%rsi
mov    -0x34(%rbp),%edx
mov    $0x1,%al
call   1030 <rhs@plt>
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x28(%rbp),%xmm0
movsd  -0x8(%rbp),%xmm1
movaps 0xe7a(%rip),%xmm2
pand   %xmm2,%xmm1
ucomisd %xmm1,%xmm0
jbe    11a3 <jacobn+0x93>
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
jmp    11b8 <jacobn+0xa8>
movsd  -0x8(%rbp),%xmm0
movaps 0xe51(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x68(%rbp),%xmm0
movsd  -0x70(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x8(%rbp),%xmm0
addsd  -0x48(%rbp),%xmm0
mov    -0x10(%rbp),%rdi
mov    -0x58(%rbp),%rsi
mov    -0x34(%rbp),%edx
mov    $0x1,%al
call   1030 <rhs@plt>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1234 <jacobn+0x124>
mov    -0x58(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x50(%rbp),%rax
movslq -0x38(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
divsd  -0x48(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    11ee <jacobn+0xde>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    135d <jacobn+0x24d>
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x28(%rbp),%xmm0
movsd  -0x60(%rbp),%xmm1
movaps 0xd8c(%rip),%xmm2
pand   %xmm2,%xmm1
ucomisd %xmm1,%xmm0
jbe    1291 <jacobn+0x181>
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
jmp    12a6 <jacobn+0x196>
movsd  -0x60(%rbp),%xmm0
movaps 0xd63(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x78(%rbp),%xmm0
movsd  -0x80(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x60(%rbp),%xmm0
addsd  -0x48(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x8(%rbp),%xmm0
mov    -0x10(%rbp),%rdi
mov    -0x58(%rbp),%rsi
mov    -0x34(%rbp),%edx
mov    $0x1,%al
call   1030 <rhs@plt>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    133d <jacobn+0x22d>
mov    -0x58(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x50(%rbp),%rax
movslq -0x3c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
divsd  -0x48(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0x34(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    12ee <jacobn+0x1de>
movsd  -0x60(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    123b <jacobn+0x12b>
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<jacobn>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x78,%rsp
mov    %rcx,%rbx
movsd  %xmm1,0x38(%rsp)
mov    %rdx,0x20(%rsp)
mov    %rsi,%r15
mov    %rdi,%r13
movaps %xmm0,0x50(%rsp)
mov    %r8d,0x4(%rsp)
movslq %r8d,%rbp
lea    (%rcx,%rbp,8),%r14
mov    %rcx,%rsi
mov    %ebp,%edx
mov    $0x1,%al
call   1030 <rhs@plt>
movapd 0xea9(%rip),%xmm0
movapd 0x50(%rsp),%xmm3
andpd  %xmm3,%xmm0
movsd  0x38(%rsp),%xmm2
movapd %xmm2,%xmm1
maxsd  %xmm0,%xmm1
mulsd  %xmm2,%xmm1
movapd %xmm1,0x10(%rsp)
movapd %xmm1,%xmm0
addsd  %xmm3,%xmm0
mov    %r13,%rdi
mov    %r14,0x8(%rsp)
mov    %r14,%rsi
mov    %ebp,%edx
mov    $0x1,%al
call   1030 <rhs@plt>
mov    %rbp,%rdi
test   %edi,%edi
jle    149d <jacobn+0x38d>
mov    0x4(%rsp),%eax
mov    %eax,%r9d
cmp    $0x1,%eax
jne    1238 <jacobn+0x128>
xor    %eax,%eax
movapd 0x10(%rsp),%xmm3
mov    %rax,%rcx
not    %rcx
test   $0x1,%r9b
je     11e1 <jacobn+0xd1>
mov    0x8(%rsp),%rdx
movsd  (%rdx,%rax,8),%xmm0
subsd  (%rbx,%rax,8),%xmm0
divsd  %xmm3,%xmm0
movsd  %xmm0,(%r15,%rax,8)
or     $0x1,%rax
add    %r9,%rcx
je     1329 <jacobn+0x219>
mov    %r9,%rcx
sub    %rax,%rcx
lea    (%r15,%rax,8),%rdx
add    $0x8,%rdx
lea    (%rbx,%rax,8),%rax
add    $0x8,%rax
movsd  -0x8(%rax,%rdi,8),%xmm0
subsd  -0x8(%rax),%xmm0
divsd  %xmm3,%xmm0
movsd  %xmm0,-0x8(%rdx)
movsd  (%rax,%rdi,8),%xmm0
subsd  (%rax),%xmm0
divsd  %xmm3,%xmm0
movsd  %xmm0,(%rdx)
add    $0x10,%rdx
add    $0x10,%rax
add    $0xfffffffffffffffe,%rcx
jne    1200 <jacobn+0xf0>
jmp    1329 <jacobn+0x219>
lea    (%r15,%r9,8),%rax
lea    (%rbx,%r9,8),%rcx
cmp    %r15,%rcx
seta   %r8b
lea    (%rdi,%r9,1),%rcx
cmp    %rbx,%rax
seta   %dl
lea    (%rbx,%rcx,8),%rcx
cmp    %r15,%rcx
seta   %cl
cmp    %rax,0x8(%rsp)
setb   %sil
xor    %eax,%eax
test   %dl,%r8b
jne    11b2 <jacobn+0xa2>
and    %sil,%cl
movapd 0x10(%rsp),%xmm3
jne    11b8 <jacobn+0xa8>
mov    %r9d,%eax
and    $0xfffffffe,%eax
movapd %xmm3,%xmm0
unpcklpd %xmm3,%xmm0
lea    -0x2(%rax),%rdx
mov    %rdx,%rcx
shr    %rcx
add    $0x1,%rcx
test   %rdx,%rdx
je     14ac <jacobn+0x39c>
mov    %rcx,%rsi
and    $0xfffffffffffffffe,%rsi
lea    (%rbx,%rdi,8),%rbp
add    $0x10,%rbp
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movupd -0x10(%rbp,%rdx,8),%xmm1
movupd (%rbx,%rdx,8),%xmm2
subpd  %xmm2,%xmm1
divpd  %xmm0,%xmm1
movupd %xmm1,(%r15,%rdx,8)
movupd 0x0(%rbp,%rdx,8),%xmm1
movupd 0x10(%rbx,%rdx,8),%xmm2
subpd  %xmm2,%xmm1
divpd  %xmm0,%xmm1
movupd %xmm1,0x10(%r15,%rdx,8)
add    $0x4,%rdx
add    $0xfffffffffffffffe,%rsi
jne    12c0 <jacobn+0x1b0>
test   $0x1,%cl
je     1320 <jacobn+0x210>
mov    0x8(%rsp),%rcx
movupd (%rcx,%rdx,8),%xmm1
movupd (%rbx,%rdx,8),%xmm2
subpd  %xmm2,%xmm1
divpd  %xmm0,%xmm1
movupd %xmm1,(%r15,%rdx,8)
cmp    %r9,%rax
jne    11b8 <jacobn+0xa8>
cmpl   $0x0,0x4(%rsp)
jle    149d <jacobn+0x38d>
mov    %r9d,%r14d
and    $0xfffffffe,%r14d
lea    (%rbx,%rdi,8),%r15
add    $0x8,%r15
mov    %rdi,%rbp
shl    $0x4,%rbp
xor    %r12d,%r12d
mov    0x20(%rsp),%rdx
mov    %r13,0x48(%rsp)
mov    %rdi,0x30(%rsp)
mov    %r9,0x28(%rsp)
jmp    138e <jacobn+0x27e>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movapd 0x60(%rsp),%xmm0
movsd  %xmm0,0x0(%r13,%r12,8)
add    $0x1,%r12
add    $0x8,%rdx
cmp    %r9,%r12
je     149d <jacobn+0x38d>
mov    %rdx,0x10(%rsp)
movsd  0x0(%r13,%r12,8),%xmm2
movapd %xmm2,%xmm0
andpd  0xc5a(%rip),%xmm0
movsd  0x38(%rsp),%xmm1
movapd %xmm1,%xmm3
maxsd  %xmm0,%xmm3
mulsd  %xmm1,%xmm3
movapd %xmm2,0x60(%rsp)
movapd %xmm2,%xmm0
movsd  %xmm3,0x40(%rsp)
addsd  %xmm3,%xmm0
movsd  %xmm0,0x0(%r13,%r12,8)
movapd 0x50(%rsp),%xmm0
mov    %r13,%rdi
mov    0x8(%rsp),%rsi
mov    0x4(%rsp),%r13d
mov    %r13d,%edx
mov    $0x1,%al
call   1030 <rhs@plt>
cmp    $0x1,%r13d
jne    1410 <jacobn+0x300>
xor    %eax,%eax
mov    0x30(%rsp),%rdi
mov    0x28(%rsp),%r9
mov    0x10(%rsp),%rdx
movsd  0x40(%rsp),%xmm1
jmp    1465 <jacobn+0x355>
nop
mov    0x10(%rsp),%rdx
mov    %rdx,%rcx
xor    %eax,%eax
mov    0x30(%rsp),%rdi
mov    0x28(%rsp),%r9
movsd  0x40(%rsp),%xmm1
nopw   0x0(%rax,%rax,1)
movsd  -0x8(%r15,%rax,8),%xmm0
subsd  (%rbx,%rax,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rcx)
movsd  (%r15,%rax,8),%xmm0
subsd  0x8(%rbx,%rax,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rcx,%rdi,8)
add    $0x2,%rax
add    %rbp,%rcx
cmp    %rax,%r14
jne    1430 <jacobn+0x320>
test   $0x1,%r9b
mov    0x48(%rsp),%r13
je     1370 <jacobn+0x260>
mov    0x8(%rsp),%rcx
movsd  (%rcx,%rax,8),%xmm0
subsd  (%rbx,%rax,8),%xmm0
divsd  %xmm1,%xmm0
imul   %rdi,%rax
add    %r12,%rax
mov    0x20(%rsp),%rcx
movsd  %xmm0,(%rcx,%rax,8)
jmp    1370 <jacobn+0x260>
add    $0x78,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %edx,%edx
test   $0x1,%cl
jne    1303 <jacobn+0x1f3>
jmp    1320 <jacobn+0x210>

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

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