#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*/ add_matvecmul(int, int, float const *, float const *, int,
                                 float *);
void compute_distances_1_nonpacked(int d, int nb, const float *a,
                                   const float *b, int ldb, float *dist2) {
  long i, j;
  double sum_d2, sum_c2;
  const float *dl;
  sum_c2 = 0;
  for (j = 0; j < d; j++) {
    sum_c2 += a[j] * a[j];
  }
  for (i = 0; i < nb; i++) {
    sum_d2 = 0.0;
    dl = b + ldb * i;
    for (j = 0; j < d; j++) {
      sum_d2 += dl[j] * dl[j];
    }
    dist2[i] = sum_d2 + sum_c2;
  }
  add_matvecmul(d, nb, a, b, ldb, dist2);
}


x86 assembly (clang, -O0):
<compute_distances_1_nonpacked>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
mov    %r9,-0x28(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
movslq -0x4(%rbp),%rcx
cmp    %rcx,%rax
jge    1188 <compute_distances_1_nonpacked+0x78>
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    113e <compute_distances_1_nonpacked+0x2e>
movq   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%rax
movslq -0x8(%rbp),%rcx
cmp    %rcx,%rax
jge    123f <compute_distances_1_nonpacked+0x12f>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
imul   -0x30(%rbp),%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x50(%rbp)
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
movslq -0x4(%rbp),%rcx
cmp    %rcx,%rax
jge    1213 <compute_distances_1_nonpacked+0x103>
mov    -0x50(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x50(%rbp),%rax
mov    -0x38(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
addsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    11c9 <compute_distances_1_nonpacked+0xb9>
movsd  -0x40(%rbp),%xmm0
addsd  -0x48(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x28(%rbp),%rax
mov    -0x30(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    1190 <compute_distances_1_nonpacked+0x80>
mov    -0x4(%rbp),%edi
mov    -0x8(%rbp),%esi
mov    -0x10(%rbp),%rdx
mov    -0x18(%rbp),%rcx
mov    -0x1c(%rbp),%r8d
mov    -0x28(%rbp),%r9
call   1030 <add_matvecmul@plt>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<compute_distances_1_nonpacked>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %r8d,-0x4(%rsp)
movslq %edi,%r15
test   %r15d,%r15d
jle    113e <compute_distances_1_nonpacked+0x2e>
lea    -0x1(%r15),%rbp
mov    %edi,%eax
and    $0x3,%eax
cmp    $0x3,%rbp
jae    1147 <compute_distances_1_nonpacked+0x37>
xorpd  %xmm0,%xmm0
xor    %ebx,%ebx
jmp    11b3 <compute_distances_1_nonpacked+0xa3>
xorpd  %xmm0,%xmm0
jmp    11db <compute_distances_1_nonpacked+0xcb>
mov    %r15,%rbp
and    $0xfffffffffffffffc,%rbp
xorpd  %xmm0,%xmm0
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  (%rdx,%rbx,4),%xmm1
movss  0x4(%rdx,%rbx,4),%xmm2
mulss  %xmm1,%xmm1
cvtss2sd %xmm1,%xmm1
addsd  %xmm0,%xmm1
mulss  %xmm2,%xmm2
xorps  %xmm0,%xmm0
cvtss2sd %xmm2,%xmm0
addsd  %xmm1,%xmm0
movss  0x8(%rdx,%rbx,4),%xmm1
mulss  %xmm1,%xmm1
cvtss2sd %xmm1,%xmm1
addsd  %xmm0,%xmm1
movss  0xc(%rdx,%rbx,4),%xmm0
mulss  %xmm0,%xmm0
cvtss2sd %xmm0,%xmm0
addsd  %xmm1,%xmm0
add    $0x4,%rbx
cmp    %rbx,%rbp
jne    1160 <compute_distances_1_nonpacked+0x50>
test   %rax,%rax
je     11db <compute_distances_1_nonpacked+0xcb>
lea    (%rdx,%rbx,4),%rbp
xor    %ebx,%ebx
xchg   %ax,%ax
movss  0x0(%rbp,%rbx,4),%xmm1
mulss  %xmm1,%xmm1
cvtss2sd %xmm1,%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rbx
cmp    %rbx,%rax
jne    11c0 <compute_distances_1_nonpacked+0xb0>
test   %esi,%esi
jle    12e1 <compute_distances_1_nonpacked+0x1d1>
movslq %esi,%r10
movslq -0x4(%rsp),%r11
lea    -0x1(%r15),%r8
mov    %edi,%ebx
and    $0x3,%ebx
and    $0xfffffffffffffffc,%r15
lea    0xc(%rcx),%rax
shl    $0x2,%r11
xor    %r12d,%r12d
mov    %rcx,%r13
jmp    1231 <compute_distances_1_nonpacked+0x121>
nopl   0x0(%rax,%rax,1)
addsd  %xmm0,%xmm1
cvtsd2ss %xmm1,%xmm1
movss  %xmm1,(%r9,%r12,4)
add    $0x1,%r12
add    %r11,%rax
add    %r11,%r13
cmp    %r10,%r12
je     12e1 <compute_distances_1_nonpacked+0x1d1>
xorpd  %xmm1,%xmm1
test   %edi,%edi
jle    1210 <compute_distances_1_nonpacked+0x100>
xor    %ebp,%ebp
cmp    $0x3,%r8
jb     12a3 <compute_distances_1_nonpacked+0x193>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movss  -0xc(%rax,%rbp,4),%xmm2
movss  -0x8(%rax,%rbp,4),%xmm3
mulss  %xmm2,%xmm2
cvtss2sd %xmm2,%xmm2
addsd  %xmm1,%xmm2
mulss  %xmm3,%xmm3
xorps  %xmm1,%xmm1
cvtss2sd %xmm3,%xmm1
addsd  %xmm2,%xmm1
movss  -0x4(%rax,%rbp,4),%xmm2
mulss  %xmm2,%xmm2
cvtss2sd %xmm2,%xmm2
addsd  %xmm1,%xmm2
movss  (%rax,%rbp,4),%xmm1
mulss  %xmm1,%xmm1
cvtss2sd %xmm1,%xmm1
addsd  %xmm2,%xmm1
add    $0x4,%rbp
cmp    %rbp,%r15
jne    1250 <compute_distances_1_nonpacked+0x140>
test   %rbx,%rbx
je     1210 <compute_distances_1_nonpacked+0x100>
lea    0x0(,%rbp,4),%rbp
add    %r13,%rbp
xor    %r14d,%r14d
nopw   0x0(%rax,%rax,1)
movss  0x0(%rbp,%r14,4),%xmm2
mulss  %xmm2,%xmm2
cvtss2sd %xmm2,%xmm2
addsd  %xmm2,%xmm1
add    $0x1,%r14
cmp    %r14,%rbx
jne    12c0 <compute_distances_1_nonpacked+0x1b0>
jmp    1210 <compute_distances_1_nonpacked+0x100>
mov    -0x4(%rsp),%r8d
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <add_matvecmul@plt>

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

LLVM IR (-O2):
define void @compute_distances_1_nonpacked(i32 %0, i32 %1, float* %2, float* %3, i32 %4, float* %5) {
%7 = sext i32 %0 to i64
%8 = icmp sgt i32 %0, 0
br i1 %8, label %9, label %32
9: ; preds = %6
%10 = add nsw i64 %7, -1
%11 = and i64 %7, 3
%12 = icmp ult i64 %10, 3
br i1 %12, label %15, label %13
13: ; preds = %9
%14 = and i64 %7, -4
br label %43
15: ; preds = %43, %9
%16 = phi double [ undef, %9 ], [ %69, %43 ]
%17 = phi double [ 0.000000e+00, %9 ], [ %69, %43 ]
%18 = phi i64 [ 0, %9 ], [ %70, %43 ]
%19 = icmp eq i64 %11, 0
br i1 %19, label %32, label %20
20: ; preds = %15, %20
%21 = phi double [ %28, %20 ], [ %17, %15 ]
%22 = phi i64 [ %29, %20 ], [ %18, %15 ]
%23 = phi i64 [ %30, %20 ], [ 0, %15 ]
%24 = getelementptr inbounds float, float* %2, i64 %22
%25 = load float, float* %24
%26 = fmul float %25, %25
%27 = fpext float %26 to double
%28 = fadd double %21, %27
%29 = add nuw nsw i64 %22, 1
%30 = add i64 %23, 1
%31 = icmp eq i64 %30, %11
br i1 %31, label %32, label %20
32: ; preds = %15, %20, %6
%33 = phi double [ 0.000000e+00, %6 ], [ %16, %15 ], [ %28, %20 ]
%34 = sext i32 %1 to i64
%35 = icmp sgt i32 %1, 0
br i1 %35, label %36, label %131
36: ; preds = %32
%37 = sext i32 %4 to i64
%38 = add nsw i64 %7, -1
%39 = and i64 %7, 3
%40 = icmp ult i64 %38, 3
%41 = and i64 %7, -4
%42 = icmp eq i64 %39, 0
br label %73
43: ; preds = %43, %13
%44 = phi double [ 0.000000e+00, %13 ], [ %69, %43 ]
%45 = phi i64 [ 0, %13 ], [ %70, %43 ]
%46 = phi i64 [ 0, %13 ], [ %71, %43 ]
%47 = getelementptr inbounds float, float* %2, i64 %45
%48 = load float, float* %47
%49 = fmul float %48, %48
%50 = fpext float %49 to double
%51 = fadd double %44, %50
%52 = or i64 %45, 1
%53 = getelementptr inbounds float, float* %2, i64 %52
%54 = load float, float* %53
%55 = fmul float %54, %54
%56 = fpext float %55 to double
%57 = fadd double %51, %56
%58 = or i64 %45, 2
%59 = getelementptr inbounds float, float* %2, i64 %58
%60 = load float, float* %59
%61 = fmul float %60, %60
%62 = fpext float %61 to double
%63 = fadd double %57, %62
%64 = or i64 %45, 3
%65 = getelementptr inbounds float, float* %2, i64 %64
%66 = load float, float* %65
%67 = fmul float %66, %66
%68 = fpext float %67 to double
%69 = fadd double %63, %68
%70 = add nuw nsw i64 %45, 4
%71 = add i64 %46, 4
%72 = icmp eq i64 %71, %14
br i1 %72, label %15, label %43
73: ; preds = %36, %124
%74 = phi i64 [ 0, %36 ], [ %129, %124 ]
%75 = mul nsw i64 %74, %37
%76 = getelementptr inbounds float, float* %3, i64 %75
br i1 %8, label %77, label %124
77: ; preds = %73
br i1 %40, label %108, label %78
78: ; preds = %77, %78
%79 = phi double [ %104, %78 ], [ 0.000000e+00, %77 ]
%80 = phi i64 [ %105, %78 ], [ 0, %77 ]
%81 = phi i64 [ %106, %78 ], [ 0, %77 ]
%82 = getelementptr inbounds float, float* %76, i64 %80
%83 = load float, float* %82
%84 = fmul float %83, %83
%85 = fpext float %84 to double
%86 = fadd double %79, %85
%87 = or i64 %80, 1
%88 = getelementptr inbounds float, float* %76, i64 %87
%89 = load float, float* %88
%90 = fmul float %89, %89
%91 = fpext float %90 to double
%92 = fadd double %86, %91
%93 = or i64 %80, 2
%94 = getelementptr inbounds float, float* %76, i64 %93
%95 = load float, float* %94
%96 = fmul float %95, %95
%97 = fpext float %96 to double
%98 = fadd double %92, %97
%99 = or i64 %80, 3
%100 = getelementptr inbounds float, float* %76, i64 %99
%101 = load float, float* %100
%102 = fmul float %101, %101
%103 = fpext float %102 to double
%104 = fadd double %98, %103
%105 = add nuw nsw i64 %80, 4
%106 = add i64 %81, 4
%107 = icmp eq i64 %106, %41
br i1 %107, label %108, label %78
108: ; preds = %78, %77
%109 = phi double [ undef, %77 ], [ %104, %78 ]
%110 = phi double [ 0.000000e+00, %77 ], [ %104, %78 ]
%111 = phi i64 [ 0, %77 ], [ %105, %78 ]
br i1 %42, label %124, label %112
112: ; preds = %108, %112
%113 = phi double [ %120, %112 ], [ %110, %108 ]
%114 = phi i64 [ %121, %112 ], [ %111, %108 ]
%115 = phi i64 [ %122, %112 ], [ 0, %108 ]
%116 = getelementptr inbounds float, float* %76, i64 %114
%117 = load float, float* %116
%118 = fmul float %117, %117
%119 = fpext float %118 to double
%120 = fadd double %113, %119
%121 = add nuw nsw i64 %114, 1
%122 = add i64 %115, 1
%123 = icmp eq i64 %122, %39
br i1 %123, label %124, label %112
124: ; preds = %108, %112, %73
%125 = phi double [ 0.000000e+00, %73 ], [ %109, %108 ], [ %120, %112 ]
%126 = fadd double %33, %125
%127 = fptrunc double %126 to float
%128 = getelementptr inbounds float, float* %5, i64 %74
store float %127, float* %128
%129 = add nuw nsw i64 %74, 1
%130 = icmp eq i64 %129, %34
br i1 %130, label %131, label %73
131: ; preds = %124, %32
%132 = tail call i32 @add_matvecmul(i32 %0, i32 %1, float* %2, float* %3, i32 %4, float* %5)
ret void
}