#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 */
double cos(double);
double exp(double);
double fabs(double);
double pow(double, double);
double sin(double);
double sqrt(double);
void cerrf(double in_real, double in_imag, double *out_real, double *out_imag) {
  int n, nc, nu;
  double constant = 1.12837916709551;
  double xLim = 5.33;
  double yLim = 4.29;
  double h, q, Saux, Sx, Sy, Tn, Tx, Ty, Wx, Wy, xh, xl, x, yh, y;
  double Rx[33];
  double Ry[33];
  x = fabs(in_real);
  y = fabs(in_imag);
  if (y < yLim && x < xLim) {
    q = (1.0 - y / yLim) * sqrt(1.0 - (x / xLim) * (x / xLim));
    h = 1.0 / (3.2 * q);
    nc = 7 + (int)(23.0 * q);
    xl = pow(h, (double)(1 - nc));
    xh = y + 0.5 / h;
    yh = x;
    nu = 10 + (int)(21.0 * q);
    Rx[nu] = 0.;
    Ry[nu] = 0.;
    for (n = nu; n > 0; n--) {
      Tx = xh + n * Rx[n];
      Ty = yh - n * Ry[n];
      Tn = Tx * Tx + Ty * Ty;
      Rx[n - 1] = 0.5 * Tx / Tn;
      Ry[n - 1] = 0.5 * Ty / Tn;
    }
    Sx = 0.;
    Sy = 0.;
    for (n = nc; n > 0; n--) {
      Saux = Sx + xl;
      Sx = Rx[n - 1] * Saux - Ry[n - 1] * Sy;
      Sy = Rx[n - 1] * Sy + Ry[n - 1] * Saux;
      xl = h * xl;
    };
    Wx = constant * Sx;
    Wy = constant * Sy;
  } else {
    xh = y;
    yh = x;
    Rx[0] = 0.;
    Ry[0] = 0.;
    for (n = 9; n > 0; n--) {
      Tx = xh + n * Rx[0];
      Ty = yh - n * Ry[0];
      Tn = Tx * Tx + Ty * Ty;
      Rx[0] = 0.5 * Tx / Tn;
      Ry[0] = 0.5 * Ty / Tn;
    };
    Wx = constant * Rx[0];
    Wy = constant * Ry[0];
  }
  if (y == 0.) {
    Wx = exp(-x * x);
  }
  if (in_imag < 0.) {
    Wx = 2.0 * exp(y * y - x * x) * cos(2.0 * x * y) - Wx;
    Wy = -2.0 * exp(y * y - x * x) * sin(2.0 * x * y) - Wy;
    if (in_real > 0.) {
      Wy = -Wy;
    }
  } else if (in_real < 0.) {
    Wy = -Wy;
  }
  *out_real = Wx;
  *out_imag = Wy;
}


