#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 */
int /*<<< orphan*/ fail(char *, int, int);
int scalar_product(float *ab, float *a, float *b, int na, int nb) {
  if (na != nb)
    fail("bad scalar product (%d %d)", na, nb);
  *ab = 0;
  for (int i = 0; i < na; i++)
    *ab += a[i] * b[i];
  return 1;
}


x86 assembly (clang, -O0):
<scalar_product>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x20(%rbp),%eax
je     1149 <scalar_product+0x39>
mov    -0x1c(%rbp),%esi
mov    -0x20(%rbp),%edx
lea    0xebc(%rip),%rdi
call   1030 <fail@plt>
mov    -0x8(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    11a3 <scalar_product+0x93>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    115b <scalar_product+0x4b>
mov    $0x1,%eax
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<scalar_product>:
push   %rbp
push   %r15
push   %r14
push   %rbx
push   %rax
mov    %ecx,%ebp
mov    %rdx,%r15
mov    %rsi,%r14
mov    %rdi,%rbx
cmp    %r8d,%ecx
je     1138 <scalar_product+0x28>
lea    0xed2(%rip),%rdi
mov    %ebp,%esi
mov    %r8d,%edx
call   1030 <fail@plt>
movl   $0x0,(%rbx)
test   %ebp,%ebp
jle    11fb <scalar_product+0xeb>
mov    %ebp,%ecx
lea    -0x1(%rcx),%rdx
mov    %ecx,%eax
and    $0x3,%eax
cmp    $0x3,%rdx
jae    115e <scalar_product+0x4e>
xorps  %xmm0,%xmm0
xor    %edx,%edx
jmp    11cf <scalar_product+0xbf>
and    $0xfffffffc,%ecx
xorps  %xmm0,%xmm0
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
movss  (%r14,%rdx,4),%xmm1
mulss  (%r15,%rdx,4),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%rbx)
movss  0x4(%r14,%rdx,4),%xmm0
mulss  0x4(%r15,%rdx,4),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rbx)
movss  0x8(%r14,%rdx,4),%xmm1
mulss  0x8(%r15,%rdx,4),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%rbx)
movss  0xc(%r14,%rdx,4),%xmm0
mulss  0xc(%r15,%rdx,4),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rbx)
add    $0x4,%rdx
cmp    %rdx,%rcx
jne    1170 <scalar_product+0x60>
test   %rax,%rax
je     11fb <scalar_product+0xeb>
lea    (%r15,%rdx,4),%rcx
lea    (%r14,%rdx,4),%rdx
xor    %esi,%esi
xchg   %ax,%ax
movss  (%rdx,%rsi,4),%xmm1
mulss  (%rcx,%rsi,4),%xmm1
addss  %xmm1,%xmm0
movss  %xmm0,(%rbx)
add    $0x1,%rsi
cmp    %rsi,%rax
jne    11e0 <scalar_product+0xd0>
mov    $0x1,%eax
add    $0x8,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@.str = [27 x i8] c"bad scalar product (%d %d)\00"
define i32 @scalar_product(float* %0, float* %1, float* %2, i32 %3, i32 %4) {
%6 = alloca float*
%7 = alloca float*
%8 = alloca float*
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
store float* %0, float** %6
store float* %1, float** %7
store float* %2, float** %8
store i32 %3, i32* %9
store i32 %4, i32* %10
%12 = load i32, i32* %9
%13 = load i32, i32* %10
%14 = icmp ne i32 %12, %13
br i1 %14, label %15, label %19
15: ; preds = %5
%16 = load i32, i32* %9
%17 = load i32, i32* %10
%18 = call i32 @fail(i8* getelementptr inbounds ([27 x i8], [27 x i8]* @.str, i64 0, i64 0), i32 %16, i32 %17)
br label %19
19: ; preds = %15, %5
%20 = load float*, float** %6
store float 0.000000e+00, float* %20
store i32 0, i32* %11
br label %21
21: ; preds = %39, %19
%22 = load i32, i32* %11
%23 = load i32, i32* %9
%24 = icmp slt i32 %22, %23
br i1 %24, label %25, label %42
25: ; preds = %21
%26 = load float*, float** %7
%27 = load i32, i32* %11
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds float, float* %26, i64 %28
%30 = load float, float* %29
%31 = load float*, float** %8
%32 = load i32, i32* %11
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds float, float* %31, i64 %33
%35 = load float, float* %34
%36 = load float*, float** %6
%37 = load float, float* %36
%38 = call float @llvm.fmuladd.f32(float %30, float %35, float %37)
store float %38, float* %36
br label %39
39: ; preds = %25
%40 = load i32, i32* %11
%41 = add nsw i32 %40, 1
store i32 %41, i32* %11
br label %21
42: ; preds = %21
ret i32 1
}

