#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 binomial(int, int);
double exp(double);
double fact2(int);
double pow(int, int);
double sqrt(double);
double three_center_1D(double xi, int ai, double alphai, double xj, int aj,
                       double alphaj, double xk, int ak, double alphak) {
  double gamma, dx, px, xpi, xpj, xpk, intgl;
  int q, r, s, n;
  gamma = alphai + alphaj + alphak;
  dx = exp(-alphai * alphaj * pow(xi - xj, 2) / gamma) *
       exp(-alphai * alphak * pow(xi - xk, 2) / gamma) *
       exp(-alphaj * alphak * pow(xj - xk, 2) / gamma);
  px = (alphai * xi + alphaj * xj + alphak * xk) / gamma;
  xpi = px - xi;
  xpj = px - xj;
  xpk = px - xk;
  intgl = 0;
  for (q = 0; q < ai + 1; q++) {
    for (r = 0; r < aj + 1; r++) {
      for (s = 0; s < ak + 1; s++) {
        n = (q + r + s) / 2;
        if ((q + r + s) % 2 == 0) {
          intgl += binomial(ai, q) * binomial(aj, r) * binomial(ak, s) *
                   pow(xpi, ai - q) * pow(xpj, aj - r) * pow(xpk, ak - s) *
                   fact2(2 * n - 1) / pow(2 * gamma, n) *
                   sqrt(3.14159265358979323846 / gamma);
        }
      }
    }
  }
  return dx * intgl;
}


