#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 exp(int);
int PMatJC69like(double P[], double t, int n) {
  int i;
  double pii = 1. / n + (1. - 1. / n) * exp(-n / (n - 1.) * t),
         pij = (1. - pii) / (n - 1.);
  for (i = 0; i < n * n; i++)
    P[i] = pij;
  for (i = 0; i < n; i++)
    P[i * n + i] = pii;
  return (0);
}


x86 assembly (clang, -O0):
<PMatJC69like>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
cvtsi2sdl -0x14(%rbp),%xmm1
movsd  0xecf(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
cvtsi2sdl -0x14(%rbp),%xmm0
movsd  0xeb9(%rip),%xmm1
divsd  %xmm0,%xmm1
movsd  0xead(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
xor    %eax,%eax
sub    -0x14(%rbp),%eax
cvtsi2sd %eax,%xmm0
cvtsi2sdl -0x14(%rbp),%xmm1
movsd  0xe8e(%rip),%xmm2
subsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mulsd  -0x10(%rbp),%xmm0
cvttsd2si %xmm0,%edi
call   1030 <exp@plt>
movsd  -0x38(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm1
cvtsi2sd %eax,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  0xe55(%rip),%xmm0
subsd  -0x20(%rbp),%xmm0
cvtsi2sdl -0x14(%rbp),%xmm1
movsd  0xe43(%rip),%xmm2
subsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
mov    -0x14(%rbp),%ecx
imul   -0x14(%rbp),%ecx
cmp    %ecx,%eax
jge    1203 <PMatJC69like+0xf3>
movsd  -0x28(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    11d1 <PMatJC69like+0xc1>
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    123f <PMatJC69like+0x12f>
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x18(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    120a <PMatJC69like+0xfa>
xor    %eax,%eax
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<PMatJC69like>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %esi,%r14d
mov    %rdi,%r15
cvtsi2sd %esi,%xmm2
movsd  0xed0(%rip),%xmm3
movapd %xmm3,%xmm1
divsd  %xmm2,%xmm1
movsd  %xmm1,0x10(%rsp)
subsd  %xmm1,%xmm3
movsd  %xmm3,0x8(%rsp)
mov    %esi,%eax
neg    %eax
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
addsd  0xead(%rip),%xmm2
movsd  %xmm2,(%rsp)
divsd  %xmm2,%xmm1
mulsd  %xmm0,%xmm1
cvttsd2si %xmm1,%edi
call   1030 <exp@plt>
movsd  0xe87(%rip),%xmm2
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mulsd  0x8(%rsp),%xmm0
addsd  0x10(%rsp),%xmm0
mov    %r14d,%ecx
imul   %r14d,%ecx
test   %ecx,%ecx
je     12cf <PMatJC69like+0x1bf>
subsd  %xmm0,%xmm2
divsd  (%rsp),%xmm2
mov    %ecx,%eax
cmp    $0x4,%ecx
jae    11b2 <PMatJC69like+0xa2>
xor    %ecx,%ecx
jmp    12c0 <PMatJC69like+0x1b0>
mov    %eax,%ecx
and    $0xfffffffc,%ecx
movapd %xmm2,%xmm1
unpcklpd %xmm2,%xmm1
lea    -0x4(%rcx),%rdi
mov    %rdi,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
mov    %esi,%edx
and    $0x7,%edx
cmp    $0x1c,%rdi
jae    11e0 <PMatJC69like+0xd0>
xor    %edi,%edi
jmp    1285 <PMatJC69like+0x175>
and    $0xfffffffffffffff8,%rsi
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
movupd %xmm1,(%r15,%rdi,8)
movupd %xmm1,0x10(%r15,%rdi,8)
movupd %xmm1,0x20(%r15,%rdi,8)
movupd %xmm1,0x30(%r15,%rdi,8)
movupd %xmm1,0x40(%r15,%rdi,8)
movupd %xmm1,0x50(%r15,%rdi,8)
movupd %xmm1,0x60(%r15,%rdi,8)
movupd %xmm1,0x70(%r15,%rdi,8)
movupd %xmm1,0x80(%r15,%rdi,8)
movupd %xmm1,0x90(%r15,%rdi,8)
movupd %xmm1,0xa0(%r15,%rdi,8)
movupd %xmm1,0xb0(%r15,%rdi,8)
movupd %xmm1,0xc0(%r15,%rdi,8)
movupd %xmm1,0xd0(%r15,%rdi,8)
movupd %xmm1,0xe0(%r15,%rdi,8)
movupd %xmm1,0xf0(%r15,%rdi,8)
add    $0x20,%rdi
add    $0xfffffffffffffff8,%rsi
jne    11f0 <PMatJC69like+0xe0>
test   %rdx,%rdx
je     12b4 <PMatJC69like+0x1a4>
lea    (%r15,%rdi,8),%rsi
add    $0x10,%rsi
shl    $0x5,%rdx
xor    %edi,%edi
nopl   0x0(%rax,%rax,1)
movupd %xmm1,-0x10(%rsi,%rdi,1)
movupd %xmm1,(%rsi,%rdi,1)
add    $0x20,%rdi
cmp    %rdi,%rdx
jne    12a0 <PMatJC69like+0x190>
cmp    %rax,%rcx
je     12cf <PMatJC69like+0x1bf>
nopl   0x0(%rax)
movsd  %xmm2,(%r15,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    12c0 <PMatJC69like+0x1b0>
test   %r14d,%r14d
jle    13b9 <PMatJC69like+0x2a9>
mov    %r14d,%r8d
lea    -0x1(%r8),%rax
mov    %r8d,%r10d
and    $0x3,%r10d
xor    %ecx,%ecx
cmp    $0x3,%rax
jb     138e <PMatJC69like+0x27e>
lea    (%r14,%r14,1),%esi
lea    0x0(,%r14,4),%r9d
lea    (%r14,%r14,2),%edx
and    $0xfffffffc,%r8d
neg    %r8
mov    $0x2,%r11d
mov    $0x1,%r12d
mov    $0x3,%r13d
xor    %ebx,%ebx
mov    %r14d,%eax
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movslq %esi,%rsi
cltq
movslq %ecx,%rcx
lea    (%rdi,%rcx,1),%rbp
movsd  %xmm0,(%r15,%rbp,8)
lea    (%r12,%rax,1),%rbp
movsd  %xmm0,(%r15,%rbp,8)
lea    (%r11,%rsi,1),%rbp
movslq %edx,%rdx
movsd  %xmm0,(%r15,%rbp,8)
lea    (%rdx,%r13,1),%rbp
movsd  %xmm0,(%r15,%rbp,8)
add    $0x4,%rdi
add    %r9d,%ecx
add    $0x4,%r11
add    %r9d,%esi
add    $0x4,%r12
add    %r9d,%eax
add    $0x4,%r13
add    %r9d,%edx
add    $0xfffffffffffffffc,%rbx
cmp    %rbx,%r8
jne    1330 <PMatJC69like+0x220>
neg    %rbx
mov    %rbx,%rcx
test   %r10,%r10
je     13b9 <PMatJC69like+0x2a9>
mov    %r14d,%eax
imul   %ecx,%eax
nopl   0x0(%rax)
cltq
lea    (%rcx,%rax,1),%rdx
movsd  %xmm0,(%r15,%rdx,8)
add    $0x1,%rcx
add    %r14d,%eax
add    $0xffffffffffffffff,%r10
jne    13a0 <PMatJC69like+0x290>
xor    %eax,%eax
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define i32 @PMatJC69like(double* %0, double %1, i32 %2) {
%4 = alloca double*
%5 = alloca double
%6 = alloca i32
%7 = alloca i32
%8 = alloca double
%9 = alloca double
store double* %0, double** %4
store double %1, double* %5
store i32 %2, i32* %6
%10 = load i32, i32* %6
%11 = sitofp i32 %10 to double
%12 = fdiv double 1.000000e+00, %11
%13 = load i32, i32* %6
%14 = sitofp i32 %13 to double
%15 = fdiv double 1.000000e+00, %14
%16 = fsub double 1.000000e+00, %15
%17 = load i32, i32* %6
%18 = sub nsw i32 0, %17
%19 = sitofp i32 %18 to double
%20 = load i32, i32* %6
%21 = sitofp i32 %20 to double
%22 = fsub double %21, 1.000000e+00
%23 = fdiv double %19, %22
%24 = load double, double* %5
%25 = fmul double %23, %24
%26 = fptosi double %25 to i32
%27 = call i32 @exp(i32 %26)
%28 = sitofp i32 %27 to double
%29 = call double @llvm.fmuladd.f64(double %16, double %28, double %12)
store double %29, double* %8
%30 = load double, double* %8
%31 = fsub double 1.000000e+00, %30
%32 = load i32, i32* %6
%33 = sitofp i32 %32 to double
%34 = fsub double %33, 1.000000e+00
%35 = fdiv double %31, %34
store double %35, double* %9
store i32 0, i32* %7
br label %36
36: ; preds = %48, %3
%37 = load i32, i32* %7
%38 = load i32, i32* %6
%39 = load i32, i32* %6
%40 = mul nsw i32 %38, %39
%41 = icmp slt i32 %37, %40
br i1 %41, label %42, label %51
42: ; preds = %36
%43 = load double, double* %9
%44 = load double*, double** %4
%45 = load i32, i32* %7
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds double, double* %44, i64 %46
store double %43, double* %47
br label %48
48: ; preds = %42
%49 = load i32, i32* %7
%50 = add nsw i32 %49, 1
store i32 %50, i32* %7
br label %36
51: ; preds = %36
store i32 0, i32* %7
br label %52
52: ; preds = %66, %51
%53 = load i32, i32* %7
%54 = load i32, i32* %6
%55 = icmp slt i32 %53, %54
br i1 %55, label %56, label %69
56: ; preds = %52
%57 = load double, double* %8
%58 = load double*, double** %4
%59 = load i32, i32* %7
%60 = load i32, i32* %6
%61 = mul nsw i32 %59, %60
%62 = load i32, i32* %7
%63 = add nsw i32 %61, %62
%64 = sext i32 %63 to i64
%65 = getelementptr inbounds double, double* %58, i64 %64
store double %57, double* %65
br label %66
66: ; preds = %56
%67 = load i32, i32* %7
%68 = add nsw i32 %67, 1
store i32 %68, i32* %7
br label %52
69: ; preds = %52
ret i32 0
}

LLVM IR (-O2):
define i32 @PMatJC69like(double* nocapture writeonly %0, double %1, i32 %2) {
%4 = sitofp i32 %2 to double
%5 = fdiv double 1.000000e+00, %4
%6 = fsub double 1.000000e+00, %5
%7 = sub nsw i32 0, %2
%8 = sitofp i32 %7 to double
%9 = fadd double %4, -1.000000e+00
%10 = fdiv double %8, %9
%11 = fmul double %10, %1
%12 = fptosi double %11 to i32
%13 = tail call i32 @exp(i32 %12)
%14 = sitofp i32 %13 to double
%15 = tail call double @llvm.fmuladd.f64(double %6, double %14, double %5)
%16 = fsub double 1.000000e+00, %15
%17 = fdiv double %16, %9
%18 = mul i32 %2, %2
%19 = icmp eq i32 %18, 0
br i1 %19, label %98, label %20
20: ; preds = %3
%21 = zext i32 %18 to i64
%22 = icmp ult i32 %18, 4
br i1 %22, label %96, label %23
23: ; preds = %20
%24 = and i64 %21, 4294967292
%25 = insertelement <2 x double> poison, double %17, i64 0
%26 = shufflevector <2 x double> %25, <2 x double> poison, <2 x i32> zeroinitializer
%27 = insertelement <2 x double> poison, double %17, i64 0
%28 = shufflevector <2 x double> %27, <2 x double> poison, <2 x i32> zeroinitializer
%29 = add nsw i64 %24, -4
%30 = lshr exact i64 %29, 2
%31 = add nuw nsw i64 %30, 1
%32 = and i64 %31, 7
%33 = icmp ult i64 %29, 28
br i1 %33, label %81, label %34
34: ; preds = %23
%35 = and i64 %31, 9223372036854775800
br label %36
36: ; preds = %36, %34
%37 = phi i64 [ 0, %34 ], [ %78, %36 ]
%38 = phi i64 [ 0, %34 ], [ %79, %36 ]
%39 = getelementptr inbounds double, double* %0, i64 %37
%40 = bitcast double* %39 to <2 x double>*
store <2 x double> %26, <2 x double>* %40
%41 = getelementptr inbounds double, double* %39, i64 2
%42 = bitcast double* %41 to <2 x double>*
store <2 x double> %28, <2 x double>* %42
%43 = or i64 %37, 4
%44 = getelementptr inbounds double, double* %0, i64 %43
%45 = bitcast double* %44 to <2 x double>*
store <2 x double> %26, <2 x double>* %45
%46 = getelementptr inbounds double, double* %44, i64 2
%47 = bitcast double* %46 to <2 x double>*
store <2 x double> %28, <2 x double>* %47
%48 = or i64 %37, 8
%49 = getelementptr inbounds double, double* %0, i64 %48
%50 = bitcast double* %49 to <2 x double>*
store <2 x double> %26, <2 x double>* %50
%51 = getelementptr inbounds double, double* %49, i64 2
%52 = bitcast double* %51 to <2 x double>*
store <2 x double> %28, <2 x double>* %52
%53 = or i64 %37, 12
%54 = getelementptr inbounds double, double* %0, i64 %53
%55 = bitcast double* %54 to <2 x double>*
store <2 x double> %26, <2 x double>* %55
%56 = getelementptr inbounds double, double* %54, i64 2
%57 = bitcast double* %56 to <2 x double>*
store <2 x double> %28, <2 x double>* %57
%58 = or i64 %37, 16
%59 = getelementptr inbounds double, double* %0, i64 %58
%60 = bitcast double* %59 to <2 x double>*
store <2 x double> %26, <2 x double>* %60
%61 = getelementptr inbounds double, double* %59, i64 2
%62 = bitcast double* %61 to <2 x double>*
store <2 x double> %28, <2 x double>* %62
%63 = or i64 %37, 20
%64 = getelementptr inbounds double, double* %0, i64 %63
%65 = bitcast double* %64 to <2 x double>*
store <2 x double> %26, <2 x double>* %65
%66 = getelementptr inbounds double, double* %64, i64 2
%67 = bitcast double* %66 to <2 x double>*
store <2 x double> %28, <2 x double>* %67
%68 = or i64 %37, 24
%69 = getelementptr inbounds double, double* %0, i64 %68
%70 = bitcast double* %69 to <2 x double>*
store <2 x double> %26, <2 x double>* %70
%71 = getelementptr inbounds double, double* %69, i64 2
%72 = bitcast double* %71 to <2 x double>*
store <2 x double> %28, <2 x double>* %72
%73 = or i64 %37, 28
%74 = getelementptr inbounds double, double* %0, i64 %73
%75 = bitcast double* %74 to <2 x double>*
store <2 x double> %26, <2 x double>* %75
%76 = getelementptr inbounds double, double* %74, i64 2
%77 = bitcast double* %76 to <2 x double>*
store <2 x double> %28, <2 x double>* %77
%78 = add nuw i64 %37, 32
%79 = add i64 %38, 8
%80 = icmp eq i64 %79, %35
br i1 %80, label %81, label %36
81: ; preds = %36, %23
%82 = phi i64 [ 0, %23 ], [ %78, %36 ]
%83 = icmp eq i64 %32, 0
br i1 %83, label %94, label %84
84: ; preds = %81, %84
%85 = phi i64 [ %91, %84 ], [ %82, %81 ]
%86 = phi i64 [ %92, %84 ], [ 0, %81 ]
%87 = getelementptr inbounds double, double* %0, i64 %85
%88 = bitcast double* %87 to <2 x double>*
store <2 x double> %26, <2 x double>* %88
%89 = getelementptr inbounds double, double* %87, i64 2
%90 = bitcast double* %89 to <2 x double>*
store <2 x double> %28, <2 x double>* %90
%91 = add nuw i64 %85, 4
%92 = add i64 %86, 1
%93 = icmp eq i64 %92, %32
br i1 %93, label %94, label %84
94: ; preds = %84, %81
%95 = icmp eq i64 %24, %21
br i1 %95, label %98, label %96
96: ; preds = %20, %94
%97 = phi i64 [ 0, %20 ], [ %24, %94 ]
br label %107
98: ; preds = %107, %94, %3
%99 = icmp sgt i32 %2, 0
br i1 %99, label %100, label %155
100: ; preds = %98
%101 = zext i32 %2 to i64
%102 = add nsw i64 %101, -1
%103 = and i64 %101, 3
%104 = icmp ult i64 %102, 3
br i1 %104, label %141, label %105
105: ; preds = %100
%106 = and i64 %101, 4294967292
br label %112
107: ; preds = %96, %107
%108 = phi i64 [ %110, %107 ], [ %97, %96 ]
%109 = getelementptr inbounds double, double* %0, i64 %108
store double %17, double* %109
%110 = add nuw nsw i64 %108, 1
%111 = icmp eq i64 %110, %21
br i1 %111, label %98, label %107
112: ; preds = %112, %105
%113 = phi i64 [ 0, %105 ], [ %138, %112 ]
%114 = phi i64 [ 0, %105 ], [ %139, %112 ]
%115 = trunc i64 %113 to i32
%116 = mul nsw i32 %115, %2
%117 = sext i32 %116 to i64
%118 = add nsw i64 %113, %117
%119 = getelementptr inbounds double, double* %0, i64 %118
store double %15, double* %119
%120 = or i64 %113, 1
%121 = trunc i64 %120 to i32
%122 = mul nsw i32 %121, %2
%123 = sext i32 %122 to i64
%124 = add nsw i64 %120, %123
%125 = getelementptr inbounds double, double* %0, i64 %124
store double %15, double* %125
%126 = or i64 %113, 2
%127 = trunc i64 %126 to i32
%128 = mul nsw i32 %127, %2
%129 = sext i32 %128 to i64
%130 = add nsw i64 %126, %129
%131 = getelementptr inbounds double, double* %0, i64 %130
store double %15, double* %131
%132 = or i64 %113, 3
%133 = trunc i64 %132 to i32
%134 = mul nsw i32 %133, %2
%135 = sext i32 %134 to i64
%136 = add nsw i64 %132, %135
%137 = getelementptr inbounds double, double* %0, i64 %136
store double %15, double* %137
%138 = add nuw nsw i64 %113, 4
%139 = add i64 %114, 4
%140 = icmp eq i64 %139, %106
br i1 %140, label %141, label %112
141: ; preds = %112, %100
%142 = phi i64 [ 0, %100 ], [ %138, %112 ]
%143 = icmp eq i64 %103, 0
br i1 %143, label %155, label %144
144: ; preds = %141, %144
%145 = phi i64 [ %152, %144 ], [ %142, %141 ]
%146 = phi i64 [ %153, %144 ], [ 0, %141 ]
%147 = trunc i64 %145 to i32
%148 = mul nsw i32 %147, %2
%149 = sext i32 %148 to i64
%150 = add nsw i64 %145, %149
%151 = getelementptr inbounds double, double* %0, i64 %150
store double %15, double* %151
%152 = add nuw nsw i64 %145, 1
%153 = add i64 %146, 1
%154 = icmp eq i64 %153, %103
br i1 %154, label %155, label %144
155: ; preds = %141, %144, %98
ret i32 0
}