#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 /*<<< orphan*/ temp2;
typedef int /*<<< orphan*/ temp1;

/* Variables and functions */
int /*<<< orphan*/ gb__memcpy_4byte(float **, float **, int);
void gb_float44_mul(float (*out)[4], float (*mat1)[4], float (*mat2)[4]) {
  int i, j;
  float temp1[4][4], temp2[4][4];
  if (mat1 == out) {
    gb__memcpy_4byte(temp1, mat1, sizeof(temp1));
    mat1 = temp1;
  }
  if (mat2 == out) {
    gb__memcpy_4byte(temp2, mat2, sizeof(temp2));
    mat2 = temp2;
  }
  for (j = 0; j < 4; j++) {
    for (i = 0; i < 4; i++) {
      out[j][i] = mat1[0][i] * mat2[j][0] + mat1[1][i] * mat2[j][1] +
                  mat1[2][i] * mat2[j][2] + mat1[3][i] * mat2[j][3];
    }
  }
}


x86 assembly (clang, -O0):
<gb_float44_mul>:
push   %rbp
mov    %rsp,%rbp
sub    $0xa0,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x10(%rbp),%rax
cmp    -0x8(%rbp),%rax
jne    114f <gb_float44_mul+0x3f>
lea    -0x60(%rbp),%rdi
mov    -0x10(%rbp),%rsi
mov    $0x40,%edx
call   1030 <gb__memcpy_4byte@plt>
lea    -0x60(%rbp),%rax
mov    %rax,-0x10(%rbp)
mov    -0x18(%rbp),%rax
cmp    -0x8(%rbp),%rax
jne    117d <gb_float44_mul+0x6d>
lea    -0xa0(%rbp),%rdi
mov    -0x18(%rbp),%rsi
mov    $0x40,%edx
call   1030 <gb__memcpy_4byte@plt>
lea    -0xa0(%rbp),%rax
mov    %rax,-0x18(%rbp)
movl   $0x0,-0x20(%rbp)
cmpl   $0x4,-0x20(%rbp)
jge    1276 <gb_float44_mul+0x166>
movl   $0x0,-0x1c(%rbp)
cmpl   $0x4,-0x1c(%rbp)
jge    1263 <gb_float44_mul+0x153>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
movss  (%rax),%xmm2
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  0x10(%rax,%rcx,4),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
mulss  0x4(%rax),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  0x20(%rax,%rcx,4),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
movss  0x8(%rax),%xmm2
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  0x30(%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
movss  0xc(%rax),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
movslq -0x1c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1195 <gb_float44_mul+0x85>
jmp    1268 <gb_float44_mul+0x158>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1184 <gb_float44_mul+0x74>
add    $0xa0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<gb_float44_mul>:
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x88,%rsp
mov    %rdx,%r15
mov    %rsi,%rbx
mov    %rdi,%r14
cmp    %rdi,%rsi
je     1263 <gb_float44_mul+0x153>
cmp    %r14,%r15
je     1284 <gb_float44_mul+0x174>
xor    %eax,%eax
nopl   0x0(%rax,%rax,1)
movss  (%rbx),%xmm0
movss  0x10(%rbx),%xmm1
mulss  0x4(%r15,%rax,1),%xmm1
mulss  (%r15,%rax,1),%xmm0
addss  %xmm1,%xmm0
movss  0x20(%rbx),%xmm1
mulss  0x8(%r15,%rax,1),%xmm1
addss  %xmm0,%xmm1
movss  0x30(%rbx),%xmm0
mulss  0xc(%r15,%rax,1),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%r14,%rax,1)
movss  0x4(%rbx),%xmm0
movss  0x14(%rbx),%xmm1
mulss  0x4(%r15,%rax,1),%xmm1
mulss  (%r15,%rax,1),%xmm0
addss  %xmm1,%xmm0
movss  0x24(%rbx),%xmm1
mulss  0x8(%r15,%rax,1),%xmm1
addss  %xmm0,%xmm1
movss  0x34(%rbx),%xmm0
mulss  0xc(%r15,%rax,1),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,0x4(%r14,%rax,1)
movss  0x8(%rbx),%xmm0
movss  0x18(%rbx),%xmm1
mulss  0x4(%r15,%rax,1),%xmm1
mulss  (%r15,%rax,1),%xmm0
addss  %xmm1,%xmm0
movss  0x28(%rbx),%xmm1
mulss  0x8(%r15,%rax,1),%xmm1
addss  %xmm0,%xmm1
movss  0x38(%rbx),%xmm0
mulss  0xc(%r15,%rax,1),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,0x8(%r14,%rax,1)
movss  0xc(%rbx),%xmm0
movss  0x1c(%rbx),%xmm1
mulss  0x4(%r15,%rax,1),%xmm1
mulss  (%r15,%rax,1),%xmm0
addss  %xmm1,%xmm0
movss  0x2c(%rbx),%xmm1
mulss  0x8(%r15,%rax,1),%xmm1
addss  %xmm0,%xmm1
movss  0x3c(%rbx),%xmm0
mulss  0xc(%r15,%rax,1),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,0xc(%r14,%rax,1)
add    $0x10,%rax
cmp    $0x40,%rax
jne    1140 <gb_float44_mul+0x30>
add    $0x88,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret
lea    0x40(%rsp),%r12
mov    %r12,%rdi
mov    %rbx,%rsi
mov    $0x40,%edx
call   1030 <gb__memcpy_4byte@plt>
mov    %r12,%rbx
cmp    %r14,%r15
jne    1139 <gb_float44_mul+0x29>
mov    %rsp,%r15
mov    %r15,%rdi
mov    %r14,%rsi
mov    $0x40,%edx
call   1030 <gb__memcpy_4byte@plt>
jmp    1139 <gb_float44_mul+0x29>

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

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