#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 distanceUB_l2_m128_blurp(__m128 *p1, __m128 *p2, int len, float dist) {
  float ret = 0;
  float dist2 = dist * dist;
  __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));
    if (i % 4 == 3) {
      _mm_store_ps(fsum, sum);
      float tmpsum = fsum[0] + fsum[1] + fsum[2] + fsum[3];
      if (tmpsum > dist2) {
        return tmpsum;
      }
    }
  }
  _mm_store_ps(fsum, sum);
  ret = fsum[0] + fsum[1] + fsum[2] + fsum[3];
  return sqrt(ret);
}


x86 assembly (clang, -O0):
<distanceUB_l2_m128_blurp>:
push   %rbp
mov    %rsp,%rbp
sub    $0x90,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
movss  %xmm0,-0x20(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x20(%rbp),%xmm0
mulss  -0x20(%rbp),%xmm0
movss  %xmm0,-0x28(%rbp)
lea    -0x38(%rbp),%rdi
xor    %esi,%esi
mov    $0x10,%edx
call   1030 <memset@plt>
movl   $0x0,-0x54(%rbp)
movl   $0x0,-0x54(%rbp)
mov    -0x54(%rbp),%eax
mov    %eax,-0x90(%rbp)
mov    -0x1c(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x90(%rbp),%eax
cmp    %ecx,%eax
jge    12dc <distanceUB_l2_m128_blurp+0x17c>
mov    -0x10(%rbp),%rcx
movslq -0x54(%rbp),%rax
shl    $0x4,%rax
add    %rax,%rcx
mov    -0x18(%rbp),%rax
movslq -0x54(%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,-0x78(%rbp)
mov    %rdx,-0x70(%rbp)
mov    -0x78(%rbp),%rax
mov    %rax,-0x68(%rbp)
mov    -0x70(%rbp),%rax
mov    %rax,-0x60(%rbp)
mov    -0x68(%rbp),%rdi
mov    -0x60(%rbp),%rsi
mov    -0x68(%rbp),%rdx
mov    -0x60(%rbp),%rcx
call   1040 <_mm_mul_ps@plt>
mov    %eax,%edx
mov    -0x38(%rbp),%rdi
mov    -0x30(%rbp),%rsi
call   1070 <_mm_add_ps@plt>
mov    %rax,-0x88(%rbp)
mov    %rdx,-0x80(%rbp)
mov    -0x88(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x80(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x54(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
cmp    $0x3,%edx
jne    12c9 <distanceUB_l2_m128_blurp+0x169>
lea    -0x50(%rbp),%rdi
mov    -0x38(%rbp),%rsi
mov    -0x30(%rbp),%rdx
call   1060 <_mm_store_ps@plt>
movss  -0x50(%rbp),%xmm0
addss  -0x4c(%rbp),%xmm0
addss  -0x48(%rbp),%xmm0
addss  -0x44(%rbp),%xmm0
movss  %xmm0,-0x8c(%rbp)
ff
movss  -0x8c(%rbp),%xmm0
ff
ucomiss -0x28(%rbp),%xmm0
jbe    12c4 <distanceUB_l2_m128_blurp+0x164>
movss  -0x8c(%rbp),%xmm0
ff
movss  %xmm0,-0x4(%rbp)
jmp    1315 <distanceUB_l2_m128_blurp+0x1b5>
jmp    12c9 <distanceUB_l2_m128_blurp+0x169>
jmp    12ce <distanceUB_l2_m128_blurp+0x16e>
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    11b0 <distanceUB_l2_m128_blurp+0x50>
lea    -0x50(%rbp),%rdi
mov    -0x38(%rbp),%rsi
mov    -0x30(%rbp),%rdx
call   1060 <_mm_store_ps@plt>
movss  -0x50(%rbp),%xmm0
addss  -0x4c(%rbp),%xmm0
addss  -0x48(%rbp),%xmm0
addss  -0x44(%rbp),%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x24(%rbp),%xmm0
call   1080 <sqrt@plt>
movss  %xmm0,-0x4(%rbp)
movss  -0x4(%rbp),%xmm0
add    $0x90,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<distanceUB_l2_m128_blurp>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %rsi,0x20(%rsp)
cmp    $0x4,%edx
jl     121c <distanceUB_l2_m128_blurp+0xcc>
mov    %edx,%r15d
mov    %rdi,%r12
mulss  %xmm0,%xmm0
movss  %xmm0,0x1c(%rsp)
shr    $0x2,%r15d
mov    $0x8,%ebp
xor    %r14d,%r14d
xor    %ebx,%ebx
xor    %r13d,%r13d
jmp    11a1 <distanceUB_l2_m128_blurp+0x51>
nop
add    $0x1,%r14
add    $0x10,%rbp
cmp    %r14,%r15
je     1221 <distanceUB_l2_m128_blurp+0xd1>
mov    -0x8(%r12,%rbp,1),%rdi
mov    (%r12,%rbp,1),%rsi
mov    0x20(%rsp),%rax
mov    -0x8(%rax,%rbp,1),%rdx
mov    (%rax,%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
mov    %r14d,%eax
not    %eax
test   $0x3,%al
jne    1190 <distanceUB_l2_m128_blurp+0x40>
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
ucomiss 0x1c(%rsp),%xmm0
jbe    1190 <distanceUB_l2_m128_blurp+0x40>
jmp    1259 <distanceUB_l2_m128_blurp+0x109>
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     1254 <distanceUB_l2_m128_blurp+0x104>
sqrtss %xmm0,%xmm0
jmp    1259 <distanceUB_l2_m128_blurp+0x109>
call   1070 <sqrt@plt>
add    $0x28,%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 @distanceUB_l2_m128_blurp(%struct.TYPE_12__* %0, %struct.TYPE_12__* %1, i32 %2, float %3) {
%5 = alloca float
%6 = alloca %struct.TYPE_12__*
%7 = alloca %struct.TYPE_12__*
%8 = alloca i32
%9 = alloca float
%10 = alloca float
%11 = alloca float
%12 = alloca %struct.TYPE_12__
%13 = alloca [4 x float]
%14 = alloca i32
%15 = alloca %struct.TYPE_12__
%16 = alloca %struct.TYPE_12__
%17 = alloca %struct.TYPE_12__
%18 = alloca float
store %struct.TYPE_12__* %0, %struct.TYPE_12__** %6
store %struct.TYPE_12__* %1, %struct.TYPE_12__** %7
store i32 %2, i32* %8
store float %3, float* %9
store float 0.000000e+00, float* %10
%19 = load float, float* %9
%20 = load float, float* %9
%21 = fmul float %19, %20
store float %21, float* %11
%22 = bitcast %struct.TYPE_12__* %12 to i8*
store i32 0, i32* %14
store i32 0, i32* %14
br label %23
23: ; preds = %108, %4
%24 = load i32, i32* %14
%25 = load i32, i32* %8
%26 = sdiv i32 %25, 4
%27 = icmp slt i32 %24, %26
br i1 %27, label %28, label %111
28: ; preds = %23
%29 = load %struct.TYPE_12__*, %struct.TYPE_12__** %6
%30 = load i32, i32* %14
%31 = sext i32 %30 to i64
%32 = getelementptr inbounds %struct.TYPE_12__, %struct.TYPE_12__* %29, i64 %31
%33 = load %struct.TYPE_12__*, %struct.TYPE_12__** %7
%34 = load i32, i32* %14
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds %struct.TYPE_12__, %struct.TYPE_12__* %33, i64 %35
%37 = bitcast %struct.TYPE_12__* %32 to { i64, i64 }*
%38 = getelementptr inbounds { i64, i64 }, { i64, i64 }* %37, i32 0, i32 0
%39 = load i64, i64* %38
%40 = getelementptr inbounds { i64, i64 }, { i64, i64 }* %37, i32 0, i32 1
%41 = load i64, i64* %40
%42 = bitcast %struct.TYPE_12__* %36 to { i64, i64 }*
%43 = getelementptr inbounds { i64, i64 }, { i64, i64 }* %42, i32 0, i32 0
%44 = load i64, i64* %43
%45 = getelementptr inbounds { i64, i64 }, { i64, i64 }* %42, i32 0, i32 1
%46 = load i64, i64* %45
%47 = call { i64, i64 } @_mm_sub_ps(i64 %39, i64 %41, i64 %44, i64 %46)
%48 = bitcast %struct.TYPE_12__* %16 to { i64, i64 }*
%49 = getelementptr inbounds { i64, i64 }, { i64, i64 }* %48, i32 0, i32 0
%50 = extractvalue { i64, i64 } %47, 0
store i64 %50, i64* %49
%51 = getelementptr inbounds { i64, i64 }, { i64, i64 }* %48, i32 0, i32 1
%52 = extractvalue { i64, i64 } %47, 1
store i64 %52, i64* %51
%53 = bitcast %struct.TYPE_12__* %15 to i8*
%54 = bitcast %struct.TYPE_12__* %16 to i8*
%55 = bitcast %struct.TYPE_12__* %15 to { i64, i64 }*
%56 = getelementptr inbounds { i64, i64 }, { i64, i64 }* %55, i32 0, i32 0
%57 = load i64, i64* %56
%58 = getelementptr inbounds { i64, i64 }, { i64, i64 }* %55, i32 0, i32 1
%59 = load i64, i64* %58
%60 = bitcast %struct.TYPE_12__* %15 to { i64, i64 }*
%61 = getelementptr inbounds { i64, i64 }, { i64, i64 }* %60, i32 0, i32 0
%62 = load i64, i64* %61
%63 = getelementptr inbounds { i64, i64 }, { i64, i64 }* %60, i32 0, i32 1
%64 = load i64, i64* %63
%65 = call i32 @_mm_mul_ps(i64 %57, i64 %59, i64 %62, i64 %64)
%66 = bitcast %struct.TYPE_12__* %12 to { i64, i64 }*
%67 = getelementptr inbounds { i64, i64 }, { i64, i64 }* %66, i32 0, i32 0
%68 = load i64, i64* %67
%69 = getelementptr inbounds { i64, i64 }, { i64, i64 }* %66, i32 0, i32 1
%70 = load i64, i64* %69
%71 = call { i64, i64 } @_mm_add_ps(i64 %68, i64 %70, i32 %65)
%72 = bitcast %struct.TYPE_12__* %17 to { i64, i64 }*
%73 = getelementptr inbounds { i64, i64 }, { i64, i64 }* %72, i32 0, i32 0
%74 = extractvalue { i64, i64 } %71, 0
store i64 %74, i64* %73
%75 = getelementptr inbounds { i64, i64 }, { i64, i64 }* %72, i32 0, i32 1
%76 = extractvalue { i64, i64 } %71, 1
store i64 %76, i64* %75
%77 = bitcast %struct.TYPE_12__* %12 to i8*
%78 = bitcast %struct.TYPE_12__* %17 to i8*
%79 = load i32, i32* %14
%80 = srem i32 %79, 4
%81 = icmp eq i32 %80, 3
br i1 %81, label %82, label %107
82: ; preds = %28
%83 = getelementptr inbounds [4 x float], [4 x float]* %13, i64 0, i64 0
%84 = bitcast %struct.TYPE_12__* %12 to { i64, i64 }*
%85 = getelementptr inbounds { i64, i64 }, { i64, i64 }* %84, i32 0, i32 0
%86 = load i64, i64* %85
%87 = getelementptr inbounds { i64, i64 }, { i64, i64 }* %84, i32 0, i32 1
%88 = load i64, i64* %87
%89 = call i32 @_mm_store_ps(float* %83, i64 %86, i64 %88)
%90 = getelementptr inbounds [4 x float], [4 x float]* %13, i64 0, i64 0
%91 = load float, float* %90
%92 = getelementptr inbounds [4 x float], [4 x float]* %13, i64 0, i64 1
%93 = load float, float* %92
%94 = fadd float %91, %93
%95 = getelementptr inbounds [4 x float], [4 x float]* %13, i64 0, i64 2
%96 = load float, float* %95
%97 = fadd float %94, %96
%98 = getelementptr inbounds [4 x float], [4 x float]* %13, i64 0, i64 3
%99 = load float, float* %98
%100 = fadd float %97, %99
store float %100, float* %18
%101 = load float, float* %18
%102 = load float, float* %11
%103 = fcmp ogt float %101, %102
br i1 %103, label %104, label %106
104: ; preds = %82
%105 = load float, float* %18
store float %105, float* %5
br label %132
106: ; preds = %82
br label %107
107: ; preds = %106, %28
br label %108
108: ; preds = %107
%109 = load i32, i32* %14
%110 = add nsw i32 %109, 1
store i32 %110, i32* %14
br label %23
111: ; preds = %23
%112 = getelementptr inbounds [4 x float], [4 x float]* %13, i64 0, i64 0
%113 = bitcast %struct.TYPE_12__* %12 to { i64, i64 }*
%114 = getelementptr inbounds { i64, i64 }, { i64, i64 }* %113, i32 0, i32 0
%115 = load i64, i64* %114
%116 = getelementptr inbounds { i64, i64 }, { i64, i64 }* %113, i32 0, i32 1
%117 = load i64, i64* %116
%118 = call i32 @_mm_store_ps(float* %112, i64 %115, i64 %117)
%119 = getelementptr inbounds [4 x float], [4 x float]* %13, i64 0, i64 0
%120 = load float, float* %119
%121 = getelementptr inbounds [4 x float], [4 x float]* %13, i64 0, i64 1
%122 = load float, float* %121
%123 = fadd float %120, %122
%124 = getelementptr inbounds [4 x float], [4 x float]* %13, i64 0, i64 2
%125 = load float, float* %124
%126 = fadd float %123, %125
%127 = getelementptr inbounds [4 x float], [4 x float]* %13, i64 0, i64 3
%128 = load float, float* %127
%129 = fadd float %126, %128
store float %129, float* %10
%130 = load float, float* %10
%131 = call float @sqrt(float %130)
store float %131, float* %5
br label %132
132: ; preds = %111, %104
%133 = load float, float* %5
ret float %133
}

LLVM IR (-O2):
%struct.TYPE_12__ = type { i32, i32, i32, i32 }
define float @distanceUB_l2_m128_blurp(%struct.TYPE_12__* nocapture %0, %struct.TYPE_12__* nocapture %1, i32 %2, float %3) {
%5 = alloca [4 x float]
%6 = fmul float %3, %3
%7 = bitcast [4 x float]* %5 to i8*
%8 = icmp sgt i32 %2, 3
br i1 %8, label %9, label %57
9: ; preds = %4
%10 = lshr i32 %2, 2
%11 = getelementptr inbounds [4 x float], [4 x float]* %5, i64 0, i64 0
%12 = getelementptr inbounds [4 x float], [4 x float]* %5, i64 0, i64 1
%13 = getelementptr inbounds [4 x float], [4 x float]* %5, i64 0, i64 2
%14 = getelementptr inbounds [4 x float], [4 x float]* %5, i64 0, i64 3
%15 = zext i32 %10 to i64
br label %16
16: ; preds = %9, %51
%17 = phi i64 [ 0, %9 ], [ %52, %51 ]
%18 = phi i64 [ 0, %9 ], [ %37, %51 ]
%19 = phi i64 [ 0, %9 ], [ %38, %51 ]
%20 = getelementptr inbounds %struct.TYPE_12__, %struct.TYPE_12__* %0, i64 %17
%21 = getelementptr inbounds %struct.TYPE_12__, %struct.TYPE_12__* %1, i64 %17
%22 = bitcast %struct.TYPE_12__* %20 to i64*
%23 = load i64, i64* %22
%24 = getelementptr inbounds %struct.TYPE_12__, %struct.TYPE_12__* %0, i64 %17, i32 2
%25 = bitcast i32* %24 to i64*
%26 = load i64, i64* %25
%27 = bitcast %struct.TYPE_12__* %21 to i64*
%28 = load i64, i64* %27
%29 = getelementptr inbounds %struct.TYPE_12__, %struct.TYPE_12__* %1, i64 %17, i32 2
%30 = bitcast i32* %29 to i64*
%31 = load i64, i64* %30
%32 = call { i64, i64 } @_mm_sub_ps(i64 %23, i64 %26, i64 %28, i64 %31)
%33 = extractvalue { i64, i64 } %32, 0
%34 = extractvalue { i64, i64 } %32, 1
%35 = call i32 @_mm_mul_ps(i64 %33, i64 %34, i64 %33, i64 %34)
%36 = call { i64, i64 } @_mm_add_ps(i64 %18, i64 %19, i32 %35)
%37 = extractvalue { i64, i64 } %36, 0
%38 = extractvalue { i64, i64 } %36, 1
%39 = and i64 %17, 3
%40 = icmp eq i64 %39, 3
br i1 %40, label %41, label %51
41: ; preds = %16
%42 = call i32 @_mm_store_ps(float* %11, i64 %37, i64 %38)
%43 = load float, float* %11
%44 = load float, float* %12
%45 = fadd float %43, %44
%46 = load float, float* %13
%47 = fadd float %45, %46
%48 = load float, float* %14
%49 = fadd float %47, %48
%50 = fcmp ogt float %49, %6
br i1 %50, label %73, label %51
51: ; preds = %16, %41
%52 = add nuw nsw i64 %17, 1
%53 = icmp eq i64 %52, %15
br i1 %53, label %54, label %16
54: ; preds = %51
%55 = extractvalue { i64, i64 } %36, 0
%56 = extractvalue { i64, i64 } %36, 1
br label %57
57: ; preds = %54, %4
%58 = phi i64 [ 0, %4 ], [ %56, %54 ]
%59 = phi i64 [ 0, %4 ], [ %55, %54 ]
%60 = getelementptr inbounds [4 x float], [4 x float]* %5, i64 0, i64 0
%61 = call i32 @_mm_store_ps(float* %60, i64 %59, i64 %58)
%62 = load float, float* %60
%63 = getelementptr inbounds [4 x float], [4 x float]* %5, i64 0, i64 1
%64 = load float, float* %63
%65 = fadd float %62, %64
%66 = getelementptr inbounds [4 x float], [4 x float]* %5, i64 0, i64 2
%67 = load float, float* %66
%68 = fadd float %65, %67
%69 = getelementptr inbounds [4 x float], [4 x float]* %5, i64 0, i64 3
%70 = load float, float* %69
%71 = fadd float %68, %70
%72 = call float @sqrt(float %71)
br label %73
73: ; preds = %41, %57
%74 = phi float [ %72, %57 ], [ %49, %41 ]
ret float %74
}