#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*/ DoDiagonalOfQ(double *, int const);
int /*<<< orphan*/ assert(int);
void MakeQ_From_S_stdfreq(double *mat, const double *pi, const int n) {
  int i, j;
  assert(NULL != mat);
  assert(NULL != pi);
  assert(n > 0);
  for (i = 0; i < n; i++)
    for (j = 0; j < n; j++)
      mat[i * n + j] *= pi[j];
  DoDiagonalOfQ(mat, n);
}


x86 assembly (clang, -O0):
<MakeQ_From_S_stdfreq>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
xor    %eax,%eax
cmp    -0x8(%rbp),%rax
setne  %al
and    $0x1,%al
movzbl %al,%edi
call   1030 <assert@plt>
xor    %eax,%eax
cmp    -0x10(%rbp),%rax
setne  %al
and    $0x1,%al
movzbl %al,%edi
call   1030 <assert@plt>
cmpl   $0x0,-0x14(%rbp)
setg   %al
and    $0x1,%al
movzbl %al,%edi
call   1030 <assert@plt>
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11d9 <MakeQ_From_S_stdfreq+0xb9>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11c6 <MakeQ_From_S_stdfreq+0xa6>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1184 <MakeQ_From_S_stdfreq+0x64>
jmp    11cb <MakeQ_From_S_stdfreq+0xab>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1171 <MakeQ_From_S_stdfreq+0x51>
mov    -0x8(%rbp),%rdi
mov    -0x14(%rbp),%esi
call   1040 <DoDiagonalOfQ@plt>
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<MakeQ_From_S_stdfreq>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %edx,%r14d
mov    %rsi,%rbp
mov    %rdi,%r13
xor    %edi,%edi
test   %r13,%r13
setne  %dil
call   1030 <assert@plt>
xor    %edi,%edi
test   %rbp,%rbp
setne  %dil
call   1030 <assert@plt>
xor    %edi,%edi
test   %r14d,%r14d
setg   %dil
call   1030 <assert@plt>
test   %r14d,%r14d
jle    1358 <MakeQ_From_S_stdfreq+0x238>
mov    %r14d,%eax
lea    0x0(,%rax,8),%r8
add    %rbp,%r8
mov    %eax,%r15d
and    $0xfffffffc,%r15d
lea    -0x4(%r15),%rcx
mov    %rcx,0x10(%rsp)
mov    %rcx,%rdx
shr    $0x2,%rdx
add    $0x1,%rdx
mov    %rdx,%rcx
mov    %rdx,0x8(%rsp)
and    $0xfffffffffffffffe,%rdx
mov    %rdx,(%rsp)
mov    %rax,%r12
neg    %r12
lea    0x30(%r13),%rdi
lea    0x0(,%rax,8),%r10
lea    0x8(%r13),%rdx
xor    %r11d,%r11d
jmp    11d3 <MakeQ_From_S_stdfreq+0xb3>
xchg   %ax,%ax
add    $0x1,%r11
add    %r10,%rdi
add    %r10,%rdx
cmp    %rax,%r11
je     1358 <MakeQ_From_S_stdfreq+0x238>
mov    %r11,%rsi
imul   %rax,%rsi
cmp    $0x4,%r14d
jb     1208 <MakeQ_From_S_stdfreq+0xe8>
lea    0x0(,%rsi,8),%rcx
add    %r13,%rcx
cmp    %r8,%rcx
jae    1270 <MakeQ_From_S_stdfreq+0x150>
lea    (%rsi,%rax,1),%rcx
lea    0x0(,%rcx,8),%rcx
add    %r13,%rcx
cmp    %rbp,%rcx
jbe    1270 <MakeQ_From_S_stdfreq+0x150>
xor    %ebx,%ebx
mov    %rbx,%rcx
test   $0x1,%al
je     122f <MakeQ_From_S_stdfreq+0x10f>
movsd  0x0(%rbp,%rbx,8),%xmm0
add    %rbx,%rsi
mulsd  0x0(%r13,%rsi,8),%xmm0
movsd  %xmm0,0x0(%r13,%rsi,8)
mov    %rbx,%rcx
or     $0x1,%rcx
not    %rbx
cmp    %r12,%rbx
je     11c0 <MakeQ_From_S_stdfreq+0xa0>
nopw   0x0(%rax,%rax,1)
movsd  0x0(%rbp,%rcx,8),%xmm0
mulsd  -0x8(%rdx,%rcx,8),%xmm0
movsd  %xmm0,-0x8(%rdx,%rcx,8)
movsd  0x8(%rbp,%rcx,8),%xmm0
mulsd  (%rdx,%rcx,8),%xmm0
movsd  %xmm0,(%rdx,%rcx,8)
add    $0x2,%rcx
cmp    %rcx,%rax
jne    1240 <MakeQ_From_S_stdfreq+0x120>
jmp    11c0 <MakeQ_From_S_stdfreq+0xa0>
mov    %r8,%rcx
cmpq   $0x0,0x10(%rsp)
je     134c <MakeQ_From_S_stdfreq+0x22c>
mov    (%rsp),%r9
xor    %r8d,%r8d
cs nopw 0x0(%rax,%rax,1)
movupd 0x0(%rbp,%r8,8),%xmm0
movupd 0x10(%rbp,%r8,8),%xmm1
movupd -0x30(%rdi,%r8,8),%xmm2
mulpd  %xmm0,%xmm2
movupd -0x20(%rdi,%r8,8),%xmm0
mulpd  %xmm1,%xmm0
movupd -0x10(%rdi,%r8,8),%xmm1
movupd (%rdi,%r8,8),%xmm3
movupd %xmm2,-0x30(%rdi,%r8,8)
movupd %xmm0,-0x20(%rdi,%r8,8)
movupd 0x20(%rbp,%r8,8),%xmm0
mulpd  %xmm1,%xmm0
movupd 0x30(%rbp,%r8,8),%xmm1
mulpd  %xmm3,%xmm1
movupd %xmm0,-0x10(%rdi,%r8,8)
movupd %xmm1,(%rdi,%r8,8)
add    $0x8,%r8
add    $0xfffffffffffffffe,%r9
jne    1290 <MakeQ_From_S_stdfreq+0x170>
testb  $0x1,0x8(%rsp)
je     1338 <MakeQ_From_S_stdfreq+0x218>
movupd 0x0(%rbp,%r8,8),%xmm0
movupd 0x10(%rbp,%r8,8),%xmm1
add    %rsi,%r8
movupd 0x0(%r13,%r8,8),%xmm2
mulpd  %xmm0,%xmm2
movupd 0x10(%r13,%r8,8),%xmm0
mulpd  %xmm1,%xmm0
movupd %xmm2,0x0(%r13,%r8,8)
movupd %xmm0,0x10(%r13,%r8,8)
mov    %r15,%rbx
cmp    %rax,%r15
mov    %rcx,%r8
je     11c0 <MakeQ_From_S_stdfreq+0xa0>
jmp    120a <MakeQ_From_S_stdfreq+0xea>
xor    %r8d,%r8d
testb  $0x1,0x8(%rsp)
jne    1303 <MakeQ_From_S_stdfreq+0x1e3>
jmp    1338 <MakeQ_From_S_stdfreq+0x218>
mov    %r13,%rdi
mov    %r14d,%esi
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1040 <DoDiagonalOfQ@plt>

