#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 sin(int);
double const sqrt(int);
double tan(double);
void SubGussFs(double *X, double *Y, double B, double L, int nCenterLongi) {
  int nzonenum;
  if (nCenterLongi == 0) {
    nzonenum = (int)L / 6 + 1;
    nCenterLongi = nzonenum * 6 - 3;
  } else
    nzonenum = (int)nCenterLongi / 6 + 1;
  double rB = B / 180 * 3.14159265358979323846;
  double rL = (L - nCenterLongi) / 180 * 3.14159265358979323846;
  const double a = 6378245.00;
  const double b = 6356863.50;
  double sqre1 = (a * a - b * b) / (a * a);
  double sinb = sin(rB);
  double cosb = cos(rB);
  double M = a * (1 - sqre1) / (1 - sqre1 * sinb * sinb) /
             sqrt(1 - sqre1 * sinb * sinb);
  double N = a / sqrt(1 - sqre1 * sinb * sinb);
  double sqrita = N / M - 1;
  double s = a * (1 - sqre1) *
             (1.00505117739 * rB - 0.00506237764 / 2 * sin(2 * rB) +
              0.0000106245 / 4 * sin(4 * rB) - 0.00000002081 / 6 * sin(6 * rB));
  double tanb = tan(rB);
  *X = s + rL * rL * N / 2 * sinb * cosb +
       rL * rL * rL * rL * N / 24 * sinb * cosb * cosb * cosb *
           (5 - tanb * tanb + 9 * sqrita * sqrita + 4 * sqrita);
  *Y = rL * N * cosb +
       rL * rL * rL * N / 6 * cosb * cosb * cosb * (1 - tanb * tanb + sqrita) +
       rL * rL * rL * rL * rL * N / 120 * cosb * cosb * cosb * cosb * cosb *
           (5 - 18 * tanb * tanb + tanb * tanb * tanb * tanb);
  *Y = *Y + 500000 + nzonenum * 1.0e+6;
}


