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

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ DiscreteGamma(double *, double *, double, double, int,
                                 int /*<<< orphan*/);
double InverseCDFNormal(int);
double LBinormal(double, double, double);
int /*<<< orphan*/ printf(char *, int, int, double);
int AutodGamma(double M[], double freqK[], double rK[], double *rho1,
               double alpha, double rho, int K) {
  int i, j, i1, i2;
  double *point = freqK;
  double x, y, large = 20, v1;
  for (i = 0; i < K - 1; i++)
    point[i] = InverseCDFNormal((i + 1.0) / K);
  for (i = 0; i < K; i++) {
    for (j = 0; j < K; j++) {
      x = (i < K - 1 ? point[i] : large);
      y = (j < K - 1 ? point[j] : large);
      M[i * K + j] = LBinormal(-(x), -(y), rho);
    }
  }
  for (i1 = 0; i1 < 2 * K - 1; i1++) {
    for (i2 = 0; i2 < K * K; i2++) {
      i = i2 / K;
      j = i2 % K;
      if (i + j != 2 * (K - 1) - i1)
        continue;
      y = 0;
      if (i > 0)
        y -= M[(i - 1) * K + j];
      if (j > 0)
        y -= M[i * K + (j - 1)];
      if (i > 0 && j > 0)
        y += M[(i - 1) * K + (j - 1)];
      M[i * K + j] = (M[i * K + j] + y) * K;
      if (M[i * K + j] < 0)
        printf("M(%d,%d) =%12.8f<0\n", i + 1, j + 1, M[i * K + j]);
    }
  }
  DiscreteGamma(freqK, rK, alpha, alpha, K, 0);
  for (i = 0, v1 = *rho1 = 0; i < K; i++) {
    v1 += rK[i] * rK[i] * freqK[i];
    for (j = 0; j < K; j++)
      *rho1 += freqK[i] * M[i * K + j] * rK[i] * rK[j];
  }
  v1 -= 1;
  *rho1 = (*rho1 - 1) / v1;
  return (0);
}


