#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*/ free_ivector(unsigned long *, int /*<<< orphan*/);
int /*<<< orphan*/ genmul(int, double *, long, unsigned long *);
unsigned long *ivector(int /*<<< orphan*/, long);
int /*<<< orphan*/ mexPrintf(char *, unsigned long);
void mdis_rndc(int n, int nevents, long ncat, double *p, double **array) {
  int i, j;
  unsigned long *a;
  a = ivector(0, ncat);
  for (i = 0; i < n; i++) {
    genmul(nevents, p, ncat, a);
    mexPrintf(" a = %ld \n", a[0]);
    for (j = 0; j < ncat; j++)
      array[i][j] = (double)a[j];
  }
  free_ivector(a, 0);
}


x86 assembly (clang, -O0):
<mdis_rndc>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    -0x10(%rbp),%rsi
xor    %edi,%edi
call   1050 <ivector@plt>
mov    %rax,-0x30(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    121d <mdis_rndc+0xdd>
mov    -0x8(%rbp),%edi
mov    -0x18(%rbp),%rsi
mov    -0x10(%rbp),%rdx
mov    -0x30(%rbp),%rcx
call   1060 <genmul@plt>
mov    -0x30(%rbp),%rax
mov    (%rax),%rsi
lea    0xe82(%rip),%rdi
call   1030 <mexPrintf@plt>
movl   $0x0,-0x28(%rbp)
movslq -0x28(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    120a <mdis_rndc+0xca>
mov    -0x30(%rbp),%rax
movslq -0x28(%rbp),%rcx
movq   (%rax,%rcx,8),%xmm0
movaps 0xe34(%rip),%xmm1
punpckldq %xmm1,%xmm0
movapd 0xe38(%rip),%xmm1
subpd  %xmm1,%xmm0
movaps %xmm0,%xmm1
unpckhpd %xmm0,%xmm0
addsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    11aa <mdis_rndc+0x6a>
jmp    120f <mdis_rndc+0xcf>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1170 <mdis_rndc+0x30>
mov    -0x30(%rbp),%rdi
xor    %esi,%esi
call   1040 <free_ivector@plt>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<mdis_rndc>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %r8,0x20(%rsp)
mov    %rcx,%r15
mov    %rdx,%rbx
mov    %esi,%r12d
mov    %edi,%r13d
xor    %edi,%edi
mov    %rdx,%rsi
call   1050 <ivector@plt>
mov    %rax,%rbp
test   %r13d,%r13d
jle    139c <mdis_rndc+0x25c>
mov    %r13d,%r13d
lea    -0x4(%rbx),%rax
mov    %rax,0x18(%rsp)
mov    %rax,%rcx
shr    $0x2,%rcx
add    $0x1,%rcx
mov    %rbx,%rax
and    $0xfffffffffffffffc,%rax
mov    %rax,0x8(%rsp)
mov    %rcx,%rax
mov    %rcx,0x10(%rsp)
and    $0xfffffffffffffffe,%rcx
mov    %rcx,(%rsp)
xor    %r14d,%r14d
jmp    11bd <mdis_rndc+0x7d>
nopl   (%rax)
add    $0x1,%r14
cmp    %r13,%r14
je     139c <mdis_rndc+0x25c>
mov    %r12d,%edi
mov    %r15,%rsi
mov    %rbx,%rdx
mov    %rbp,%rcx
call   1060 <genmul@plt>
mov    0x0(%rbp),%rsi
lea    0xe87(%rip),%rdi
call   1030 <mexPrintf@plt>
movapd 0xe49(%rip),%xmm8
movdqa 0xe31(%rip),%xmm7
movdqa 0xe19(%rip),%xmm6
movdqa 0xe01(%rip),%xmm5
movapd 0xe49(%rip),%xmm4
movapd 0xe31(%rip),%xmm3
test   %rbx,%rbx
jle    11b0 <mdis_rndc+0x70>
mov    0x20(%rsp),%rax
mov    (%rax,%r14,8),%rax
cmp    $0x4,%rbx
jae    1230 <mdis_rndc+0xf0>
xor    %ecx,%ecx
jmp    1370 <mdis_rndc+0x230>
nopw   0x0(%rax,%rax,1)
cmpq   $0x4,0x18(%rsp)
jae    123f <mdis_rndc+0xff>
xor    %ecx,%ecx
jmp    1305 <mdis_rndc+0x1c5>
mov    (%rsp),%rdx
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nop
movdqu 0x0(%rbp,%rcx,8),%xmm0
movdqu 0x10(%rbp,%rcx,8),%xmm1
movdqa %xmm0,%xmm2
pand   %xmm5,%xmm2
por    %xmm6,%xmm2
psrlq  $0x20,%xmm0
por    %xmm7,%xmm0
subpd  %xmm8,%xmm0
addpd  %xmm2,%xmm0
movdqa %xmm1,%xmm2
pand   %xmm5,%xmm2
por    %xmm6,%xmm2
psrlq  $0x20,%xmm1
por    %xmm7,%xmm1
subpd  %xmm8,%xmm1
addpd  %xmm2,%xmm1
movupd %xmm0,(%rax,%rcx,8)
movupd %xmm1,0x10(%rax,%rcx,8)
movdqu 0x20(%rbp,%rcx,8),%xmm0
movdqu 0x30(%rbp,%rcx,8),%xmm1
movdqa %xmm0,%xmm2
pand   %xmm5,%xmm2
por    %xmm6,%xmm2
psrlq  $0x20,%xmm0
por    %xmm7,%xmm0
subpd  %xmm8,%xmm0
addpd  %xmm2,%xmm0
movdqa %xmm1,%xmm2
pand   %xmm5,%xmm2
por    %xmm6,%xmm2
psrlq  $0x20,%xmm1
por    %xmm7,%xmm1
subpd  %xmm8,%xmm1
addpd  %xmm2,%xmm1
movupd %xmm0,0x20(%rax,%rcx,8)
movupd %xmm1,0x30(%rax,%rcx,8)
add    $0x8,%rcx
add    $0xfffffffffffffffe,%rdx
jne    1250 <mdis_rndc+0x110>
testb  $0x1,0x10(%rsp)
je     135f <mdis_rndc+0x21f>
movdqu 0x0(%rbp,%rcx,8),%xmm0
movdqu 0x10(%rbp,%rcx,8),%xmm1
movdqa %xmm0,%xmm2
pand   %xmm5,%xmm2
por    %xmm6,%xmm2
psrlq  $0x20,%xmm0
por    %xmm7,%xmm0
subpd  %xmm8,%xmm0
addpd  %xmm2,%xmm0
movdqa %xmm1,%xmm2
pand   %xmm5,%xmm2
por    %xmm6,%xmm2
psrlq  $0x20,%xmm1
por    %xmm7,%xmm1
subpd  %xmm8,%xmm1
addpd  %xmm2,%xmm1
movupd %xmm0,(%rax,%rcx,8)
movupd %xmm1,0x10(%rax,%rcx,8)
mov    0x8(%rsp),%rdx
mov    %rdx,%rcx
cmp    %rbx,%rdx
je     11b0 <mdis_rndc+0x70>
movsd  0x0(%rbp,%rcx,8),%xmm0
unpcklps %xmm3,%xmm0
subpd  %xmm4,%xmm0
movapd %xmm0,%xmm1
unpckhpd %xmm0,%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rax,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%rbx
jne    1370 <mdis_rndc+0x230>
jmp    11b0 <mdis_rndc+0x70>
mov    %rbp,%rdi
xor    %esi,%esi
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1040 <free_ivector@plt>

LLVM IR (-O0):
@.str = [11 x i8] c" a = %ld \0A\00"
define void @mdis_rndc(i32 %0, i32 %1, i64 %2, double* %3, double** %4) {
%6 = alloca i32
%7 = alloca i32
%8 = alloca i64
%9 = alloca double*
%10 = alloca double**
%11 = alloca i32
%12 = alloca i32
%13 = alloca i64*
store i32 %0, i32* %6
store i32 %1, i32* %7
store i64 %2, i64* %8
store double* %3, double** %9
store double** %4, double*** %10
%14 = load i64, i64* %8
%15 = call i64* @ivector(i32 0, i64 %14)
store i64* %15, i64** %13
store i32 0, i32* %11
br label %16
16: ; preds = %54, %5
%17 = load i32, i32* %11
%18 = load i32, i32* %6
%19 = icmp slt i32 %17, %18
br i1 %19, label %20, label %57
20: ; preds = %16
%21 = load i32, i32* %7
%22 = load double*, double** %9
%23 = load i64, i64* %8
%24 = load i64*, i64** %13
%25 = call i32 @genmul(i32 %21, double* %22, i64 %23, i64* %24)
%26 = load i64*, i64** %13
%27 = getelementptr inbounds i64, i64* %26, i64 0
%28 = load i64, i64* %27
%29 = call i32 @mexPrintf(i8* getelementptr inbounds ([11 x i8], [11 x i8]* @.str, i64 0, i64 0), i64 %28)
store i32 0, i32* %12
br label %30
30: ; preds = %50, %20
%31 = load i32, i32* %12
%32 = sext i32 %31 to i64
%33 = load i64, i64* %8
%34 = icmp slt i64 %32, %33
br i1 %34, label %35, label %53
35: ; preds = %30
%36 = load i64*, i64** %13
%37 = load i32, i32* %12
%38 = sext i32 %37 to i64
%39 = getelementptr inbounds i64, i64* %36, i64 %38
%40 = load i64, i64* %39
%41 = uitofp i64 %40 to double
%42 = load double**, double*** %10
%43 = load i32, i32* %11
%44 = sext i32 %43 to i64
%45 = getelementptr inbounds double*, double** %42, i64 %44
%46 = load double*, double** %45
%47 = load i32, i32* %12
%48 = sext i32 %47 to i64
%49 = getelementptr inbounds double, double* %46, i64 %48
store double %41, double* %49
br label %50
50: ; preds = %35
%51 = load i32, i32* %12
%52 = add nsw i32 %51, 1
store i32 %52, i32* %12
br label %30
53: ; preds = %30
br label %54
54: ; preds = %53
%55 = load i32, i32* %11
%56 = add nsw i32 %55, 1
store i32 %56, i32* %11
br label %16
57: ; preds = %16
%58 = load i64*, i64** %13
%59 = call i32 @free_ivector(i64* %58, i32 0)
ret void
}

LLVM IR (-O2):
@.str = [11 x i8] c" a = %ld \0A\00"
define void @mdis_rndc(i32 %0, i32 %1, i64 %2, double* %3, double** nocapture %4) {
%6 = tail call i64* @ivector(i32 0, i64 %2)
%7 = icmp sgt i32 %0, 0
br i1 %7, label %8, label %90
8: ; preds = %5
%9 = icmp sgt i64 %2, 0
%10 = zext i32 %0 to i64
%11 = add i64 %2, -4
%12 = lshr i64 %11, 2
%13 = add nuw nsw i64 %12, 1
%14 = icmp ult i64 %2, 4
%15 = and i64 %2, -4
%16 = and i64 %13, 1
%17 = icmp ult i64 %11, 4
%18 = and i64 %13, 9223372036854775806
%19 = icmp eq i64 %16, 0
%20 = icmp eq i64 %15, %2
br label %21
21: ; preds = %8, %87
%22 = phi i64 [ 0, %8 ], [ %88, %87 ]
%23 = tail call i32 @genmul(i32 %1, double* %3, i64 %2, i64* %6)
%24 = load i64, i64* %6
%25 = tail call i32 @mexPrintf(i8* getelementptr inbounds ([11 x i8], [11 x i8]* @.str, i64 0, i64 0), i64 %24)
br i1 %9, label %26, label %87
26: ; preds = %21
%27 = getelementptr inbounds double*, double** %4, i64 %22
%28 = load double*, double** %27
br i1 %14, label %77, label %29
29: ; preds = %26
br i1 %17, label %61, label %30
30: ; preds = %29, %30
%31 = phi i64 [ %58, %30 ], [ 0, %29 ]
%32 = phi i64 [ %59, %30 ], [ 0, %29 ]
%33 = getelementptr inbounds i64, i64* %6, i64 %31
%34 = bitcast i64* %33 to <2 x i64>*
%35 = load <2 x i64>, <2 x i64>* %34
%36 = getelementptr inbounds i64, i64* %33, i64 2
%37 = bitcast i64* %36 to <2 x i64>*
%38 = load <2 x i64>, <2 x i64>* %37
%39 = uitofp <2 x i64> %35 to <2 x double>
%40 = uitofp <2 x i64> %38 to <2 x double>
%41 = getelementptr inbounds double, double* %28, i64 %31
%42 = bitcast double* %41 to <2 x double>*
store <2 x double> %39, <2 x double>* %42
%43 = getelementptr inbounds double, double* %41, i64 2
%44 = bitcast double* %43 to <2 x double>*
store <2 x double> %40, <2 x double>* %44
%45 = or i64 %31, 4
%46 = getelementptr inbounds i64, i64* %6, i64 %45
%47 = bitcast i64* %46 to <2 x i64>*
%48 = load <2 x i64>, <2 x i64>* %47
%49 = getelementptr inbounds i64, i64* %46, i64 2
%50 = bitcast i64* %49 to <2 x i64>*
%51 = load <2 x i64>, <2 x i64>* %50
%52 = uitofp <2 x i64> %48 to <2 x double>
%53 = uitofp <2 x i64> %51 to <2 x double>
%54 = getelementptr inbounds double, double* %28, i64 %45
%55 = bitcast double* %54 to <2 x double>*
store <2 x double> %52, <2 x double>* %55
%56 = getelementptr inbounds double, double* %54, i64 2
%57 = bitcast double* %56 to <2 x double>*
store <2 x double> %53, <2 x double>* %57
%58 = add nuw i64 %31, 8
%59 = add i64 %32, 2
%60 = icmp eq i64 %59, %18
br i1 %60, label %61, label %30
61: ; preds = %30, %29
%62 = phi i64 [ 0, %29 ], [ %58, %30 ]
br i1 %19, label %76, label %63
63: ; preds = %61
%64 = getelementptr inbounds i64, i64* %6, i64 %62
%65 = bitcast i64* %64 to <2 x i64>*
%66 = load <2 x i64>, <2 x i64>* %65
%67 = getelementptr inbounds i64, i64* %64, i64 2
%68 = bitcast i64* %67 to <2 x i64>*
%69 = load <2 x i64>, <2 x i64>* %68
%70 = uitofp <2 x i64> %66 to <2 x double>
%71 = uitofp <2 x i64> %69 to <2 x double>
%72 = getelementptr inbounds double, double* %28, i64 %62
%73 = bitcast double* %72 to <2 x double>*
store <2 x double> %70, <2 x double>* %73
%74 = getelementptr inbounds double, double* %72, i64 2
%75 = bitcast double* %74 to <2 x double>*
store <2 x double> %71, <2 x double>* %75
br label %76
76: ; preds = %61, %63
br i1 %20, label %87, label %77
77: ; preds = %26, %76
%78 = phi i64 [ 0, %26 ], [ %15, %76 ]
br label %79
79: ; preds = %77, %79
%80 = phi i64 [ %85, %79 ], [ %78, %77 ]
%81 = getelementptr inbounds i64, i64* %6, i64 %80
%82 = load i64, i64* %81
%83 = uitofp i64 %82 to double
%84 = getelementptr inbounds double, double* %28, i64 %80
store double %83, double* %84
%85 = add nuw nsw i64 %80, 1
%86 = icmp eq i64 %85, %2
br i1 %86, label %87, label %79
87: ; preds = %79, %76, %21
%88 = add nuw nsw i64 %22, 1
%89 = icmp eq i64 %88, %10
br i1 %89, label %90, label %21
90: ; preds = %87, %5
%91 = tail call i32 @free_ivector(i64* %6, i32 0)
ret void
}