void sample_P_dp(int N, int M, int *cla_t, double alpha, double *P_dp) {
  int *M_arr;
  M_arr = (int *)malloc(M * sizeof(int));
  for (int m = 0; m < M; m++) {
    M_arr[m] = 0;
  }
  double *V_arr, sum;
  V_arr = (double *)malloc(M * sizeof(double));
  for (int m = 0; m <= M - 1; m++) {
    for (int n = 0; n < N; n++) {
      if (cla_t[n] == m) {
        M_arr[m] += 1;
      }
    }
  }
  for (int i = 0; i < M - 1; i++) {
    sum = 0.0;
    for (int j = i + 1; j < M; j++) {
      sum += M_arr[j];
    }
    V_arr[i] = rbeta(alpha / M + M_arr[i], alpha / M * (M - i - 1) + sum);
  }
  P_dp[0] = V_arr[0];
  for (int i = 1; i < M - 1; i++) {
    sum = 0;
    for (int j = 0; j < i; j++) {
      sum += P_dp[j];
    }
    P_dp[i] = (1 - sum) * V_arr[i];
  }
  sum = 0;
  for (int i = 0; i < M - 1; i++) {
    sum += P_dp[i];
  }
  P_dp[M - 1] = 1 - sum;
  free(M_arr);
  free(V_arr);
}


