#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 */
typedef int GLdouble;

/* Variables and functions */
void __gluMultMatricesd(const GLdouble a[16], const GLdouble b[16],
                        GLdouble r[16]) {
  int i, j;
  for (i = 0; i < 4; i++) {
    for (j = 0; j < 4; j++) {
      r[i * 4 + j] = a[i * 4 + 0] * b[0 * 4 + j] + a[i * 4 + 1] * b[1 * 4 + j] +
                     a[i * 4 + 2] * b[2 * 4 + j] + a[i * 4 + 3] * b[3 * 4 + j];
    }
  }
}


x86 assembly (clang, -O0):
<__gluMultMatricesd>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
cmpl   $0x4,-0x1c(%rbp)
jge    11fc <__gluMultMatricesd+0xfc>
movl   $0x0,-0x20(%rbp)
cmpl   $0x4,-0x20(%rbp)
jge    11e9 <__gluMultMatricesd+0xe9>
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
shl    $0x2,%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%edx
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
imul   (%rax,%rcx,4),%edx
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
shl    $0x2,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
mov    -0x10(%rbp),%rcx
mov    -0x20(%rbp),%esi
add    $0x4,%esi
movslq %esi,%rsi
imul   (%rcx,%rsi,4),%eax
add    %eax,%edx
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
shl    $0x2,%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
mov    -0x10(%rbp),%rcx
mov    -0x20(%rbp),%esi
add    $0x8,%esi
movslq %esi,%rsi
imul   (%rcx,%rsi,4),%eax
add    %eax,%edx
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
shl    $0x2,%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
mov    -0x10(%rbp),%rcx
mov    -0x20(%rbp),%esi
add    $0xc,%esi
movslq %esi,%rsi
imul   (%rcx,%rsi,4),%eax
add    %eax,%edx
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
shl    $0x2,%ecx
add    -0x20(%rbp),%ecx
movslq %ecx,%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1128 <__gluMultMatricesd+0x28>
jmp    11ee <__gluMultMatricesd+0xee>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1117 <__gluMultMatricesd+0x17>
pop    %rbp
ret

x86 assembly (clang, -O2):
<__gluMultMatricesd>:
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    (%rsi),%r8d
mov    0x10(%rsi),%r9d
imul   (%rdi,%rax,1),%r8d
imul   0x4(%rdi,%rax,1),%r9d
add    %r8d,%r9d
mov    0x20(%rsi),%r8d
imul   0x8(%rdi,%rax,1),%r8d
add    %r9d,%r8d
mov    0x30(%rsi),%ecx
imul   0xc(%rdi,%rax,1),%ecx
add    %r8d,%ecx
mov    %ecx,(%rdx,%rax,1)
mov    0x4(%rsi),%r8d
mov    0x14(%rsi),%r9d
imul   (%rdi,%rax,1),%r8d
imul   0x4(%rdi,%rax,1),%r9d
add    %r8d,%r9d
mov    0x24(%rsi),%r8d
imul   0x8(%rdi,%rax,1),%r8d
add    %r9d,%r8d
mov    0x34(%rsi),%ecx
imul   0xc(%rdi,%rax,1),%ecx
add    %r8d,%ecx
mov    %ecx,0x4(%rdx,%rax,1)
mov    0x8(%rsi),%r8d
mov    0x18(%rsi),%r9d
imul   (%rdi,%rax,1),%r8d
imul   0x4(%rdi,%rax,1),%r9d
add    %r8d,%r9d
mov    0x28(%rsi),%r8d
imul   0x8(%rdi,%rax,1),%r8d
add    %r9d,%r8d
mov    0x38(%rsi),%ecx
imul   0xc(%rdi,%rax,1),%ecx
add    %r8d,%ecx
mov    %ecx,0x8(%rdx,%rax,1)
mov    0xc(%rsi),%r8d
mov    0x1c(%rsi),%r9d
imul   (%rdi,%rax,1),%r8d
imul   0x4(%rdi,%rax,1),%r9d
add    %r8d,%r9d
mov    0x2c(%rsi),%r8d
imul   0x8(%rdi,%rax,1),%r8d
add    %r9d,%r8d
mov    0x3c(%rsi),%ecx
imul   0xc(%rdi,%rax,1),%ecx
add    %r8d,%ecx
mov    %ecx,0xc(%rdx,%rax,1)
add    $0x10,%rax
cmp    $0x40,%rax
jne    1110 <__gluMultMatricesd+0x10>
ret

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

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