#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 exp(double);
int fabs(double);
int karlin(int low, int high, double *pr, double *lambda_p, double *H_p) {
  int i, range, nit;
  double up, new, sum, av, beta, ftemp;
  double lambda;
  double *ptr1;
  range = high - low;
  sum = 0;
  ptr1 = pr;
  for (i = low; i <= high; i++)
    sum += i * (*ptr1++);
  if (sum >= 0.0) {
    return 0;
  }
  up = 0.5;
  do {
    up *= 2.0;
    ptr1 = pr;
    beta = exp(up);
    ftemp = exp(up * (low - 1));
    sum = 0.0;
    for (i = 0; i <= range; ++i)
      sum += *ptr1++ * (ftemp *= beta);
  } while (sum < 1.0);
  lambda = 0.0;
  nit = 0;
  while (nit++ < 25) {
    new = (lambda + up) / 2.0;
    beta = exp(new);
    ftemp = exp(new *(low - 1));
    ptr1 = pr;
    sum = 0.0;
    for (i = 0; i <= range; ++i)
      sum += *ptr1++ * (ftemp *= beta);
    if (sum > 1.0 + 1e-6)
      up = new;
    else {
      if (fabs(lambda - new) < 1e-6)
        goto done;
      lambda = new;
    }
  }
  if (lambda <= 1e-10) {
    lambda = -1.0;
    return 0;
  }
done:
  *lambda_p = lambda;
  ptr1 = pr;
  ftemp = exp(lambda * (low - 1));
  for (av = 0.0, i = low; i <= high; ++i)
    av += *ptr1++ * i * (ftemp *= beta);
  *H_p = lambda * av;
  return 1;
}


