void leesuto_gradient(double *pars, double *r, double *grad) {
  double x, y, z, _r, _r2, _r4, ax, ay, az;
  double v_h2, x0, x2, x22;
  double x20, x21, x7, x1;
  double x10, x13, x15, x16, x17;
  double e_b2 = 1 - pow(pars[3] / pars[2], 2);
  double e_c2 = 1 - pow(pars[4] / pars[2], 2);
  v_h2 = pars[0] * pars[0] /
         (log(2.) - 0.5 + (log(2.) - 0.75) * e_b2 + (log(2.) - 0.75) * e_c2);
  x = pars[5] * r[0] + pars[6] * r[1] + pars[7] * r[2];
  y = pars[8] * r[0] + pars[9] * r[1] + pars[10] * r[2];
  z = pars[11] * r[0] + pars[12] * r[1] + pars[13] * r[2];
  _r2 = x * x + y * y + z * z;
  _r = sqrt(_r2);
  _r4 = _r2 * _r2;
  x0 = _r + pars[1];
  x1 = x0 * x0;
  x2 = v_h2 / (12. * _r4 * _r2 * _r * x1);
  x10 = log(x0 / pars[1]);
  x13 = _r * 3. * pars[1];
  x15 = x13 - _r2;
  x16 = x15 + 6. * (pars[1] * pars[1]);
  x17 = 6. * pars[1] * x0 * (_r * x16 - x0 * x10 * 6. * (pars[1] * pars[1]));
  x20 = x0 * _r2;
  x21 = 2. * _r * x0;
  x7 = e_b2 * y * y + e_c2 * z * z;
  x22 = -12. * _r4 * _r * pars[1] * x0 + 12. * _r4 * pars[1] * x1 * x10 +
        3. * pars[1] * x7 *
            (x16 * _r2 - 18. * x1 * x10 * (pars[1] * pars[1]) +
             x20 * (2. * _r - 3. * pars[1]) +
             x21 * (x15 + 9. * (pars[1] * pars[1]))) -
        x20 * (e_b2 + e_c2) *
            (-6. * _r * pars[1] * (_r2 - (pars[1] * pars[1])) +
             6. * pars[1] * x0 * x10 * (_r2 - 3. * (pars[1] * pars[1])) +
             x20 * (-4. * _r + 3. * pars[1]) +
             x21 * (-x13 + 2. * _r2 + 6. * (pars[1] * pars[1])));
  ax = x2 * x * (x17 * x7 + x22);
  ay = x2 * y * (x17 * (x7 - _r2 * e_b2) + x22);
  az = x2 * z * (x17 * (x7 - _r2 * e_c2) + x22);
  grad[0] = pars[5] * ax + pars[8] * ay + pars[11] * az;
  grad[1] = pars[6] * ax + pars[9] * ay + pars[12] * az;
  grad[2] = pars[7] * ax + pars[10] * ay + pars[13] * az;
}


