#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 fabsf(float const);
float sqrtf(float);
float geom_norm3f(const float v[3]) {
  float a[3] = {fabsf(v[0]), fabsf(v[1]), fabsf(v[2])};
  float w = a[0];
  if (a[1] > w) {
    w = a[1];
  }
  if (a[2] > w) {
    w = a[2];
  }
  if (0 == w) {
    return a[0] + a[1] + a[2];
  } else {
    a[0] /= w;
    a[1] /= w;
    a[2] /= w;
    w *= sqrtf(a[0] * a[0] + a[1] * a[1] + a[2] * a[2]);
    return w;
  }
}


x86 assembly (clang, -O0):
<geom_norm3f>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
mov    %rdi,-0x10(%rbp)
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm0
movaps 0xed5(%rip),%xmm1
pand   %xmm1,%xmm0
movss  %xmm0,-0x1c(%rbp)
mov    -0x10(%rbp),%rax
movss  0x4(%rax),%xmm0
movaps 0xebc(%rip),%xmm1
pand   %xmm1,%xmm0
movss  %xmm0,-0x18(%rbp)
mov    -0x10(%rbp),%rax
movss  0x8(%rax),%xmm0
movaps 0xea3(%rip),%xmm1
pand   %xmm1,%xmm0
movss  %xmm0,-0x14(%rbp)
movss  -0x1c(%rbp),%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x18(%rbp),%xmm0
ucomiss -0x20(%rbp),%xmm0
jbe    1189 <geom_norm3f+0x79>
movss  -0x18(%rbp),%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x14(%rbp),%xmm0
ucomiss -0x20(%rbp),%xmm0
jbe    11a2 <geom_norm3f+0x92>
movss  -0x14(%rbp),%xmm0
movss  %xmm0,-0x20(%rbp)
xorps  %xmm0,%xmm0
ucomiss -0x20(%rbp),%xmm0
jne    11ce <geom_norm3f+0xbe>
jp     11ce <geom_norm3f+0xbe>
movss  -0x1c(%rbp),%xmm0
addss  -0x18(%rbp),%xmm0
addss  -0x14(%rbp),%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    124e <geom_norm3f+0x13e>
movss  -0x20(%rbp),%xmm1
movss  -0x1c(%rbp),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,-0x1c(%rbp)
movss  -0x20(%rbp),%xmm1
movss  -0x18(%rbp),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,-0x18(%rbp)
movss  -0x20(%rbp),%xmm1
movss  -0x14(%rbp),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,-0x14(%rbp)
movss  -0x1c(%rbp),%xmm1
movss  -0x1c(%rbp),%xmm2
movss  -0x18(%rbp),%xmm0
mulss  -0x18(%rbp),%xmm0
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
movss  -0x14(%rbp),%xmm0
movss  -0x14(%rbp),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
call   1030 <sqrtf@plt>
mulss  -0x20(%rbp),%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x20(%rbp),%xmm0
movss  %xmm0,-0x4(%rbp)
movss  -0x4(%rbp),%xmm0
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<geom_norm3f>:
sub    $0x18,%rsp
movq   (%rdi),%xmm1
movdqa 0xee0(%rip),%xmm2
pand   %xmm2,%xmm1
movd   0x8(%rdi),%xmm0
pand   %xmm2,%xmm0
pshufd $0x55,%xmm1,%xmm2
movdqa %xmm2,%xmm3
maxss  %xmm1,%xmm3
movdqa %xmm0,%xmm4
maxss  %xmm3,%xmm4
xorps  %xmm3,%xmm3
ucomiss %xmm3,%xmm4
jne    1159 <geom_norm3f+0x49>
jp     1159 <geom_norm3f+0x49>
addss  %xmm2,%xmm1
addss  %xmm0,%xmm1
movaps %xmm1,%xmm4
jmp    1195 <geom_norm3f+0x85>
divss  %xmm4,%xmm1
divss  %xmm4,%xmm2
divss  %xmm4,%xmm0
mulss  %xmm2,%xmm2
mulss  %xmm1,%xmm1
addss  %xmm2,%xmm1
mulss  %xmm0,%xmm0
addss  %xmm1,%xmm0
ucomiss %xmm3,%xmm0
jb     1184 <geom_norm3f+0x74>
sqrtss %xmm0,%xmm0
jmp    1191 <geom_norm3f+0x81>
movaps %xmm4,(%rsp)
call   1030 <sqrtf@plt>
movaps (%rsp),%xmm4
mulss  %xmm0,%xmm4
movaps %xmm4,%xmm0
add    $0x18,%rsp
ret