LLVM IR (-O0):
define void @MakeQ_From_S_stdfreq(double* %0, double* %1, i32 %2) {
%4 = alloca double*
%5 = alloca double*
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
store double* %0, double** %4
store double* %1, double** %5
store i32 %2, i32* %6
%9 = load double*, double** %4
%10 = icmp ne double* null, %9
%11 = zext i1 %10 to i32
%12 = call i32 @assert(i32 %11)
%13 = load double*, double** %5
%14 = icmp ne double* null, %13
%15 = zext i1 %14 to i32
%16 = call i32 @assert(i32 %15)
%17 = load i32, i32* %6
%18 = icmp sgt i32 %17, 0
%19 = zext i1 %18 to i32
%20 = call i32 @assert(i32 %19)
store i32 0, i32* %7
br label %21
21: ; preds = %50, %3
%22 = load i32, i32* %7
%23 = load i32, i32* %6
%24 = icmp slt i32 %22, %23
br i1 %24, label %25, label %53
25: ; preds = %21
store i32 0, i32* %8
br label %26
26: ; preds = %46, %25
%27 = load i32, i32* %8
%28 = load i32, i32* %6
%29 = icmp slt i32 %27, %28
br i1 %29, label %30, label %49
30: ; preds = %26
%31 = load double*, double** %5
%32 = load i32, i32* %8
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds double, double* %31, i64 %33
%35 = load double, double* %34
%36 = load double*, double** %4
%37 = load i32, i32* %7
%38 = load i32, i32* %6
%39 = mul nsw i32 %37, %38
%40 = load i32, i32* %8
%41 = add nsw i32 %39, %40
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds double, double* %36, i64 %42
%44 = load double, double* %43
%45 = fmul double %44, %35
store double %45, double* %43
br label %46
46: ; preds = %30
%47 = load i32, i32* %8
%48 = add nsw i32 %47, 1
store i32 %48, i32* %8
br label %26
49: ; preds = %26
br label %50
50: ; preds = %49
%51 = load i32, i32* %7
%52 = add nsw i32 %51, 1
store i32 %52, i32* %7
br label %21
53: ; preds = %21
%54 = load double*, double** %4
%55 = load i32, i32* %6
%56 = call i32 @DoDiagonalOfQ(double* %54, i32 %55)
ret void
}

