void derivatives(f1, f2, n, Q, x1, x2, gradient, hessian) int n;
double f1, f2, Q, *x1, *x2;
double *gradient, **hessian;
{
  int i, j;
  double exp1, exp2, dexp, expo;
  double f, df, omega;
  double *dadp, *dbdp, *dqdp, d2qdp2;
  double tau_e, tau_s, a, b, Q_omega;
  double *dvector();
  void free_dvector();
  dadp = dvector(1, n);
  dbdp = dvector(1, n);
  dqdp = dvector(1, n);
  exp1 = log10(f1);
  exp2 = log10(f2);
  dexp = (exp2 - exp1) / 100.0;
  for (i = 1; i <= n; i++) {
    gradient[i] = 0.0;
    for (j = 1; j <= i; j++) {
      hessian[j][i] = 0.0;
      hessian[j][i] = hessian[i][j];
    }
  }
  for (expo = exp1; expo <= exp2; expo += dexp) {
    f = pow(10.0, expo);
    df = pow(10.0, expo + dexp) - f;
    omega = 6.28318530717958 * f;
    a = (double)(1 - n);
    b = 0.0;
    for (i = 1; i <= n; i++) {
      tau_e = x1[i] + x2[i];
      tau_s = x2[i];
      a += (1.0 + omega * omega * tau_e * tau_s) /
           (1.0 + omega * omega * tau_s * tau_s);
      b += omega * (tau_e - tau_s) / (1.0 + omega * omega * tau_s * tau_s);
      dadp[i] = omega * omega * tau_s / (1.0 + omega * omega * tau_s * tau_s);
      dbdp[i] = omega / (1.0 + omega * omega * tau_s * tau_s);
    }
    Q_omega = a / b;
    for (i = 1; i <= n; i++) {
      dqdp[i] = (dbdp[i] - (b / a) * dadp[i]) / a;
      gradient[i] += 2.0 * (1.0 / Q_omega - 1.0 / Q) * dqdp[i] * df / (f2 - f1);
      for (j = 1; j <= i; j++) {
        d2qdp2 = -(dadp[i] * dbdp[j] + dbdp[i] * dadp[j] -
                   2.0 * (b / a) * dadp[i] * dadp[j]) /
                 (a * a);
        hessian[i][j] += (2.0 * dqdp[i] * dqdp[j] +
                          2.0 * (1.0 / Q_omega - 1.0 / Q) * d2qdp2) *
                         df / (f2 - f1);
        hessian[j][i] = hessian[i][j];
      }
    }
  }
  free_dvector(dadp, 1, n);
  free_dvector(dbdp, 1, n);
  free_dvector(dqdp, 1, n);
}