LLVM IR (-O0):
define float @geom_norm3f(float* %0) {
%2 = alloca float
%3 = alloca float*
%4 = alloca [3 x float]
%5 = alloca float
store float* %0, float** %3
%6 = getelementptr inbounds [3 x float], [3 x float]* %4, i64 0, i64 0
%7 = load float*, float** %3
%8 = getelementptr inbounds float, float* %7, i64 0
%9 = load float, float* %8
%10 = call float @llvm.fabs.f32(float %9)
store float %10, float* %6
%11 = getelementptr inbounds float, float* %6, i64 1
%12 = load float*, float** %3
%13 = getelementptr inbounds float, float* %12, i64 1
%14 = load float, float* %13
%15 = call float @llvm.fabs.f32(float %14)
store float %15, float* %11
%16 = getelementptr inbounds float, float* %11, i64 1
%17 = load float*, float** %3
%18 = getelementptr inbounds float, float* %17, i64 2
%19 = load float, float* %18
%20 = call float @llvm.fabs.f32(float %19)
store float %20, float* %16
%21 = getelementptr inbounds [3 x float], [3 x float]* %4, i64 0, i64 0
%22 = load float, float* %21
store float %22, float* %5
%23 = getelementptr inbounds [3 x float], [3 x float]* %4, i64 0, i64 1
%24 = load float, float* %23
%25 = load float, float* %5
%26 = fcmp ogt float %24, %25
br i1 %26, label %27, label %30
27: ; preds = %1
%28 = getelementptr inbounds [3 x float], [3 x float]* %4, i64 0, i64 1
%29 = load float, float* %28
store float %29, float* %5
br label %30
30: ; preds = %27, %1
%31 = getelementptr inbounds [3 x float], [3 x float]* %4, i64 0, i64 2
%32 = load float, float* %31
%33 = load float, float* %5
%34 = fcmp ogt float %32, %33
br i1 %34, label %35, label %38
35: ; preds = %30
%36 = getelementptr inbounds [3 x float], [3 x float]* %4, i64 0, i64 2
%37 = load float, float* %36
store float %37, float* %5
br label %38
38: ; preds = %35, %30
%39 = load float, float* %5
%40 = fcmp oeq float 0.000000e+00, %39
br i1 %40, label %41, label %50
41: ; preds = %38
%42 = getelementptr inbounds [3 x float], [3 x float]* %4, i64 0, i64 0
%43 = load float, float* %42
%44 = getelementptr inbounds [3 x float], [3 x float]* %4, i64 0, i64 1
%45 = load float, float* %44
%46 = fadd float %43, %45
%47 = getelementptr inbounds [3 x float], [3 x float]* %4, i64 0, i64 2
%48 = load float, float* %47
%49 = fadd float %46, %48
store float %49, float* %2
br label %82
50: ; preds = %38
%51 = load float, float* %5
%52 = getelementptr inbounds [3 x float], [3 x float]* %4, i64 0, i64 0
%53 = load float, float* %52
%54 = fdiv float %53, %51
store float %54, float* %52
%55 = load float, float* %5
%56 = getelementptr inbounds [3 x float], [3 x float]* %4, i64 0, i64 1
%57 = load float, float* %56
%58 = fdiv float %57, %55
store float %58, float* %56
%59 = load float, float* %5
%60 = getelementptr inbounds [3 x float], [3 x float]* %4, i64 0, i64 2
%61 = load float, float* %60
%62 = fdiv float %61, %59
store float %62, float* %60
%63 = getelementptr inbounds [3 x float], [3 x float]* %4, i64 0, i64 0
%64 = load float, float* %63
%65 = getelementptr inbounds [3 x float], [3 x float]* %4, i64 0, i64 0
%66 = load float, float* %65
%67 = getelementptr inbounds [3 x float], [3 x float]* %4, i64 0, i64 1
%68 = load float, float* %67
%69 = getelementptr inbounds [3 x float], [3 x float]* %4, i64 0, i64 1
%70 = load float, float* %69
%71 = fmul float %68, %70
%72 = call float @llvm.fmuladd.f32(float %64, float %66, float %71)
%73 = getelementptr inbounds [3 x float], [3 x float]* %4, i64 0, i64 2
%74 = load float, float* %73
%75 = getelementptr inbounds [3 x float], [3 x float]* %4, i64 0, i64 2
%76 = load float, float* %75
%77 = call float @llvm.fmuladd.f32(float %74, float %76, float %72)
%78 = call float @sqrtf(float %77)
%79 = load float, float* %5
%80 = fmul float %79, %78
store float %80, float* %5
%81 = load float, float* %5
store float %81, float* %2
br label %82
82: ; preds = %50, %41
%83 = load float, float* %2
ret float %83
}

LLVM IR (-O2):
define float @geom_norm3f(float* nocapture %0) {
%2 = bitcast float* %0 to <2 x float>*
%3 = load <2 x float>, <2 x float>* %2
%4 = call <2 x float> @llvm.fabs.v2f32(<2 x float> %3)
%5 = getelementptr inbounds float, float* %0, i64 2
%6 = load float, float* %5
%7 = tail call float @llvm.fabs.f32(float %6)
%8 = extractelement <2 x float> %4, i64 0
%9 = extractelement <2 x float> %4, i64 1
%10 = fcmp ogt float %9, %8
%11 = select i1 %10, float %9, float %8
%12 = fcmp ogt float %7, %11
%13 = select i1 %12, float %7, float %11
%14 = fcmp oeq float %13, 0.000000e+00
br i1 %14, label %15, label %18
15: ; preds = %1
%16 = fadd float %8, %9
%17 = fadd float %16, %7
br label %27
18: ; preds = %1
%19 = fdiv float %8, %13
%20 = fdiv float %9, %13
%21 = fdiv float %7, %13
%22 = fmul float %20, %20
%23 = tail call float @llvm.fmuladd.f32(float %19, float %19, float %22)
%24 = tail call float @llvm.fmuladd.f32(float %21, float %21, float %23)
%25 = tail call float @sqrtf(float %24)
%26 = fmul float %13, %25
br label %27
27: ; preds = %18, %15
%28 = phi float [ %17, %15 ], [ %26, %18 ]
ret float %28
}