#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 */
size_t mpos(int, int);
void D_eq9(int N, int M, double *C, double *D) {
  for (int i = 0; i < M; i++) {
    for (int j = i; j < M; j++) {
      double s = 0.0;
      for (int k = 0; k < N; k++) {
        s += C[k * M + i] * C[k * M + j];
      }
      D[mpos(i, j)] = s;
    }
  }
  return;
}


x86 assembly (clang, -O0):
<D_eq9>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    1201 <D_eq9+0xf1>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    11ee <D_eq9+0xde>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11b2 <D_eq9+0xa2>
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    -0x20(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x28(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    115a <D_eq9+0x4a>
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x18(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x1c(%rbp),%edi
mov    -0x20(%rbp),%esi
call   1030 <mpos@plt>
movsd  -0x40(%rbp),%xmm0
mov    %rax,%rcx
mov    -0x38(%rbp),%rax
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    113f <D_eq9+0x2f>
jmp    11f3 <D_eq9+0xe3>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    112d <D_eq9+0x1d>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<D_eq9>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
mov    %rcx,0x38(%rsp)
mov    %rdx,0x20(%rsp)
test   %esi,%esi
jle    125c <D_eq9+0x14c>
mov    %esi,%r13d
mov    %edi,%eax
mov    %rax,0x30(%rsp)
and    $0xfffffffe,%eax
mov    %rax,0x28(%rsp)
mov    0x20(%rsp),%rbp
lea    0x0(,%r13,8),%r14
add    %rbp,%r14
mov    %r13,0x10(%rsp)
shl    $0x4,%r13
xor    %ebx,%ebx
mov    %edi,0x1c(%rsp)
jmp    118c <D_eq9+0x7c>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    0x8(%rsp),%rbx
add    $0x1,%rbx
add    $0x8,%r14
add    $0x8,%rbp
cmp    0x10(%rsp),%rbx
je     125c <D_eq9+0x14c>
mov    %rbp,%r12
mov    %r14,%r15
mov    %rbx,0x8(%rsp)
jmp    11d9 <D_eq9+0xc9>
nopl   0x0(%rax)
mov    0x8(%rsp),%rdi
movsd  %xmm1,0x40(%rsp)
mov    %ebx,%esi
call   1030 <mpos@plt>
mov    0x38(%rsp),%rcx
movsd  0x40(%rsp),%xmm0
movsd  %xmm0,(%rcx,%rax,8)
add    $0x1,%rbx
add    $0x8,%r15
add    $0x8,%r12
cmp    0x10(%rsp),%rbx
mov    0x1c(%rsp),%edi
je     1170 <D_eq9+0x60>
xorpd  %xmm1,%xmm1
test   %edi,%edi
jle    11a0 <D_eq9+0x90>
cmp    $0x1,%edi
jne    11f0 <D_eq9+0xe0>
xor    %eax,%eax
mov    0x8(%rsp),%rdi
jmp    122c <D_eq9+0x11c>
nop
xor    %ecx,%ecx
xor    %eax,%eax
mov    0x28(%rsp),%rdx
mov    0x8(%rsp),%rdi
xchg   %ax,%ax
movsd  0x0(%rbp,%rcx,1),%xmm0
mulsd  (%r12,%rcx,1),%xmm0
addsd  %xmm1,%xmm0
movsd  (%r14,%rcx,1),%xmm1
mulsd  (%r15,%rcx,1),%xmm1
addsd  %xmm0,%xmm1
add    $0x2,%rax
add    %r13,%rcx
cmp    %rax,%rdx
jne    1200 <D_eq9+0xf0>
testb  $0x1,0x30(%rsp)
je     11a5 <D_eq9+0x95>
imul   0x10(%rsp),%rax
lea    (%rax,%rdi,1),%rcx
mov    0x20(%rsp),%rdx
movsd  (%rdx,%rcx,8),%xmm0
add    %rbx,%rax
mulsd  (%rdx,%rax,8),%xmm0
addsd  %xmm0,%xmm1
jmp    11a5 <D_eq9+0x95>
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @D_eq9(i32 %0, i32 %1, double* nocapture %2, double* nocapture writeonly %3) {
%5 = icmp sgt i32 %1, 0
br i1 %5, label %6, label %18
6: ; preds = %4
%7 = icmp sgt i32 %0, 0
%8 = zext i32 %1 to i64
%9 = zext i32 %1 to i64
%10 = zext i32 %0 to i64
%11 = and i64 %10, 1
%12 = icmp eq i32 %0, 1
%13 = and i64 %10, 4294967294
%14 = icmp eq i64 %11, 0
br label %15
15: ; preds = %22, %6
%16 = phi i64 [ 0, %6 ], [ %23, %22 ]
%17 = trunc i64 %16 to i32
br label %19
18: ; preds = %22, %4
ret void
19: ; preds = %15, %38
%20 = phi i64 [ %16, %15 ], [ %43, %38 ]
br i1 %7, label %21, label %38
21: ; preds = %19
br i1 %12, label %25, label %45
22: ; preds = %38
%23 = add nuw nsw i64 %16, 1
%24 = icmp eq i64 %23, %9
br i1 %24, label %18, label %15
25: ; preds = %45, %21
%26 = phi double [ undef, %21 ], [ %65, %45 ]
%27 = phi i64 [ 0, %21 ], [ %66, %45 ]
%28 = phi double [ 0.000000e+00, %21 ], [ %65, %45 ]
br i1 %14, label %38, label %29
29: ; preds = %25
%30 = mul nsw i64 %27, %8
%31 = add nuw nsw i64 %30, %16
%32 = getelementptr inbounds double, double* %2, i64 %31
%33 = load double, double* %32
%34 = add nuw nsw i64 %30, %20
%35 = getelementptr inbounds double, double* %2, i64 %34
%36 = load double, double* %35
%37 = tail call double @llvm.fmuladd.f64(double %33, double %36, double %28)
br label %38
38: ; preds = %29, %25, %19
%39 = phi double [ 0.000000e+00, %19 ], [ %26, %25 ], [ %37, %29 ]
%40 = trunc i64 %20 to i32
%41 = tail call i64 @mpos(i32 %17, i32 %40)
%42 = getelementptr inbounds double, double* %3, i64 %41
store double %39, double* %42
%43 = add nuw nsw i64 %20, 1
%44 = icmp eq i64 %43, %9
br i1 %44, label %22, label %19
45: ; preds = %21, %45
%46 = phi i64 [ %66, %45 ], [ 0, %21 ]
%47 = phi double [ %65, %45 ], [ 0.000000e+00, %21 ]
%48 = phi i64 [ %67, %45 ], [ 0, %21 ]
%49 = mul nsw i64 %46, %8
%50 = add nuw nsw i64 %49, %16
%51 = getelementptr inbounds double, double* %2, i64 %50
%52 = load double, double* %51
%53 = add nuw nsw i64 %49, %20
%54 = getelementptr inbounds double, double* %2, i64 %53
%55 = load double, double* %54
%56 = tail call double @llvm.fmuladd.f64(double %52, double %55, double %47)
%57 = or i64 %46, 1
%58 = mul nsw i64 %57, %8
%59 = add nuw nsw i64 %58, %16
%60 = getelementptr inbounds double, double* %2, i64 %59
%61 = load double, double* %60
%62 = add nuw nsw i64 %58, %20
%63 = getelementptr inbounds double, double* %2, i64 %62
%64 = load double, double* %63
%65 = tail call double @llvm.fmuladd.f64(double %61, double %64, double %56)
%66 = add nuw nsw i64 %46, 2
%67 = add i64 %48, 2
%68 = icmp eq i64 %67, %13
br i1 %68, label %25, label %45
}