x86 assembly (clang, -O0):
<derivatives>:
push   %rbp
mov    %rsp,%rbp
sub    $0xd0,%rsp
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
mov    %edi,-0x14(%rbp)
movsd  %xmm2,-0x20(%rbp)
mov    %rsi,-0x28(%rbp)
mov    %rdx,-0x30(%rbp)
mov    %rcx,-0x38(%rbp)
mov    %r8,-0x40(%rbp)
mov    -0x14(%rbp),%esi
mov    $0x1,%edi
mov    $0x0,%al
call   1030 <dvector@plt>
mov    %rax,-0x88(%rbp)
mov    -0x14(%rbp),%esi
mov    $0x1,%edi
mov    $0x0,%al
call   1030 <dvector@plt>
mov    %rax,-0x90(%rbp)
mov    -0x14(%rbp),%esi
mov    $0x1,%edi
mov    $0x0,%al
call   1030 <dvector@plt>
mov    %rax,-0x98(%rbp)
movsd  -0x8(%rbp),%xmm0
call   1060 <log10@plt>
movsd  %xmm0,-0x50(%rbp)
movsd  -0x10(%rbp),%xmm0
call   1060 <log10@plt>
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm0
subsd  -0x50(%rbp),%xmm0
movsd  0xe21(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movl   $0x1,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     1281 <derivatives+0x141>
mov    -0x38(%rbp),%rax
movslq -0x44(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x1,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x44(%rbp),%eax
jg     126e <derivatives+0x12e>
mov    -0x40(%rbp),%rax
movslq -0x48(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x44(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%rax
movslq -0x44(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%rbp),%rax
movslq -0x48(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    1212 <derivatives+0xd2>
jmp    1273 <derivatives+0x133>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    11ef <derivatives+0xaf>
movsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm1
movsd  -0x58(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jb     1785 <derivatives+0x645>
movsd  -0x68(%rbp),%xmm1
movsd  0xd64(%rip),%xmm0
call   1050 <pow@plt>
movsd  %xmm0,-0x70(%rbp)
movsd  -0x68(%rbp),%xmm1
addsd  -0x60(%rbp),%xmm1
movsd  0xd48(%rip),%xmm0
call   1050 <pow@plt>
subsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  0xd29(%rip),%xmm0
mulsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    $0x1,%eax
sub    -0x14(%rbp),%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
movl   $0x1,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     14b2 <derivatives+0x372>
mov    -0x28(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
movslq -0x44(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
mov    -0x30(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0x80(%rbp),%xmm0
mulsd  -0x80(%rbp),%xmm0
mulsd  -0xa8(%rbp),%xmm0
ff
movsd  -0xb0(%rbp),%xmm1
ff
mulsd  %xmm1,%xmm0
movsd  0xca8(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  -0x80(%rbp),%xmm1
mulsd  -0x80(%rbp),%xmm1
mulsd  -0xb0(%rbp),%xmm1
ff
movsd  -0xb0(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm1
movsd  0xc7e(%rip),%xmm2
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
addsd  -0xb8(%rbp),%xmm0
ff
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0x80(%rbp),%xmm0
movsd  -0xa8(%rbp),%xmm1
ff
subsd  -0xb0(%rbp),%xmm1
ff
mulsd  %xmm1,%xmm0
movsd  -0x80(%rbp),%xmm1
mulsd  -0x80(%rbp),%xmm1
mulsd  -0xb0(%rbp),%xmm1
ff
movsd  -0xb0(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm1
movsd  0xc27(%rip),%xmm2
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
addsd  -0xc0(%rbp),%xmm0
ff
movsd  %xmm0,-0xc0(%rbp)
ff
movsd  -0x80(%rbp),%xmm0
mulsd  -0x80(%rbp),%xmm0
mulsd  -0xb0(%rbp),%xmm0
ff
movsd  -0x80(%rbp),%xmm1
mulsd  -0x80(%rbp),%xmm1
mulsd  -0xb0(%rbp),%xmm1
ff
movsd  -0xb0(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm1
movsd  0xbd7(%rip),%xmm2
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x88(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x80(%rbp),%xmm0
movsd  -0x80(%rbp),%xmm1
mulsd  -0x80(%rbp),%xmm1
mulsd  -0xb0(%rbp),%xmm1
ff
movsd  -0xb0(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm1
movsd  0xb94(%rip),%xmm2
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x90(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    130f <derivatives+0x1cf>
movsd  -0xb8(%rbp),%xmm0
ff
divsd  -0xc0(%rbp),%xmm0
ff
movsd  %xmm0,-0xc8(%rbp)
ff
movl   $0x1,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     176c <derivatives+0x62c>
mov    -0x90(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0xc0(%rbp),%xmm0
ff
divsd  -0xb8(%rbp),%xmm0
ff
mov    -0x88(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
divsd  -0xb8(%rbp),%xmm0
ff
mov    -0x98(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  0xad4(%rip),%xmm1
divsd  -0xc8(%rbp),%xmm1
ff
movsd  0xac4(%rip),%xmm0
divsd  -0x20(%rbp),%xmm0
subsd  %xmm0,%xmm1
movsd  0xaab(%rip),%xmm0
mulsd  %xmm1,%xmm0
mov    -0x98(%rbp),%rax
movslq -0x44(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mulsd  -0x78(%rbp),%xmm0
movsd  -0x10(%rbp),%xmm1
subsd  -0x8(%rbp),%xmm1
divsd  %xmm1,%xmm0
mov    -0x38(%rbp),%rax
movslq -0x44(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x1,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x44(%rbp),%eax
jg     1759 <derivatives+0x619>
mov    -0x88(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x90(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x90(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x88(%rbp),%rax
movslq -0x48(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0xc0(%rbp),%xmm2
ff
divsd  -0xb8(%rbp),%xmm2
ff
movsd  0x9ff(%rip),%xmm0
mulsd  %xmm2,%xmm0
mov    -0x88(%rbp),%rax
movslq -0x44(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x88(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  -0xb8(%rbp),%xmm1
ff
mulsd  -0xb8(%rbp),%xmm1
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
mov    -0x98(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  0x976(%rip),%xmm0
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x98(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  0x961(%rip),%xmm3
divsd  -0xc8(%rbp),%xmm3
ff
movsd  0x951(%rip),%xmm1
divsd  -0x20(%rbp),%xmm1
subsd  %xmm1,%xmm3
movsd  0x938(%rip),%xmm1
mulsd  %xmm3,%xmm1
mulsd  -0xa0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mulsd  -0x78(%rbp),%xmm0
movsd  -0x10(%rbp),%xmm1
subsd  -0x8(%rbp),%xmm1
divsd  %xmm1,%xmm0
mov    -0x40(%rbp),%rax
movslq -0x44(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x48(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%rax
movslq -0x44(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%rbp),%rax
movslq -0x48(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    15ad <derivatives+0x46d>
jmp    175e <derivatives+0x61e>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    14d1 <derivatives+0x391>
jmp    1771 <derivatives+0x631>
movsd  -0x60(%rbp),%xmm0
addsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
jmp    128b <derivatives+0x14b>
mov    -0x88(%rbp),%rdi
mov    -0x14(%rbp),%edx
mov    $0x1,%esi
mov    $0x0,%al
call   1040 <free_dvector@plt>
mov    -0x90(%rbp),%rdi
mov    -0x14(%rbp),%edx
mov    $0x1,%esi
mov    $0x0,%al
call   1040 <free_dvector@plt>
mov    -0x98(%rbp),%rdi
mov    -0x14(%rbp),%edx
mov    $0x1,%esi
mov    $0x0,%al
call   1040 <free_dvector@plt>
add    $0xd0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<derivatives>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
mov    %r8,%rbx
mov    %rcx,%r15
mov    %rdx,0x58(%rsp)
mov    %rsi,0x50(%rsp)
movsd  %xmm2,0x30(%rsp)
mov    %edi,%r13d
movsd  %xmm1,0x18(%rsp)
movsd  %xmm0,0x20(%rsp)
mov    $0x1,%edi
mov    %r13d,%esi
xor    %eax,%eax
call   1030 <dvector@plt>
mov    %rax,%rbp
mov    $0x1,%edi
mov    %r13d,%esi
xor    %eax,%eax
call   1030 <dvector@plt>
mov    %rax,%r14
mov    $0x1,%edi
mov    %r13d,%esi
xor    %eax,%eax
call   1030 <dvector@plt>
mov    %rax,%r12
movsd  0x20(%rsp),%xmm0
call   1060 <log10@plt>
movsd  %xmm0,0x10(%rsp)
movsd  0x18(%rsp),%xmm0
call   1060 <log10@plt>
movsd  %xmm0,0x28(%rsp)
mov    %r13,0x8(%rsp)
test   %r13d,%r13d
jle    1295 <derivatives+0x155>
mov    0x8(%rsp),%rax
lea    0x1(%rax),%r8d
mov    $0x1,%ecx
xor    %r9d,%r9d
jmp    11fd <derivatives+0xbd>
nopl   0x0(%rax)
add    $0x1,%rcx
cmp    %r8,%rcx
je     1295 <derivatives+0x155>
mov    %r9,%rax
add    $0x1,%r9
movq   $0x0,(%r15,%rcx,8)
mov    (%rbx,%rcx,8),%rsi
mov    $0x1,%edi
test   %rax,%rax
je     1270 <derivatives+0x130>
mov    %rcx,%rax
and    $0xfffffffffffffffe,%rax
neg    %rax
mov    $0x1,%edi
nopl   0x0(%rax)
mov    (%rbx,%rdi,8),%rdx
movq   $0x0,(%rdx,%rcx,8)
movsd  (%rsi,%rdi,8),%xmm0
movsd  %xmm0,(%rdx,%rcx,8)
mov    0x8(%rbx,%rdi,8),%rdx
movq   $0x0,(%rdx,%rcx,8)
movsd  0x8(%rsi,%rdi,8),%xmm0
movsd  %xmm0,(%rdx,%rcx,8)
lea    (%rax,%rdi,1),%rdx
add    $0x2,%rdx
add    $0x2,%rdi
cmp    $0x1,%rdx
jne    1230 <derivatives+0xf0>
test   $0x1,%r9b
je     11f0 <derivatives+0xb0>
mov    (%rbx,%rdi,8),%rax
movq   $0x0,(%rax,%rcx,8)
movsd  (%rsi,%rdi,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    11f0 <derivatives+0xb0>
movsd  0x10(%rsp),%xmm2
movsd  0x28(%rsp),%xmm0
ucomisd %xmm2,%xmm0
jae    12ed <derivatives+0x1ad>
mov    %rbp,%rdi
mov    $0x1,%esi
mov    0x8(%rsp),%rbx
mov    %ebx,%edx
xor    %eax,%eax
call   1040 <free_dvector@plt>
mov    %r14,%rdi
mov    $0x1,%esi
mov    %ebx,%edx
xor    %eax,%eax
call   1040 <free_dvector@plt>
mov    %r12,%rdi
mov    $0x1,%esi
mov    %ebx,%edx
xor    %eax,%eax
add    $0x68,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1040 <free_dvector@plt>
movsd  0x28(%rsp),%xmm0
subsd  %xmm2,%xmm0
divsd  0xd01(%rip),%xmm0
movsd  %xmm0,0x60(%rsp)
mov    $0x1,%eax
mov    0x8(%rsp),%rcx
sub    %ecx,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0xce8(%rip),%xmm4
divsd  0x30(%rsp),%xmm4
movsd  %xmm4,0x48(%rsp)
movsd  0x18(%rsp),%xmm1
subsd  0x20(%rsp),%xmm1
movsd  %xmm1,0x18(%rsp)
lea    0x1(%rcx),%r13d
xorpd  %xmm1,%xmm1
unpcklpd %xmm0,%xmm1
movapd %xmm1,0x30(%rsp)
jmp    1376 <derivatives+0x236>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  0x10(%rsp),%xmm2
movsd  0x28(%rsp),%xmm0
ucomisd %xmm2,%xmm0
jb     12a7 <derivatives+0x167>
movsd  %xmm2,0x10(%rsp)
movsd  0xc8c(%rip),%xmm0
movapd %xmm2,%xmm1
call   1050 <pow@plt>
movsd  %xmm0,0x20(%rsp)
movsd  0x10(%rsp),%xmm1
addsd  0x60(%rsp),%xmm1
movsd  0xc69(%rip),%xmm0
movsd  %xmm1,0x10(%rsp)
call   1050 <pow@plt>
movsd  0x20(%rsp),%xmm1
movsd  0xc47(%rip),%xmm11
movapd %xmm1,%xmm8
mulsd  0xc49(%rip),%xmm8
cmpl   $0x0,0x8(%rsp)
movsd  0x18(%rsp),%xmm10
jle    1360 <derivatives+0x220>
subsd  %xmm1,%xmm0
movapd %xmm8,%xmm1
mulsd  %xmm8,%xmm1
mov    $0x1,%eax
movapd 0x30(%rsp),%xmm2
mov    0x58(%rsp),%rcx
mov    0x50(%rsp),%rdx
movsd  (%rcx,%rax,8),%xmm5
movsd  (%rdx,%rax,8),%xmm4
addsd  %xmm5,%xmm4
movapd %xmm1,%xmm6
mulsd  %xmm4,%xmm6
mulsd  %xmm5,%xmm6
addsd  %xmm11,%xmm6
movapd %xmm1,%xmm3
mulsd  %xmm5,%xmm3
movapd %xmm3,%xmm7
mulsd  %xmm5,%xmm7
addsd  %xmm11,%xmm7
subsd  %xmm5,%xmm4
mulsd  %xmm8,%xmm4
unpcklpd %xmm6,%xmm4
divsd  %xmm7,%xmm3
movapd %xmm8,%xmm5
divsd  %xmm7,%xmm5
unpcklpd %xmm7,%xmm7
divpd  %xmm7,%xmm4
addpd  %xmm4,%xmm2
movsd  %xmm3,0x0(%rbp,%rax,8)
movsd  %xmm5,(%r14,%rax,8)
add    $0x1,%rax
cmp    %rax,%r13
jne    1400 <derivatives+0x2c0>
cmpl   $0x0,0x8(%rsp)
jle    1360 <derivatives+0x220>
movapd %xmm2,%xmm9
unpckhpd %xmm2,%xmm9
movapd %xmm9,%xmm5
divsd  %xmm2,%xmm5
divsd  %xmm9,%xmm2
movapd %xmm2,%xmm8
xorpd  0xb90(%rip),%xmm8
movapd %xmm11,%xmm4
divsd  %xmm5,%xmm4
subsd  0x48(%rsp),%xmm4
addsd  %xmm4,%xmm4
mulsd  0xb65(%rip),%xmm2
movapd %xmm9,%xmm5
mulsd  %xmm9,%xmm5
mov    $0x1,%eax
mov    $0x8,%ecx
nop
movsd  0x0(%rbp,%rax,8),%xmm6
mulsd  %xmm8,%xmm6
addsd  (%r14,%rax,8),%xmm6
divsd  %xmm9,%xmm6
movsd  %xmm6,(%r12,%rax,8)
mulsd  %xmm4,%xmm6
mulsd  %xmm0,%xmm6
divsd  %xmm10,%xmm6
addsd  (%r15,%rax,8),%xmm6
movsd  %xmm6,(%r15,%rax,8)
mov    (%rbx,%rax,8),%rdx
xor    %esi,%esi
nopl   0x0(%rax,%rax,1)
movsd  0x0(%rbp,%rax,8),%xmm6
movsd  0x8(%rbp,%rsi,1),%xmm7
movsd  (%r14,%rax,8),%xmm3
movapd %xmm2,%xmm1
mulsd  %xmm6,%xmm1
mulsd  0x8(%r14,%rsi,1),%xmm6
mulsd  %xmm7,%xmm3
addsd  %xmm3,%xmm6
mulsd  %xmm7,%xmm1
addsd  %xmm6,%xmm1
divsd  %xmm5,%xmm1
movsd  (%r12,%rax,8),%xmm3
addsd  %xmm3,%xmm3
mulsd  %xmm4,%xmm1
mulsd  0x8(%r12,%rsi,1),%xmm3
subsd  %xmm1,%xmm3
mulsd  %xmm0,%xmm3
divsd  %xmm10,%xmm3
addsd  0x8(%rdx,%rsi,1),%xmm3
movsd  %xmm3,0x8(%rdx,%rsi,1)
mov    0x8(%rbx,%rsi,1),%rdi
movsd  %xmm3,(%rdi,%rax,8)
add    $0x8,%rsi
cmp    %rsi,%rcx
jne    1510 <derivatives+0x3d0>
add    $0x1,%rax
add    $0x8,%rcx
cmp    %r13,%rax
jne    14d0 <derivatives+0x390>
jmp    1360 <derivatives+0x220>

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

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