LLVM IR (-O2):
define void @MakeQ_From_S_stdfreq(double* %0, double* %1, i32 %2) {
%4 = icmp ne double* %0, null
%5 = zext i1 %4 to i32
%6 = tail call i32 @assert(i32 %5)
%7 = icmp ne double* %1, null
%8 = zext i1 %7 to i32
%9 = tail call i32 @assert(i32 %8)
%10 = icmp sgt i32 %2, 0
%11 = zext i1 %10 to i32
%12 = tail call i32 @assert(i32 %11)
br i1 %10, label %13, label %139
13: ; preds = %3
%14 = zext i32 %2 to i64
%15 = zext i32 %2 to i64
%16 = getelementptr double, double* %1, i64 %15
%17 = and i64 %15, 4294967292
%18 = add nsw i64 %17, -4
%19 = lshr exact i64 %18, 2
%20 = add nuw nsw i64 %19, 1
%21 = icmp ult i32 %2, 4
%22 = and i64 %15, 4294967292
%23 = and i64 %20, 1
%24 = icmp eq i64 %18, 0
%25 = and i64 %20, 9223372036854775806
%26 = icmp eq i64 %23, 0
%27 = icmp eq i64 %22, %15
%28 = and i64 %15, 1
%29 = icmp eq i64 %28, 0
%30 = sub nsw i64 0, %15
br label %31
31: ; preds = %136, %13
%32 = phi i64 [ 0, %13 ], [ %137, %136 ]
%33 = mul i64 %32, %15
%34 = mul nsw i64 %32, %14
br i1 %21, label %105, label %35
35: ; preds = %31
%36 = add i64 %33, %15
%37 = getelementptr double, double* %0, i64 %36
%38 = getelementptr double, double* %0, i64 %33
%39 = icmp ult double* %38, %16
%40 = icmp ugt double* %37, %1
%41 = and i1 %39, %40
br i1 %41, label %105, label %42
42: ; preds = %35
br i1 %24, label %84, label %43
43: ; preds = %42, %43
%44 = phi i64 [ %81, %43 ], [ 0, %42 ]
%45 = phi i64 [ %82, %43 ], [ 0, %42 ]
%46 = getelementptr inbounds double, double* %1, i64 %44
%47 = bitcast double* %46 to <2 x double>*
%48 = load <2 x double>, <2 x double>* %47
%49 = getelementptr inbounds double, double* %46, i64 2
%50 = bitcast double* %49 to <2 x double>*
%51 = load <2 x double>, <2 x double>* %50
%52 = add nuw nsw i64 %44, %34
%53 = getelementptr inbounds double, double* %0, i64 %52
%54 = bitcast double* %53 to <2 x double>*
%55 = load <2 x double>, <2 x double>* %54
%56 = getelementptr inbounds double, double* %53, i64 2
%57 = bitcast double* %56 to <2 x double>*
%58 = load <2 x double>, <2 x double>* %57
%59 = fmul <2 x double> %48, %55
%60 = fmul <2 x double> %51, %58
%61 = bitcast double* %53 to <2 x double>*
store <2 x double> %59, <2 x double>* %61
%62 = bitcast double* %56 to <2 x double>*
store <2 x double> %60, <2 x double>* %62
%63 = or i64 %44, 4
%64 = getelementptr inbounds double, double* %1, i64 %63
%65 = bitcast double* %64 to <2 x double>*
%66 = load <2 x double>, <2 x double>* %65
%67 = getelementptr inbounds double, double* %64, i64 2
%68 = bitcast double* %67 to <2 x double>*
%69 = load <2 x double>, <2 x double>* %68
%70 = add nuw nsw i64 %63, %34
%71 = getelementptr inbounds double, double* %0, i64 %70
%72 = bitcast double* %71 to <2 x double>*
%73 = load <2 x double>, <2 x double>* %72
%74 = getelementptr inbounds double, double* %71, i64 2
%75 = bitcast double* %74 to <2 x double>*
%76 = load <2 x double>, <2 x double>* %75
%77 = fmul <2 x double> %66, %73
%78 = fmul <2 x double> %69, %76
%79 = bitcast double* %71 to <2 x double>*
store <2 x double> %77, <2 x double>* %79
%80 = bitcast double* %74 to <2 x double>*
store <2 x double> %78, <2 x double>* %80
%81 = add nuw i64 %44, 8
%82 = add i64 %45, 2
%83 = icmp eq i64 %82, %25
br i1 %83, label %84, label %43
84: ; preds = %43, %42
%85 = phi i64 [ 0, %42 ], [ %81, %43 ]
br i1 %26, label %104, label %86
86: ; preds = %84
%87 = getelementptr inbounds double, double* %1, i64 %85
%88 = bitcast double* %87 to <2 x double>*
%89 = load <2 x double>, <2 x double>* %88
%90 = getelementptr inbounds double, double* %87, i64 2
%91 = bitcast double* %90 to <2 x double>*
%92 = load <2 x double>, <2 x double>* %91
%93 = add nuw nsw i64 %85, %34
%94 = getelementptr inbounds double, double* %0, i64 %93
%95 = bitcast double* %94 to <2 x double>*
%96 = load <2 x double>, <2 x double>* %95
%97 = getelementptr inbounds double, double* %94, i64 2
%98 = bitcast double* %97 to <2 x double>*
%99 = load <2 x double>, <2 x double>* %98
%100 = fmul <2 x double> %89, %96
%101 = fmul <2 x double> %92, %99
%102 = bitcast double* %94 to <2 x double>*
store <2 x double> %100, <2 x double>* %102
%103 = bitcast double* %97 to <2 x double>*
store <2 x double> %101, <2 x double>* %103
br label %104
104: ; preds = %84, %86
br i1 %27, label %136, label %105
105: ; preds = %35, %31, %104
%106 = phi i64 [ 0, %35 ], [ 0, %31 ], [ %22, %104 ]
%107 = xor i64 %106, -1
br i1 %29, label %116, label %108
108: ; preds = %105
%109 = getelementptr inbounds double, double* %1, i64 %106
%110 = load double, double* %109
%111 = add nuw nsw i64 %106, %34
%112 = getelementptr inbounds double, double* %0, i64 %111
%113 = load double, double* %112
%114 = fmul double %110, %113
store double %114, double* %112
%115 = or i64 %106, 1
br label %116
116: ; preds = %108, %105
%117 = phi i64 [ %106, %105 ], [ %115, %108 ]
%118 = icmp eq i64 %107, %30
br i1 %118, label %136, label %119
119: ; preds = %116, %119
%120 = phi i64 [ %134, %119 ], [ %117, %116 ]
%121 = getelementptr inbounds double, double* %1, i64 %120
%122 = load double, double* %121
%123 = add nuw nsw i64 %120, %34
%124 = getelementptr inbounds double, double* %0, i64 %123
%125 = load double, double* %124
%126 = fmul double %122, %125
store double %126, double* %124
%127 = add nuw nsw i64 %120, 1
%128 = getelementptr inbounds double, double* %1, i64 %127
%129 = load double, double* %128
%130 = add nuw nsw i64 %127, %34
%131 = getelementptr inbounds double, double* %0, i64 %130
%132 = load double, double* %131
%133 = fmul double %129, %132
store double %133, double* %131
%134 = add nuw nsw i64 %120, 2
%135 = icmp eq i64 %134, %15
br i1 %135, label %136, label %119
136: ; preds = %116, %119, %104
%137 = add nuw nsw i64 %32, 1
%138 = icmp eq i64 %137, %15
br i1 %138, label %139, label %31
139: ; preds = %136, %3
%140 = tail call i32 @DoDiagonalOfQ(double* %0, i32 %2)
ret void
}