LLVM IR (-O2):
@.str = [27 x i8] c"bad scalar product (%d %d)\00"
define i32 @scalar_product(float* nocapture %0, float* nocapture %1, float* nocapture %2, i32 %3, i32 %4) {
%6 = icmp eq i32 %3, %4
br i1 %6, label %9, label %7
7: ; preds = %5
%8 = tail call i32 @fail(i8* getelementptr inbounds ([27 x i8], [27 x i8]* @.str, i64 0, i64 0), i32 %3, i32 %4)
br label %9
9: ; preds = %7, %5
store float 0.000000e+00, float* %0
%10 = icmp sgt i32 %3, 0
br i1 %10, label %11, label %34
11: ; preds = %9
%12 = zext i32 %3 to i64
%13 = add nsw i64 %12, -1
%14 = and i64 %12, 3
%15 = icmp ult i64 %13, 3
br i1 %15, label %18, label %16
16: ; preds = %11
%17 = and i64 %12, 4294967292
br label %35
18: ; preds = %35, %11
%19 = phi float [ 0.000000e+00, %11 ], [ %61, %35 ]
%20 = phi i64 [ 0, %11 ], [ %62, %35 ]
%21 = icmp eq i64 %14, 0
br i1 %21, label %34, label %22
22: ; preds = %18, %22
%23 = phi float [ %30, %22 ], [ %19, %18 ]
%24 = phi i64 [ %31, %22 ], [ %20, %18 ]
%25 = phi i64 [ %32, %22 ], [ 0, %18 ]
%26 = getelementptr inbounds float, float* %1, i64 %24
%27 = load float, float* %26
%28 = getelementptr inbounds float, float* %2, i64 %24
%29 = load float, float* %28
%30 = tail call float @llvm.fmuladd.f32(float %27, float %29, float %23)
store float %30, float* %0
%31 = add nuw nsw i64 %24, 1
%32 = add i64 %25, 1
%33 = icmp eq i64 %32, %14
br i1 %33, label %34, label %22
34: ; preds = %18, %22, %9
ret i32 1
35: ; preds = %35, %16
%36 = phi float [ 0.000000e+00, %16 ], [ %61, %35 ]
%37 = phi i64 [ 0, %16 ], [ %62, %35 ]
%38 = phi i64 [ 0, %16 ], [ %63, %35 ]
%39 = getelementptr inbounds float, float* %1, i64 %37
%40 = load float, float* %39
%41 = getelementptr inbounds float, float* %2, i64 %37
%42 = load float, float* %41
%43 = tail call float @llvm.fmuladd.f32(float %40, float %42, float %36)
store float %43, float* %0
%44 = or i64 %37, 1
%45 = getelementptr inbounds float, float* %1, i64 %44
%46 = load float, float* %45
%47 = getelementptr inbounds float, float* %2, i64 %44
%48 = load float, float* %47
%49 = tail call float @llvm.fmuladd.f32(float %46, float %48, float %43)
store float %49, float* %0
%50 = or i64 %37, 2
%51 = getelementptr inbounds float, float* %1, i64 %50
%52 = load float, float* %51
%53 = getelementptr inbounds float, float* %2, i64 %50
%54 = load float, float* %53
%55 = tail call float @llvm.fmuladd.f32(float %52, float %54, float %49)
store float %55, float* %0
%56 = or i64 %37, 3
%57 = getelementptr inbounds float, float* %1, i64 %56
%58 = load float, float* %57
%59 = getelementptr inbounds float, float* %2, i64 %56
%60 = load float, float* %59
%61 = tail call float @llvm.fmuladd.f32(float %58, float %60, float %55)
store float %61, float* %0
%62 = add nuw nsw i64 %37, 4
%63 = add i64 %38, 4
%64 = icmp eq i64 %63, %17
br i1 %64, label %18, label %35
}