x86 assembly (clang, -O0):
<leesuto_gradient>:
push   %rbp
mov    %rsp,%rbp
sub    $0xf0,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x18(%rax),%xmm0
mov    -0x8(%rbp),%rax
divsd  0x10(%rax),%xmm0
movsd  0xee7(%rip),%xmm1
call   1040 <pow@plt>
movaps %xmm0,%xmm1
movsd  0xeef(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0xd0(%rbp)
ff
mov    -0x8(%rbp),%rax
movsd  0x20(%rax),%xmm0
mov    -0x8(%rbp),%rax
divsd  0x10(%rax),%xmm0
movsd  0xeb1(%rip),%xmm1
call   1040 <pow@plt>
movaps %xmm0,%xmm1
movsd  0xeb9(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x8(%rbp),%rax
mulsd  (%rax),%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
movsd  0xe75(%rip),%xmm0
call   1030 <log@plt>
movsd  0xe78(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0xf0(%rbp)
ff
movsd  0xe54(%rip),%xmm0
call   1030 <log@plt>
movsd  -0xf0(%rbp),%xmm1
movsd  0xe47(%rip),%xmm2
subsd  %xmm2,%xmm0
movsd  -0xd0(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xe8(%rbp)
ff
movsd  0xe1b(%rip),%xmm0
call   1030 <log@plt>
movsd  -0xe8(%rbp),%xmm2
ff
movaps %xmm0,%xmm1
movsd  -0xe0(%rbp),%xmm0
ff
movsd  0xe03(%rip),%xmm3
subsd  %xmm3,%xmm1
movsd  -0xd8(%rbp),%xmm3
ff
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x28(%rax),%xmm1
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm2
mov    -0x8(%rbp),%rax
movsd  0x30(%rax),%xmm0
mov    -0x10(%rbp),%rax
mulsd  0x8(%rax),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x8(%rbp),%rax
movsd  0x38(%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  0x10(%rax),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x40(%rax),%xmm1
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm2
mov    -0x8(%rbp),%rax
movsd  0x48(%rax),%xmm0
mov    -0x10(%rbp),%rax
mulsd  0x8(%rax),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x8(%rbp),%rax
movsd  0x50(%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  0x10(%rax),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x58(%rax),%xmm1
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm2
mov    -0x8(%rbp),%rax
movsd  0x60(%rax),%xmm0
mov    -0x10(%rbp),%rax
mulsd  0x8(%rax),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x8(%rbp),%rax
movsd  0x68(%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  0x10(%rax),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x20(%rbp),%xmm1
movsd  -0x20(%rbp),%xmm2
movsd  -0x28(%rbp),%xmm0
mulsd  -0x28(%rbp),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x30(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
call   1050 <sqrt@plt>
movsd  %xmm0,-0x38(%rbp)
movsd  -0x40(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
addsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x70(%rbp),%xmm0
mulsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
movsd  0xc6d(%rip),%xmm1
mulsd  -0x48(%rbp),%xmm1
mulsd  -0x40(%rbp),%xmm1
mulsd  -0x38(%rbp),%xmm1
mulsd  -0xa0(%rbp),%xmm1
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x70(%rbp),%xmm0
mov    -0x8(%rbp),%rax
divsd  0x8(%rax),%xmm0
call   1030 <log@plt>
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  0xc02(%rip),%xmm0
mulsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mulsd  0x8(%rax),%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0xb0(%rbp),%xmm0
ff
subsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0xb8(%rbp),%xmm1
ff
mov    -0x8(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x8(%rbp),%rax
mulsd  0x8(%rax),%xmm0
movsd  0xba5(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
mov    -0x8(%rbp),%rax
movsd  0xb89(%rip),%xmm0
mulsd  0x8(%rax),%xmm0
mulsd  -0x70(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm1
movsd  -0xc0(%rbp),%xmm3
ff
movsd  -0x70(%rbp),%xmm2
mulsd  -0xa8(%rbp),%xmm2
ff
movsd  0xb5d(%rip),%xmm4
mulsd  %xmm4,%xmm2
mov    -0x8(%rbp),%rax
movsd  0x8(%rax),%xmm4
mov    -0x8(%rbp),%rax
mulsd  0x8(%rax),%xmm4
mulsd  %xmm4,%xmm2
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  0xb46(%rip),%xmm0
mulsd  -0x38(%rbp),%xmm0
mulsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0xd0(%rbp),%xmm0
ff
mulsd  -0x28(%rbp),%xmm0
movsd  -0x28(%rbp),%xmm2
movsd  -0xd8(%rbp),%xmm1
ff
mulsd  -0x30(%rbp),%xmm1
mulsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  0xaf0(%rip),%xmm0
mulsd  -0x48(%rbp),%xmm0
mulsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mulsd  0x8(%rax),%xmm0
movsd  -0x70(%rbp),%xmm2
movsd  0xac8(%rip),%xmm1
mulsd  -0x48(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mulsd  0x8(%rax),%xmm1
mulsd  -0xa0(%rbp),%xmm1
ff
mulsd  -0xa8(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movsd  0xa6e(%rip),%xmm1
mulsd  0x8(%rax),%xmm1
mulsd  -0x98(%rbp),%xmm1
ff
movsd  -0xc0(%rbp),%xmm2
ff
movsd  -0x40(%rbp),%xmm4
movsd  0xa6c(%rip),%xmm3
mulsd  -0xa0(%rbp),%xmm3
ff
mulsd  -0xa8(%rbp),%xmm3
ff
mov    -0x8(%rbp),%rax
movsd  0x8(%rax),%xmm5
mov    -0x8(%rbp),%rax
mulsd  0x8(%rax),%xmm5
mulsd  %xmm5,%xmm3
movq   %xmm3,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm3
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
movsd  -0x88(%rbp),%xmm3
ff
movsd  -0x38(%rbp),%xmm4
mov    -0x8(%rbp),%rax
movsd  0x9ee(%rip),%xmm5
mulsd  0x8(%rax),%xmm5
movq   %xmm5,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm5
addsd  %xmm4,%xmm4
addsd  %xmm5,%xmm4
mulsd  %xmm4,%xmm3
addsd  %xmm2,%xmm3
movsd  -0x90(%rbp),%xmm2
ff
movsd  -0xb8(%rbp),%xmm5
ff
mov    -0x8(%rbp),%rax
movsd  0x8(%rax),%xmm4
mov    -0x8(%rbp),%rax
mulsd  0x8(%rax),%xmm4
movsd  0x9b0(%rip),%xmm6
mulsd  %xmm6,%xmm4
addsd  %xmm5,%xmm4
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x88(%rbp),%xmm0
ff
movsd  -0xd0(%rbp),%xmm2
ff
addsd  -0xd8(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
movsd  0x96c(%rip),%xmm2
mulsd  -0x38(%rbp),%xmm2
mov    -0x8(%rbp),%rax
mulsd  0x8(%rax),%xmm2
movsd  -0x40(%rbp),%xmm3
mov    -0x8(%rbp),%rax
movsd  0x8(%rax),%xmm4
mov    -0x8(%rbp),%rax
movsd  0x8(%rax),%xmm5
movq   %xmm4,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm4
mulsd  %xmm5,%xmm4
addsd  %xmm3,%xmm4
mov    -0x8(%rbp),%rax
movsd  0x8fc(%rip),%xmm3
mulsd  0x8(%rax),%xmm3
mulsd  -0x70(%rbp),%xmm3
mulsd  -0xa8(%rbp),%xmm3
ff
movsd  -0x40(%rbp),%xmm6
mov    -0x8(%rbp),%rax
movsd  0x8(%rax),%xmm5
mov    -0x8(%rbp),%rax
mulsd  0x8(%rax),%xmm5
movsd  0x8e3(%rip),%xmm7
mulsd  %xmm7,%xmm5
addsd  %xmm6,%xmm5
mulsd  %xmm5,%xmm3
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
movsd  -0x88(%rbp),%xmm3
ff
movsd  -0x38(%rbp),%xmm4
mov    -0x8(%rbp),%rax
movsd  0x8ae(%rip),%xmm5
mulsd  0x8(%rax),%xmm5
movsd  0x899(%rip),%xmm6
mulsd  %xmm6,%xmm4
addsd  %xmm5,%xmm4
mulsd  %xmm4,%xmm3
addsd  %xmm2,%xmm3
movsd  -0x90(%rbp),%xmm2
ff
movsd  -0xb0(%rbp),%xmm4
ff
movq   %xmm4,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm4
movsd  -0x40(%rbp),%xmm5
addsd  %xmm5,%xmm5
addsd  %xmm4,%xmm5
mov    -0x8(%rbp),%rax
movsd  0x8(%rax),%xmm4
mov    -0x8(%rbp),%rax
mulsd  0x8(%rax),%xmm4
movsd  0x833(%rip),%xmm6
mulsd  %xmm6,%xmm4
addsd  %xmm5,%xmm4
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x78(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  -0xc8(%rbp),%xmm1
ff
movsd  -0x98(%rbp),%xmm3
ff
movsd  -0x80(%rbp),%xmm2
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x78(%rbp),%xmm0
mulsd  -0x28(%rbp),%xmm0
movsd  -0xc8(%rbp),%xmm1
ff
movsd  -0x98(%rbp),%xmm2
ff
movsd  -0x40(%rbp),%xmm3
movsd  -0xd0(%rbp),%xmm4
ff
movq   %xmm3,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm3
mulsd  %xmm4,%xmm3
addsd  %xmm2,%xmm3
movsd  -0x80(%rbp),%xmm2
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x78(%rbp),%xmm0
mulsd  -0x30(%rbp),%xmm0
movsd  -0xc8(%rbp),%xmm1
ff
movsd  -0x98(%rbp),%xmm2
ff
movsd  -0x40(%rbp),%xmm3
movsd  -0xd8(%rbp),%xmm4
ff
movq   %xmm3,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm3
mulsd  %xmm4,%xmm3
addsd  %xmm2,%xmm3
movsd  -0x80(%rbp),%xmm2
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x28(%rax),%xmm1
movsd  -0x50(%rbp),%xmm2
mov    -0x8(%rbp),%rax
movsd  0x40(%rax),%xmm0
mulsd  -0x58(%rbp),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x8(%rbp),%rax
movsd  0x58(%rax),%xmm0
movsd  -0x60(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x8(%rbp),%rax
movsd  0x30(%rax),%xmm1
movsd  -0x50(%rbp),%xmm2
mov    -0x8(%rbp),%rax
movsd  0x48(%rax),%xmm0
mulsd  -0x58(%rbp),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x8(%rbp),%rax
movsd  0x60(%rax),%xmm0
movsd  -0x60(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,0x8(%rax)
mov    -0x8(%rbp),%rax
movsd  0x38(%rax),%xmm1
movsd  -0x50(%rbp),%xmm2
mov    -0x8(%rbp),%rax
movsd  0x50(%rax),%xmm0
mulsd  -0x58(%rbp),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x8(%rbp),%rax
movsd  0x68(%rax),%xmm0
movsd  -0x60(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,0x10(%rax)
add    $0xf0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<leesuto_gradient>:
push   %r14
push   %rbx
sub    $0x78,%rsp
mov    %rdx,%r14
mov    %rdi,%rbx
movsd  0x18(%rdi),%xmm0
movsd  (%rdi),%xmm2
movsd  0x10(%rdi),%xmm1
divsd  %xmm1,%xmm0
mulsd  %xmm0,%xmm0
movsd  0xeb5(%rip),%xmm3
movapd %xmm3,%xmm4
subsd  %xmm0,%xmm4
movsd  0x20(%rdi),%xmm0
divsd  %xmm1,%xmm0
mulsd  %xmm0,%xmm0
subsd  %xmm0,%xmm3
mulsd  %xmm2,%xmm2
movsd  0xe98(%rip),%xmm0
movsd  %xmm4,0x70(%rsp)
movapd %xmm4,%xmm1
mulsd  %xmm0,%xmm1
addsd  0xe8a(%rip),%xmm1
movsd  %xmm3,0x28(%rsp)
mulsd  %xmm3,%xmm0
addsd  %xmm1,%xmm0
divsd  %xmm0,%xmm2
movsd  %xmm2,0x18(%rsp)
movsd  (%rsi),%xmm0
movsd  0x8(%rsi),%xmm1
movsd  0x30(%rdi),%xmm2
mulsd  %xmm1,%xmm2
movsd  0x28(%rdi),%xmm3
mulsd  %xmm0,%xmm3
addsd  %xmm2,%xmm3
movsd  0x10(%rsi),%xmm4
movsd  0x38(%rdi),%xmm5
mulsd  %xmm4,%xmm5
addsd  %xmm3,%xmm5
movsd  0x48(%rdi),%xmm2
mulsd  %xmm1,%xmm2
movsd  0x40(%rdi),%xmm3
mulsd  %xmm0,%xmm3
addsd  %xmm2,%xmm3
movsd  0x50(%rdi),%xmm2
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
mulsd  0x60(%rdi),%xmm1
mulsd  0x58(%rdi),%xmm0
addsd  %xmm1,%xmm0
mulsd  0x68(%rdi),%xmm4
addsd  %xmm0,%xmm4
movsd  %xmm2,0x60(%rsp)
movapd %xmm2,%xmm0
mulsd  %xmm2,%xmm0
movsd  %xmm5,0x68(%rsp)
movapd %xmm5,%xmm1
mulsd  %xmm5,%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm4,0x20(%rsp)
mulsd  %xmm4,%xmm4
addsd  %xmm1,%xmm4
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm4
movsd  %xmm4,0x8(%rsp)
jb     1250 <leesuto_gradient+0x130>
xorps  %xmm0,%xmm0
sqrtsd %xmm4,%xmm0
jmp    125f <leesuto_gradient+0x13f>
movapd %xmm4,%xmm0
call   1040 <sqrt@plt>
movsd  0x8(%rsp),%xmm4
movapd %xmm4,%xmm3
mulsd  %xmm4,%xmm3
movsd  %xmm3,0x10(%rsp)
movsd  0x8(%rbx),%xmm1
movapd %xmm0,%xmm2
addsd  %xmm1,%xmm2
movapd %xmm2,%xmm6
movapd %xmm2,%xmm7
movsd  %xmm2,0x48(%rsp)
mulsd  %xmm2,%xmm6
movsd  %xmm6,0x50(%rsp)
movsd  0xd7e(%rip),%xmm5
mulsd  %xmm3,%xmm5
movsd  %xmm5,0x58(%rsp)
movapd %xmm0,%xmm2
movsd  %xmm0,0x38(%rsp)
movapd %xmm4,%xmm0
mulsd  %xmm5,%xmm0
mulsd  %xmm2,%xmm0
mulsd  %xmm6,%xmm0
movsd  0x18(%rsp),%xmm2
divsd  %xmm0,%xmm2
movsd  %xmm2,0x18(%rsp)
movapd %xmm7,%xmm0
divsd  %xmm1,%xmm0
call   1030 <log@plt>
movapd %xmm0,%xmm8
movsd  0x38(%rsp),%xmm12
movapd %xmm12,%xmm3
mulsd  0xd2c(%rip),%xmm3
movsd  0x8(%rbx),%xmm0
mulsd  %xmm0,%xmm3
movapd %xmm0,%xmm1
movsd  0x8(%rsp),%xmm0
movapd %xmm0,%xmm5
addsd  %xmm0,%xmm5
subsd  %xmm3,%xmm5
movapd %xmm3,%xmm10
subsd  %xmm0,%xmm10
movapd %xmm1,%xmm0
movapd %xmm1,%xmm6
mulsd  %xmm1,%xmm6
movsd  0xcf7(%rip),%xmm7
movapd %xmm6,%xmm1
mulsd  %xmm7,%xmm1
addsd  %xmm1,%xmm5
addsd  %xmm10,%xmm1
movsd  0x48(%rsp),%xmm9
movapd %xmm9,%xmm3
movsd  %xmm8,0x40(%rsp)
mulsd  %xmm8,%xmm3
movsd  0xccd(%rip),%xmm13
mulsd  %xmm13,%xmm3
mulsd  %xmm6,%xmm3
movapd %xmm12,%xmm11
mulsd  %xmm1,%xmm11
addsd  %xmm3,%xmm11
movsd  0x70(%rsp),%xmm15
movapd %xmm15,%xmm3
movsd  0x60(%rsp),%xmm14
mulsd  %xmm14,%xmm3
movsd  0x28(%rsp),%xmm4
movsd  0x20(%rsp),%xmm2
mulsd  %xmm2,%xmm4
mulsd  %xmm2,%xmm4
mulsd  %xmm14,%xmm3
addsd  %xmm4,%xmm3
movsd  0x10(%rsp),%xmm2
mulsd  0xc7a(%rip),%xmm2
mulsd  %xmm12,%xmm2
movapd %xmm0,%xmm4
movsd  %xmm0,0x30(%rsp)
mulsd  %xmm0,%xmm2
movsd  0x58(%rsp),%xmm0
mulsd  %xmm4,%xmm0
movsd  0x50(%rsp),%xmm4
mulsd  %xmm4,%xmm0
mulsd  %xmm8,%xmm0
mulsd  %xmm9,%xmm2
mulsd  0xc49(%rip),%xmm4
addsd  %xmm0,%xmm2
movsd  %xmm2,0x10(%rsp)
mulsd  %xmm8,%xmm4
mulsd  %xmm6,%xmm4
movsd  0x8(%rsp),%xmm8
mulsd  %xmm8,%xmm1
addsd  %xmm4,%xmm1
movsd  0x30(%rsp),%xmm0
mulsd  %xmm0,%xmm7
movapd %xmm9,%xmm2
mulsd  %xmm9,%xmm7
movapd %xmm12,%xmm4
movapd %xmm12,%xmm9
addsd  %xmm12,%xmm4
movapd %xmm4,%xmm12
mulsd  %xmm2,%xmm12
mulsd  %xmm8,%xmm2
movsd  0xbcb(%rip),%xmm8
mulsd  %xmm0,%xmm8
subsd  %xmm8,%xmm4
mulsd  %xmm2,%xmm4
addsd  %xmm1,%xmm4
movsd  0xbd9(%rip),%xmm1
mulsd  %xmm6,%xmm1
addsd  %xmm10,%xmm1
mulsd  %xmm12,%xmm1
addsd  %xmm4,%xmm1
mulsd  %xmm9,%xmm13
mulsd  0xbc9(%rip),%xmm9
addsd  %xmm8,%xmm9
mulsd  %xmm3,%xmm8
mulsd  %xmm8,%xmm1
addsd  0x10(%rsp),%xmm1
mulsd  %xmm0,%xmm13
movsd  0x8(%rsp),%xmm0
movapd %xmm0,%xmm4
subsd  %xmm6,%xmm4
mulsd  %xmm13,%xmm4
mulsd  %xmm7,%xmm11
mulsd  0xb87(%rip),%xmm6
mulsd  0x40(%rsp),%xmm7
addsd  %xmm0,%xmm6
movapd %xmm0,%xmm8
mulsd  %xmm7,%xmm6
addsd  %xmm4,%xmm6
mulsd  %xmm2,%xmm9
addsd  %xmm6,%xmm9
mulsd  %xmm12,%xmm5
addsd  %xmm9,%xmm5
movapd %xmm15,%xmm0
movsd  0x28(%rsp),%xmm6
addsd  %xmm6,%xmm0
mulsd  %xmm2,%xmm0
mulsd  %xmm5,%xmm0
subsd  %xmm0,%xmm1
movsd  0x68(%rsp),%xmm2
movsd  0x18(%rsp),%xmm4
mulsd  %xmm4,%xmm2
movapd %xmm11,%xmm0
mulsd  %xmm3,%xmm0
addsd  %xmm1,%xmm0
mulsd  %xmm2,%xmm0
mulsd  %xmm8,%xmm15
movapd %xmm3,%xmm2
subsd  %xmm15,%xmm2
mulsd  %xmm4,%xmm14
mulsd  %xmm11,%xmm2
addsd  %xmm1,%xmm2
mulsd  %xmm14,%xmm2
mulsd  0x20(%rsp),%xmm4
mulsd  %xmm6,%xmm8
subsd  %xmm8,%xmm3
mulsd  %xmm11,%xmm3
addsd  %xmm1,%xmm3
mulsd  %xmm4,%xmm3
movsd  0x40(%rbx),%xmm1
mulsd  %xmm2,%xmm1
movsd  0x28(%rbx),%xmm4
mulsd  %xmm0,%xmm4
addsd  %xmm1,%xmm4
movsd  0x58(%rbx),%xmm1
mulsd  %xmm3,%xmm1
addsd  %xmm4,%xmm1
movsd  %xmm1,(%r14)
movsd  0x48(%rbx),%xmm1
mulsd  %xmm2,%xmm1
movsd  0x30(%rbx),%xmm4
mulsd  %xmm0,%xmm4
addsd  %xmm1,%xmm4
movsd  0x60(%rbx),%xmm1
mulsd  %xmm3,%xmm1
addsd  %xmm4,%xmm1
movsd  %xmm1,0x8(%r14)
mulsd  0x50(%rbx),%xmm2
mulsd  0x38(%rbx),%xmm0
addsd  %xmm2,%xmm0
mulsd  0x68(%rbx),%xmm3
addsd  %xmm0,%xmm3
movsd  %xmm3,0x10(%r14)
add    $0x78,%rsp
pop    %rbx
pop    %r14
ret

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

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