x86 assembly (clang, -O0):
<sample_P_dp>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movslq -0x8(%rbp),%rdi
shl    $0x2,%rdi
call   1050 <malloc@plt>
mov    %rax,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    118c <sample_P_dp+0x5c>
mov    -0x28(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movl   $0x0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1163 <sample_P_dp+0x33>
movslq -0x8(%rbp),%rdi
shl    $0x3,%rdi
call   1050 <malloc@plt>
mov    %rax,-0x38(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
mov    -0x8(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     1213 <sample_P_dp+0xe3>
movl   $0x0,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1200 <sample_P_dp+0xd0>
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
cmp    -0x44(%rbp),%eax
jne    11ed <sample_P_dp+0xbd>
mov    -0x28(%rbp),%rax
movslq -0x44(%rbp),%rcx
mov    (%rax,%rcx,4),%edx
add    $0x1,%edx
mov    %edx,(%rax,%rcx,4)
jmp    11f2 <sample_P_dp+0xc2>
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    11bc <sample_P_dp+0x8c>
jmp    1205 <sample_P_dp+0xd5>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    11a4 <sample_P_dp+0x74>
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
mov    -0x8(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    12da <sample_P_dp+0x1aa>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
mov    -0x50(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    126d <sample_P_dp+0x13d>
mov    -0x28(%rbp),%rax
movslq -0x50(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
addsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    123c <sample_P_dp+0x10c>
movsd  -0x18(%rbp),%xmm0
cvtsi2sdl -0x8(%rbp),%xmm1
divsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movslq -0x4c(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm1
addsd  %xmm1,%xmm0
movsd  -0x18(%rbp),%xmm1
cvtsi2sdl -0x8(%rbp),%xmm2
divsd  %xmm2,%xmm1
mov    -0x8(%rbp),%eax
sub    -0x4c(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm3
movsd  -0x40(%rbp),%xmm2
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
mov    $0x2,%al
call   1040 <rbeta@plt>
cvtsi2sd %eax,%xmm0
mov    -0x38(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    121a <sample_P_dp+0xea>
mov    -0x38(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x1,-0x54(%rbp)
mov    -0x54(%rbp),%eax
mov    -0x8(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1377 <sample_P_dp+0x247>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x58(%rbp)
mov    -0x58(%rbp),%eax
cmp    -0x54(%rbp),%eax
jge    1342 <sample_P_dp+0x212>
mov    -0x20(%rbp),%rax
movslq -0x58(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    1311 <sample_P_dp+0x1e1>
movsd  0xcb6(%rip),%xmm0
subsd  -0x40(%rbp),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x54(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    12f1 <sample_P_dp+0x1c1>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
mov    -0x8(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    13bc <sample_P_dp+0x28c>
mov    -0x20(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x5c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
jmp    1386 <sample_P_dp+0x256>
movsd  0xc3c(%rip),%xmm0
subsd  -0x40(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x8(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rdi
call   1030 <free@plt>
mov    -0x38(%rbp),%rdi
call   1030 <free@plt>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<sample_P_dp>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %rcx,%r13
movsd  %xmm0,0x20(%rsp)
mov    %rdx,%rbx
mov    %esi,%r15d
mov    %edi,%ebp
movslq %esi,%r12
lea    0x0(,%r12,4),%rdi
call   1060 <malloc@plt>
mov    %rax,%r14
test   %r12d,%r12d
jle    12fa <sample_P_dp+0x1ba>
mov    %r15,0x18(%rsp)
mov    %r12,0x28(%rsp)
mov    %r12d,%r12d
lea    0x0(,%r12,4),%rdx
xor    %r15d,%r15d
mov    %r14,%rdi
xor    %esi,%esi
call   1050 <memset@plt>
mov    0x28(%rsp),%rax
lea    0x0(,%rax,8),%rdi
call   1060 <malloc@plt>
mov    %rax,0x10(%rsp)
mov    %ebp,%eax
mov    %eax,%ecx
and    $0xfffffffe,%ecx
jmp    11c9 <sample_P_dp+0x89>
nopl   (%rax)
add    $0x1,%r15
cmp    %r12,%r15
je     1216 <sample_P_dp+0xd6>
test   %ebp,%ebp
jle    11c0 <sample_P_dp+0x80>
xor    %edx,%edx
cmp    $0x1,%ebp
jne    11f9 <sample_P_dp+0xb9>
test   $0x1,%al
je     11c0 <sample_P_dp+0x80>
mov    (%rbx,%rdx,4),%edx
cmp    %rdx,%r15
jne    11c0 <sample_P_dp+0x80>
addl   $0x1,(%r14,%r15,4)
jmp    11c0 <sample_P_dp+0x80>
nopw   0x0(%rax,%rax,1)
add    $0x2,%rdx
cmp    %rdx,%rcx
je     11d4 <sample_P_dp+0x94>
mov    (%rbx,%rdx,4),%esi
cmp    %rsi,%r15
jne    1206 <sample_P_dp+0xc6>
addl   $0x1,(%r14,%r15,4)
mov    0x4(%rbx,%rdx,4),%esi
cmp    %rsi,%r15
jne    11f0 <sample_P_dp+0xb0>
addl   $0x1,(%r14,%r15,4)
jmp    11f0 <sample_P_dp+0xb0>
mov    0x18(%rsp),%rcx
lea    -0x1(%rcx),%eax
xorpd  %xmm0,%xmm0
cmp    $0x2,%ecx
jl     131a <sample_P_dp+0x1da>
xorps  %xmm0,%xmm0
cvtsi2sd %ecx,%xmm0
movsd  0x20(%rsp),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x20(%rsp)
mov    %eax,0xc(%rsp)
mov    %eax,%r8d
xor    %r15d,%r15d
mov    %r14,%rbx
mov    %r8,0x30(%rsp)
jmp    12b9 <sample_P_dp+0x179>
cs nopw 0x0(%rax,%rax,1)
xorps  %xmm0,%xmm0
cvtsi2sdl (%r14,%r15,4),%xmm0
movsd  0x20(%rsp),%xmm3
addsd  %xmm3,%xmm0
mov    %r15d,%eax
not    %eax
add    0x18(%rsp),%eax
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
mov    $0x2,%al
call   1040 <rbeta@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mov    0x10(%rsp),%rax
movsd  %xmm0,(%rax,%r15,8)
add    $0x4,%rbx
add    $0xffffffffffffffff,%r12
mov    %rbp,%r15
mov    0x30(%rsp),%r8
cmp    %r8,%rbp
je     1324 <sample_P_dp+0x1e4>
lea    0x1(%r15),%rbp
xorpd  %xmm2,%xmm2
cmp    0x28(%rsp),%rbp
jge    1260 <sample_P_dp+0x120>
xorpd  %xmm2,%xmm2
mov    $0x1,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
xorps  %xmm0,%xmm0
cvtsi2sdl (%rbx,%rax,4),%xmm0
addsd  %xmm0,%xmm2
add    $0x1,%rax
cmp    %rax,%r12
jne    12e0 <sample_P_dp+0x1a0>
jmp    1260 <sample_P_dp+0x120>
lea    0x0(,%r12,8),%rdi
call   1060 <malloc@plt>
mov    %rax,%rbx
add    $0xffffffff,%r12d
xorpd  %xmm0,%xmm0
mov    %r12d,%eax
jmp    14d2 <sample_P_dp+0x392>
mov    0x10(%rsp),%rbx
jmp    14d2 <sample_P_dp+0x392>
mov    0x10(%rsp),%rbx
movsd  (%rbx),%xmm0
movsd  %xmm0,0x0(%r13)
mov    0x18(%rsp),%rbp
cmp    $0x3,%ebp
jl     1423 <sample_P_dp+0x2e3>
mov    $0x1,%eax
xor    %ecx,%ecx
movsd  0xcb0(%rip),%xmm0
jmp    1385 <sample_P_dp+0x245>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movapd %xmm0,%xmm2
subsd  %xmm1,%xmm2
mulsd  (%rbx,%rax,8),%xmm2
movsd  %xmm2,0x0(%r13,%rax,8)
add    $0x1,%rax
add    $0x1,%rcx
cmp    %r8,%rax
je     1423 <sample_P_dp+0x2e3>
cmp    $0x7,%rcx
jae    13a0 <sample_P_dp+0x260>
xorpd  %xmm1,%xmm1
xor    %edx,%edx
jmp    13f1 <sample_P_dp+0x2b1>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %rax,%rsi
and    $0xfffffffffffffff8,%rsi
xorpd  %xmm1,%xmm1
xor    %edx,%edx
nopl   (%rax)
addsd  0x0(%r13,%rdx,8),%xmm1
addsd  0x8(%r13,%rdx,8),%xmm1
addsd  0x10(%r13,%rdx,8),%xmm1
addsd  0x18(%r13,%rdx,8),%xmm1
addsd  0x20(%r13,%rdx,8),%xmm1
addsd  0x28(%r13,%rdx,8),%xmm1
addsd  0x30(%r13,%rdx,8),%xmm1
addsd  0x38(%r13,%rdx,8),%xmm1
add    $0x8,%rdx
cmp    %rdx,%rsi
jne    13b0 <sample_P_dp+0x270>
test   $0x7,%al
je     1360 <sample_P_dp+0x220>
mov    %eax,%esi
and    $0x7,%esi
lea    0x0(,%rdx,8),%rdx
add    %r13,%rdx
xor    %edi,%edi
nopl   0x0(%rax,%rax,1)
addsd  (%rdx,%rdi,8),%xmm1
add    $0x1,%rdi
cmp    %rdi,%rsi
jne    1410 <sample_P_dp+0x2d0>
jmp    1360 <sample_P_dp+0x220>
xorpd  %xmm0,%xmm0
cmp    $0x2,%ebp
jl     14ce <sample_P_dp+0x38e>
lea    -0x1(%r8),%rcx
mov    %r8d,%eax
and    $0x7,%eax
cmp    $0x7,%rcx
jae    1448 <sample_P_dp+0x308>
xorpd  %xmm0,%xmm0
xor    %ecx,%ecx
jmp    14a1 <sample_P_dp+0x361>
and    $0xfffffff8,%r8d
xorpd  %xmm0,%xmm0
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
addsd  0x0(%r13,%rcx,8),%xmm0
addsd  0x8(%r13,%rcx,8),%xmm0
addsd  0x10(%r13,%rcx,8),%xmm0
addsd  0x18(%r13,%rcx,8),%xmm0
addsd  0x20(%r13,%rcx,8),%xmm0
addsd  0x28(%r13,%rcx,8),%xmm0
addsd  0x30(%r13,%rcx,8),%xmm0
addsd  0x38(%r13,%rcx,8),%xmm0
add    $0x8,%rcx
cmp    %rcx,%r8
jne    1460 <sample_P_dp+0x320>
test   %rax,%rax
je     14ce <sample_P_dp+0x38e>
lea    0x0(,%rcx,8),%rcx
add    %r13,%rcx
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
addsd  (%rcx,%rdx,8),%xmm0
add    $0x1,%rdx
cmp    %rdx,%rax
jne    14c0 <sample_P_dp+0x380>
mov    0xc(%rsp),%eax
movsd  0xb26(%rip),%xmm1
subsd  %xmm0,%xmm1
cltq
movsd  %xmm1,0x0(%r13,%rax,8)
mov    %r14,%rdi
call   1030 <free@plt>
mov    %rbx,%rdi
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

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

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