void kernel_mvt(int n, double x1[2000 + 0], double x2[2000 + 0],
                double y_1[2000 + 0], double y_2[2000 + 0],
                double A[2000 + 0][2000 + 0]) {
  int i, j;
#pragma scop
  for (i = 0; i < n; i++)
    for (j = 0; j < n; j++)
      x1[i] = x1[i] + A[i][j] * y_1[j];
  for (i = 0; i < n; i++)
    for (j = 0; j < n; j++)
      x2[i] = x2[i] + A[j][i] * y_2[j];
#pragma endscop
}


x86 assembly (clang, -O0):
<kernel_mvt>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11ac <kernel_mvt+0xac>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1199 <kernel_mvt+0x99>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x3e80,%rcx,%rcx
add    %rcx,%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1135 <kernel_mvt+0x35>
jmp    119e <kernel_mvt+0x9e>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1122 <kernel_mvt+0x22>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    123d <kernel_mvt+0x13d>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    122a <kernel_mvt+0x12a>
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x30(%rbp),%rax
movslq -0x38(%rbp),%rcx
imul   $0x3e80,%rcx,%rcx
add    %rcx,%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    11c6 <kernel_mvt+0xc6>
jmp    122f <kernel_mvt+0x12f>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    11b3 <kernel_mvt+0xb3>
pop    %rbp
ret