x86 assembly (clang, -O0):
<cerrf>:
push   %rbp
mov    %rsp,%rbp
sub    $0x300,%rsp
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
mov    %rdi,-0x18(%rbp)
mov    %rsi,-0x20(%rbp)
movsd  0xeab(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  0xe96(%rip),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  0xe81(%rip),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x8(%rbp),%xmm0
movaps 0xe60(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0x10(%rbp),%xmm0
movaps 0xe48(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
movsd  -0xc0(%rbp),%xmm1
ff
movsd  -0x48(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1525 <cerrf+0x3d5>
movsd  -0xb0(%rbp),%xmm1
ff
movsd  -0x40(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1525 <cerrf+0x3d5>
movsd  -0xc0(%rbp),%xmm1
ff
divsd  -0x48(%rbp),%xmm1
movsd  0xe39(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x2e8(%rbp)
ff
movsd  -0xb0(%rbp),%xmm0
ff
divsd  -0x40(%rbp),%xmm0
movsd  -0xb0(%rbp),%xmm1
ff
divsd  -0x40(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm1,%xmm0
movsd  0xdf0(%rip),%xmm1
addsd  %xmm1,%xmm0
call   1070 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0x2e8(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  0xdd3(%rip),%xmm1
mulsd  -0x58(%rbp),%xmm1
movsd  0xdbe(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  0xda5(%rip),%xmm0
mulsd  -0x58(%rbp),%xmm0
cvttsd2si %xmm0,%eax
add    $0x7,%eax
mov    %eax,-0x28(%rbp)
movsd  -0x50(%rbp),%xmm0
mov    $0x1,%eax
sub    -0x28(%rbp),%eax
cvtsi2sd %eax,%xmm1
call   1060 <pow@plt>
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0xc0(%rbp),%xmm0
ff
movsd  0xd58(%rip),%xmm1
divsd  -0x50(%rbp),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0xb0(%rbp),%xmm0
ff
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  0xd37(%rip),%xmm0
mulsd  -0x58(%rbp),%xmm0
cvttsd2si %xmm0,%eax
add    $0xa,%eax
mov    %eax,-0x2c(%rbp)
movslq -0x2c(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x1d0(%rbp,%rax,8)
ff ff
movslq -0x2c(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x2e0(%rbp,%rax,8)
ff ff
mov    -0x2c(%rbp),%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
jle    141e <cerrf+0x2ce>
movsd  -0xa0(%rbp),%xmm1
ff
cvtsi2sdl -0x24(%rbp),%xmm0
movslq -0x24(%rbp),%rax
movsd  -0x1d0(%rbp,%rax,8),%xmm2
ff ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0xb8(%rbp),%xmm1
ff
cvtsi2sdl -0x24(%rbp),%xmm0
movslq -0x24(%rbp),%rax
movsd  -0x2e0(%rbp,%rax,8),%xmm2
ff ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x80(%rbp),%xmm0
movsd  -0x80(%rbp),%xmm2
movsd  -0x88(%rbp),%xmm1
ff
mulsd  -0x88(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  0xc59(%rip),%xmm0
mulsd  -0x80(%rbp),%xmm0
divsd  -0x78(%rbp),%xmm0
mov    -0x24(%rbp),%eax
sub    $0x1,%eax
cltq
movsd  %xmm0,-0x1d0(%rbp,%rax,8)
ff ff
movsd  0xc36(%rip),%xmm0
mulsd  -0x88(%rbp),%xmm0
ff
divsd  -0x78(%rbp),%xmm0
mov    -0x24(%rbp),%eax
sub    $0x1,%eax
cltq
movsd  %xmm0,-0x2e0(%rbp,%rax,8)
ff ff
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
jmp    132e <cerrf+0x1de>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x68(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x28(%rbp),%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
jle    14fc <cerrf+0x3ac>
movsd  -0x68(%rbp),%xmm0
addsd  -0xa8(%rbp),%xmm0
ff
movsd  %xmm0,-0x60(%rbp)
mov    -0x24(%rbp),%eax
sub    $0x1,%eax
cltq
movsd  -0x1d0(%rbp,%rax,8),%xmm0
ff ff
movsd  -0x60(%rbp),%xmm2
mov    -0x24(%rbp),%eax
sub    $0x1,%eax
cltq
movsd  -0x2e0(%rbp,%rax,8),%xmm1
ff ff
mulsd  -0x70(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x24(%rbp),%eax
sub    $0x1,%eax
cltq
movsd  -0x1d0(%rbp,%rax,8),%xmm0
ff ff
movsd  -0x70(%rbp),%xmm2
mov    -0x24(%rbp),%eax
sub    $0x1,%eax
cltq
movsd  -0x2e0(%rbp,%rax,8),%xmm1
ff ff
mulsd  -0x60(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x50(%rbp),%xmm0
mulsd  -0xa8(%rbp),%xmm0
ff
movsd  %xmm0,-0xa8(%rbp)
ff
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
jmp    1434 <cerrf+0x2e4>
movsd  -0x38(%rbp),%xmm0
mulsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x38(%rbp),%xmm0
mulsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
jmp    1660 <cerrf+0x510>
movsd  -0xc0(%rbp),%xmm0
ff
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0xb0(%rbp),%xmm0
ff
movsd  %xmm0,-0xb8(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x1d0(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x2e0(%rbp)
ff
movl   $0x9,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
jle    1636 <cerrf+0x4e6>
movsd  -0xa0(%rbp),%xmm1
ff
cvtsi2sdl -0x24(%rbp),%xmm0
movsd  -0x1d0(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0xb8(%rbp),%xmm1
ff
cvtsi2sdl -0x24(%rbp),%xmm0
movsd  -0x2e0(%rbp),%xmm2
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x80(%rbp),%xmm0
movsd  -0x80(%rbp),%xmm2
movsd  -0x88(%rbp),%xmm1
ff
mulsd  -0x88(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  0xa2f(%rip),%xmm0
mulsd  -0x80(%rbp),%xmm0
divsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x1d0(%rbp)
ff
movsd  0xa15(%rip),%xmm0
mulsd  -0x88(%rbp),%xmm0
ff
divsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x2e0(%rbp)
ff
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
jmp    1562 <cerrf+0x412>
movsd  -0x38(%rbp),%xmm0
mulsd  -0x1d0(%rbp),%xmm0
ff
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x38(%rbp),%xmm0
mulsd  -0x2e0(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0xc0(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    16af <cerrf+0x55f>
jp     16af <cerrf+0x55f>
movsd  -0xb0(%rbp),%xmm0
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  -0xb0(%rbp),%xmm0
ff
call   1040 <exp@plt>
movsd  %xmm0,-0x90(%rbp)
ff
xorps  %xmm0,%xmm0
ucomisd -0x10(%rbp),%xmm0
jbe    185f <cerrf+0x70f>
movsd  -0xc0(%rbp),%xmm0
ff
movsd  -0xc0(%rbp),%xmm2
ff
movsd  -0xb0(%rbp),%xmm1
ff
mulsd  -0xb0(%rbp),%xmm1
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1040 <exp@plt>
movaps %xmm0,%xmm1
movsd  0x944(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x2f8(%rbp)
ff
movsd  0x930(%rip),%xmm0
mulsd  -0xb0(%rbp),%xmm0
ff
mulsd  -0xc0(%rbp),%xmm0
ff
call   1050 <cos@plt>
movaps %xmm0,%xmm2
movsd  -0x2f8(%rbp),%xmm0
ff
movsd  -0x90(%rbp),%xmm1
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0xc0(%rbp),%xmm0
ff
movsd  -0xc0(%rbp),%xmm2
ff
movsd  -0xb0(%rbp),%xmm1
ff
mulsd  -0xb0(%rbp),%xmm1
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1040 <exp@plt>
movaps %xmm0,%xmm1
movsd  0x89a(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x2f0(%rbp)
ff
movsd  0x87e(%rip),%xmm0
mulsd  -0xb0(%rbp),%xmm0
ff
mulsd  -0xc0(%rbp),%xmm0
ff
call   1030 <sin@plt>
movaps %xmm0,%xmm2
movsd  -0x2f0(%rbp),%xmm0
ff
movsd  -0x98(%rbp),%xmm1
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x8(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    185a <cerrf+0x70a>
movsd  -0x98(%rbp),%xmm0
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
jmp    1899 <cerrf+0x749>
xorps  %xmm0,%xmm0
ucomisd -0x8(%rbp),%xmm0
jbe    1894 <cerrf+0x744>
movsd  -0x98(%rbp),%xmm0
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
jmp    1899 <cerrf+0x749>
movsd  -0x90(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x98(%rbp),%xmm0
ff
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
add    $0x300,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<cerrf>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x288,%rsp
mov    %rsi,%r14
mov    %rdi,%r15
movapd %xmm0,%xmm2
movapd 0xe8d(%rip),%xmm5
movapd %xmm0,%xmm12
andpd  %xmm5,%xmm12
andpd  %xmm1,%xmm5
movsd  0xec7(%rip),%xmm0
ucomisd %xmm5,%xmm0
movapd %xmm2,0x50(%rsp)
movapd %xmm12,0x30(%rsp)
movapd %xmm1,0x10(%rsp)
movapd %xmm5,0x20(%rsp)
jbe    11ff <cerrf+0xaf>
movsd  0xea8(%rip),%xmm0
ucomisd %xmm12,%xmm0
jbe    11ff <cerrf+0xaf>
movapd %xmm5,%xmm6
divsd  0xf05(%rip),%xmm6
movsd  0xf05(%rip),%xmm4
addsd  %xmm4,%xmm6
movapd %xmm12,%xmm3
divsd  0xe7c(%rip),%xmm3
mulsd  %xmm3,%xmm3
movapd %xmm4,%xmm0
subsd  %xmm3,%xmm0
xorpd  %xmm3,%xmm3
ucomisd %xmm3,%xmm0
jb     1426 <cerrf+0x2d6>
sqrtsd %xmm0,%xmm0
jmp    143d <cerrf+0x2ed>
movq   $0x0,0x180(%rsp)
movq   $0x0,0x70(%rsp)
movapd %xmm2,%xmm0
mulsd  %xmm2,%xmm0
movapd %xmm1,%xmm2
mulsd  %xmm1,%xmm2
addsd  %xmm0,%xmm2
movsd  0xe30(%rip),%xmm0
movapd %xmm5,%xmm7
mulsd  %xmm0,%xmm7
divsd  %xmm2,%xmm7
movapd %xmm12,%xmm4
mulsd  %xmm0,%xmm4
divsd  %xmm2,%xmm4
mulsd  0xe17(%rip),%xmm7
addsd  %xmm5,%xmm7
mulsd  0xe13(%rip),%xmm4
addsd  %xmm12,%xmm4
movapd %xmm4,%xmm2
mulsd  %xmm4,%xmm2
movapd %xmm7,%xmm3
mulsd  %xmm7,%xmm3
addsd  %xmm2,%xmm3
mulsd  %xmm0,%xmm7
divsd  %xmm3,%xmm7
mulsd  %xmm0,%xmm4
divsd  %xmm3,%xmm4
mulsd  0xdea(%rip),%xmm7
addsd  %xmm5,%xmm7
mulsd  0xde6(%rip),%xmm4
addsd  %xmm12,%xmm4
movapd %xmm4,%xmm2
mulsd  %xmm4,%xmm2
movapd %xmm7,%xmm3
mulsd  %xmm7,%xmm3
addsd  %xmm2,%xmm3
mulsd  %xmm0,%xmm7
divsd  %xmm3,%xmm7
mulsd  %xmm0,%xmm4
mulsd  0xdc1(%rip),%xmm7
divsd  %xmm3,%xmm4
addsd  %xmm5,%xmm7
mulsd  0xdb9(%rip),%xmm4
addsd  %xmm12,%xmm4
movapd %xmm4,%xmm2
mulsd  %xmm4,%xmm2
movapd %xmm7,%xmm3
mulsd  %xmm7,%xmm3
addsd  %xmm2,%xmm3
mulsd  %xmm0,%xmm7
divsd  %xmm3,%xmm7
mulsd  %xmm0,%xmm4
divsd  %xmm3,%xmm4
mulsd  0xd90(%rip),%xmm7
mulsd  0xd90(%rip),%xmm4
addsd  %xmm5,%xmm7
addsd  %xmm12,%xmm4
movapd %xmm4,%xmm2
mulsd  %xmm4,%xmm2
movapd %xmm7,%xmm3
mulsd  %xmm7,%xmm3
addsd  %xmm2,%xmm3
mulsd  %xmm0,%xmm7
divsd  %xmm3,%xmm7
mulsd  %xmm0,%xmm4
divsd  %xmm3,%xmm4
mulsd  0xd63(%rip),%xmm7
addsd  %xmm5,%xmm7
mulsd  0xd5f(%rip),%xmm4
addsd  %xmm12,%xmm4
movapd %xmm4,%xmm2
mulsd  %xmm4,%xmm2
movapd %xmm7,%xmm3
mulsd  %xmm7,%xmm3
addsd  %xmm2,%xmm3
mulsd  %xmm0,%xmm7
divsd  %xmm3,%xmm7
mulsd  %xmm0,%xmm4
divsd  %xmm3,%xmm4
mulsd  0xd36(%rip),%xmm7
addsd  %xmm5,%xmm7
mulsd  0xd32(%rip),%xmm4
addsd  %xmm12,%xmm4
movapd %xmm4,%xmm2
mulsd  %xmm4,%xmm2
movapd %xmm7,%xmm3
mulsd  %xmm7,%xmm3
addsd  %xmm2,%xmm3
mulsd  %xmm0,%xmm7
divsd  %xmm3,%xmm7
mulsd  %xmm0,%xmm4
divsd  %xmm3,%xmm4
addsd  %xmm7,%xmm7
addsd  %xmm5,%xmm7
addsd  %xmm4,%xmm4
movapd %xmm12,%xmm2
subsd  %xmm4,%xmm2
movapd %xmm2,%xmm4
mulsd  %xmm2,%xmm4
movapd %xmm7,%xmm3
mulsd  %xmm7,%xmm3
addsd  %xmm4,%xmm3
mulsd  %xmm0,%xmm7
divsd  %xmm3,%xmm7
mulsd  %xmm0,%xmm2
divsd  %xmm3,%xmm2
addsd  %xmm5,%xmm7
movapd %xmm12,%xmm0
subsd  %xmm2,%xmm0
movapd %xmm7,%xmm2
unpcklpd %xmm0,%xmm7
mulsd  %xmm0,%xmm0
mulsd  %xmm2,%xmm2
addsd  %xmm0,%xmm2
mulpd  0xbf7(%rip),%xmm7
unpcklpd %xmm2,%xmm2
divpd  %xmm2,%xmm7
jmp    182c <cerrf+0x6dc>
movsd  %xmm6,(%rsp)
call   1070 <sqrt@plt>
movsd  (%rsp),%xmm6
movsd  0xc93(%rip),%xmm4
mulsd  %xmm0,%xmm6
movsd  %xmm6,(%rsp)
movsd  0xc8a(%rip),%xmm0
mulsd  %xmm6,%xmm0
divsd  %xmm0,%xmm4
movsd  %xmm4,0x40(%rsp)
movsd  0xc7c(%rip),%xmm0
mulsd  %xmm6,%xmm0
cvttsd2si %xmm0,%r12d
mov    $0xfffffffa,%eax
sub    %r12d,%eax
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
movapd %xmm4,%xmm0
call   1060 <pow@plt>
movsd  0x40(%rsp),%xmm13
movsd  (%rsp),%xmm1
mulsd  0xc4f(%rip),%xmm1
cvttsd2si %xmm1,%ecx
movslq %ecx,%rax
movq   $0x0,0x1d0(%rsp,%rax,8)
movq   $0x0,0xc0(%rsp,%rax,8)
cmp    $0xfffffff7,%ecx
movapd 0x30(%rsp),%xmm12
jl     179e <cerrf+0x64e>
movsd  0xb90(%rip),%xmm1
movapd %xmm1,%xmm2
divsd  %xmm13,%xmm2
movapd 0x20(%rsp),%xmm3
addsd  %xmm3,%xmm2
lea    0xa(%rcx),%esi
mov    %esi,%r10d
test   %esi,%esi
mov    $0x1,%eax
mov    $0x1,%edx
cmovle %esi,%edx
mov    %ecx,%r9d
sub    %edx,%r9d
add    $0xa,%r9d
je     1708 <cerrf+0x5b8>
test   %esi,%esi
cmovle %esi,%eax
neg    %eax
lea    (%rcx,%rax,1),%edx
add    $0xa,%edx
lea    0x9(%rcx),%eax
lea    (%rsp,%rax,8),%rbx
add    $0x180,%rbx
shl    $0x3,%rdx
mov    %rbx,%rdi
sub    %rdx,%rdi
cmp    %rbx,%rdi
ja     1708 <cerrf+0x5b8>
lea    (%rsp,%r10,8),%rbx
add    $0x180,%rbx
mov    %rbx,%rdi
sub    %rdx,%rdi
cmp    %rbx,%rdi
ja     1708 <cerrf+0x5b8>
lea    (%rsp,%rax,8),%rbx
add    $0x70,%rbx
mov    %rbx,%rdi
sub    %rdx,%rdi
cmp    %rbx,%rdi
ja     1708 <cerrf+0x5b8>
lea    (%rsp,%r10,8),%rbx
add    $0x70,%rbx
mov    %rbx,%rdi
sub    %rdx,%rdi
cmp    %rbx,%rdi
ja     1708 <cerrf+0x5b8>
test   %esi,%esi
mov    $0x1,%edx
cmovle %esi,%edx
neg    %edx
add    %edx,%ecx
add    $0xa,%ecx
mov    %rax,%rdx
sub    %rcx,%rdx
lea    (%rsp,%rdx,8),%r11
add    $0x180,%r11
lea    (%rsp,%rax,8),%r8
add    $0x188,%r8
mov    %r10,%rbx
sub    %rcx,%rbx
lea    (%rsp,%rbx,8),%rbp
add    $0x180,%rbp
lea    0x188(%rsp,%r10,8),%rdi
lea    0x70(%rsp,%rdx,8),%rdx
lea    (%rsp,%rax,8),%r13
add    $0x78,%r13
lea    (%rsp,%rbx,8),%rbx
add    $0x70,%rbx
lea    (%rsp,%r10,8),%rax
add    $0x78,%rax
cmp    %rdi,%r11
setb   %dil
cmp    %r8,%rbp
setb   %cl
cmp    %rax,%rdx
setb   %al
cmp    %r13,%rbx
setb   %dl
test   %cl,%dil
jne    1708 <cerrf+0x5b8>
and    %dl,%al
jne    1708 <cerrf+0x5b8>
add    $0x1,%r9
mov    %r9,%r8
and    $0xfffffffffffffffe,%r8
mov    %r10,%rcx
sub    %r8,%rcx
movd   %esi,%xmm3
pshufd $0x50,%xmm3,%xmm3
paddd  0x9f3(%rip),%xmm3
movapd %xmm2,%xmm9
unpcklpd %xmm2,%xmm9
movapd %xmm12,%xmm10
unpcklpd %xmm12,%xmm10
mov    %r8,%rsi
neg    %rsi
lea    (%rsp,%r10,8),%rdi
add    $0x178,%rdi
lea    (%rsp,%r10,8),%rax
add    $0x68,%rax
mov    $0xffffffff,%edx
add    %r10,%rdx
xor    %ebx,%ebx
movapd 0x9a3(%rip),%xmm11
movdqa 0x9ba(%rip),%xmm8
cs nopw 0x0(%rax,%rax,1)
cvtdq2pd %xmm3,%xmm4
movupd (%rdi,%rbx,8),%xmm7
shufpd $0x1,%xmm7,%xmm7
mulpd  %xmm4,%xmm7
movupd (%rax,%rbx,8),%xmm5
shufpd $0x1,%xmm5,%xmm5
mulpd  %xmm4,%xmm5
movapd %xmm10,%xmm4
addpd  %xmm9,%xmm7
subpd  %xmm5,%xmm4
movapd %xmm4,%xmm5
mulpd  %xmm4,%xmm5
movapd %xmm7,%xmm6
mulpd  %xmm7,%xmm6
addpd  %xmm5,%xmm6
mulpd  %xmm11,%xmm7
divpd  %xmm6,%xmm7
lea    (%rdx,%rbx,1),%ebp
shufpd $0x1,%xmm7,%xmm7
mulpd  %xmm11,%xmm4
divpd  %xmm6,%xmm4
shufpd $0x1,%xmm4,%xmm4
movupd %xmm7,0x178(%rsp,%rbp,8)
movupd %xmm4,0x68(%rsp,%rbp,8)
paddd  %xmm8,%xmm3
add    $0xfffffffffffffffe,%rbx
cmp    %rbx,%rsi
jne    1680 <cerrf+0x530>
cmp    %r8,%r9
jne    170b <cerrf+0x5bb>
jmp    179e <cerrf+0x64e>
mov    %r10,%rcx
lea    (%rsp,%rcx,8),%rax
add    $0x70,%rax
lea    0x180(%rsp,%rcx,8),%rdx
lea    -0x1(%rcx),%edi
lea    (%rsp,%rdi,8),%rsi
add    $0x70,%rsi
lea    (%rsp,%rdi,8),%rdi
add    $0x180,%rdi
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
xorps  %xmm3,%xmm3
cvtsi2sd %ecx,%xmm3
movsd  (%rdx,%rbx,1),%xmm4
mulsd  %xmm3,%xmm4
addsd  %xmm2,%xmm4
mulsd  (%rax,%rbx,1),%xmm3
movapd %xmm12,%xmm5
subsd  %xmm3,%xmm5
movapd %xmm5,%xmm3
mulsd  %xmm5,%xmm3
movapd %xmm4,%xmm6
mulsd  %xmm4,%xmm6
addsd  %xmm3,%xmm6
mulsd  %xmm1,%xmm4
divsd  %xmm6,%xmm4
movsd  %xmm4,(%rdi,%rbx,1)
mulsd  %xmm1,%xmm5
divsd  %xmm6,%xmm5
movsd  %xmm5,(%rsi,%rbx,1)
lea    -0x1(%rcx),%ebp
add    $0xfffffffffffffff8,%rbx
cmp    $0x1,%ecx
mov    %ebp,%ecx
jg     1740 <cerrf+0x5f0>
xorpd  %xmm7,%xmm7
cmp    $0xfffffffa,%r12d
movapd 0x10(%rsp),%xmm1
jl     182c <cerrf+0x6dc>
lea    0x8(%r12),%eax
add    $0x6,%r12d
lea    (%rsp,%r12,8),%rcx
add    $0x70,%rcx
lea    (%rsp,%r12,8),%rdx
add    $0x180,%rdx
xorpd  %xmm7,%xmm7
xor    %esi,%esi
movapd 0x868(%rip),%xmm6
nopl   0x0(%rax,%rax,1)
movapd %xmm0,%xmm2
addsd  %xmm7,%xmm2
movsd  (%rcx,%rsi,8),%xmm3
movapd %xmm3,%xmm4
xorpd  %xmm6,%xmm4
movapd %xmm7,%xmm5
shufpd $0x1,%xmm2,%xmm5
unpcklpd %xmm3,%xmm4
mulpd  %xmm5,%xmm4
movsd  (%rdx,%rsi,8),%xmm2
unpcklpd %xmm2,%xmm2
addsd  %xmm0,%xmm7
mulpd  %xmm2,%xmm7
addpd  %xmm4,%xmm7
mulsd  %xmm13,%xmm0
add    $0xffffffffffffffff,%rsi
lea    (%rax,%rsi,1),%edi
cmp    $0x1,%edi
jg     17e0 <cerrf+0x690>
movapd %xmm7,%xmm3
unpckhpd %xmm7,%xmm3
movsd  0x8b4(%rip),%xmm0
xorpd  %xmm4,%xmm4
ucomisd %xmm4,%xmm1
jne    1881 <cerrf+0x731>
jp     1881 <cerrf+0x731>
movapd 0x7f0(%rip),%xmm0
xorpd  %xmm12,%xmm0
mulsd  %xmm12,%xmm0
movapd %xmm3,(%rsp)
call   1040 <exp@plt>
xorpd  %xmm4,%xmm4
movapd (%rsp),%xmm3
movapd 0x10(%rsp),%xmm1
movapd %xmm0,%xmm7
movsd  0x871(%rip),%xmm0
jmp    1885 <cerrf+0x735>
mulsd  %xmm0,%xmm7
movapd 0x50(%rsp),%xmm2
mulsd  %xmm0,%xmm3
ucomisd %xmm1,%xmm4
jbe    194d <cerrf+0x7fd>
movapd %xmm2,%xmm0
mulsd  %xmm2,%xmm0
mulsd  %xmm1,%xmm1
subsd  %xmm0,%xmm1
movapd %xmm1,0x10(%rsp)
movapd %xmm1,%xmm0
movapd %xmm7,0x40(%rsp)
movapd %xmm3,(%rsp)
call   1040 <exp@plt>
addsd  %xmm0,%xmm0
movsd  %xmm0,0x68(%rsp)
movapd 0x30(%rsp),%xmm0
addsd  %xmm0,%xmm0
mulsd  0x20(%rsp),%xmm0
movapd %xmm0,0x30(%rsp)
call   1050 <cos@plt>
mulsd  0x68(%rsp),%xmm0
subsd  0x40(%rsp),%xmm0
movsd  %xmm0,0x20(%rsp)
movapd 0x10(%rsp),%xmm0
call   1040 <exp@plt>
mulsd  0x7eb(%rip),%xmm0
movsd  %xmm0,0x10(%rsp)
movapd 0x30(%rsp),%xmm0
call   1030 <sin@plt>
mulsd  0x10(%rsp),%xmm0
subsd  (%rsp),%xmm0
xorpd  %xmm1,%xmm1
movapd 0x50(%rsp),%xmm2
ucomisd %xmm1,%xmm2
jbe    1941 <cerrf+0x7f1>
xorpd  0x6ff(%rip),%xmm0
movapd %xmm0,%xmm3
movsd  0x20(%rsp),%xmm7
jmp    195f <cerrf+0x80f>
xorpd  %xmm0,%xmm0
ucomisd %xmm2,%xmm0
jbe    195f <cerrf+0x80f>
xorpd  0x6e1(%rip),%xmm3
movsd  %xmm7,(%r15)
movsd  %xmm3,(%r14)
add    $0x288,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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