#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 */
double sum_square_matrix(const double *matrix, size_t size) {
  size_t i_j = 0;
  double sum = 0;
#pragma omp parallel for private(i_j) reduction(+ : sum)
  for (i_j = 0; i_j < size * size; i_j++) {
    sum += matrix[i_j];
  }
  return sum;
}


x86 assembly (clang, -O0):
<sum_square_matrix>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movq   $0x0,-0x18(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
movq   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x10(%rbp),%rcx
imul   -0x10(%rbp),%rcx
cmp    %rcx,%rax
jae    1162 <sum_square_matrix+0x62>
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x18(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x18(%rbp)
jmp    1124 <sum_square_matrix+0x24>
movsd  -0x20(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<sum_square_matrix>:
imul   %rsi,%rsi
test   %rsi,%rsi
je     1120 <sum_square_matrix+0x20>
lea    -0x1(%rsi),%rcx
mov    %esi,%eax
and    $0x7,%eax
cmp    $0x7,%rcx
jae    1124 <sum_square_matrix+0x24>
xorpd  %xmm0,%xmm0
xor    %ecx,%ecx
jmp    1168 <sum_square_matrix+0x68>
xorps  %xmm0,%xmm0
ret
and    $0xfffffffffffffff8,%rsi
xorpd  %xmm0,%xmm0
xor    %ecx,%ecx
xchg   %ax,%ax
addsd  (%rdi,%rcx,8),%xmm0
addsd  0x8(%rdi,%rcx,8),%xmm0
addsd  0x10(%rdi,%rcx,8),%xmm0
addsd  0x18(%rdi,%rcx,8),%xmm0
addsd  0x20(%rdi,%rcx,8),%xmm0
addsd  0x28(%rdi,%rcx,8),%xmm0
addsd  0x30(%rdi,%rcx,8),%xmm0
addsd  0x38(%rdi,%rcx,8),%xmm0
add    $0x8,%rcx
cmp    %rcx,%rsi
jne    1130 <sum_square_matrix+0x30>
test   %rax,%rax
je     118e <sum_square_matrix+0x8e>
lea    (%rdi,%rcx,8),%rcx
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
addsd  (%rcx,%rdx,8),%xmm0
add    $0x1,%rdx
cmp    %rdx,%rax
jne    1180 <sum_square_matrix+0x80>
ret

LLVM IR (-O0):
define double @sum_square_matrix(double* %0, i64 %1) {
%3 = alloca double*
%4 = alloca i64
%5 = alloca i64
%6 = alloca double
store double* %0, double** %3
store i64 %1, i64* %4
store i64 0, i64* %5
store double 0.000000e+00, double* %6
store i64 0, i64* %5
br label %7
7: ; preds = %20, %2
%8 = load i64, i64* %5
%9 = load i64, i64* %4
%10 = load i64, i64* %4
%11 = mul i64 %9, %10
%12 = icmp ult i64 %8, %11
br i1 %12, label %13, label %23
13: ; preds = %7
%14 = load double*, double** %3
%15 = load i64, i64* %5
%16 = getelementptr inbounds double, double* %14, i64 %15
%17 = load double, double* %16
%18 = load double, double* %6
%19 = fadd double %18, %17
store double %19, double* %6
br label %20
20: ; preds = %13
%21 = load i64, i64* %5
%22 = add i64 %21, 1
store i64 %22, i64* %5
br label %7
23: ; preds = %7
%24 = load double, double* %6
ret double %24
}

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