#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 */
double sqrt(float);
void sqrt_vecf(int n, float *source, float *target) {
  int i;
  double d;
  float v;
  for (i = 0; i < n; i++) {
    if ((v = source[i]) >= 0.0) {
      d = sqrt(v);
      target[i] = (float)d;
    }
  }
}


x86 assembly (clang, -O0):
<sqrt_vecf>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1191 <sqrt_vecf+0x81>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x2c(%rbp)
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     117e <sqrt_vecf+0x6e>
movss  -0x2c(%rbp),%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    1183 <sqrt_vecf+0x73>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    112a <sqrt_vecf+0x1a>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<sqrt_vecf>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %edi,%edi
jle    1150 <sqrt_vecf+0x40>
mov    %rdx,%r14
mov    %rsi,%r15
mov    %edi,%r12d
cmp    $0x1,%edi
jne    115a <sqrt_vecf+0x4a>
xor    %ebx,%ebx
test   $0x1,%r12b
je     1150 <sqrt_vecf+0x40>
movss  (%r15,%rbx,4),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jb     1150 <sqrt_vecf+0x40>
call   1030 <sqrt@plt>
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,(%r14,%rbx,4)
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret
mov    %r12d,%r13d
and    $0xfffffffe,%r13d
xor    %ebx,%ebx
xorps  %xmm1,%xmm1
jmp    1179 <sqrt_vecf+0x69>
nopl   0x0(%rax,%rax,1)
add    $0x2,%rbx
cmp    %rbx,%r13
je     112d <sqrt_vecf+0x1d>
movss  (%r15,%rbx,4),%xmm0
ucomiss %xmm1,%xmm0
jb     1196 <sqrt_vecf+0x86>
call   1030 <sqrt@plt>
xorps  %xmm1,%xmm1
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,(%r14,%rbx,4)
movss  0x4(%r15,%rbx,4),%xmm0
ucomiss %xmm1,%xmm0
jb     1170 <sqrt_vecf+0x60>
call   1030 <sqrt@plt>
xorps  %xmm1,%xmm1
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0x4(%r14,%rbx,4)
jmp    1170 <sqrt_vecf+0x60>

LLVM IR (-O0):
define void @sqrt_vecf(i32 %0, float* %1, float* %2) {
%4 = alloca i32
%5 = alloca float*
%6 = alloca float*
%7 = alloca i32
%8 = alloca double
%9 = alloca float
store i32 %0, i32* %4
store float* %1, float** %5
store float* %2, float** %6
store i32 0, i32* %7
br label %10
10: ; preds = %32, %3
%11 = load i32, i32* %7
%12 = load i32, i32* %4
%13 = icmp slt i32 %11, %12
br i1 %13, label %14, label %35
14: ; preds = %10
%15 = load float*, float** %5
%16 = load i32, i32* %7
%17 = sext i32 %16 to i64
%18 = getelementptr inbounds float, float* %15, i64 %17
%19 = load float, float* %18
store float %19, float* %9
%20 = fpext float %19 to double
%21 = fcmp oge double %20, 0.000000e+00
br i1 %21, label %22, label %31
22: ; preds = %14
%23 = load float, float* %9
%24 = call double @sqrt(float %23)
store double %24, double* %8
%25 = load double, double* %8
%26 = fptrunc double %25 to float
%27 = load float*, float** %6
%28 = load i32, i32* %7
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds float, float* %27, i64 %29
store float %26, float* %30
br label %31
31: ; preds = %22, %14
br label %32
32: ; preds = %31
%33 = load i32, i32* %7
%34 = add nsw i32 %33, 1
store i32 %34, i32* %7
br label %10
35: ; preds = %10
ret void
}

LLVM IR (-O2):
define void @sqrt_vecf(i32 %0, float* nocapture %1, float* nocapture writeonly %2) {
%4 = icmp sgt i32 %0, 0
br i1 %4, label %5, label %45
5: ; preds = %3
%6 = zext i32 %0 to i64
%7 = and i64 %6, 1
%8 = icmp eq i32 %0, 1
br i1 %8, label %34, label %9
9: ; preds = %5
%10 = and i64 %6, 4294967294
br label %11
11: ; preds = %30, %9
%12 = phi i64 [ 0, %9 ], [ %31, %30 ]
%13 = phi i64 [ 0, %9 ], [ %32, %30 ]
%14 = getelementptr inbounds float, float* %1, i64 %12
%15 = load float, float* %14
%16 = fcmp ult float %15, 0.000000e+00
br i1 %16, label %21, label %17
17: ; preds = %11
%18 = tail call double @sqrt(float %15)
%19 = fptrunc double %18 to float
%20 = getelementptr inbounds float, float* %2, i64 %12
store float %19, float* %20
br label %21
21: ; preds = %11, %17
%22 = or i64 %12, 1
%23 = getelementptr inbounds float, float* %1, i64 %22
%24 = load float, float* %23
%25 = fcmp ult float %24, 0.000000e+00
br i1 %25, label %30, label %26
26: ; preds = %21
%27 = tail call double @sqrt(float %24)
%28 = fptrunc double %27 to float
%29 = getelementptr inbounds float, float* %2, i64 %22
store float %28, float* %29
br label %30
30: ; preds = %26, %21
%31 = add nuw nsw i64 %12, 2
%32 = add i64 %13, 2
%33 = icmp eq i64 %32, %10
br i1 %33, label %34, label %11
34: ; preds = %30, %5
%35 = phi i64 [ 0, %5 ], [ %31, %30 ]
%36 = icmp eq i64 %7, 0
br i1 %36, label %45, label %37
37: ; preds = %34
%38 = getelementptr inbounds float, float* %1, i64 %35
%39 = load float, float* %38
%40 = fcmp ult float %39, 0.000000e+00
br i1 %40, label %45, label %41
41: ; preds = %37
%42 = tail call double @sqrt(float %39)
%43 = fptrunc double %42 to float
%44 = getelementptr inbounds float, float* %2, i64 %35
store float %43, float* %44
br label %45
45: ; preds = %34, %41, %37, %3
ret void
}