#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*/ LUDecompose_det2(double *, int);
double *dvector(int);
double exp(double const);
int /*<<< orphan*/ free(double *);
double log(int /*<<< orphan*/);
void inverse_transform_msimplex2(double *x, const double *y, int K,
                                 double *lp) {
  double sum = 0;
  double sum_one = 1.0;
  int dim = K - 1;
  int k = 0;
  for (; k < dim; k++) {
    x[k] = exp(y[k]);
    sum_one += x[k];
  }
  x[k] = 1.0 / (sum_one);
  for (int i = 0; i < dim; i++) {
    x[i] /= x[k];
  }
  double inv_sum_one_squared = 1.0 / (sum_one * sum_one);
  double *m = dvector(dim * dim);
  k = 0;
  for (int i = 0; i < dim; i++) {
    double temp = 1;
    int j = 0;
    for (; j < i; j++)
      temp += x[j];
    j++;
    for (; j < dim; j++)
      temp += x[j];
    m[dim * i + i] = x[i] * temp / inv_sum_one_squared;
    for (int j = i + 1; j < dim; j++) {
      m[dim * i + j] = -exp(y[i] + y[j]) / inv_sum_one_squared;
      m[dim * j + i] = m[dim * i + j];
    }
  }
  *lp = log(LUDecompose_det2(m, dim));
  free(m);
}


