#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int fabs(double);
double pow(double, double);
double stub1(double *);
double stub2(double *);
double stub3(double *);
double stub4(double *);
double stub5(double *);
double stub6(double *);
void Yanggradient(int n, double x[], int need[], double f0, double g[],
                  double (*fun)(double x[]), double space[], int central,
                  double *lowbound, double *upbound) {
  int i, j;
  double *x0 = space, *x1 = space + n, eh0, eh01, eh;
  eh0 = eh01 = 1.e-8;
  if (central) {
    for (i = 0; i < n; i++) {
      if (need[i]) {
        for (j = 0; j < n; j++)
          x0[j] = x1[j] = x[j];
        eh = pow(eh01 * (fabs(x[i]) + 1), 0.67);
        x0[i] -= eh;
        x1[i] += eh;
        if (x0[i] < lowbound[i]) {
          x1[i] += eh;
          g[i] = ((*fun)(x1)-f0) / (eh * 2.0);
        } else if (x1[i] > upbound[i]) {
          x0[i] -= eh;
          g[i] = (f0 - (*fun)(x0)) / (eh * 2.0);
        } else
          g[i] = ((*fun)(x1) - (*fun)(x0)) / (eh * 2.0);
      }
    }
  } else {
    for (i = 0; i < n; i++) {
      if (need[i]) {
        for (j = 0; j < n; j++)
          x1[j] = x[j];
        eh = 2.0 * pow(eh0 * (fabs(x[i]) + 1), 0.67);
        if (x1[i] + eh > upbound[i]) {
          x1[i] -= eh;
          g[i] = (f0 - (*fun)(x1)) / eh;
        } else {
          x1[i] += eh;
          g[i] = ((*fun)(x1)-f0) / eh;
        }
      }
    }
  }
}