x86 assembly (clang, -O0):
<SubGussFs>:
push   %rbp
mov    %rsp,%rbp
sub    $0xb0,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
movsd  %xmm1,-0x20(%rbp)
mov    %edx,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
jne    118c <SubGussFs+0x4c>
cvttsd2si -0x20(%rbp),%eax
mov    $0x6,%ecx
cltd
idiv   %ecx
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
imul   $0x6,-0x28(%rbp),%eax
sub    $0x3,%eax
mov    %eax,-0x24(%rbp)
jmp    119d <SubGussFs+0x5d>
mov    -0x24(%rbp),%eax
mov    $0x6,%ecx
cltd
idiv   %ecx
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  0xeee(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  0xeda(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x20(%rbp),%xmm0
cvtsi2sdl -0x24(%rbp),%xmm1
subsd  %xmm1,%xmm0
movsd  0xec3(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  0xeaf(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  0xe86(%rip),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  0xe89(%rip),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  0xe74(%rip),%xmm0
movsd  %xmm0,-0x50(%rbp)
cvttsd2si -0x30(%rbp),%edi
call   1030 <sin@plt>
movsd  %xmm0,-0x58(%rbp)
movsd  -0x30(%rbp),%xmm0
call   1050 <cos@plt>
movsd  %xmm0,-0x60(%rbp)
movsd  0xdf1(%rip),%xmm1
subsd  -0x50(%rbp),%xmm1
movsd  0xe34(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  -0x50(%rbp),%xmm1
mulsd  -0x58(%rbp),%xmm1
movsd  -0x58(%rbp),%xmm2
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm1
movsd  0xdae(%rip),%xmm2
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0x50(%rbp),%xmm0
mulsd  -0x58(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm1,%xmm0
movsd  0xd6c(%rip),%xmm1
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1060 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0xb0(%rbp),%xmm0
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x50(%rbp),%xmm0
mulsd  -0x58(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm1,%xmm0
movsd  0xd19(%rip),%xmm1
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1060 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  0xd51(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x70(%rbp),%xmm0
divsd  -0x68(%rbp),%xmm0
movsd  0xce6(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  0xcd5(%rip),%xmm1
subsd  -0x50(%rbp),%xmm1
movsd  0xd18(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  0xccf(%rip),%xmm0
mulsd  -0x30(%rbp),%xmm0
cvttsd2si %xmm0,%edi
call   1030 <sin@plt>
movaps %xmm0,%xmm2
movsd  -0xa8(%rbp),%xmm0
ff
movsd  0xcce(%rip),%xmm1
mulsd  %xmm2,%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
movsd  0xca3(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  0xc5b(%rip),%xmm0
mulsd  -0x30(%rbp),%xmm0
cvttsd2si %xmm0,%edi
call   1030 <sin@plt>
movsd  -0xa0(%rbp),%xmm1
ff
movsd  0xc65(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  0xc1d(%rip),%xmm0
mulsd  -0x30(%rbp),%xmm0
cvttsd2si %xmm0,%edi
call   1030 <sin@plt>
movsd  -0x98(%rbp),%xmm2
ff
movaps %xmm0,%xmm1
movsd  -0x90(%rbp),%xmm0
ff
movsd  0xc1c(%rip),%xmm3
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x30(%rbp),%xmm0
call   1040 <tan@plt>
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x80(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm1
mulsd  -0x38(%rbp),%xmm1
mulsd  -0x70(%rbp),%xmm1
movsd  0xbd5(%rip),%xmm2
divsd  %xmm2,%xmm1
mulsd  -0x58(%rbp),%xmm1
movsd  -0x60(%rbp),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x38(%rbp),%xmm0
mulsd  -0x38(%rbp),%xmm0
mulsd  -0x38(%rbp),%xmm0
mulsd  -0x38(%rbp),%xmm0
mulsd  -0x70(%rbp),%xmm0
movsd  0xb96(%rip),%xmm2
divsd  %xmm2,%xmm0
mulsd  -0x58(%rbp),%xmm0
mulsd  -0x60(%rbp),%xmm0
mulsd  -0x60(%rbp),%xmm0
mulsd  -0x60(%rbp),%xmm0
movsd  -0x88(%rbp),%xmm2
ff
movsd  -0x88(%rbp),%xmm3
ff
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm3,%xmm2
movsd  0xb13(%rip),%xmm3
addsd  %xmm3,%xmm2
movsd  0xb37(%rip),%xmm3
mulsd  -0x78(%rbp),%xmm3
movsd  -0x78(%rbp),%xmm4
mulsd  %xmm4,%xmm3
addsd  %xmm2,%xmm3
movsd  -0x78(%rbp),%xmm2
movsd  0xb10(%rip),%xmm4
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x38(%rbp),%xmm1
mulsd  -0x70(%rbp),%xmm1
movsd  -0x60(%rbp),%xmm2
movsd  -0x38(%rbp),%xmm0
mulsd  -0x38(%rbp),%xmm0
mulsd  -0x38(%rbp),%xmm0
mulsd  -0x70(%rbp),%xmm0
movsd  0xac5(%rip),%xmm3
divsd  %xmm3,%xmm0
mulsd  -0x60(%rbp),%xmm0
mulsd  -0x60(%rbp),%xmm0
mulsd  -0x60(%rbp),%xmm0
movsd  -0x88(%rbp),%xmm3
ff
movsd  -0x88(%rbp),%xmm4
ff
movq   %xmm3,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm3
mulsd  %xmm4,%xmm3
movsd  0xa77(%rip),%xmm4
addsd  %xmm4,%xmm3
addsd  -0x78(%rbp),%xmm3
mulsd  %xmm3,%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x38(%rbp),%xmm0
mulsd  -0x38(%rbp),%xmm0
mulsd  -0x38(%rbp),%xmm0
mulsd  -0x38(%rbp),%xmm0
mulsd  -0x38(%rbp),%xmm0
mulsd  -0x70(%rbp),%xmm0
movsd  0xa34(%rip),%xmm2
divsd  %xmm2,%xmm0
mulsd  -0x60(%rbp),%xmm0
mulsd  -0x60(%rbp),%xmm0
mulsd  -0x60(%rbp),%xmm0
mulsd  -0x60(%rbp),%xmm0
mulsd  -0x60(%rbp),%xmm0
movsd  0xa07(%rip),%xmm3
mulsd  -0x88(%rbp),%xmm3
ff
movsd  -0x88(%rbp),%xmm2
ff
movq   %xmm3,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm3
mulsd  %xmm2,%xmm3
movsd  0x9cc(%rip),%xmm2
addsd  %xmm2,%xmm3
movsd  -0x88(%rbp),%xmm2
ff
mulsd  -0x88(%rbp),%xmm2
ff
mulsd  -0x88(%rbp),%xmm2
ff
movsd  -0x88(%rbp),%xmm4
ff
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x10(%rbp),%rax
movsd  0x97c(%rip),%xmm1
addsd  (%rax),%xmm1
cvtsi2sdl -0x28(%rbp),%xmm0
movsd  0x963(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,(%rax)
add    $0xb0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<SubGussFs>:
push   %rbp
push   %r15
push   %r14
push   %rbx
sub    $0x38,%rsp
mov    %rsi,%r14
mov    %rdi,%r15
test   %edx,%edx
je     1171 <SubGussFs+0x31>
movslq %edx,%rax
imul   $0x2aaaaaab,%rax,%rax
mov    %rax,%rcx
shr    $0x20,%rcx
shr    $0x3f,%rax
lea    (%rcx,%rax,1),%ebx
add    $0x1,%ebx
jmp    1198 <SubGussFs+0x58>
cvttsd2si %xmm1,%eax
cltq
imul   $0x2aaaaaab,%rax,%rax
mov    %rax,%rcx
shr    $0x20,%rcx
shr    $0x3f,%rax
lea    (%rcx,%rax,1),%ebx
add    $0x1,%ebx
lea    (%rbx,%rbx,2),%eax
lea    (%rax,%rax,1),%edx
add    $0xfffffffd,%edx
movsd  0xe60(%rip),%xmm4
divsd  %xmm4,%xmm0
movsd  0xe5c(%rip),%xmm3
mulsd  %xmm3,%xmm0
movsd  %xmm0,0x10(%rsp)
cvtsi2sd %edx,%xmm2
subsd  %xmm2,%xmm1
divsd  %xmm4,%xmm1
mulsd  %xmm3,%xmm1
movsd  %xmm1,0x28(%rsp)
cvttsd2si %xmm0,%edi
call   1030 <sin@plt>
movsd  %xmm0,0x18(%rsp)
movsd  0x10(%rsp),%xmm0
call   1050 <cos@plt>
movsd  %xmm0,0x30(%rsp)
movsd  0xe1c(%rip),%xmm0
movsd  0x18(%rsp),%xmm1
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
addsd  0xe0e(%rip),%xmm0
movsd  0xe0e(%rip),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rsp)
cvttsd2si %xmm0,%ebp
mov    %ebp,%edi
call   1060 <sqrt@plt>
movsd  0x8(%rsp),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rsp)
mov    %ebp,%edi
call   1060 <sqrt@plt>
movsd  0xde2(%rip),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x20(%rsp)
movapd %xmm1,%xmm0
divsd  0x8(%rsp),%xmm0
addsd  0xdce(%rip),%xmm0
movsd  %xmm0,0x8(%rsp)
movsd  0x10(%rsp),%xmm0
addsd  %xmm0,%xmm0
cvttsd2si %xmm0,%edi
call   1030 <sin@plt>
mulsd  0xdb5(%rip),%xmm0
movsd  0xdb5(%rip),%xmm2
movsd  0x10(%rsp),%xmm1
mulsd  %xmm1,%xmm2
addsd  %xmm0,%xmm2
movsd  %xmm2,(%rsp)
movapd %xmm1,%xmm0
mulsd  0xd9e(%rip),%xmm0
cvttsd2si %xmm0,%edi
call   1030 <sin@plt>
mulsd  0xd95(%rip),%xmm0
addsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
movsd  0x10(%rsp),%xmm0
mulsd  0xd85(%rip),%xmm0
cvttsd2si %xmm0,%edi
call   1030 <sin@plt>
mulsd  0xd7c(%rip),%xmm0
addsd  (%rsp),%xmm0
mulsd  0xd2f(%rip),%xmm0
movsd  %xmm0,(%rsp)
movsd  0x10(%rsp),%xmm0
call   1040 <tan@plt>
movsd  0x28(%rsp),%xmm9
movapd %xmm9,%xmm4
mulsd  %xmm9,%xmm4
movapd %xmm9,%xmm13
mulsd  %xmm4,%xmm13
movsd  0x20(%rsp),%xmm12
mulsd  %xmm12,%xmm4
mulsd  0xd38(%rip),%xmm4
movsd  0x18(%rsp),%xmm3
mulsd  %xmm3,%xmm4
movsd  0x30(%rsp),%xmm10
mulsd  %xmm10,%xmm4
addsd  (%rsp),%xmm4
movapd %xmm9,%xmm1
mulsd  %xmm13,%xmm1
movapd %xmm12,%xmm6
mulsd  %xmm1,%xmm6
divsd  0xd0a(%rip),%xmm6
mulsd  %xmm3,%xmm6
mulsd  %xmm10,%xmm6
mulsd  %xmm10,%xmm6
mulsd  %xmm10,%xmm6
movapd 0xd2f(%rip),%xmm5
xorpd  %xmm0,%xmm5
mulsd  %xmm0,%xmm5
movsd  0xce6(%rip),%xmm8
movapd %xmm5,%xmm11
addsd  %xmm8,%xmm11
movsd  0xcdc(%rip),%xmm3
movsd  0x8(%rsp),%xmm7
mulsd  %xmm7,%xmm3
mulsd  %xmm7,%xmm3
addsd  %xmm11,%xmm3
movsd  0xc89(%rip),%xmm2
mulsd  %xmm7,%xmm2
addsd  %xmm3,%xmm2
mulsd  %xmm6,%xmm2
addsd  %xmm4,%xmm2
movsd  %xmm2,(%r15)
movapd %xmm9,%xmm4
mulsd  %xmm12,%xmm4
mulsd  %xmm12,%xmm13
divsd  0xc6c(%rip),%xmm13
mulsd  %xmm10,%xmm13
mulsd  %xmm10,%xmm13
mulsd  %xmm10,%xmm13
addsd  0xc15(%rip),%xmm5
addsd  %xmm7,%xmm5
mulsd  %xmm13,%xmm5
mulsd  %xmm10,%xmm4
mulsd  %xmm9,%xmm1
mulsd  %xmm12,%xmm1
divsd  0xc65(%rip),%xmm1
addsd  %xmm5,%xmm4
mulsd  %xmm10,%xmm1
mulsd  %xmm10,%xmm1
mulsd  %xmm10,%xmm1
mulsd  %xmm10,%xmm1
mulsd  %xmm10,%xmm1
movsd  0xc48(%rip),%xmm2
mulsd  %xmm0,%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm8,%xmm2
movapd %xmm0,%xmm3
mulsd  %xmm0,%xmm3
mulsd  %xmm0,%xmm3
mulsd  %xmm0,%xmm3
addsd  %xmm2,%xmm3
mulsd  %xmm1,%xmm3
addsd  %xmm4,%xmm3
addsd  0xc1f(%rip),%xmm3
cvtsi2sd %ebx,%xmm0
mulsd  0xc1b(%rip),%xmm0
addsd  %xmm3,%xmm0
movsd  %xmm0,(%r14)
add    $0x38,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @SubGussFs(double* %0, double* %1, double %2, double %3, i32 %4) {
%6 = alloca double*
%7 = alloca double*
%8 = alloca double
%9 = alloca double
%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
store double* %0, double** %6
store double* %1, double** %7
store double %2, double* %8
store double %3, double* %9
store i32 %4, i32* %10
%24 = load i32, i32* %10
%25 = icmp eq i32 %24, 0
br i1 %25, label %26, label %34
26: ; preds = %5
%27 = load double, double* %9
%28 = fptosi double %27 to i32
%29 = sdiv i32 %28, 6
%30 = add nsw i32 %29, 1
store i32 %30, i32* %11
%31 = load i32, i32* %11
%32 = mul nsw i32 %31, 6
%33 = sub nsw i32 %32, 3
store i32 %33, i32* %10
br label %38
34: ; preds = %5
%35 = load i32, i32* %10
%36 = sdiv i32 %35, 6
%37 = add nsw i32 %36, 1
store i32 %37, i32* %11
br label %38
38: ; preds = %34, %26
%39 = load double, double* %8
%40 = fdiv double %39, 1.800000e+02
%41 = fmul double %40, 0x400921FB54442D18
store double %41, double* %12
%42 = load double, double* %9
%43 = load i32, i32* %10
%44 = sitofp i32 %43 to double
%45 = fsub double %42, %44
%46 = fdiv double %45, 1.800000e+02
%47 = fmul double %46, 0x400921FB54442D18
store double %47, double* %13
store double 0x415854C140000000, double* %14
store double 0x41583FDFE0000000, double* %15
store double 0x3F7B6A6750BF3640, double* %16
%48 = load double, double* %12
%49 = fptosi double %48 to i32
%50 = call double @sin(i32 %49)
store double %50, double* %17
%51 = load double, double* %12
%52 = call double @cos(double %51)
store double %52, double* %18
%53 = load double, double* %16
%54 = fsub double 1.000000e+00, %53
%55 = fmul double 0x415854C140000000, %54
%56 = load double, double* %16
%57 = load double, double* %17
%58 = fmul double %56, %57
%59 = load double, double* %17
%60 = fneg double %58
%61 = call double @llvm.fmuladd.f64(double %60, double %59, double 1.000000e+00)
%62 = fdiv double %55, %61
%63 = load double, double* %16
%64 = load double, double* %17
%65 = fmul double %63, %64
%66 = load double, double* %17
%67 = fneg double %65
%68 = call double @llvm.fmuladd.f64(double %67, double %66, double 1.000000e+00)
%69 = fptosi double %68 to i32
%70 = call double @sqrt(i32 %69)
%71 = fdiv double %62, %70
store double %71, double* %19
%72 = load double, double* %16
%73 = load double, double* %17
%74 = fmul double %72, %73
%75 = load double, double* %17
%76 = fneg double %74
%77 = call double @llvm.fmuladd.f64(double %76, double %75, double 1.000000e+00)
%78 = fptosi double %77 to i32
%79 = call double @sqrt(i32 %78)
%80 = fdiv double 0x415854C140000000, %79
store double %80, double* %20
%81 = load double, double* %20
%82 = load double, double* %19
%83 = fdiv double %81, %82
%84 = fsub double %83, 1.000000e+00
store double %84, double* %21
%85 = load double, double* %16
%86 = fsub double 1.000000e+00, %85
%87 = fmul double 0x415854C140000000, %86
%88 = load double, double* %12
%89 = load double, double* %12
%90 = fmul double 2.000000e+00, %89
%91 = fptosi double %90 to i32
%92 = call double @sin(i32 %91)
%93 = fmul double 0x3F64BC49A675F8BB, %92
%94 = fneg double %93
%95 = call double @llvm.fmuladd.f64(double 0x3FF014B08B1B243A, double %88, double %94)
%96 = load double, double* %12
%97 = fmul double 4.000000e+00, %96
%98 = fptosi double %97 to i32
%99 = call double @sin(i32 %98)
%100 = call double @llvm.fmuladd.f64(double 0x3EC647FC2941FEE4, double %99, double %95)
%101 = load double, double* %12
%102 = fmul double 6.000000e+00, %101
%103 = fptosi double %102 to i32
%104 = call double @sin(i32 %103)
%105 = call double @llvm.fmuladd.f64(double 0xBE2DCAF216A4A1C5, double %104, double %100)
%106 = fmul double %87, %105
store double %106, double* %22
%107 = load double, double* %12
%108 = call double @tan(double %107)
store double %108, double* %23
%109 = load double, double* %22
%110 = load double, double* %13
%111 = load double, double* %13
%112 = fmul double %110, %111
%113 = load double, double* %20
%114 = fmul double %112, %113
%115 = fdiv double %114, 2.000000e+00
%116 = load double, double* %17
%117 = fmul double %115, %116
%118 = load double, double* %18
%119 = call double @llvm.fmuladd.f64(double %117, double %118, double %109)
%120 = load double, double* %13
%121 = load double, double* %13
%122 = fmul double %120, %121
%123 = load double, double* %13
%124 = fmul double %122, %123
%125 = load double, double* %13
%126 = fmul double %124, %125
%127 = load double, double* %20
%128 = fmul double %126, %127
%129 = fdiv double %128, 2.400000e+01
%130 = load double, double* %17
%131 = fmul double %129, %130
%132 = load double, double* %18
%133 = fmul double %131, %132
%134 = load double, double* %18
%135 = fmul double %133, %134
%136 = load double, double* %18
%137 = fmul double %135, %136
%138 = load double, double* %23
%139 = load double, double* %23
%140 = fneg double %138
%141 = call double @llvm.fmuladd.f64(double %140, double %139, double 5.000000e+00)
%142 = load double, double* %21
%143 = fmul double 9.000000e+00, %142
%144 = load double, double* %21
%145 = call double @llvm.fmuladd.f64(double %143, double %144, double %141)
%146 = load double, double* %21
%147 = call double @llvm.fmuladd.f64(double 4.000000e+00, double %146, double %145)
%148 = call double @llvm.fmuladd.f64(double %137, double %147, double %119)
%149 = load double*, double** %6
store double %148, double* %149
%150 = load double, double* %13
%151 = load double, double* %20
%152 = fmul double %150, %151
%153 = load double, double* %18
%154 = load double, double* %13
%155 = load double, double* %13
%156 = fmul double %154, %155
%157 = load double, double* %13
%158 = fmul double %156, %157
%159 = load double, double* %20
%160 = fmul double %158, %159
%161 = fdiv double %160, 6.000000e+00
%162 = load double, double* %18
%163 = fmul double %161, %162
%164 = load double, double* %18
%165 = fmul double %163, %164
%166 = load double, double* %18
%167 = fmul double %165, %166
%168 = load double, double* %23
%169 = load double, double* %23
%170 = fneg double %168
%171 = call double @llvm.fmuladd.f64(double %170, double %169, double 1.000000e+00)
%172 = load double, double* %21
%173 = fadd double %171, %172
%174 = fmul double %167, %173
%175 = call double @llvm.fmuladd.f64(double %152, double %153, double %174)
%176 = load double, double* %13
%177 = load double, double* %13
%178 = fmul double %176, %177
%179 = load double, double* %13
%180 = fmul double %178, %179
%181 = load double, double* %13
%182 = fmul double %180, %181
%183 = load double, double* %13
%184 = fmul double %182, %183
%185 = load double, double* %20
%186 = fmul double %184, %185
%187 = fdiv double %186, 1.200000e+02
%188 = load double, double* %18
%189 = fmul double %187, %188
%190 = load double, double* %18
%191 = fmul double %189, %190
%192 = load double, double* %18
%193 = fmul double %191, %192
%194 = load double, double* %18
%195 = fmul double %193, %194
%196 = load double, double* %18
%197 = fmul double %195, %196
%198 = load double, double* %23
%199 = fmul double 1.800000e+01, %198
%200 = load double, double* %23
%201 = fneg double %199
%202 = call double @llvm.fmuladd.f64(double %201, double %200, double 5.000000e+00)
%203 = load double, double* %23
%204 = load double, double* %23
%205 = fmul double %203, %204
%206 = load double, double* %23
%207 = fmul double %205, %206
%208 = load double, double* %23
%209 = call double @llvm.fmuladd.f64(double %207, double %208, double %202)
%210 = call double @llvm.fmuladd.f64(double %197, double %209, double %175)
%211 = load double*, double** %7
store double %210, double* %211
%212 = load double*, double** %7
%213 = load double, double* %212
%214 = fadd double %213, 5.000000e+05
%215 = load i32, i32* %11
%216 = sitofp i32 %215 to double
%217 = call double @llvm.fmuladd.f64(double %216, double 1.000000e+06, double %214)
%218 = load double*, double** %7
store double %217, double* %218
ret void
}

LLVM IR (-O2):
define void @SubGussFs(double* nocapture writeonly %0, double* nocapture writeonly %1, double %2, double %3, i32 %4) {
%6 = icmp eq i32 %4, 0
br i1 %6, label %7, label %13
7: ; preds = %5
%8 = fptosi double %3 to i32
%9 = sdiv i32 %8, 6
%10 = add nsw i32 %9, 1
%11 = mul nsw i32 %10, 6
%12 = add nsw i32 %11, -3
br label %16
13: ; preds = %5
%14 = sdiv i32 %4, 6
%15 = add nsw i32 %14, 1
br label %16
16: ; preds = %13, %7
%17 = phi i32 [ %10, %7 ], [ %15, %13 ]
%18 = phi i32 [ %12, %7 ], [ %4, %13 ]
%19 = fdiv double %2, 1.800000e+02
%20 = fmul double %19, 0x400921FB54442D18
%21 = sitofp i32 %18 to double
%22 = fsub double %3, %21
%23 = fdiv double %22, 1.800000e+02
%24 = fmul double %23, 0x400921FB54442D18
%25 = fptosi double %20 to i32
%26 = tail call double @sin(i32 %25)
%27 = tail call double @cos(double %20)
%28 = fmul double %26, 0xBF7B6A6750BF3640
%29 = tail call double @llvm.fmuladd.f64(double %28, double %26, double 1.000000e+00)
%30 = fdiv double 0x41582B106B475168, %29
%31 = fptosi double %29 to i32
%32 = tail call double @sqrt(i32 %31)
%33 = fdiv double %30, %32
%34 = tail call double @sqrt(i32 %31)
%35 = fdiv double 0x415854C140000000, %34
%36 = fdiv double %35, %33
%37 = fadd double %36, -1.000000e+00
%38 = fmul double %20, 2.000000e+00
%39 = fptosi double %38 to i32
%40 = tail call double @sin(i32 %39)
%41 = fmul double %40, 0xBF64BC49A675F8BB
%42 = tail call double @llvm.fmuladd.f64(double %20, double 0x3FF014B08B1B243A, double %41)
%43 = fmul double %20, 4.000000e+00
%44 = fptosi double %43 to i32
%45 = tail call double @sin(i32 %44)
%46 = tail call double @llvm.fmuladd.f64(double %45, double 0x3EC647FC2941FEE4, double %42)
%47 = fmul double %20, 6.000000e+00
%48 = fptosi double %47 to i32
%49 = tail call double @sin(i32 %48)
%50 = tail call double @llvm.fmuladd.f64(double %49, double 0xBE2DCAF216A4A1C5, double %46)
%51 = fmul double %50, 0x41582B106B475168
%52 = tail call double @tan(double %20)
%53 = fmul double %24, %24
%54 = fmul double %53, %35
%55 = fmul double %54, 5.000000e-01
%56 = fmul double %26, %55
%57 = tail call double @llvm.fmuladd.f64(double %56, double %27, double %51)
%58 = fmul double %24, %53
%59 = fmul double %24, %58
%60 = fmul double %35, %59
%61 = fdiv double %60, 2.400000e+01
%62 = fmul double %26, %61
%63 = fmul double %27, %62
%64 = fmul double %27, %63
%65 = fmul double %27, %64
%66 = fneg double %52
%67 = tail call double @llvm.fmuladd.f64(double %66, double %52, double 5.000000e+00)
%68 = fmul double %37, 9.000000e+00
%69 = tail call double @llvm.fmuladd.f64(double %68, double %37, double %67)
%70 = tail call double @llvm.fmuladd.f64(double %37, double 4.000000e+00, double %69)
%71 = tail call double @llvm.fmuladd.f64(double %65, double %70, double %57)
store double %71, double* %0
%72 = fmul double %24, %35
%73 = fmul double %35, %58
%74 = fdiv double %73, 6.000000e+00
%75 = fmul double %27, %74
%76 = fmul double %27, %75
%77 = fmul double %27, %76
%78 = tail call double @llvm.fmuladd.f64(double %66, double %52, double 1.000000e+00)
%79 = fadd double %37, %78
%80 = fmul double %79, %77
%81 = tail call double @llvm.fmuladd.f64(double %72, double %27, double %80)
%82 = fmul double %24, %59
%83 = fmul double %35, %82
%84 = fdiv double %83, 1.200000e+02
%85 = fmul double %27, %84
%86 = fmul double %27, %85
%87 = fmul double %27, %86
%88 = fmul double %27, %87
%89 = fmul double %27, %88
%90 = fmul double %52, -1.800000e+01
%91 = tail call double @llvm.fmuladd.f64(double %90, double %52, double 5.000000e+00)
%92 = fmul double %52, %52
%93 = fmul double %52, %92
%94 = tail call double @llvm.fmuladd.f64(double %93, double %52, double %91)
%95 = tail call double @llvm.fmuladd.f64(double %89, double %94, double %81)
%96 = fadd double %95, 5.000000e+05
%97 = sitofp i32 %17 to double
%98 = tail call double @llvm.fmuladd.f64(double %97, double 1.000000e+06, double %96)
store double %98, double* %1
ret void
}