x86 assembly (clang, -O0):
<inverse_transform_msimplex2>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  0xe89(%rip),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x14(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x34(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    11dc <inverse_transform_msimplex2+0x8c>
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1040 <exp@plt>
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    118c <inverse_transform_msimplex2+0x3c>
movsd  0xe1c(%rip),%xmm0
divsd  -0x30(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    123a <inverse_transform_msimplex2+0xea>
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    11fd <inverse_transform_msimplex2+0xad>
movsd  -0x30(%rbp),%xmm1
mulsd  -0x30(%rbp),%xmm1
movsd  0xdb4(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x34(%rbp),%edi
imul   -0x34(%rbp),%edi
call   1050 <dvector@plt>
mov    %rax,-0x50(%rbp)
movl   $0x0,-0x38(%rbp)
movl   $0x0,-0x54(%rbp)
mov    -0x54(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    13de <inverse_transform_msimplex2+0x28e>
movsd  0xd79(%rip),%xmm0
movsd  %xmm0,-0x60(%rbp)
movl   $0x0,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x54(%rbp),%eax
jge    12c4 <inverse_transform_msimplex2+0x174>
mov    -0x8(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    1293 <inverse_transform_msimplex2+0x143>
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    12fe <inverse_transform_msimplex2+0x1ae>
mov    -0x8(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    12cd <inverse_transform_msimplex2+0x17d>
mov    -0x8(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x60(%rbp),%xmm0
divsd  -0x48(%rbp),%xmm0
mov    -0x50(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x54(%rbp),%ecx
add    -0x54(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x68(%rbp)
mov    -0x68(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    13cb <inverse_transform_msimplex2+0x27b>
mov    -0x10(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x68(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
call   1040 <exp@plt>
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
divsd  -0x48(%rbp),%xmm0
mov    -0x50(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x54(%rbp),%ecx
add    -0x68(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x50(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x54(%rbp),%ecx
add    -0x68(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x50(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x68(%rbp),%ecx
add    -0x54(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x68(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x68(%rbp)
jmp    1334 <inverse_transform_msimplex2+0x1e4>
jmp    13d0 <inverse_transform_msimplex2+0x280>
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    1273 <inverse_transform_msimplex2+0x123>
mov    -0x50(%rbp),%rdi
mov    -0x34(%rbp),%esi
call   1070 <LUDecompose_det2@plt>
mov    %eax,%edi
call   1060 <log@plt>
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x50(%rbp),%rdi
call   1030 <free@plt>
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<inverse_transform_msimplex2>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x78,%rsp
mov    %rcx,%r15
mov    %rsi,%r13
mov    %rdi,%rbp
lea    -0x1(%rdx),%r12d
mov    %rdx,0x10(%rsp)
cmp    $0x1,%edx
mov    %rdi,0x58(%rsp)
mov    %rcx,0x40(%rsp)
jle    1271 <inverse_transform_msimplex2+0x121>
mov    %r12d,%r14d
movsd  0xe72(%rip),%xmm2
xor    %ebx,%ebx
movsd  %xmm2,0x8(%rsp)
movsd  0x0(%r13,%rbx,8),%xmm0
call   1040 <exp@plt>
movsd  0x8(%rsp),%xmm2
movsd  %xmm0,0x0(%rbp,%rbx,8)
addsd  %xmm0,%xmm2
add    $0x1,%rbx
cmp    %rbx,%r14
jne    1190 <inverse_transform_msimplex2+0x40>
movsd  0xe3d(%rip),%xmm0
divsd  %xmm2,%xmm0
movsd  %xmm0,0x0(%rbp,%r14,8)
cmpl   $0x2,0x10(%rsp)
movsd  %xmm2,0x8(%rsp)
jl     128d <inverse_transform_msimplex2+0x13d>
movsd  0x0(%rbp),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x0(%rbp)
cmp    $0x1,%r12d
je     128d <inverse_transform_msimplex2+0x13d>
lea    -0x1(%r14),%rax
mov    $0x1,%ecx
cmp    $0x2,%r12d
je     1258 <inverse_transform_msimplex2+0x108>
mov    %rax,%rdx
and    $0xfffffffffffffffe,%rdx
neg    %rdx
mov    $0x1,%ecx
cs nopw 0x0(%rax,%rax,1)
nop
movsd  0x0(%rbp,%rcx,8),%xmm0
divsd  0x0(%rbp,%r14,8),%xmm0
movsd  0x8(%rbp,%rcx,8),%xmm1
movsd  %xmm0,0x0(%rbp,%rcx,8)
divsd  0x0(%rbp,%r14,8),%xmm1
movsd  %xmm1,0x8(%rbp,%rcx,8)
lea    (%rdx,%rcx,1),%rsi
add    $0x2,%rsi
add    $0x2,%rcx
cmp    $0x1,%rsi
jne    1220 <inverse_transform_msimplex2+0xd0>
test   $0x1,%al
je     128d <inverse_transform_msimplex2+0x13d>
movsd  0x0(%rbp,%rcx,8),%xmm0
divsd  0x0(%rbp,%r14,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rcx,8)
jmp    128d <inverse_transform_msimplex2+0x13d>
movabs $0x3ff0000000000000,%rax
mov    %rax,0x0(%rbp)
movsd  0xd79(%rip),%xmm0
movsd  %xmm0,0x8(%rsp)
mov    %r12d,%edi
imul   %r12d,%edi
call   1050 <dvector@plt>
mov    %rax,0x18(%rsp)
cmpl   $0x2,0x10(%rsp)
jl     1559 <inverse_transform_msimplex2+0x409>
movsd  0x8(%rsp),%xmm2
mulsd  %xmm2,%xmm2
movsd  0xd45(%rip),%xmm1
divsd  %xmm2,%xmm1
mov    %r12d,%r8d
mov    0x10(%rsp),%rax
lea    0x6(%rax),%ecx
mov    %ecx,0x28(%rsp)
add    $0xfffffffd,%eax
mov    %eax,0x24(%rsp)
mov    0x18(%rsp),%r14
lea    (%r14,%r8,8),%rax
mov    %rax,0x30(%rsp)
lea    0x0(,%r8,8),%r9
lea    0x8(,%r8,8),%rax
mov    %rax,0x48(%rsp)
mov    %r14,%rbx
add    $0x8,%rbx
lea    0x8(%r13),%r15
lea    -0x1(%r8),%rdi
mov    $0xffffffff,%r10d
xor    %edx,%edx
mov    %r13,0x38(%rsp)
mov    %r12d,0x2c(%rsp)
movsd  %xmm1,0x70(%rsp)
mov    %r8,0x50(%rsp)
jmp    1366 <inverse_transform_msimplex2+0x216>
cs nopw 0x0(%rax,%rax,1)
mov    0x68(%rsp),%r10
add    $0x1,%r10d
mov    0x48(%rsp),%rax
add    %rax,0x30(%rsp)
add    %rax,%rbx
add    $0x8,%r15
add    $0xffffffffffffffff,%rdi
mov    0x50(%rsp),%r8
mov    0x60(%rsp),%rdx
cmp    %r8,%rdx
mov    0x2c(%rsp),%r12d
je     1559 <inverse_transform_msimplex2+0x409>
movsd  0xc92(%rip),%xmm0
mov    $0x0,%r14d
test   %rdx,%rdx
mov    0x58(%rsp),%rbp
je     142e <inverse_transform_msimplex2+0x2de>
lea    -0x1(%rdx),%rax
cmp    $0x7,%rax
jae    13a0 <inverse_transform_msimplex2+0x250>
xor    %eax,%eax
movsd  0xc6a(%rip),%xmm0
jmp    13f9 <inverse_transform_msimplex2+0x2a9>
nopl   0x0(%rax,%rax,1)
mov    %rdx,%rcx
and    $0xfffffffffffffff8,%rcx
xor    %eax,%eax
movsd  0xc4f(%rip),%xmm0
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
addsd  0x0(%rbp,%rax,8),%xmm0
addsd  0x8(%rbp,%rax,8),%xmm0
addsd  0x10(%rbp,%rax,8),%xmm0
addsd  0x18(%rbp,%rax,8),%xmm0
addsd  0x20(%rbp,%rax,8),%xmm0
addsd  0x28(%rbp,%rax,8),%xmm0
addsd  0x30(%rbp,%rax,8),%xmm0
addsd  0x38(%rbp,%rax,8),%xmm0
add    $0x8,%rax
cmp    %rax,%rcx
jne    13c0 <inverse_transform_msimplex2+0x270>
lea    0x1(%r10),%r14d
test   $0x7,%dl
je     142e <inverse_transform_msimplex2+0x2de>
mov    %edx,%ecx
and    $0x7,%ecx
lea    0x0(,%rax,8),%rax
add    %rbp,%rax
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
addsd  (%rax,%rsi,8),%xmm0
add    $0x1,%rsi
cmp    %rsi,%rcx
jne    1420 <inverse_transform_msimplex2+0x2d0>
lea    0x1(%r14),%eax
cmp    %r12d,%eax
mov    %r10,0x68(%rsp)
jge    14b9 <inverse_transform_msimplex2+0x369>
mov    %r14d,%eax
add    $0x1,%rax
mov    0x28(%rsp),%esi
sub    %r14d,%esi
mov    0x24(%rsp),%r11d
sub    %r14d,%r11d
and    $0x7,%esi
je     146f <inverse_transform_msimplex2+0x31f>
nopw   0x0(%rax,%rax,1)
addsd  0x0(%rbp,%rax,8),%xmm0
add    $0x1,%rax
add    $0xffffffff,%esi
jne    1460 <inverse_transform_msimplex2+0x310>
cmp    $0x7,%r11d
jb     14b9 <inverse_transform_msimplex2+0x369>
cs nopw 0x0(%rax,%rax,1)
nop
addsd  0x0(%rbp,%rax,8),%xmm0
addsd  0x8(%rbp,%rax,8),%xmm0
addsd  0x10(%rbp,%rax,8),%xmm0
addsd  0x18(%rbp,%rax,8),%xmm0
addsd  0x20(%rbp,%rax,8),%xmm0
addsd  0x28(%rbp,%rax,8),%xmm0
addsd  0x30(%rbp,%rax,8),%xmm0
addsd  0x38(%rbp,%rax,8),%xmm0
add    $0x8,%rax
cmp    %eax,%r12d
jg     1480 <inverse_transform_msimplex2+0x330>
mulsd  0x0(%rbp,%rdx,8),%xmm0
divsd  %xmm1,%xmm0
mov    %edx,%eax
imul   0x10(%rsp),%eax
cltq
mov    0x18(%rsp),%rcx
movsd  %xmm0,(%rcx,%rax,8)
lea    0x1(%rdx),%rax
mov    %rax,0x60(%rsp)
cmp    %r8,%rax
mov    0x38(%rsp),%rax
jae    1330 <inverse_transform_msimplex2+0x1e0>
xor    %r13d,%r13d
mov    0x30(%rsp),%r12
mov    %rdx,0x8(%rsp)
nopw   0x0(%rax,%rax,1)
movsd  (%rax,%rdx,8),%xmm0
addsd  (%r15,%r13,8),%xmm0
mov    %r9,%r14
mov    %rdi,%rbp
call   1040 <exp@plt>
mov    0x38(%rsp),%rax
mov    0x8(%rsp),%rdx
movapd 0xae8(%rip),%xmm2
mov    %rbp,%rdi
mov    %r14,%r9
movsd  0x70(%rsp),%xmm1
xorpd  %xmm2,%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rbx,%r13,8)
movsd  %xmm0,(%r12)
add    %r14,%r12
add    $0x1,%r13
cmp    %r13,%rbp
jne    1500 <inverse_transform_msimplex2+0x3b0>
jmp    1330 <inverse_transform_msimplex2+0x1e0>
mov    0x18(%rsp),%rbx
mov    %rbx,%rdi
mov    %r12d,%esi
call   1070 <LUDecompose_det2@plt>
mov    %eax,%edi
call   1060 <log@plt>
mov    0x40(%rsp),%rax
movsd  %xmm0,(%rax)
mov    %rbx,%rdi
add    $0x78,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

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

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