x86 assembly (clang, -O2):
<kernel_mvt>:
push   %r15
push   %r14
push   %r12
push   %rbx
test   %edi,%edi
jle    125b <kernel_mvt+0x15b>
mov    %edi,%r11d
lea    -0x1(%r11),%r10
mov    %r11d,%r14d
and    $0xfffffffe,%r14d
lea    0x8(%r9),%r15
xor    %r12d,%r12d
jmp    1140 <kernel_mvt+0x40>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%r12
add    $0x3e80,%r15
cmp    %r11,%r12
je     11ab <kernel_mvt+0xab>
movsd  (%rsi,%r12,8),%xmm0
xor    %ebx,%ebx
test   %r10,%r10
je     1185 <kernel_mvt+0x85>
nopl   (%rax)
movsd  -0x8(%r15,%rbx,8),%xmm1
mulsd  (%rcx,%rbx,8),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rsi,%r12,8)
movsd  (%r15,%rbx,8),%xmm0
mulsd  0x8(%rcx,%rbx,8),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rsi,%r12,8)
add    $0x2,%rbx
cmp    %rbx,%r14
jne    1150 <kernel_mvt+0x50>
test   $0x1,%r11b
je     1130 <kernel_mvt+0x30>
imul   $0x3e80,%r12,%rax
add    %r9,%rax
movsd  (%rax,%rbx,8),%xmm1
mulsd  (%rcx,%rbx,8),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rsi,%r12,8)
jmp    1130 <kernel_mvt+0x30>
test   %edi,%edi
jle    125b <kernel_mvt+0x15b>
mov    %r11d,%ecx
and    $0xfffffffe,%ecx
xor    %esi,%esi
mov    %r9,%rax
jmp    11d1 <kernel_mvt+0xd1>
add    $0x1,%rsi
add    $0x8,%rax
cmp    %r11,%rsi
je     125b <kernel_mvt+0x15b>
movsd  (%rdx,%rsi,8),%xmm0
test   %r10,%r10
je     1230 <kernel_mvt+0x130>
mov    %rax,%rbx
xor    %edi,%edi
movsd  (%rbx),%xmm1
mulsd  (%r8,%rdi,8),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%rsi,8)
movsd  0x3e80(%rbx),%xmm0
mulsd  0x8(%r8,%rdi,8),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rdx,%rsi,8)
add    $0x2,%rdi
add    $0x7d00,%rbx
cmp    %rdi,%rcx
jne    11e0 <kernel_mvt+0xe0>
test   $0x1,%r11b
je     11c0 <kernel_mvt+0xc0>
jmp    1238 <kernel_mvt+0x138>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
xor    %edi,%edi
test   $0x1,%r11b
je     11c0 <kernel_mvt+0xc0>
imul   $0x3e80,%rdi,%rbx
add    %r9,%rbx
movsd  (%rbx,%rsi,8),%xmm1
mulsd  (%r8,%rdi,8),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%rsi,8)
jmp    11c0 <kernel_mvt+0xc0>
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define void @kernel_mvt(i32 %0, double* %1, double* %2, double* %3, double* %4, [2000 x double]* %5) {
%7 = alloca i32
%8 = alloca double*
%9 = alloca double*
%10 = alloca double*
%11 = alloca double*
%12 = alloca [2000 x double]*
%13 = alloca i32
%14 = alloca i32
store i32 %0, i32* %7
store double* %1, double** %8
store double* %2, double** %9
store double* %3, double** %10
store double* %4, double** %11
store [2000 x double]* %5, [2000 x double]** %12
store i32 0, i32* %13
br label %15
15: ; preds = %52, %6
%16 = load i32, i32* %13
%17 = load i32, i32* %7
%18 = icmp slt i32 %16, %17
br i1 %18, label %19, label %55
19: ; preds = %15
store i32 0, i32* %14
br label %20
20: ; preds = %48, %19
%21 = load i32, i32* %14
%22 = load i32, i32* %7
%23 = icmp slt i32 %21, %22
br i1 %23, label %24, label %51
24: ; preds = %20
%25 = load double*, double** %8
%26 = load i32, i32* %13
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds double, double* %25, i64 %27
%29 = load double, double* %28
%30 = load [2000 x double]*, [2000 x double]** %12
%31 = load i32, i32* %13
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds [2000 x double], [2000 x double]* %30, i64 %32
%34 = load i32, i32* %14
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds [2000 x double], [2000 x double]* %33, i64 0, i64 %35
%37 = load double, double* %36
%38 = load double*, double** %10
%39 = load i32, i32* %14
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds double, double* %38, i64 %40
%42 = load double, double* %41
%43 = call double @llvm.fmuladd.f64(double %37, double %42, double %29)
%44 = load double*, double** %8
%45 = load i32, i32* %13
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds double, double* %44, i64 %46
store double %43, double* %47
br label %48
48: ; preds = %24
%49 = load i32, i32* %14
%50 = add nsw i32 %49, 1
store i32 %50, i32* %14
br label %20
51: ; preds = %20
br label %52
52: ; preds = %51
%53 = load i32, i32* %13
%54 = add nsw i32 %53, 1
store i32 %54, i32* %13
br label %15
55: ; preds = %15
store i32 0, i32* %13
br label %56
56: ; preds = %93, %55
%57 = load i32, i32* %13
%58 = load i32, i32* %7
%59 = icmp slt i32 %57, %58
br i1 %59, label %60, label %96
60: ; preds = %56
store i32 0, i32* %14
br label %61
61: ; preds = %89, %60
%62 = load i32, i32* %14
%63 = load i32, i32* %7
%64 = icmp slt i32 %62, %63
br i1 %64, label %65, label %92
65: ; preds = %61
%66 = load double*, double** %9
%67 = load i32, i32* %13
%68 = sext i32 %67 to i64
%69 = getelementptr inbounds double, double* %66, i64 %68
%70 = load double, double* %69
%71 = load [2000 x double]*, [2000 x double]** %12
%72 = load i32, i32* %14
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds [2000 x double], [2000 x double]* %71, i64 %73
%75 = load i32, i32* %13
%76 = sext i32 %75 to i64
%77 = getelementptr inbounds [2000 x double], [2000 x double]* %74, i64 0, i64 %76
%78 = load double, double* %77
%79 = load double*, double** %11
%80 = load i32, i32* %14
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds double, double* %79, i64 %81
%83 = load double, double* %82
%84 = call double @llvm.fmuladd.f64(double %78, double %83, double %70)
%85 = load double*, double** %9
%86 = load i32, i32* %13
%87 = sext i32 %86 to i64
%88 = getelementptr inbounds double, double* %85, i64 %87
store double %84, double* %88
br label %89
89: ; preds = %65
%90 = load i32, i32* %14
%91 = add nsw i32 %90, 1
store i32 %91, i32* %14
br label %61
92: ; preds = %61
br label %93
93: ; preds = %92
%94 = load i32, i32* %13
%95 = add nsw i32 %94, 1
store i32 %95, i32* %13
br label %56
96: ; preds = %56
ret void
}

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