x86 assembly (clang, -O0):
<karlin>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %edi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0xc(%rbp),%eax
sub    -0x8(%rbp),%eax
mov    %eax,-0x30(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x18(%rbp),%rax
mov    %rax,-0x78(%rbp)
mov    -0x8(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     11a0 <karlin+0x80>
cvtsi2sdl -0x2c(%rbp),%xmm0
mov    -0x78(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x78(%rbp)
movsd  (%rax),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    115c <karlin+0x3c>
movsd  -0x50(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     11be <karlin+0x9e>
movl   $0x0,-0x4(%rbp)
jmp    1495 <karlin+0x375>
movsd  0xe3a(%rip),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  0xe35(%rip),%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x18(%rbp),%rax
mov    %rax,-0x78(%rbp)
movsd  -0x40(%rbp),%xmm0
call   1040 <exp@plt>
movsd  %xmm0,-0x60(%rbp)
movsd  -0x40(%rbp),%xmm0
mov    -0x8(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm1,%xmm0
call   1040 <exp@plt>
movsd  %xmm0,-0x68(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x30(%rbp),%eax
jg     126e <karlin+0x14e>
mov    -0x78(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x78(%rbp)
movsd  (%rax),%xmm0
movsd  -0x60(%rbp),%xmm2
mulsd  -0x68(%rbp),%xmm2
movsd  %xmm2,-0x68(%rbp)
movsd  -0x50(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1220 <karlin+0x100>
jmp    1273 <karlin+0x153>
movsd  0xd95(%rip),%xmm0
ucomisd -0x50(%rbp),%xmm0
ja     11cb <karlin+0xab>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x70(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    %eax,%ecx
add    $0x1,%ecx
mov    %ecx,-0x34(%rbp)
cmp    $0x19,%eax
jge    13b4 <karlin+0x294>
movsd  -0x70(%rbp),%xmm0
addsd  -0x40(%rbp),%xmm0
movsd  0xd4d(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
call   1040 <exp@plt>
movsd  %xmm0,-0x60(%rbp)
movsd  -0x48(%rbp),%xmm0
mov    -0x8(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm1,%xmm0
call   1040 <exp@plt>
movsd  %xmm0,-0x68(%rbp)
mov    -0x18(%rbp),%rax
mov    %rax,-0x78(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x30(%rbp),%eax
jg     1355 <karlin+0x235>
mov    -0x78(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x78(%rbp)
movsd  (%rax),%xmm0
movsd  -0x60(%rbp),%xmm2
mulsd  -0x68(%rbp),%xmm2
movsd  %xmm2,-0x68(%rbp)
movsd  -0x50(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1307 <karlin+0x1e7>
movsd  -0x50(%rbp),%xmm0
movsd  0xcc6(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    137b <karlin+0x25b>
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
jmp    13af <karlin+0x28f>
movsd  -0x70(%rbp),%xmm0
subsd  -0x48(%rbp),%xmm0
call   1030 <fabs@plt>
cvtsi2sd %eax,%xmm1
movsd  0xc9a(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    13a5 <karlin+0x285>
jmp    13e5 <karlin+0x2c5>
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
jmp    1295 <karlin+0x175>
movsd  0xc5c(%rip),%xmm0
ucomisd -0x70(%rbp),%xmm0
jb     13e0 <karlin+0x2c0>
movsd  0xc51(%rip),%xmm0
movsd  %xmm0,-0x70(%rbp)
movl   $0x0,-0x4(%rbp)
jmp    1495 <karlin+0x375>
jmp    13e5 <karlin+0x2c5>
movsd  -0x70(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x18(%rbp),%rax
mov    %rax,-0x78(%rbp)
movsd  -0x70(%rbp),%xmm0
mov    -0x8(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm1,%xmm0
call   1040 <exp@plt>
movsd  %xmm0,-0x68(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x8(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     147c <karlin+0x35c>
mov    -0x78(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x78(%rbp)
movsd  (%rax),%xmm0
cvtsi2sdl -0x2c(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x60(%rbp),%xmm2
mulsd  -0x68(%rbp),%xmm2
movsd  %xmm2,-0x68(%rbp)
movsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1425 <karlin+0x305>
movsd  -0x70(%rbp),%xmm0
mulsd  -0x58(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x1,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<karlin>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
mov    %edi,%ebx
xor    %eax,%eax
mov    %esi,%edi
sub    %ebx,%edi
jl     1578 <karlin+0x458>
mov    %rdx,%r15
mov    %esi,%r13d
lea    0x1(%rdi),%esi
mov    %esi,0x4(%rsp)
and    $0x3,%esi
mov    %rcx,0x40(%rsp)
xorpd  %xmm0,%xmm0
mov    %rdx,%rcx
mov    %ebx,%edx
je     117b <karlin+0x5b>
nopl   0x0(%rax)
xorps  %xmm1,%xmm1
cvtsi2sd %edx,%xmm1
mulsd  (%rcx),%xmm1
add    $0x8,%rcx
addsd  %xmm1,%xmm0
add    $0x1,%edx
add    $0xffffffff,%esi
jne    1160 <karlin+0x40>
mov    %rdi,0x20(%rsp)
cmp    $0x3,%edi
jb     11e4 <karlin+0xc4>
cs nopw 0x0(%rax,%rax,1)
nop
xorps  %xmm1,%xmm1
cvtsi2sd %edx,%xmm1
mulsd  (%rcx),%xmm1
addsd  %xmm0,%xmm1
lea    0x1(%rdx),%esi
xorps  %xmm0,%xmm0
cvtsi2sd %esi,%xmm0
mulsd  0x8(%rcx),%xmm0
addsd  %xmm1,%xmm0
lea    0x2(%rdx),%esi
xorps  %xmm1,%xmm1
cvtsi2sd %esi,%xmm1
mulsd  0x10(%rcx),%xmm1
addsd  %xmm0,%xmm1
lea    0x3(%rdx),%esi
xorps  %xmm0,%xmm0
cvtsi2sd %esi,%xmm0
mulsd  0x18(%rcx),%xmm0
addsd  %xmm1,%xmm0
add    $0x4,%edx
add    $0x20,%rcx
cmp    %r13d,%esi
jne    1190 <karlin+0x70>
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jae    1578 <karlin+0x458>
mov    %r8,0x38(%rsp)
lea    -0x1(%rbx),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,0x28(%rsp)
mov    0x4(%rsp),%r14d
and    $0xfffffffc,%r14d
mov    %r13d,%eax
sub    %bl,%al
add    $0x1,%al
movzbl %al,%ebp
and    $0x3,%ebp
movsd  0xddb(%rip),%xmm0
movsd  %xmm0,0x18(%rsp)
mov    0x20(%rsp),%r12
jmp    1252 <karlin+0x132>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  0xdc0(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    132e <karlin+0x20e>
movsd  0x18(%rsp),%xmm0
addsd  %xmm0,%xmm0
movsd  %xmm0,0x18(%rsp)
call   1040 <exp@plt>
movsd  %xmm0,0x8(%rsp)
movsd  0x18(%rsp),%xmm0
mulsd  0x28(%rsp),%xmm0
call   1040 <exp@plt>
movsd  0x8(%rsp),%xmm4
xorpd  %xmm1,%xmm1
cmp    %ebx,%r13d
jl     1240 <karlin+0x120>
xorpd  %xmm1,%xmm1
cmp    $0x3,%r12d
jae    12a0 <karlin+0x180>
mov    %r15,%rax
jmp    12fc <karlin+0x1dc>
nopl   0x0(%rax)
xorpd  %xmm1,%xmm1
mov    %r14d,%ecx
mov    %r15,%rax
nopw   0x0(%rax,%rax,1)
mulsd  %xmm4,%xmm0
movapd %xmm4,%xmm2
mulsd  %xmm0,%xmm2
mulsd  (%rax),%xmm0
addsd  %xmm1,%xmm0
movapd %xmm4,%xmm3
mulsd  %xmm2,%xmm3
mulsd  0x8(%rax),%xmm2
addsd  %xmm0,%xmm2
movapd %xmm4,%xmm0
mulsd  %xmm3,%xmm0
mulsd  0x10(%rax),%xmm3
addsd  %xmm2,%xmm3
movsd  0x18(%rax),%xmm1
mulsd  %xmm0,%xmm1
add    $0x20,%rax
addsd  %xmm3,%xmm1
add    $0xfffffffc,%ecx
jne    12b0 <karlin+0x190>
testb  $0x3,0x4(%rsp)
je     1240 <karlin+0x120>
xor    %ecx,%ecx
nopl   0x0(%rax)
mulsd  %xmm4,%xmm0
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm2,%xmm1
add    $0x1,%rcx
cmp    %ecx,%ebp
jne    1310 <karlin+0x1f0>
jmp    1240 <karlin+0x120>
mov    0x4(%rsp),%ebp
and    $0xfffffffc,%ebp
mov    %r13d,%eax
sub    %bl,%al
add    $0x1,%al
movzbl %al,%r12d
and    $0x3,%r12d
xorpd  %xmm2,%xmm2
mov    $0x1,%r14d
jmp    1381 <karlin+0x261>
movsd  0x10(%rsp),%xmm0
subsd  0x8(%rsp),%xmm0
call   1030 <fabs@plt>
movsd  0x8(%rsp),%xmm1
test   %eax,%eax
jle    14a9 <karlin+0x389>
add    $0x1,%r14d
movapd %xmm1,%xmm2
cmp    $0x1a,%r14d
je     148f <karlin+0x36f>
movsd  0x18(%rsp),%xmm0
movsd  %xmm2,0x10(%rsp)
addsd  %xmm2,%xmm0
mulsd  0xc67(%rip),%xmm0
movsd  %xmm0,0x8(%rsp)
call   1040 <exp@plt>
movsd  %xmm0,0x30(%rsp)
movsd  0x8(%rsp),%xmm0
mulsd  0x28(%rsp),%xmm0
call   1040 <exp@plt>
cmp    %ebx,%r13d
jl     1350 <karlin+0x230>
xorpd  %xmm1,%xmm1
cmpl   $0x3,0x20(%rsp)
movsd  0x30(%rsp),%xmm4
jae    13e0 <karlin+0x2c0>
mov    %r15,%rax
jmp    143c <karlin+0x31c>
cs nopw 0x0(%rax,%rax,1)
xorpd  %xmm1,%xmm1
mov    %ebp,%ecx
mov    %r15,%rax
nopl   0x0(%rax)
mulsd  %xmm4,%xmm0
movapd %xmm4,%xmm2
mulsd  %xmm0,%xmm2
mulsd  (%rax),%xmm0
addsd  %xmm1,%xmm0
movapd %xmm4,%xmm3
mulsd  %xmm2,%xmm3
mulsd  0x8(%rax),%xmm2
addsd  %xmm0,%xmm2
movapd %xmm4,%xmm0
mulsd  %xmm3,%xmm0
mulsd  0x10(%rax),%xmm3
addsd  %xmm2,%xmm3
movsd  0x18(%rax),%xmm1
mulsd  %xmm0,%xmm1
add    $0x20,%rax
addsd  %xmm3,%xmm1
add    $0xfffffffc,%ecx
jne    13f0 <karlin+0x2d0>
testb  $0x3,0x4(%rsp)
je     146a <karlin+0x34a>
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nop
mulsd  %xmm4,%xmm0
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm2,%xmm1
add    $0x1,%rcx
cmp    %ecx,%r12d
jne    1450 <karlin+0x330>
ucomisd 0xb9e(%rip),%xmm1
jbe    1350 <karlin+0x230>
movsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x18(%rsp)
movsd  0x10(%rsp),%xmm1
jmp    136f <karlin+0x24f>
xor    %eax,%eax
movsd  0xb7f(%rip),%xmm0
ucomisd %xmm1,%xmm0
movsd  %xmm1,0x10(%rsp)
jae    1578 <karlin+0x458>
mov    0x40(%rsp),%rax
movsd  0x10(%rsp),%xmm1
movsd  %xmm1,(%rax)
movsd  0x28(%rsp),%xmm0
mulsd  %xmm1,%xmm0
call   1040 <exp@plt>
xorpd  %xmm1,%xmm1
cmp    %ebx,%r13d
movsd  0x30(%rsp),%xmm3
jl     1560 <karlin+0x440>
testb  $0x1,0x4(%rsp)
jne    14e7 <karlin+0x3c7>
xorpd  %xmm1,%xmm1
jmp    150a <karlin+0x3ea>
xorps  %xmm2,%xmm2
cvtsi2sd %ebx,%xmm2
mulsd  (%r15),%xmm2
add    $0x8,%r15
mulsd  %xmm3,%xmm0
mulsd  %xmm0,%xmm2
xorpd  %xmm1,%xmm1
addsd  %xmm2,%xmm1
add    $0x1,%ebx
mov    0x20(%rsp),%rax
test   %eax,%eax
je     1560 <karlin+0x440>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
xorps  %xmm2,%xmm2
cvtsi2sd %ebx,%xmm2
mulsd  (%r15),%xmm2
mulsd  %xmm3,%xmm0
mulsd  %xmm0,%xmm2
addsd  %xmm1,%xmm2
lea    0x1(%rbx),%eax
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
mulsd  0x8(%r15),%xmm1
mulsd  %xmm3,%xmm0
mulsd  %xmm0,%xmm1
addsd  %xmm2,%xmm1
add    $0x10,%r15
add    $0x2,%ebx
cmp    %r13d,%eax
jne    1520 <karlin+0x400>
movsd  0x10(%rsp),%xmm0
mulsd  %xmm1,%xmm0
mov    0x38(%rsp),%rax
movsd  %xmm0,(%rax)
mov    $0x1,%eax
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define i32 @karlin(i32 %0, i32 %1, double* nocapture %2, double* nocapture writeonly %3, double* nocapture writeonly %4) {
%6 = icmp slt i32 %1, %0
br i1 %6, label %231, label %7
7: ; preds = %5
%8 = add i32 %1, 1
%9 = sub i32 %8, %0
%10 = sub i32 %1, %0
%11 = and i32 %9, 3
%12 = icmp eq i32 %11, 0
br i1 %12, label %25, label %13
13: ; preds = %7, %13
%14 = phi double* [ %19, %13 ], [ %2, %7 ]
%15 = phi double [ %21, %13 ], [ 0.000000e+00, %7 ]
%16 = phi i32 [ %22, %13 ], [ %0, %7 ]
%17 = phi i32 [ %23, %13 ], [ 0, %7 ]
%18 = sitofp i32 %16 to double
%19 = getelementptr inbounds double, double* %14, i64 1
%20 = load double, double* %14
%21 = tail call double @llvm.fmuladd.f64(double %18, double %20, double %15)
%22 = add i32 %16, 1
%23 = add i32 %17, 1
%24 = icmp eq i32 %23, %11
br i1 %24, label %25, label %13
25: ; preds = %13, %7
%26 = phi double [ undef, %7 ], [ %21, %13 ]
%27 = phi double* [ %2, %7 ], [ %19, %13 ]
%28 = phi double [ 0.000000e+00, %7 ], [ %21, %13 ]
%29 = phi i32 [ %0, %7 ], [ %22, %13 ]
%30 = icmp ult i32 %10, 3
br i1 %30, label %56, label %31
31: ; preds = %25, %31
%32 = phi double* [ %51, %31 ], [ %27, %25 ]
%33 = phi double [ %53, %31 ], [ %28, %25 ]
%34 = phi i32 [ %54, %31 ], [ %29, %25 ]
%35 = sitofp i32 %34 to double
%36 = getelementptr inbounds double, double* %32, i64 1
%37 = load double, double* %32
%38 = tail call double @llvm.fmuladd.f64(double %35, double %37, double %33)
%39 = add i32 %34, 1
%40 = sitofp i32 %39 to double
%41 = getelementptr inbounds double, double* %32, i64 2
%42 = load double, double* %36
%43 = tail call double @llvm.fmuladd.f64(double %40, double %42, double %38)
%44 = add i32 %34, 2
%45 = sitofp i32 %44 to double
%46 = getelementptr inbounds double, double* %32, i64 3
%47 = load double, double* %41
%48 = tail call double @llvm.fmuladd.f64(double %45, double %47, double %43)
%49 = add i32 %34, 3
%50 = sitofp i32 %49 to double
%51 = getelementptr inbounds double, double* %32, i64 4
%52 = load double, double* %46
%53 = tail call double @llvm.fmuladd.f64(double %50, double %52, double %48)
%54 = add i32 %34, 4
%55 = icmp eq i32 %49, %1
br i1 %55, label %56, label %31
56: ; preds = %31, %25
%57 = phi double [ %26, %25 ], [ %53, %31 ]
%58 = fcmp ult double %57, 0.000000e+00
br i1 %58, label %59, label %231
59: ; preds = %56
%60 = add nsw i32 %0, -1
%61 = sitofp i32 %60 to double
%62 = icmp slt i32 %1, %0
%63 = add i32 %1, 1
%64 = sub i32 %63, %0
%65 = and i32 %64, 3
%66 = icmp ult i32 %10, 3
%67 = and i32 %64, -4
%68 = icmp eq i32 %65, 0
br label %69
69: ; preds = %59, %115
%70 = phi double [ %71, %115 ], [ 5.000000e-01, %59 ]
%71 = fmul double %70, 2.000000e+00
%72 = tail call double @exp(double %71)
%73 = fmul double %71, %61
%74 = tail call double @exp(double %73)
br i1 %62, label %115, label %75
75: ; preds = %69
br i1 %66, label %99, label %76
76: ; preds = %75, %76
%77 = phi double* [ %93, %76 ], [ %2, %75 ]
%78 = phi double [ %95, %76 ], [ %74, %75 ]
%79 = phi double [ %96, %76 ], [ 0.000000e+00, %75 ]
%80 = phi i32 [ %97, %76 ], [ 0, %75 ]
%81 = getelementptr inbounds double, double* %77, i64 1
%82 = load double, double* %77
%83 = fmul double %72, %78
%84 = tail call double @llvm.fmuladd.f64(double %82, double %83, double %79)
%85 = getelementptr inbounds double, double* %77, i64 2
%86 = load double, double* %81
%87 = fmul double %72, %83
%88 = tail call double @llvm.fmuladd.f64(double %86, double %87, double %84)
%89 = getelementptr inbounds double, double* %77, i64 3
%90 = load double, double* %85
%91 = fmul double %72, %87
%92 = tail call double @llvm.fmuladd.f64(double %90, double %91, double %88)
%93 = getelementptr inbounds double, double* %77, i64 4
%94 = load double, double* %89
%95 = fmul double %72, %91
%96 = tail call double @llvm.fmuladd.f64(double %94, double %95, double %92)
%97 = add i32 %80, 4
%98 = icmp eq i32 %97, %67
br i1 %98, label %99, label %76
99: ; preds = %76, %75
%100 = phi double [ undef, %75 ], [ %96, %76 ]
%101 = phi double* [ %2, %75 ], [ %93, %76 ]
%102 = phi double [ %74, %75 ], [ %95, %76 ]
%103 = phi double [ 0.000000e+00, %75 ], [ %96, %76 ]
br i1 %68, label %115, label %104
104: ; preds = %99, %104
%105 = phi double* [ %109, %104 ], [ %101, %99 ]
%106 = phi double [ %111, %104 ], [ %102, %99 ]
%107 = phi double [ %112, %104 ], [ %103, %99 ]
%108 = phi i32 [ %113, %104 ], [ 0, %99 ]
%109 = getelementptr inbounds double, double* %105, i64 1
%110 = load double, double* %105
%111 = fmul double %72, %106
%112 = tail call double @llvm.fmuladd.f64(double %110, double %111, double %107)
%113 = add i32 %108, 1
%114 = icmp eq i32 %113, %65
br i1 %114, label %115, label %104
115: ; preds = %99, %104, %69
%116 = phi double [ 0.000000e+00, %69 ], [ %100, %99 ], [ %112, %104 ]
%117 = fcmp olt double %116, 1.000000e+00
br i1 %117, label %69, label %118
118: ; preds = %115
%119 = and i32 %64, 3
%120 = icmp ult i32 %10, 3
%121 = and i32 %64, -4
%122 = icmp eq i32 %119, 0
br label %123
123: ; preds = %118, %179
%124 = phi i32 [ %182, %179 ], [ 1, %118 ]
%125 = phi double [ %181, %179 ], [ 0.000000e+00, %118 ]
%126 = phi double [ %180, %179 ], [ %71, %118 ]
%127 = fadd double %126, %125
%128 = fmul double %127, 5.000000e-01
%129 = tail call double @exp(double %128)
%130 = fmul double %128, %61
%131 = tail call double @exp(double %130)
br i1 %62, label %175, label %132
132: ; preds = %123
br i1 %120, label %156, label %133
133: ; preds = %132, %133
%134 = phi double* [ %150, %133 ], [ %2, %132 ]
%135 = phi double [ %152, %133 ], [ %131, %132 ]
%136 = phi double [ %153, %133 ], [ 0.000000e+00, %132 ]
%137 = phi i32 [ %154, %133 ], [ 0, %132 ]
%138 = getelementptr inbounds double, double* %134, i64 1
%139 = load double, double* %134
%140 = fmul double %129, %135
%141 = tail call double @llvm.fmuladd.f64(double %139, double %140, double %136)
%142 = getelementptr inbounds double, double* %134, i64 2
%143 = load double, double* %138
%144 = fmul double %129, %140
%145 = tail call double @llvm.fmuladd.f64(double %143, double %144, double %141)
%146 = getelementptr inbounds double, double* %134, i64 3
%147 = load double, double* %142
%148 = fmul double %129, %144
%149 = tail call double @llvm.fmuladd.f64(double %147, double %148, double %145)
%150 = getelementptr inbounds double, double* %134, i64 4
%151 = load double, double* %146
%152 = fmul double %129, %148
%153 = tail call double @llvm.fmuladd.f64(double %151, double %152, double %149)
%154 = add i32 %137, 4
%155 = icmp eq i32 %154, %121
br i1 %155, label %156, label %133
156: ; preds = %133, %132
%157 = phi double [ undef, %132 ], [ %153, %133 ]
%158 = phi double* [ %2, %132 ], [ %150, %133 ]
%159 = phi double [ %131, %132 ], [ %152, %133 ]
%160 = phi double [ 0.000000e+00, %132 ], [ %153, %133 ]
br i1 %122, label %172, label %161
161: ; preds = %156, %161
%162 = phi double* [ %166, %161 ], [ %158, %156 ]
%163 = phi double [ %168, %161 ], [ %159, %156 ]
%164 = phi double [ %169, %161 ], [ %160, %156 ]
%165 = phi i32 [ %170, %161 ], [ 0, %156 ]
%166 = getelementptr inbounds double, double* %162, i64 1
%167 = load double, double* %162
%168 = fmul double %129, %163
%169 = tail call double @llvm.fmuladd.f64(double %167, double %168, double %164)
%170 = add i32 %165, 1
%171 = icmp eq i32 %170, %119
br i1 %171, label %172, label %161
172: ; preds = %161, %156
%173 = phi double [ %157, %156 ], [ %169, %161 ]
%174 = fcmp ogt double %173, 0x3FF000010C6F7A0B
br i1 %174, label %179, label %175
175: ; preds = %123, %172
%176 = fsub double %125, %128
%177 = tail call i32 @fabs(double %176)
%178 = icmp slt i32 %177, 1
br i1 %178, label %186, label %179
179: ; preds = %175, %172
%180 = phi double [ %128, %172 ], [ %126, %175 ]
%181 = phi double [ %125, %172 ], [ %128, %175 ]
%182 = add nuw nsw i32 %124, 1
%183 = icmp eq i32 %182, 26
br i1 %183, label %184, label %123
184: ; preds = %179
%185 = fcmp ugt double %181, 1.000000e-10
br i1 %185, label %186, label %231
186: ; preds = %175, %184
%187 = phi double [ %181, %184 ], [ %125, %175 ]
store double %187, double* %3
%188 = fmul double %187, %61
%189 = tail call double @exp(double %188)
br i1 %6, label %228, label %190
190: ; preds = %186
%191 = and i32 %64, 1
%192 = icmp eq i32 %191, 0
br i1 %192, label %201, label %193
193: ; preds = %190
%194 = getelementptr inbounds double, double* %2, i64 1
%195 = load double, double* %2
%196 = sitofp i32 %0 to double
%197 = fmul double %195, %196
%198 = fmul double %129, %189
%199 = tail call double @llvm.fmuladd.f64(double %197, double %198, double 0.000000e+00)
%200 = add i32 %0, 1
br label %201
201: ; preds = %193, %190
%202 = phi double [ undef, %190 ], [ %199, %193 ]
%203 = phi double* [ %2, %190 ], [ %194, %193 ]
%204 = phi double [ %189, %190 ], [ %198, %193 ]
%205 = phi double [ 0.000000e+00, %190 ], [ %199, %193 ]
%206 = phi i32 [ %0, %190 ], [ %200, %193 ]
%207 = icmp eq i32 %10, 0
br i1 %207, label %228, label %208
208: ; preds = %201, %208
%209 = phi double* [ %220, %208 ], [ %203, %201 ]
%210 = phi double [ %224, %208 ], [ %204, %201 ]
%211 = phi double [ %225, %208 ], [ %205, %201 ]
%212 = phi i32 [ %226, %208 ], [ %206, %201 ]
%213 = getelementptr inbounds double, double* %209, i64 1
%214 = load double, double* %209
%215 = sitofp i32 %212 to double
%216 = fmul double %214, %215
%217 = fmul double %129, %210
%218 = tail call double @llvm.fmuladd.f64(double %216, double %217, double %211)
%219 = add i32 %212, 1
%220 = getelementptr inbounds double, double* %209, i64 2
%221 = load double, double* %213
%222 = sitofp i32 %219 to double
%223 = fmul double %221, %222
%224 = fmul double %129, %217
%225 = tail call double @llvm.fmuladd.f64(double %223, double %224, double %218)
%226 = add i32 %212, 2
%227 = icmp eq i32 %219, %1
br i1 %227, label %228, label %208
228: ; preds = %201, %208, %186
%229 = phi double [ 0.000000e+00, %186 ], [ %202, %201 ], [ %225, %208 ]
%230 = fmul double %187, %229
store double %230, double* %4
br label %231
231: ; preds = %5, %184, %56, %228
%232 = phi i32 [ 1, %228 ], [ 0, %56 ], [ 0, %184 ], [ 0, %5 ]
ret i32 %232
}