#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 */
float sqrt(double);
double unit_var(float *f, int npts, char *format) {
  int i, n;
  double v;
  for (v = n = i = 0; i < npts; i++)
    if (format[i] != 'x') {
      n++;
      v += f[i] * f[i];
    }
  v /= n;
  for (i = 0; i < npts; i++) {
    if (format[i] == 'x') {
      f[i] = 0.0;
    } else {
      f[i] /= sqrt(v);
    }
  }
  return v;
}


x86 assembly (clang, -O0):
<unit_var>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
movl   $0x0,-0x20(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    119e <unit_var+0x8e>
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsbl (%rax,%rcx,1),%eax
cmp    $0x78,%eax
je     118b <unit_var+0x7b>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
addsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
jmp    1190 <unit_var+0x80>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1139 <unit_var+0x29>
cvtsi2sdl -0x20(%rbp),%xmm1
movsd  -0x28(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1224 <unit_var+0x114>
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsbl (%rax,%rcx,1),%eax
cmp    $0x78,%eax
jne    11ee <unit_var+0xde>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    1211 <unit_var+0x101>
movsd  -0x28(%rbp),%xmm0
call   1030 <sqrt@plt>
movaps %xmm0,%xmm1
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    1216 <unit_var+0x106>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    11b8 <unit_var+0xa8>
movsd  -0x28(%rbp),%xmm0
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<unit_var>:
push   %r15
push   %r14
push   %r12
push   %rbx
push   %rax
test   %esi,%esi
jle    11b4 <unit_var+0xa4>
mov    %rdx,%r12
mov    %rdi,%r14
mov    %esi,%r15d
cmp    $0x1,%esi
jne    11c8 <unit_var+0xb8>
xor    %ecx,%ecx
xorpd  %xmm0,%xmm0
xor    %eax,%eax
test   $0x1,%r15b
je     115c <unit_var+0x4c>
cmpb   $0x78,(%r12,%rcx,1)
je     115c <unit_var+0x4c>
add    $0x1,%eax
movss  (%r14,%rcx,4),%xmm1
mulss  %xmm1,%xmm1
cvtss2sd %xmm1,%xmm1
addsd  %xmm1,%xmm0
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
test   %esi,%esi
jle    11bc <unit_var+0xac>
xor    %ebx,%ebx
movsd  %xmm0,(%rsp)
jmp    11a3 <unit_var+0x93>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
call   1030 <sqrt@plt>
movss  (%r14,%rbx,4),%xmm1
divss  %xmm0,%xmm1
movss  %xmm1,(%r14,%rbx,4)
movsd  (%rsp),%xmm0
add    $0x1,%rbx
cmp    %rbx,%r15
je     11bc <unit_var+0xac>
cmpb   $0x78,(%r12,%rbx,1)
jne    1180 <unit_var+0x70>
movl   $0x0,(%r14,%rbx,4)
jmp    119a <unit_var+0x8a>
movsd  0xe44(%rip),%xmm0
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret
mov    %r15d,%edx
and    $0xfffffffe,%edx
xor    %ecx,%ecx
xorpd  %xmm0,%xmm0
xor    %eax,%eax
jmp    11ed <unit_var+0xdd>
nopl   0x0(%rax,%rax,1)
add    $0x2,%rcx
cmp    %rcx,%rdx
je     113a <unit_var+0x2a>
cmpb   $0x78,(%r12,%rcx,1)
je     1209 <unit_var+0xf9>
add    $0x1,%eax
movss  (%r14,%rcx,4),%xmm1
mulss  %xmm1,%xmm1
cvtss2sd %xmm1,%xmm1
addsd  %xmm1,%xmm0
cmpb   $0x78,0x1(%r12,%rcx,1)
je     11e0 <unit_var+0xd0>
add    $0x1,%eax
movss  0x4(%r14,%rcx,4),%xmm1
mulss  %xmm1,%xmm1
cvtss2sd %xmm1,%xmm1
addsd  %xmm1,%xmm0
jmp    11e0 <unit_var+0xd0>

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

LLVM IR (-O2):
define double @unit_var(float* nocapture %0, i32 %1, i8* nocapture %2) {
%4 = icmp sgt i32 %1, 0
br i1 %4, label %5, label %86
5: ; preds = %3
%6 = zext i32 %1 to i64
%7 = and i64 %6, 1
%8 = icmp eq i32 %1, 1
br i1 %8, label %46, label %9
9: ; preds = %5
%10 = and i64 %6, 4294967294
br label %11
11: ; preds = %40, %9
%12 = phi i64 [ 0, %9 ], [ %43, %40 ]
%13 = phi double [ 0.000000e+00, %9 ], [ %42, %40 ]
%14 = phi i32 [ 0, %9 ], [ %41, %40 ]
%15 = phi i64 [ 0, %9 ], [ %44, %40 ]
%16 = getelementptr inbounds i8, i8* %2, i64 %12
%17 = load i8, i8* %16
%18 = icmp eq i8 %17, 120
br i1 %18, label %26, label %19
19: ; preds = %11
%20 = add nsw i32 %14, 1
%21 = getelementptr inbounds float, float* %0, i64 %12
%22 = load float, float* %21
%23 = fmul float %22, %22
%24 = fpext float %23 to double
%25 = fadd double %13, %24
br label %26
26: ; preds = %11, %19
%27 = phi i32 [ %20, %19 ], [ %14, %11 ]
%28 = phi double [ %25, %19 ], [ %13, %11 ]
%29 = or i64 %12, 1
%30 = getelementptr inbounds i8, i8* %2, i64 %29
%31 = load i8, i8* %30
%32 = icmp eq i8 %31, 120
br i1 %32, label %40, label %33
33: ; preds = %26
%34 = add nsw i32 %27, 1
%35 = getelementptr inbounds float, float* %0, i64 %29
%36 = load float, float* %35
%37 = fmul float %36, %36
%38 = fpext float %37 to double
%39 = fadd double %28, %38
br label %40
40: ; preds = %33, %26
%41 = phi i32 [ %34, %33 ], [ %27, %26 ]
%42 = phi double [ %39, %33 ], [ %28, %26 ]
%43 = add nuw nsw i64 %12, 2
%44 = add i64 %15, 2
%45 = icmp eq i64 %44, %10
br i1 %45, label %46, label %11
46: ; preds = %40, %5
%47 = phi i32 [ undef, %5 ], [ %41, %40 ]
%48 = phi double [ undef, %5 ], [ %42, %40 ]
%49 = phi i64 [ 0, %5 ], [ %43, %40 ]
%50 = phi double [ 0.000000e+00, %5 ], [ %42, %40 ]
%51 = phi i32 [ 0, %5 ], [ %41, %40 ]
%52 = icmp eq i64 %7, 0
br i1 %52, label %64, label %53
53: ; preds = %46
%54 = getelementptr inbounds i8, i8* %2, i64 %49
%55 = load i8, i8* %54
%56 = icmp eq i8 %55, 120
br i1 %56, label %64, label %57
57: ; preds = %53
%58 = add nsw i32 %51, 1
%59 = getelementptr inbounds float, float* %0, i64 %49
%60 = load float, float* %59
%61 = fmul float %60, %60
%62 = fpext float %61 to double
%63 = fadd double %50, %62
br label %64
64: ; preds = %53, %57, %46
%65 = phi i32 [ %47, %46 ], [ %58, %57 ], [ %51, %53 ]
%66 = phi double [ %48, %46 ], [ %63, %57 ], [ %50, %53 ]
%67 = sitofp i32 %65 to double
%68 = fdiv double %66, %67
br i1 %4, label %69, label %86
69: ; preds = %64
%70 = zext i32 %1 to i64
br label %71
71: ; preds = %69, %83
%72 = phi i64 [ 0, %69 ], [ %84, %83 ]
%73 = getelementptr inbounds i8, i8* %2, i64 %72
%74 = load i8, i8* %73
%75 = icmp eq i8 %74, 120
br i1 %75, label %76, label %78
76: ; preds = %71
%77 = getelementptr inbounds float, float* %0, i64 %72
store float 0.000000e+00, float* %77
br label %83
78: ; preds = %71
%79 = tail call float @sqrt(double %68)
%80 = getelementptr inbounds float, float* %0, i64 %72
%81 = load float, float* %80
%82 = fdiv float %81, %79
store float %82, float* %80
br label %83
83: ; preds = %76, %78
%84 = add nuw nsw i64 %72, 1
%85 = icmp eq i64 %84, %70
br i1 %85, label %86, label %71
86: ; preds = %83, %3, %64
%87 = phi double [ %68, %64 ], [ 0x7FF8000000000000, %3 ], [ %68, %83 ]
ret double %87
}