x86 assembly (clang, -O0):
<AutodGamma>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movsd  %xmm0,-0x28(%rbp)
movsd  %xmm1,-0x30(%rbp)
mov    %r8d,-0x34(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x50(%rbp)
movsd  0xe87(%rip),%xmm0
movsd  %xmm0,-0x68(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    11d4 <AutodGamma+0x94>
cvtsi2sdl -0x38(%rbp),%xmm0
movsd  0xe65(%rip),%xmm1
addsd  %xmm1,%xmm0
cvtsi2sdl -0x34(%rbp),%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1030 <InverseCDFNormal@plt>
mov    -0x50(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1185 <AutodGamma+0x45>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    12eb <AutodGamma+0x1ab>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    12d8 <AutodGamma+0x198>
mov    -0x38(%rbp),%eax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1222 <AutodGamma+0xe2>
mov    -0x50(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
jmp    122c <AutodGamma+0xec>
movsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x3c(%rbp),%eax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    125e <AutodGamma+0x11e>
mov    -0x50(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
jmp    1268 <AutodGamma+0x128>
movsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x58(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  -0x60(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
movsd  -0x30(%rbp),%xmm2
call   1060 <LBinormal@plt>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   -0x34(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    11ee <AutodGamma+0xae>
jmp    12dd <AutodGamma+0x19d>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    11db <AutodGamma+0x9b>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
mov    -0x34(%rbp),%ecx
shl    $0x1,%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    14a7 <AutodGamma+0x367>
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
mov    -0x34(%rbp),%ecx
imul   -0x34(%rbp),%ecx
cmp    %ecx,%eax
jge    1494 <AutodGamma+0x354>
mov    -0x44(%rbp),%eax
cltd
idivl  -0x34(%rbp)
mov    %eax,-0x38(%rbp)
mov    -0x44(%rbp),%eax
cltd
idivl  -0x34(%rbp)
mov    %edx,-0x3c(%rbp)
mov    -0x38(%rbp),%eax
add    -0x3c(%rbp),%eax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
shl    $0x1,%ecx
sub    -0x40(%rbp),%ecx
cmp    %ecx,%eax
je     1352 <AutodGamma+0x212>
jmp    1486 <AutodGamma+0x346>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x60(%rbp)
cmpl   $0x0,-0x38(%rbp)
jle    138b <AutodGamma+0x24b>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x34(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x60(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
cmpl   $0x0,-0x3c(%rbp)
jle    13be <AutodGamma+0x27e>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   -0x34(%rbp),%ecx
mov    -0x3c(%rbp),%edx
sub    $0x1,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x60(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
cmpl   $0x0,-0x38(%rbp)
jle    13fa <AutodGamma+0x2ba>
cmpl   $0x0,-0x3c(%rbp)
jle    13fa <AutodGamma+0x2ba>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x34(%rbp),%ecx
mov    -0x3c(%rbp),%edx
sub    $0x1,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   -0x34(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x60(%rbp),%xmm0
cvtsi2sdl -0x34(%rbp),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   -0x34(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   -0x34(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
ucomisd (%rax,%rcx,8),%xmm0
jbe    1481 <AutodGamma+0x341>
mov    -0x38(%rbp),%esi
add    $0x1,%esi
mov    -0x3c(%rbp),%edx
add    $0x1,%edx
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   -0x34(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
lea    0xb94(%rip),%rdi
call   1040 <printf@plt>
jmp    1486 <AutodGamma+0x346>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    130d <AutodGamma+0x1cd>
jmp    1499 <AutodGamma+0x359>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    12f2 <AutodGamma+0x1b2>
mov    -0x10(%rbp),%rdi
mov    -0x18(%rbp),%rsi
movsd  -0x28(%rbp),%xmm0
movsd  -0x28(%rbp),%xmm1
mov    -0x34(%rbp),%edx
xor    %ecx,%ecx
call   1050 <DiscreteGamma@plt>
movl   $0x0,-0x38(%rbp)
mov    -0x20(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    15a7 <AutodGamma+0x467>
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x70(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1594 <AutodGamma+0x454>
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   -0x34(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1529 <AutodGamma+0x3e9>
jmp    1599 <AutodGamma+0x459>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    14dd <AutodGamma+0x39d>
movsd  -0x70(%rbp),%xmm0
movsd  0xa54(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  0xa3b(%rip),%xmm1
subsd  %xmm1,%xmm0
divsd  -0x70(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
xor    %eax,%eax
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<AutodGamma>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
movsd  %xmm1,0x20(%rsp)
movsd  %xmm0,0x28(%rsp)
mov    %rcx,0x38(%rsp)
mov    %rdx,0x30(%rsp)
mov    %rsi,%r15
mov    %rdi,0x10(%rsp)
lea    -0x1(%r8),%r12d
mov    %r12d,%ebp
mov    %r8,%r14
cmp    $0x2,%r8d
jl     11bd <AutodGamma+0x7d>
xorps  %xmm0,%xmm0
cvtsi2sd %r14d,%xmm0
movsd  %xmm0,0x8(%rsp)
xor    %ebx,%ebx
nopl   0x0(%rax)
xorps  %xmm0,%xmm0
cvtsi2sd %ebx,%xmm0
addsd  0xe61(%rip),%xmm0
divsd  0x8(%rsp),%xmm0
cvttsd2si %xmm0,%edi
call   1030 <InverseCDFNormal@plt>
movsd  %xmm0,(%r15,%rbx,8)
add    $0x1,%rbx
cmp    %rbx,%rbp
jne    1190 <AutodGamma+0x50>
mov    %r15,0x8(%rsp)
mov    %r14,%rbx
test   %ebx,%ebx
jle    13ab <AutodGamma+0x26b>
mov    %r12d,0x4(%rsp)
mov    %ebx,%eax
lea    0x0(,%rax,8),%rcx
mov    %rcx,0x18(%rsp)
xor    %r12d,%r12d
movsd  0xe1c(%rip),%xmm1
mov    0x10(%rsp),%r13
mov    %rax,0x40(%rsp)
jmp    120e <AutodGamma+0xce>
nopl   0x0(%rax,%rax,1)
add    $0x1,%r12
add    0x18(%rsp),%r13
cmp    %rax,%r12
je     127e <AutodGamma+0x13e>
xor    %r15d,%r15d
jmp    1258 <AutodGamma+0x118>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movapd 0xdf8(%rip),%xmm2
xorpd  %xmm2,%xmm0
xorpd  %xmm2,%xmm1
movsd  0x20(%rsp),%xmm2
call   1060 <LBinormal@plt>
movsd  %xmm0,0x0(%r13,%r15,8)
add    $0x1,%r15
mov    0x40(%rsp),%rax
cmp    %r15,%rax
movsd  0xdb2(%rip),%xmm1
je     1200 <AutodGamma+0xc0>
movapd %xmm1,%xmm0
cmp    %rbp,%r12
jae    126c <AutodGamma+0x12c>
mov    0x8(%rsp),%rcx
movsd  (%rcx,%r12,8),%xmm0
cmp    %rbp,%r15
jae    1220 <AutodGamma+0xe0>
mov    0x8(%rsp),%rax
movsd  (%rax,%r15,8),%xmm1
jmp    1220 <AutodGamma+0xe0>
mov    %r14,%rbx
test   %ebx,%ebx
mov    0x4(%rsp),%r13d
jle    13ab <AutodGamma+0x26b>
mov    %ebx,%r8d
imul   %ebx,%r8d
add    %r13d,%r13d
xorps  %xmm1,%xmm1
cvtsi2sd %ebx,%xmm1
cmp    $0x2,%r8d
mov    $0x1,%ecx
cmovb  %ecx,%r8d
lea    (%rbx,%rbx,1),%eax
add    $0xffffffff,%eax
cmp    $0x2,%eax
cmovge %eax,%ecx
mov    %ecx,0x18(%rsp)
xor    %eax,%eax
xorpd  %xmm2,%xmm2
mov    0x10(%rsp),%r15
movsd  %xmm1,0x20(%rsp)
mov    %r13d,0x4(%rsp)
jmp    12f5 <AutodGamma+0x1b5>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    %r12d,%eax
add    $0x1,%eax
cmp    0x18(%rsp),%eax
mov    0x4(%rsp),%r13d
je     13dd <AutodGamma+0x29d>
mov    %eax,%r12d
sub    %eax,%r13d
xor    %ebp,%ebp
jmp    1308 <AutodGamma+0x1c8>
nop
add    $0x1,%ebp
cmp    %ebp,%r8d
je     12e0 <AutodGamma+0x1a0>
mov    %ebp,%eax
cltd
idiv   %ebx
lea    (%rdx,%rax,1),%ecx
cmp    %r13d,%ecx
jne    1300 <AutodGamma+0x1c0>
test   %eax,%eax
jle    1330 <AutodGamma+0x1f0>
lea    -0x1(%rax),%ecx
imul   %ebx,%ecx
add    %edx,%ecx
movslq %ecx,%rcx
xorpd  %xmm0,%xmm0
subsd  (%r15,%rcx,8),%xmm0
jmp    1334 <AutodGamma+0x1f4>
xorpd  %xmm0,%xmm0
mov    %eax,%ecx
imul   %ebx,%ecx
test   %edx,%edx
je     1364 <AutodGamma+0x224>
lea    (%rdx,%rcx,1),%esi
add    $0xffffffff,%esi
movslq %esi,%rsi
subsd  (%r15,%rsi,8),%xmm0
test   %eax,%eax
jle    1364 <AutodGamma+0x224>
lea    -0x1(%rdx),%esi
lea    -0x1(%rax),%edi
imul   %ebx,%edi
add    %esi,%edi
movslq %edi,%rsi
addsd  (%r15,%rsi,8),%xmm0
add    %edx,%ecx
movslq %ecx,%rcx
addsd  (%r15,%rcx,8),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,(%r15,%rcx,8)
ucomisd %xmm0,%xmm2
jbe    1300 <AutodGamma+0x1c0>
add    $0x1,%eax
add    $0x1,%edx
lea    0xca4(%rip),%rdi
mov    %eax,%esi
mov    %r8d,%ebx
call   1040 <printf@plt>
xorpd  %xmm2,%xmm2
movsd  0x20(%rsp),%xmm1
mov    %ebx,%r8d
mov    %r14,%rbx
jmp    1300 <AutodGamma+0x1c0>
mov    0x8(%rsp),%rdi
mov    0x30(%rsp),%rsi
movsd  0x28(%rsp),%xmm0
movaps %xmm0,%xmm1
mov    %ebx,%edx
xor    %ecx,%ecx
call   1050 <DiscreteGamma@plt>
movsd  0xc41(%rip),%xmm0
movapd %xmm0,%xmm1
mov    0x38(%rsp),%rdx
jmp    1569 <AutodGamma+0x429>
mov    0x8(%rsp),%r13
mov    %r13,%rdi
mov    0x30(%rsp),%rbp
mov    %rbp,%rsi
movsd  0x28(%rsp),%xmm0
movapd %xmm0,%xmm1
mov    %ebx,%edx
xor    %ecx,%ecx
call   1050 <DiscreteGamma@plt>
mov    0x38(%rsp),%rdx
movq   $0x0,(%rdx)
test   %r14d,%r14d
jle    155d <AutodGamma+0x41d>
mov    0x40(%rsp),%r12
lea    -0x1(%r12),%r8
lea    0x8(%r15),%rcx
lea    0x0(,%r12,8),%r9
mov    %r8,%rsi
and    $0xfffffffffffffffe,%rsi
neg    %rsi
xorpd  %xmm0,%xmm0
xor    %edi,%edi
xorpd  %xmm1,%xmm1
jmp    1471 <AutodGamma+0x331>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mulsd  %xmm2,%xmm2
mulsd  %xmm3,%xmm2
addsd  %xmm2,%xmm1
add    $0x1,%rdi
add    %r9,%rcx
cmp    %r12,%rdi
mov    0x10(%rsp),%r15
je     154b <AutodGamma+0x40b>
movapd %xmm0,%xmm4
movsd  0x0(%rbp,%rdi,8),%xmm2
movsd  0x0(%r13,%rdi,8),%xmm3
mov    %rdi,%r10
imul   %r12,%r10
movsd  (%r15,%r10,8),%xmm0
mulsd  %xmm3,%xmm0
mulsd  %xmm2,%xmm0
mulsd  0x0(%rbp),%xmm0
addsd  %xmm4,%xmm0
movsd  %xmm0,(%rdx)
cmp    $0x1,%r14d
je     1450 <AutodGamma+0x310>
mov    $0x1,%ebx
cmp    $0x2,%r14d
je     1513 <AutodGamma+0x3d3>
mov    $0x1,%ebx
nopw   0x0(%rax,%rax,1)
movsd  0x0(%r13,%rdi,8),%xmm4
mulsd  -0x8(%rcx,%rbx,8),%xmm4
mulsd  0x0(%rbp,%rdi,8),%xmm4
mulsd  0x0(%rbp,%rbx,8),%xmm4
addsd  %xmm0,%xmm4
movsd  %xmm4,(%rdx)
movsd  0x0(%r13,%rdi,8),%xmm0
mulsd  (%rcx,%rbx,8),%xmm0
mulsd  0x0(%rbp,%rdi,8),%xmm0
mulsd  0x8(%rbp,%rbx,8),%xmm0
addsd  %xmm4,%xmm0
movsd  %xmm0,(%rdx)
lea    (%rsi,%rbx,1),%rax
add    $0x2,%rax
add    $0x2,%rbx
cmp    $0x1,%rax
jne    14c0 <AutodGamma+0x380>
test   $0x1,%r8b
je     1450 <AutodGamma+0x310>
movsd  0x0(%r13,%rdi,8),%xmm4
add    %rbx,%r10
mov    0x10(%rsp),%rax
mulsd  (%rax,%r10,8),%xmm4
mulsd  0x0(%rbp,%rdi,8),%xmm4
mulsd  0x0(%rbp,%rbx,8),%xmm4
addsd  %xmm4,%xmm0
movsd  %xmm0,(%rdx)
jmp    1450 <AutodGamma+0x310>
movsd  0xabd(%rip),%xmm2
addsd  %xmm2,%xmm1
addsd  %xmm2,%xmm0
jmp    1569 <AutodGamma+0x429>
movsd  0xaab(%rip),%xmm0
movapd %xmm0,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rdx)
xor    %eax,%eax
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
@.str = [20 x i8] c"M(%d,%d) =%12.8f<0\0A\00"
define i32 @AutodGamma(double* nocapture %0, double* %1, double* %2, double* nocapture %3, double %4, double %5, i32 %6) {
%8 = add i32 %6, -1
%9 = icmp sgt i32 %6, 1
br i1 %9, label %10, label %13
10: ; preds = %7
%11 = sitofp i32 %6 to double
%12 = zext i32 %8 to i64
br label %19
13: ; preds = %19, %7
%14 = icmp sgt i32 %6, 0
br i1 %14, label %15, label %123
15: ; preds = %13
%16 = zext i32 %8 to i64
%17 = zext i32 %6 to i64
%18 = zext i32 %6 to i64
br label %30
19: ; preds = %10, %19
%20 = phi i64 [ 0, %10 ], [ %28, %19 ]
%21 = trunc i64 %20 to i32
%22 = sitofp i32 %21 to double
%23 = fadd double %22, 1.000000e+00
%24 = fdiv double %23, %11
%25 = fptosi double %24 to i32
%26 = tail call double @InverseCDFNormal(i32 %25)
%27 = getelementptr inbounds double, double* %1, i64 %20
store double %26, double* %27
%28 = add nuw nsw i64 %20, 1
%29 = icmp eq i64 %28, %12
br i1 %29, label %13, label %19
30: ; preds = %63, %15
%31 = phi i64 [ 0, %15 ], [ %64, %63 ]
%32 = icmp ult i64 %31, %16
%33 = getelementptr inbounds double, double* %1, i64 %31
%34 = mul nsw i64 %31, %17
br label %44
35: ; preds = %63
br i1 %14, label %36, label %123
36: ; preds = %35
%37 = shl nuw nsw i32 %6, 1
%38 = add i32 %37, -1
%39 = mul i32 %6, %6
%40 = shl nuw nsw i32 %8, 1
%41 = sitofp i32 %6 to double
%42 = call i32 @llvm.umax.i32(i32 %39, i32 1)
%43 = call i32 @llvm.smax.i32(i32 %38, i32 1)
br label %66
44: ; preds = %30, %54
%45 = phi i64 [ 0, %30 ], [ %61, %54 ]
br i1 %32, label %46, label %48
46: ; preds = %44
%47 = load double, double* %33
br label %48
48: ; preds = %44, %46
%49 = phi double [ %47, %46 ], [ 2.000000e+01, %44 ]
%50 = icmp ult i64 %45, %16
br i1 %50, label %51, label %54
51: ; preds = %48
%52 = getelementptr inbounds double, double* %1, i64 %45
%53 = load double, double* %52
br label %54
54: ; preds = %48, %51
%55 = phi double [ %53, %51 ], [ 2.000000e+01, %48 ]
%56 = fneg double %49
%57 = fneg double %55
%58 = tail call double @LBinormal(double %56, double %57, double %5)
%59 = add nuw nsw i64 %45, %34
%60 = getelementptr inbounds double, double* %0, i64 %59
store double %58, double* %60
%61 = add nuw nsw i64 %45, 1
%62 = icmp eq i64 %61, %18
br i1 %62, label %63, label %44
63: ; preds = %54
%64 = add nuw nsw i64 %31, 1
%65 = icmp eq i64 %64, %18
br i1 %65, label %35, label %30
66: ; preds = %120, %36
%67 = phi i32 [ 0, %36 ], [ %121, %120 ]
%68 = sub nsw i32 %40, %67
br label %69
69: ; preds = %66, %117
%70 = phi i32 [ 0, %66 ], [ %118, %117 ]
%71 = sdiv i32 %70, %6
%72 = srem i32 %70, %6
%73 = add nsw i32 %72, %71
%74 = icmp eq i32 %73, %68
br i1 %74, label %75, label %117
75: ; preds = %69
%76 = icmp sgt i32 %71, 0
br i1 %76, label %77, label %85
77: ; preds = %75
%78 = add nsw i32 %71, -1
%79 = mul nsw i32 %78, %6
%80 = add nsw i32 %79, %72
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds double, double* %0, i64 %81
%83 = load double, double* %82
%84 = fsub double 0.000000e+00, %83
br label %85
85: ; preds = %77, %75
%86 = phi double [ %84, %77 ], [ 0.000000e+00, %75 ]
%87 = icmp eq i32 %72, 0
%88 = mul nsw i32 %71, %6
br i1 %87, label %104, label %89
89: ; preds = %85
%90 = add nsw i32 %72, -1
%91 = add nsw i32 %90, %88
%92 = sext i32 %91 to i64
%93 = getelementptr inbounds double, double* %0, i64 %92
%94 = load double, double* %93
%95 = fsub double %86, %94
br i1 %76, label %96, label %104
96: ; preds = %89
%97 = add nsw i32 %71, -1
%98 = mul nsw i32 %97, %6
%99 = add nsw i32 %90, %98
%100 = sext i32 %99 to i64
%101 = getelementptr inbounds double, double* %0, i64 %100
%102 = load double, double* %101
%103 = fadd double %95, %102
br label %104
104: ; preds = %85, %96, %89
%105 = phi double [ %103, %96 ], [ %95, %89 ], [ %86, %85 ]
%106 = add nsw i32 %88, %72
%107 = sext i32 %106 to i64
%108 = getelementptr inbounds double, double* %0, i64 %107
%109 = load double, double* %108
%110 = fadd double %105, %109
%111 = fmul double %110, %41
store double %111, double* %108
%112 = fcmp olt double %111, 0.000000e+00
br i1 %112, label %113, label %117
113: ; preds = %104
%114 = add nsw i32 %71, 1
%115 = add nuw nsw i32 %72, 1
%116 = tail call i32 @printf(i8* getelementptr inbounds ([20 x i8], [20 x i8]* @.str, i64 0, i64 0), i32 %114, i32 %115, double %111)
br label %117
117: ; preds = %104, %113, %69
%118 = add nuw i32 %70, 1
%119 = icmp eq i32 %118, %42
br i1 %119, label %120, label %69
120: ; preds = %117
%121 = add nuw nsw i32 %67, 1
%122 = icmp eq i32 %121, %43
br i1 %122, label %125, label %66
123: ; preds = %35, %13
%124 = tail call i32 @DiscreteGamma(double* %1, double* %2, double %4, double %4, i32 %6, i32 0)
br label %202
125: ; preds = %120
%126 = tail call i32 @DiscreteGamma(double* %1, double* %2, double %4, double %4, i32 %6, i32 0)
store double 0.000000e+00, double* %3
br i1 %14, label %127, label %202
127: ; preds = %125
%128 = zext i32 %6 to i64
%129 = zext i32 %6 to i64
%130 = add nsw i64 %129, -1
%131 = and i64 %130, 1
%132 = icmp eq i32 %6, 2
%133 = and i64 %130, -2
%134 = icmp eq i64 %131, 0
br label %135
135: ; preds = %195, %127
%136 = phi double [ 0.000000e+00, %127 ], [ %196, %195 ]
%137 = phi i64 [ 0, %127 ], [ %197, %195 ]
%138 = phi double [ 0.000000e+00, %127 ], [ %144, %195 ]
%139 = getelementptr inbounds double, double* %2, i64 %137
%140 = load double, double* %139
%141 = fmul double %140, %140
%142 = getelementptr inbounds double, double* %1, i64 %137
%143 = load double, double* %142
%144 = tail call double @llvm.fmuladd.f64(double %141, double %143, double %138)
%145 = mul nsw i64 %137, %128
%146 = getelementptr inbounds double, double* %0, i64 %145
%147 = load double, double* %146
%148 = fmul double %143, %147
%149 = fmul double %148, %140
%150 = load double, double* %2
%151 = tail call double @llvm.fmuladd.f64(double %149, double %150, double %136)
store double %151, double* %3
switch i32 %6, label %152 [
i32 1, label %195
i32 2, label %180
]
152: ; preds = %135, %152
%153 = phi i64 [ %177, %152 ], [ 1, %135 ]
%154 = phi double [ %176, %152 ], [ %151, %135 ]
%155 = phi i64 [ %178, %152 ], [ 0, %135 ]
%156 = load double, double* %142
%157 = load double, double* %139
%158 = add nuw nsw i64 %153, %145
%159 = getelementptr inbounds double, double* %0, i64 %158
%160 = load double, double* %159
%161 = fmul double %156, %160
%162 = fmul double %161, %157
%163 = getelementptr inbounds double, double* %2, i64 %153
%164 = load double, double* %163
%165 = tail call double @llvm.fmuladd.f64(double %162, double %164, double %154)
store double %165, double* %3
%166 = add nuw nsw i64 %153, 1
%167 = load double, double* %142
%168 = load double, double* %139
%169 = add nuw nsw i64 %166, %145
%170 = getelementptr inbounds double, double* %0, i64 %169
%171 = load double, double* %170
%172 = fmul double %167, %171
%173 = fmul double %172, %168
%174 = getelementptr inbounds double, double* %2, i64 %166
%175 = load double, double* %174
%176 = tail call double @llvm.fmuladd.f64(double %173, double %175, double %165)
store double %176, double* %3
%177 = add nuw nsw i64 %153, 2
%178 = add i64 %155, 2
%179 = icmp eq i64 %178, %133
br i1 %179, label %180, label %152
180: ; preds = %152, %135
%181 = phi double [ undef, %135 ], [ %176, %152 ]
%182 = phi i64 [ 1, %135 ], [ %177, %152 ]
%183 = phi double [ %151, %135 ], [ %176, %152 ]
br i1 %134, label %195, label %184
184: ; preds = %180
%185 = load double, double* %142
%186 = load double, double* %139
%187 = add nuw nsw i64 %182, %145
%188 = getelementptr inbounds double, double* %0, i64 %187
%189 = load double, double* %188
%190 = fmul double %185, %189
%191 = fmul double %190, %186
%192 = getelementptr inbounds double, double* %2, i64 %182
%193 = load double, double* %192
%194 = tail call double @llvm.fmuladd.f64(double %191, double %193, double %183)
store double %194, double* %3
br label %195
195: ; preds = %184, %180, %135
%196 = phi double [ %151, %135 ], [ %181, %180 ], [ %194, %184 ]
%197 = add nuw nsw i64 %137, 1
%198 = icmp eq i64 %197, %129
br i1 %198, label %199, label %135
199: ; preds = %195
%200 = fadd double %144, -1.000000e+00
%201 = fadd double %196, -1.000000e+00
br label %202
202: ; preds = %123, %199, %125
%203 = phi double [ -1.000000e+00, %125 ], [ %201, %199 ], [ -1.000000e+00, %123 ]
%204 = phi double [ -1.000000e+00, %125 ], [ %200, %199 ], [ -1.000000e+00, %123 ]
%205 = fdiv double %203, %204
store double %205, double* %3
ret i32 0
}