x86 assembly (clang, -O0):
<three_center_1D>:
push   %rbp
mov    %rsp,%rbp
sub    $0x100,%rsp
movsd  %xmm0,-0x8(%rbp)
mov    %edi,-0xc(%rbp)
movsd  %xmm1,-0x18(%rbp)
movsd  %xmm2,-0x20(%rbp)
mov    %esi,-0x24(%rbp)
movsd  %xmm3,-0x30(%rbp)
movsd  %xmm4,-0x38(%rbp)
mov    %edx,-0x3c(%rbp)
movsd  %xmm5,-0x48(%rbp)
movsd  -0x18(%rbp),%xmm0
addsd  -0x30(%rbp),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x18(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0x8(%rbp),%xmm0
subsd  -0x20(%rbp),%xmm0
cvttsd2si %xmm0,%edi
mov    $0x2,%esi
call   1060 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0xb8(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
divsd  -0x50(%rbp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0x18(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0x8(%rbp),%xmm0
subsd  -0x38(%rbp),%xmm0
cvttsd2si %xmm0,%edi
mov    $0x2,%esi
call   1060 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0xb0(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
divsd  -0x50(%rbp),%xmm0
call   1050 <exp@plt>
movaps %xmm0,%xmm1
movsd  -0xa8(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x30(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x20(%rbp),%xmm0
subsd  -0x38(%rbp),%xmm0
cvttsd2si %xmm0,%edi
mov    $0x2,%esi
call   1060 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0xa0(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
divsd  -0x50(%rbp),%xmm0
call   1050 <exp@plt>
movaps %xmm0,%xmm1
movsd  -0x98(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x18(%rbp),%xmm1
movsd  -0x8(%rbp),%xmm2
movsd  -0x30(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x48(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
divsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
subsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x60(%rbp),%xmm0
subsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x60(%rbp),%xmm0
subsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x80(%rbp)
movl   $0x0,-0x84(%rbp)
mov    -0x84(%rbp),%eax
mov    -0xc(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    1594 <three_center_1D+0x444>
movl   $0x0,-0x88(%rbp)
mov    -0x88(%rbp),%eax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    157b <three_center_1D+0x42b>
movl   $0x0,-0x8c(%rbp)
mov    -0x8c(%rbp),%eax
mov    -0x3c(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    1562 <three_center_1D+0x412>
mov    -0x84(%rbp),%eax
add    -0x88(%rbp),%eax
add    -0x8c(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x90(%rbp)
mov    -0x84(%rbp),%eax
add    -0x88(%rbp),%eax
add    -0x8c(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
jne    1549 <three_center_1D+0x3f9>
mov    -0xc(%rbp),%edi
mov    -0x84(%rbp),%esi
call   1030 <binomial@plt>
movsd  %xmm0,-0xf8(%rbp)
ff
mov    -0x24(%rbp),%edi
mov    -0x88(%rbp),%esi
call   1030 <binomial@plt>
movaps %xmm0,%xmm1
movsd  -0xf8(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xf0(%rbp)
ff
mov    -0x3c(%rbp),%edi
mov    -0x8c(%rbp),%esi
call   1030 <binomial@plt>
movaps %xmm0,%xmm1
movsd  -0xf0(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xe8(%rbp)
ff
cvttsd2si -0x68(%rbp),%edi
mov    -0xc(%rbp),%esi
sub    -0x84(%rbp),%esi
call   1060 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0xe8(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
cvttsd2si -0x70(%rbp),%edi
mov    -0x24(%rbp),%esi
sub    -0x88(%rbp),%esi
call   1060 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0xe0(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
cvttsd2si -0x78(%rbp),%edi
mov    -0x3c(%rbp),%esi
sub    -0x8c(%rbp),%esi
call   1060 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0xd8(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xd0(%rbp)
ff
mov    -0x90(%rbp),%edi
shl    $0x1,%edi
sub    $0x1,%edi
call   1040 <fact2@plt>
movaps %xmm0,%xmm1
movsd  -0xd0(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
movsd  0xb19(%rip),%xmm0
mulsd  -0x50(%rbp),%xmm0
cvttsd2si %xmm0,%edi
mov    -0x90(%rbp),%esi
call   1060 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0xc8(%rbp),%xmm0
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
movsd  0xade(%rip),%xmm0
divsd  -0x50(%rbp),%xmm0
call   1070 <sqrt@plt>
movaps %xmm0,%xmm2
movsd  -0xc0(%rbp),%xmm0
ff
movsd  -0x80(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
jmp    154e <three_center_1D+0x3fe>
mov    -0x8c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x8c(%rbp)
jmp    138a <three_center_1D+0x23a>
jmp    1567 <three_center_1D+0x417>
mov    -0x88(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x88(%rbp)
jmp    136c <three_center_1D+0x21c>
jmp    1580 <three_center_1D+0x430>
mov    -0x84(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x84(%rbp)
jmp    134e <three_center_1D+0x1fe>
movsd  -0x58(%rbp),%xmm0
mulsd  -0x80(%rbp),%xmm0
add    $0x100,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<three_center_1D>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xd8,%rsp
movsd  %xmm5,0x18(%rsp)
mov    %edx,%r14d
movsd  %xmm4,0x10(%rsp)
movapd %xmm3,0x90(%rsp)
mov    %esi,0x8(%rsp)
movsd  %xmm2,0x38(%rsp)
mov    %edi,%ebx
movsd  %xmm0,0x30(%rsp)
movapd %xmm1,%xmm4
addsd  %xmm3,%xmm4
addsd  %xmm5,%xmm4
movapd %xmm4,0x70(%rsp)
movapd %xmm1,0x80(%rsp)
xorpd  0xe52(%rip),%xmm1
movapd %xmm1,0x40(%rsp)
mulsd  %xmm3,%xmm1
movapd %xmm1,0xa0(%rsp)
movapd %xmm0,%xmm1
subsd  %xmm2,%xmm1
cvttsd2si %xmm1,%edi
mov    $0x2,%esi
call   1060 <pow@plt>
mulsd  0xa0(%rsp),%xmm0
divsd  0x70(%rsp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,0xc0(%rsp)
movapd 0x40(%rsp),%xmm0
mulsd  0x18(%rsp),%xmm0
movapd %xmm0,0x40(%rsp)
movsd  0x30(%rsp),%xmm0
subsd  0x10(%rsp),%xmm0
cvttsd2si %xmm0,%edi
mov    $0x2,%esi
call   1060 <pow@plt>
mulsd  0x40(%rsp),%xmm0
divsd  0x70(%rsp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,0xb8(%rsp)
movapd 0xdbe(%rip),%xmm0
xorpd  0x90(%rsp),%xmm0
mulsd  0x18(%rsp),%xmm0
movapd %xmm0,0x40(%rsp)
movsd  0x38(%rsp),%xmm0
subsd  0x10(%rsp),%xmm0
cvttsd2si %xmm0,%edi
mov    $0x2,%esi
call   1060 <pow@plt>
mulsd  0x40(%rsp),%xmm0
divsd  0x70(%rsp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,0xb0(%rsp)
mov    %ebx,0xc(%rsp)
test   %ebx,%ebx
js     1504 <three_center_1D+0x3b4>
movsd  0x38(%rsp),%xmm5
movapd 0x90(%rsp),%xmm1
mulsd  %xmm5,%xmm1
movsd  0x30(%rsp),%xmm6
movapd 0x80(%rsp),%xmm0
mulsd  %xmm6,%xmm0
addsd  %xmm1,%xmm0
movsd  0x18(%rsp),%xmm2
movsd  0x10(%rsp),%xmm3
mulsd  %xmm3,%xmm2
addsd  %xmm0,%xmm2
movapd 0x70(%rsp),%xmm4
divsd  %xmm4,%xmm2
movapd %xmm2,%xmm0
movapd %xmm2,%xmm1
subsd  %xmm5,%xmm1
subsd  %xmm6,%xmm2
cvttsd2si %xmm2,%eax
mov    %eax,0x64(%rsp)
cvttsd2si %xmm1,%eax
mov    %eax,0x60(%rsp)
subsd  %xmm3,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,0x5c(%rsp)
movapd %xmm4,%xmm0
addsd  %xmm4,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,0x58(%rsp)
movsd  0xced(%rip),%xmm0
divsd  %xmm4,%xmm0
xor    %eax,%eax
mov    %rax,0x28(%rsp)
movsd  %xmm0,0x80(%rsp)
sqrtsd %xmm0,%xmm0
movsd  %xmm0,0xd0(%rsp)
xorpd  %xmm1,%xmm1
mov    %r14d,0x24(%rsp)
jmp    1367 <three_center_1D+0x217>
nop
mov    0x28(%rsp),%rcx
lea    0x1(%rcx),%eax
cmp    0xc(%rsp),%ecx
mov    %rax,0x28(%rsp)
je     1508 <three_center_1D+0x3b8>
cmpl   $0x0,0x8(%rsp)
js     1350 <three_center_1D+0x200>
mov    0xc(%rsp),%eax
mov    0x28(%rsp),%rbp
sub    %ebp,%eax
mov    %eax,0x68(%rsp)
xor    %r13d,%r13d
jmp    13a6 <three_center_1D+0x256>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
lea    0x1(%r13),%eax
add    $0x1,%ebp
cmp    0x8(%rsp),%r13d
mov    %eax,%r13d
mov    0x24(%rsp),%r14d
je     1350 <three_center_1D+0x200>
test   %r14d,%r14d
js     1390 <three_center_1D+0x240>
mov    0x8(%rsp),%eax
sub    %r13d,%eax
mov    %eax,0x6c(%rsp)
xor    %ebx,%ebx
jmp    140b <three_center_1D+0x2bb>
nopw   0x0(%rax,%rax,1)
movsd  0x40(%rsp),%xmm2
mulsd  0x18(%rsp),%xmm2
mulsd  0x10(%rsp),%xmm2
mulsd  0x38(%rsp),%xmm2
mulsd  0x30(%rsp),%xmm2
mulsd  0x90(%rsp),%xmm2
mulsd  0xa0(%rsp),%xmm2
divsd  %xmm1,%xmm2
mulsd  %xmm0,%xmm2
movsd  0x70(%rsp),%xmm1
addsd  %xmm2,%xmm1
add    $0x1,%ebx
add    $0xffffffff,%r14d
jae    1390 <three_center_1D+0x240>
lea    (%rbx,%rbp,1),%r15d
test   $0x1,%r15b
jne    1402 <three_center_1D+0x2b2>
movsd  %xmm1,0x70(%rsp)
mov    %r15d,%r12d
shr    %r12d
mov    0xc(%rsp),%edi
mov    0x28(%rsp),%rsi
call   1030 <binomial@plt>
movsd  %xmm0,0x40(%rsp)
mov    0x8(%rsp),%edi
mov    %r13d,%esi
call   1030 <binomial@plt>
movsd  %xmm0,0x18(%rsp)
mov    0x24(%rsp),%edi
mov    %ebx,%esi
call   1030 <binomial@plt>
movsd  %xmm0,0x10(%rsp)
mov    0x64(%rsp),%edi
mov    0x68(%rsp),%esi
call   1060 <pow@plt>
movsd  %xmm0,0x38(%rsp)
mov    0x60(%rsp),%edi
mov    0x6c(%rsp),%esi
call   1060 <pow@plt>
movsd  %xmm0,0x30(%rsp)
mov    0x5c(%rsp),%edi
mov    %r14d,%esi
call   1060 <pow@plt>
movsd  %xmm0,0x90(%rsp)
and    $0x7ffffffe,%r15d
add    $0xffffffff,%r15d
mov    %r15d,%edi
call   1040 <fact2@plt>
movsd  %xmm0,0xa0(%rsp)
mov    0x58(%rsp),%edi
mov    %r12d,%esi
call   1060 <pow@plt>
movapd %xmm0,%xmm1
movsd  0x80(%rsp),%xmm0
ucomisd 0xb48(%rip),%xmm0
movsd  0xd0(%rsp),%xmm0
jae    13c0 <three_center_1D+0x270>
movsd  0x80(%rsp),%xmm0
movsd  %xmm1,0xc8(%rsp)
call   1070 <sqrt@plt>
movsd  0xc8(%rsp),%xmm1
jmp    13c0 <three_center_1D+0x270>
xorpd  %xmm1,%xmm1
movsd  0xc0(%rsp),%xmm0
mulsd  0xb8(%rsp),%xmm0
mulsd  0xb0(%rsp),%xmm0
mulsd  %xmm1,%xmm0
add    $0xd8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define double @three_center_1D(double %0, i32 %1, double %2, double %3, i32 %4, double %5, double %6, i32 %7, double %8) {
%10 = fadd double %2, %5
%11 = fadd double %10, %8
%12 = fneg double %2
%13 = fmul double %12, %5
%14 = fsub double %0, %3
%15 = fptosi double %14 to i32
%16 = tail call double @pow(i32 %15, i32 2)
%17 = fmul double %13, %16
%18 = fdiv double %17, %11
%19 = tail call double @exp(double %18)
%20 = fmul double %12, %8
%21 = fsub double %0, %6
%22 = fptosi double %21 to i32
%23 = tail call double @pow(i32 %22, i32 2)
%24 = fmul double %20, %23
%25 = fdiv double %24, %11
%26 = tail call double @exp(double %25)
%27 = fneg double %5
%28 = fmul double %27, %8
%29 = fsub double %3, %6
%30 = fptosi double %29 to i32
%31 = tail call double @pow(i32 %30, i32 2)
%32 = fmul double %28, %31
%33 = fdiv double %32, %11
%34 = tail call double @exp(double %33)
%35 = fmul double %3, %5
%36 = tail call double @llvm.fmuladd.f64(double %2, double %0, double %35)
%37 = tail call double @llvm.fmuladd.f64(double %8, double %6, double %36)
%38 = fdiv double %37, %11
%39 = icmp slt i32 %1, 0
br i1 %39, label %103, label %40
40: ; preds = %9
%41 = fsub double %38, %6
%42 = fsub double %38, %3
%43 = fsub double %38, %0
%44 = icmp slt i32 %4, 0
%45 = icmp slt i32 %7, 0
%46 = fptosi double %43 to i32
%47 = fptosi double %42 to i32
%48 = fptosi double %41 to i32
%49 = fmul double %11, 2.000000e+00
%50 = fptosi double %49 to i32
%51 = fdiv double 0x400921FB54442D18, %11
br label %52
52: ; preds = %40, %99
%53 = phi double [ 0.000000e+00, %40 ], [ %100, %99 ]
%54 = phi i32 [ 0, %40 ], [ %101, %99 ]
br i1 %44, label %99, label %55
55: ; preds = %52
%56 = sub nsw i32 %1, %54
br label %57
57: ; preds = %55, %95
%58 = phi double [ %53, %55 ], [ %96, %95 ]
%59 = phi i32 [ 0, %55 ], [ %97, %95 ]
br i1 %45, label %95, label %60
60: ; preds = %57
%61 = add nuw nsw i32 %59, %54
%62 = sub nsw i32 %4, %59
br label %63
63: ; preds = %60, %91
%64 = phi double [ %58, %60 ], [ %92, %91 ]
%65 = phi i32 [ 0, %60 ], [ %93, %91 ]
%66 = add nuw nsw i32 %61, %65
%67 = and i32 %66, 1
%68 = icmp eq i32 %67, 0
br i1 %68, label %69, label %91
69: ; preds = %63
%70 = lshr i32 %66, 1
%71 = tail call double @binomial(i32 %1, i32 %54)
%72 = tail call double @binomial(i32 %4, i32 %59)
%73 = fmul double %71, %72
%74 = tail call double @binomial(i32 %7, i32 %65)
%75 = fmul double %73, %74
%76 = tail call double @pow(i32 %46, i32 %56)
%77 = fmul double %75, %76
%78 = tail call double @pow(i32 %47, i32 %62)
%79 = fmul double %77, %78
%80 = sub nsw i32 %7, %65
%81 = tail call double @pow(i32 %48, i32 %80)
%82 = fmul double %79, %81
%83 = and i32 %66, 2147483646
%84 = add nsw i32 %83, -1
%85 = tail call double @fact2(i32 %84)
%86 = fmul double %82, %85
%87 = tail call double @pow(i32 %50, i32 %70)
%88 = fdiv double %86, %87
%89 = tail call double @sqrt(double %51)
%90 = tail call double @llvm.fmuladd.f64(double %88, double %89, double %64)
br label %91
91: ; preds = %63, %69
%92 = phi double [ %90, %69 ], [ %64, %63 ]
%93 = add nuw i32 %65, 1
%94 = icmp eq i32 %65, %7
br i1 %94, label %95, label %63
95: ; preds = %91, %57
%96 = phi double [ %58, %57 ], [ %92, %91 ]
%97 = add nuw i32 %59, 1
%98 = icmp eq i32 %59, %4
br i1 %98, label %99, label %57
99: ; preds = %95, %52
%100 = phi double [ %53, %52 ], [ %96, %95 ]
%101 = add nuw i32 %54, 1
%102 = icmp eq i32 %54, %1
br i1 %102, label %103, label %52
103: ; preds = %99, %9
%104 = phi double [ 0.000000e+00, %9 ], [ %100, %99 ]
%105 = fmul double %19, %26
%106 = fmul double %105, %34
%107 = fmul double %106, %104
ret double %107
}