#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 */
typedef struct TYPE_12__ TYPE_1__;

/* Type definitions */
struct TYPE_12__ {
  int /*<<< orphan*/ member_3;
  int /*<<< orphan*/ member_2;
  int /*<<< orphan*/ member_1;
  int /*<<< orphan*/ member_0;
};
typedef TYPE_1__ __m128;

/* Variables and functions */
TYPE_1__ _mm_add_ps(TYPE_1__, int /*<<< orphan*/);
int /*<<< orphan*/ _mm_mul_ps(TYPE_1__, TYPE_1__);
int /*<<< orphan*/ _mm_store_ps(float *, TYPE_1__);
TYPE_1__ _mm_sub_ps(TYPE_1__, TYPE_1__);
float sqrt(float);
float distance_l2_m128(__m128 *p1, __m128 *p2, int len) {
  float ret = 0;
  __m128 sum = {0, 0, 0, 0};
  float fsum[4];
  int i = 0;
  for (i = 0; i < len / 4; i++) {
    __m128 diff;
    diff = _mm_sub_ps(p1[i], p2[i]);
    sum = _mm_add_ps(sum, _mm_mul_ps(diff, diff));
  }
  _mm_store_ps(fsum, sum);
  ret = fsum[0] + fsum[1] + fsum[2] + fsum[3];
  return sqrt(ret);
}