x86 assembly (clang, -O0):
<Yanggradient>:
push   %rbp
mov    %rsp,%rbp
sub    $0x90,%rsp
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%eax
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movsd  %xmm0,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    %r8,-0x30(%rbp)
mov    %r9,-0x38(%rbp)
mov    -0x38(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x38(%rbp),%rax
movslq -0x4(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x50(%rbp)
movsd  0xe8b(%rip),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  0xe7e(%rip),%xmm0
movsd  %xmm0,-0x58(%rbp)
cmpl   $0x0,0x10(%rbp)
je     13be <Yanggradient+0x29e>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    13b9 <Yanggradient+0x299>
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
je     13a6 <Yanggradient+0x286>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11fe <Yanggradient+0xde>
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x50(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x48(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    11bd <Yanggradient+0x9d>
movsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <fabs@plt>
movsd  -0x70(%rbp),%xmm0
add    $0x1,%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm1,%xmm0
movsd  0xdde(%rip),%xmm1
call   1040 <pow@plt>
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm1
mov    -0x48(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x68(%rbp),%xmm0
mov    -0x50(%rbp),%rax
movslq -0x3c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x48(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    12db <Yanggradient+0x1bb>
movsd  -0x68(%rbp),%xmm0
mov    -0x50(%rbp),%rax
movslq -0x3c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%rax
mov    -0x50(%rbp),%rdi
call   *%rax
subsd  -0x20(%rbp),%xmm0
movsd  0xd48(%rip),%xmm1
mulsd  -0x68(%rbp),%xmm1
divsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    13a1 <Yanggradient+0x281>
mov    -0x50(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
ucomisd (%rax,%rcx,8),%xmm0
jbe    1359 <Yanggradient+0x239>
movsd  -0x68(%rbp),%xmm1
mov    -0x48(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x30(%rbp),%rax
mov    -0x48(%rbp),%rdi
call   *%rax
movaps %xmm0,%xmm1
movsd  -0x78(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  0xcca(%rip),%xmm1
mulsd  -0x68(%rbp),%xmm1
divsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    139c <Yanggradient+0x27c>
mov    -0x30(%rbp),%rax
mov    -0x50(%rbp),%rdi
call   *%rax
movsd  %xmm0,-0x80(%rbp)
mov    -0x30(%rbp),%rax
mov    -0x48(%rbp),%rdi
call   *%rax
movaps %xmm0,%xmm1
movsd  -0x80(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  0xc82(%rip),%xmm1
mulsd  -0x68(%rbp),%xmm1
divsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    13a1 <Yanggradient+0x281>
jmp    13a6 <Yanggradient+0x286>
jmp    13ab <Yanggradient+0x28b>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1198 <Yanggradient+0x78>
jmp    1543 <Yanggradient+0x423>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    153e <Yanggradient+0x41e>
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
je     152b <Yanggradient+0x40b>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    141e <Yanggradient+0x2fe>
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x50(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    13ea <Yanggradient+0x2ca>
movsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <fabs@plt>
movsd  -0x88(%rbp),%xmm0
ff
add    $0x1,%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm1,%xmm0
movsd  0xbb8(%rip),%xmm1
call   1040 <pow@plt>
movaps %xmm0,%xmm1
movsd  0xba0(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x50(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x68(%rbp),%xmm0
mov    0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
ucomisd (%rax,%rcx,8),%xmm0
jbe    14ee <Yanggradient+0x3ce>
movsd  -0x68(%rbp),%xmm1
mov    -0x50(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0x30(%rbp),%rax
mov    -0x50(%rbp),%rdi
call   *%rax
movaps %xmm0,%xmm1
movsd  -0x90(%rbp),%xmm0
ff
subsd  %xmm1,%xmm0
divsd  -0x68(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1526 <Yanggradient+0x406>
movsd  -0x68(%rbp),%xmm0
mov    -0x50(%rbp),%rax
movslq -0x3c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%rax
mov    -0x50(%rbp),%rdi
call   *%rax
subsd  -0x20(%rbp),%xmm0
divsd  -0x68(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    152b <Yanggradient+0x40b>
jmp    1530 <Yanggradient+0x410>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    13c5 <Yanggradient+0x2a5>
jmp    1543 <Yanggradient+0x423>
add    $0x90,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<Yanggradient>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x88,%rsp
mov    %r9,%r14
mov    %r8,0x20(%rsp)
mov    %rcx,0x58(%rsp)
movsd  %xmm0,0x30(%rsp)
mov    %rdx,0x60(%rsp)
mov    %rsi,%rbx
movslq %edi,%r15
lea    (%r9,%r15,8),%rax
mov    %rax,0x8(%rsp)
cmpl   $0x0,0xc0(%rsp)
mov    %edi,0x4(%rsp)
je     14e9 <Yanggradient+0x3c9>
test   %edi,%edi
jle    1796 <Yanggradient+0x676>
mov    %edi,%ebp
lea    (%r15,%rbp,1),%rax
lea    (%r14,%rax,8),%r9
lea    (%r14,%rbp,8),%r8
mov    %ebp,%eax
and    $0xfffffffc,%eax
mov    %rax,0x48(%rsp)
add    $0xfffffffffffffffc,%rax
mov    %rax,0x40(%rsp)
mov    %rax,%rdx
shr    $0x2,%rdx
add    $0x1,%rdx
mov    0x8(%rsp),%rsi
cmp    %r8,%rsi
setb   %al
cmp    %r14,%r9
seta   %cl
lea    (%rbx,%rbp,8),%r10
and    %al,%cl
cmp    %r10,%rsi
setb   %sil
cmp    %rbx,%r9
seta   %al
and    %sil,%al
or     %cl,%al
cmp    %r14,%r10
seta   %sil
cmp    %rbx,%r8
seta   %cl
and    %sil,%cl
or     %al,%cl
mov    %cl,0x50(%rsp)
movabs $0x7ffffffffffffffc,%rcx
ff ff 7f
or     $0x2,%rcx
mov    %rdx,%rax
mov    %rdx,0x38(%rsp)
and    %rdx,%rcx
mov    %rcx,0x68(%rsp)
mov    %rbp,0x10(%rsp)
and    $0x3,%ebp
lea    (%r14,%r15,8),%r12
add    $0x30,%r12
lea    0x18(%r14),%rax
mov    %rax,0x78(%rsp)
lea    0x18(%rbx),%rax
mov    %rax,0x70(%rsp)
xor    %r13d,%r13d
mov    %rbp,0x18(%rsp)
jmp    127d <Yanggradient+0x15d>
nopw   0x0(%rax,%rax,1)
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rdi,%r13,8)
call   *0x20(%rsp)
movsd  0x28(%rsp),%xmm1
movapd %xmm0,%xmm2
movsd  0x30(%rsp),%xmm0
mov    0x4(%rsp),%edi
mov    0x18(%rsp),%rbp
subsd  %xmm0,%xmm2
addsd  %xmm1,%xmm1
divsd  %xmm1,%xmm2
mov    0x58(%rsp),%rax
movsd  %xmm2,(%rax,%r13,8)
add    $0x1,%r13
cmp    0x10(%rsp),%r13
je     1796 <Yanggradient+0x676>
mov    0x60(%rsp),%rax
cmpl   $0x0,(%rax,%r13,4)
je     126e <Yanggradient+0x14e>
cmp    $0x4,%edi
setb   %al
or     0x50(%rsp),%al
je     12b0 <Yanggradient+0x190>
xor    %eax,%eax
mov    0x8(%rsp),%rsi
jmp    135e <Yanggradient+0x23e>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
cmpq   $0x0,0x40(%rsp)
mov    0x8(%rsp),%rsi
je     14d7 <Yanggradient+0x3b7>
mov    0x68(%rsp),%rcx
xor    %eax,%eax
nopl   0x0(%rax,%rax,1)
movups (%rbx,%rax,8),%xmm0
movups 0x10(%rbx,%rax,8),%xmm1
movups %xmm0,-0x30(%r12,%rax,8)
movups %xmm1,-0x20(%r12,%rax,8)
movups %xmm0,(%r14,%rax,8)
movups %xmm1,0x10(%r14,%rax,8)
movupd 0x20(%rbx,%rax,8),%xmm0
movupd 0x30(%rbx,%rax,8),%xmm1
movupd %xmm0,-0x10(%r12,%rax,8)
movupd %xmm1,(%r12,%rax,8)
movupd %xmm0,0x20(%r14,%rax,8)
movupd %xmm1,0x30(%r14,%rax,8)
add    $0x8,%rax
add    $0xfffffffffffffffe,%rcx
jne    12d0 <Yanggradient+0x1b0>
testb  $0x1,0x38(%rsp)
je     134b <Yanggradient+0x22b>
movupd (%rbx,%rax,8),%xmm0
movupd 0x10(%rbx,%rax,8),%xmm1
movupd %xmm0,(%rsi,%rax,8)
movupd %xmm1,0x10(%rsi,%rax,8)
movupd %xmm0,(%r14,%rax,8)
movupd %xmm1,0x10(%r14,%rax,8)
mov    0x48(%rsp),%rcx
mov    %rcx,%rax
cmp    0x10(%rsp),%rcx
je     140f <Yanggradient+0x2ef>
mov    %rax,%rcx
not    %rcx
add    0x10(%rsp),%rcx
test   %rbp,%rbp
je     139a <Yanggradient+0x27a>
mov    %rbp,%rdx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  (%rbx,%rax,8),%xmm0
movsd  %xmm0,(%rsi,%rax,8)
movsd  %xmm0,(%r14,%rax,8)
add    $0x1,%rax
add    $0xffffffffffffffff,%rdx
jne    1380 <Yanggradient+0x260>
cmp    $0x3,%rcx
jb     140f <Yanggradient+0x2ef>
mov    0x10(%rsp),%rcx
sub    %rax,%rcx
mov    0x78(%rsp),%rdx
lea    (%rdx,%rax,8),%rdx
mov    0x70(%rsp),%rsi
lea    (%rsi,%rax,8),%rax
nopw   0x0(%rax,%rax,1)
movsd  -0x18(%rax),%xmm0
movsd  %xmm0,-0x18(%rdx,%r15,8)
movsd  %xmm0,-0x18(%rdx)
movsd  -0x10(%rax),%xmm0
movsd  %xmm0,-0x10(%rdx,%r15,8)
movsd  %xmm0,-0x10(%rdx)
movsd  -0x8(%rax),%xmm0
movsd  %xmm0,-0x8(%rdx,%r15,8)
movsd  %xmm0,-0x8(%rdx)
movsd  (%rax),%xmm0
movsd  %xmm0,(%rdx,%r15,8)
movsd  %xmm0,(%rdx)
add    $0x20,%rdx
add    $0x20,%rax
add    $0xfffffffffffffffc,%rcx
jne    13c0 <Yanggradient+0x2a0>
movsd  (%rbx,%r13,8),%xmm0
call   1030 <fabs@plt>
add    $0x1,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mulsd  0xbd4(%rip),%xmm0
movsd  0xbd4(%rip),%xmm1
call   1040 <pow@plt>
mov    0x8(%rsp),%rdi
movsd  (%r14,%r13,8),%xmm3
subsd  %xmm0,%xmm3
movsd  %xmm3,(%r14,%r13,8)
movsd  (%rdi,%r13,8),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rdi,%r13,8)
mov    0xc8(%rsp),%rax
movsd  (%rax,%r13,8),%xmm2
ucomisd %xmm3,%xmm2
movsd  %xmm0,0x28(%rsp)
ja     1230 <Yanggradient+0x110>
mov    0xd0(%rsp),%rax
ucomisd (%rax,%r13,8),%xmm1
jbe    14ae <Yanggradient+0x38e>
subsd  %xmm0,%xmm3
movsd  %xmm3,(%r14,%r13,8)
mov    %r14,%rdi
call   *0x20(%rsp)
movsd  0x28(%rsp),%xmm1
movsd  0x30(%rsp),%xmm2
jmp    124e <Yanggradient+0x12e>
mov    0x20(%rsp),%rbp
call   *%rbp
movsd  %xmm0,0x80(%rsp)
mov    %r14,%rdi
call   *%rbp
movsd  0x80(%rsp),%xmm2
movsd  0x28(%rsp),%xmm1
jmp    124e <Yanggradient+0x12e>
xor    %eax,%eax
testb  $0x1,0x38(%rsp)
jne    1328 <Yanggradient+0x208>
jmp    134b <Yanggradient+0x22b>
test   %edi,%edi
jle    1796 <Yanggradient+0x676>
mov    %edi,%r12d
lea    (%r15,%r12,1),%rax
lea    (%r14,%rax,8),%rax
lea    (%rbx,%r12,8),%rcx
mov    %r12d,%edx
and    $0xfffffffc,%edx
mov    %rdx,0x50(%rsp)
add    $0xfffffffffffffffc,%rdx
mov    %rdx,0x18(%rsp)
shr    $0x2,%rdx
add    $0x1,%rdx
cmp    %rcx,0x8(%rsp)
setb   %cl
cmp    %rbx,%rax
seta   %al
and    %cl,%al
mov    %al,0x28(%rsp)
mov    %edx,%eax
and    $0x3,%eax
and    $0xfffffffffffffffc,%rdx
mov    %rdx,0x38(%rsp)
mov    %r12d,%r13d
and    $0x3,%r13d
lea    (%r14,%r15,8),%rbp
add    $0x70,%rbp
lea    0x10(%rbx),%rcx
mov    %rcx,0x40(%rsp)
mov    %rax,0x48(%rsp)
mov    %rax,%r15
shl    $0x5,%r15
xor    %r14d,%r14d
jmp    15aa <Yanggradient+0x48a>
nopw   0x0(%rax,%rax,1)
movsd  %xmm1,(%rdi,%r14,8)
call   *0x20(%rsp)
movsd  0x10(%rsp),%xmm2
movsd  0x30(%rsp),%xmm1
mov    0x4(%rsp),%edi
subsd  %xmm1,%xmm0
divsd  %xmm2,%xmm0
mov    0x58(%rsp),%rax
movsd  %xmm0,(%rax,%r14,8)
add    $0x1,%r14
cmp    %r12,%r14
je     1796 <Yanggradient+0x676>
mov    0x60(%rsp),%rax
cmpl   $0x0,(%rax,%r14,4)
je     159d <Yanggradient+0x47d>
cmp    $0x4,%edi
setb   %al
or     0x28(%rsp),%al
je     15d0 <Yanggradient+0x4b0>
xor    %eax,%eax
jmp    169c <Yanggradient+0x57c>
nopl   0x0(%rax)
cmpq   $0xc,0x18(%rsp)
jae    15dc <Yanggradient+0x4bc>
xor    %eax,%eax
jmp    164d <Yanggradient+0x52d>
mov    0x38(%rsp),%rcx
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movups (%rbx,%rax,8),%xmm0
movups 0x10(%rbx,%rax,8),%xmm1
movups %xmm0,-0x70(%rbp,%rax,8)
movups %xmm1,-0x60(%rbp,%rax,8)
movups 0x20(%rbx,%rax,8),%xmm0
movups 0x30(%rbx,%rax,8),%xmm1
movups %xmm0,-0x50(%rbp,%rax,8)
movups %xmm1,-0x40(%rbp,%rax,8)
movups 0x40(%rbx,%rax,8),%xmm0
movups 0x50(%rbx,%rax,8),%xmm1
movups %xmm0,-0x30(%rbp,%rax,8)
movups %xmm1,-0x20(%rbp,%rax,8)
movupd 0x60(%rbx,%rax,8),%xmm0
movupd 0x70(%rbx,%rax,8),%xmm1
movupd %xmm0,-0x10(%rbp,%rax,8)
movupd %xmm1,0x0(%rbp,%rax,8)
add    $0x10,%rax
add    $0xfffffffffffffffc,%rcx
jne    15f0 <Yanggradient+0x4d0>
cmpq   $0x0,0x48(%rsp)
je     168f <Yanggradient+0x56f>
lea    -0x60(,%rax,8),%rcx
ff
add    %rbp,%rcx
mov    0x40(%rsp),%rdx
lea    (%rdx,%rax,8),%rax
xor    %edx,%edx
nopl   0x0(%rax,%rax,1)
movupd -0x10(%rax,%rdx,1),%xmm0
movupd (%rax,%rdx,1),%xmm1
movupd %xmm0,-0x10(%rcx,%rdx,1)
movupd %xmm1,(%rcx,%rdx,1)
add    $0x20,%rdx
cmp    %rdx,%r15
jne    1670 <Yanggradient+0x550>
mov    0x50(%rsp),%rcx
mov    %rcx,%rax
cmp    %r12,%rcx
je     1718 <Yanggradient+0x5f8>
mov    %rax,%rcx
not    %rcx
add    %r12,%rcx
test   %r13,%r13
je     16d4 <Yanggradient+0x5b4>
mov    %r13,%rdx
mov    0x8(%rsp),%rsi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%rbx,%rax,8),%xmm0
movsd  %xmm0,(%rsi,%rax,8)
add    $0x1,%rax
add    $0xffffffffffffffff,%rdx
jne    16c0 <Yanggradient+0x5a0>
cmp    $0x3,%rcx
jb     1718 <Yanggradient+0x5f8>
nopw   0x0(%rax,%rax,1)
movsd  (%rbx,%rax,8),%xmm0
movsd  %xmm0,-0x70(%rbp,%rax,8)
movsd  0x8(%rbx,%rax,8),%xmm0
movsd  %xmm0,-0x68(%rbp,%rax,8)
movsd  0x10(%rbx,%rax,8),%xmm0
movsd  %xmm0,-0x60(%rbp,%rax,8)
movsd  0x18(%rbx,%rax,8),%xmm0
movsd  %xmm0,-0x58(%rbp,%rax,8)
add    $0x4,%rax
cmp    %rax,%r12
jne    16e0 <Yanggradient+0x5c0>
movsd  (%rbx,%r14,8),%xmm0
call   1030 <fabs@plt>
add    $0x1,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mulsd  0x8cb(%rip),%xmm0
movsd  0x8cb(%rip),%xmm1
call   1040 <pow@plt>
addsd  %xmm0,%xmm0
mov    0x8(%rsp),%rdi
movsd  (%rdi,%r14,8),%xmm2
movapd %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    0xd0(%rsp),%rax
ucomisd (%rax,%r14,8),%xmm1
movsd  %xmm0,0x10(%rsp)
jbe    1570 <Yanggradient+0x450>
subsd  %xmm0,%xmm2
movsd  %xmm2,(%rdi,%r14,8)
call   *0x20(%rsp)
movsd  0x10(%rsp),%xmm2
movapd %xmm0,%xmm1
movsd  0x30(%rsp),%xmm0
jmp    1586 <Yanggradient+0x466>
add    $0x88,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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