x86 assembly (clang, -O0):
<distance_l2_m128>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x18(%rbp)
lea    -0x28(%rbp),%rdi
xor    %esi,%esi
mov    $0x10,%edx
call   1030 <memset@plt>
movl   $0x0,-0x44(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
mov    %eax,-0x7c(%rbp)
mov    -0x14(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x7c(%rbp),%eax
cmp    %ecx,%eax
jge    124d <distance_l2_m128+0xed>
mov    -0x8(%rbp),%rcx
movslq -0x44(%rbp),%rax
shl    $0x4,%rax
add    %rax,%rcx
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rdx
shl    $0x4,%rdx
add    %rdx,%rax
mov    (%rcx),%rdi
mov    0x8(%rcx),%rsi
mov    (%rax),%rdx
mov    0x8(%rax),%rcx
call   1050 <_mm_sub_ps@plt>
mov    %rax,-0x68(%rbp)
mov    %rdx,-0x60(%rbp)
mov    -0x68(%rbp),%rax
mov    %rax,-0x58(%rbp)
mov    -0x60(%rbp),%rax
mov    %rax,-0x50(%rbp)
mov    -0x58(%rbp),%rdi
mov    -0x50(%rbp),%rsi
mov    -0x58(%rbp),%rdx
mov    -0x50(%rbp),%rcx
call   1040 <_mm_mul_ps@plt>
mov    %eax,%edx
mov    -0x28(%rbp),%rdi
mov    -0x20(%rbp),%rsi
call   1070 <_mm_add_ps@plt>
mov    %rax,-0x78(%rbp)
mov    %rdx,-0x70(%rbp)
mov    -0x78(%rbp),%rax
mov    %rax,-0x28(%rbp)
mov    -0x70(%rbp),%rax
mov    %rax,-0x20(%rbp)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    119c <distance_l2_m128+0x3c>
lea    -0x40(%rbp),%rdi
mov    -0x28(%rbp),%rsi
mov    -0x20(%rbp),%rdx
call   1060 <_mm_store_ps@plt>
movss  -0x40(%rbp),%xmm0
addss  -0x3c(%rbp),%xmm0
addss  -0x38(%rbp),%xmm0
addss  -0x34(%rbp),%xmm0
movss  %xmm0,-0x18(%rbp)
movss  -0x18(%rbp),%xmm0
call   1080 <sqrt@plt>
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<distance_l2_m128>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
cmp    $0x4,%edx
jl     11c6 <distance_l2_m128+0x76>
mov    %edx,%r14d
mov    %rsi,%r15
mov    %rdi,%r12
shr    $0x2,%r14d
shl    $0x4,%r14
xor    %ebp,%ebp
xor    %ebx,%ebx
xor    %r13d,%r13d
nopl   0x0(%rax,%rax,1)
mov    (%r12,%rbp,1),%rdi
mov    0x8(%r12,%rbp,1),%rsi
mov    (%r15,%rbp,1),%rdx
mov    0x8(%r15,%rbp,1),%rcx
call   1040 <_mm_sub_ps@plt>
mov    %rdx,%rcx
mov    %rax,%rdi
mov    %rdx,%rsi
mov    %rax,%rdx
call   1030 <_mm_mul_ps@plt>
mov    %rbx,%rdi
mov    %r13,%rsi
mov    %eax,%edx
call   1060 <_mm_add_ps@plt>
mov    %rax,%rbx
mov    %rdx,%r13
add    $0x10,%rbp
cmp    %rbp,%r14
jne    1180 <distance_l2_m128+0x30>
jmp    11cb <distance_l2_m128+0x7b>
xor    %r13d,%r13d
xor    %ebx,%ebx
mov    %rsp,%rdi
mov    %rbx,%rsi
mov    %r13,%rdx
call   1050 <_mm_store_ps@plt>
movss  (%rsp),%xmm0
addss  0x4(%rsp),%xmm0
addss  0x8(%rsp),%xmm0
addss  0xc(%rsp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jb     11fe <distance_l2_m128+0xae>
sqrtss %xmm0,%xmm0
jmp    1203 <distance_l2_m128+0xb3>
call   1070 <sqrt@plt>
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
%struct.TYPE_12__ = type { i32, i32, i32, i32 }
define float @distance_l2_m128(%struct.TYPE_12__* nocapture %0, %struct.TYPE_12__* nocapture %1, i32 %2) {
%4 = alloca [4 x float]
%5 = bitcast [4 x float]* %4 to i8*
%6 = icmp sgt i32 %2, 3
br i1 %6, label %7, label %38
7: ; preds = %3
%8 = lshr i32 %2, 2
%9 = zext i32 %8 to i64
br label %10
10: ; preds = %7, %10
%11 = phi i64 [ 0, %7 ], [ %33, %10 ]
%12 = phi i64 [ 0, %7 ], [ %31, %10 ]
%13 = phi i64 [ 0, %7 ], [ %32, %10 ]
%14 = getelementptr inbounds %struct.TYPE_12__, %struct.TYPE_12__* %0, i64 %11
%15 = getelementptr inbounds %struct.TYPE_12__, %struct.TYPE_12__* %1, i64 %11
%16 = bitcast %struct.TYPE_12__* %14 to i64*
%17 = load i64, i64* %16
%18 = getelementptr inbounds %struct.TYPE_12__, %struct.TYPE_12__* %0, i64 %11, i32 2
%19 = bitcast i32* %18 to i64*
%20 = load i64, i64* %19
%21 = bitcast %struct.TYPE_12__* %15 to i64*
%22 = load i64, i64* %21
%23 = getelementptr inbounds %struct.TYPE_12__, %struct.TYPE_12__* %1, i64 %11, i32 2
%24 = bitcast i32* %23 to i64*
%25 = load i64, i64* %24
%26 = tail call { i64, i64 } @_mm_sub_ps(i64 %17, i64 %20, i64 %22, i64 %25)
%27 = extractvalue { i64, i64 } %26, 0
%28 = extractvalue { i64, i64 } %26, 1
%29 = tail call i32 @_mm_mul_ps(i64 %27, i64 %28, i64 %27, i64 %28)
%30 = tail call { i64, i64 } @_mm_add_ps(i64 %12, i64 %13, i32 %29)
%31 = extractvalue { i64, i64 } %30, 0
%32 = extractvalue { i64, i64 } %30, 1
%33 = add nuw nsw i64 %11, 1
%34 = icmp eq i64 %33, %9
br i1 %34, label %35, label %10
35: ; preds = %10
%36 = extractvalue { i64, i64 } %30, 0
%37 = extractvalue { i64, i64 } %30, 1
br label %38
38: ; preds = %35, %3
%39 = phi i64 [ 0, %3 ], [ %37, %35 ]
%40 = phi i64 [ 0, %3 ], [ %36, %35 ]
%41 = getelementptr inbounds [4 x float], [4 x float]* %4, i64 0, i64 0
%42 = call i32 @_mm_store_ps(float* %41, i64 %40, i64 %39)
%43 = load float, float* %41
%44 = getelementptr inbounds [4 x float], [4 x float]* %4, i64 0, i64 1
%45 = load float, float* %44
%46 = fadd float %43, %45
%47 = getelementptr inbounds [4 x float], [4 x float]* %4, i64 0, i64 2
%48 = load float, float* %47
%49 = fadd float %46, %48
%50 = getelementptr inbounds [4 x float], [4 x float]* %4, i64 0, i64 3
%51 = load float, float* %50
%52 = fadd float %49, %51
%53 = call float @sqrt(float %52)
ret float %53
}