float pythag(float a, float b) {
  float x = fabs(a);
  float y = fabs(b);
  if (x < y)
    do {
      typeof(x) t = (x);
      (x) = (y);
      (y) = t;
    } while (0);
  if (x != 0.0)
    return x * sqrt(1.0 + (y * y) / (x * x));
  return 0.0;
}


x86 assembly (clang, -O0):
<pythag>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
movss  %xmm0,-0x8(%rbp)
movss  %xmm1,-0xc(%rbp)
movss  -0x8(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movaps 0xece(%rip),%xmm1
pand   %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x10(%rbp)
movss  -0xc(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movaps 0xeb1(%rip),%xmm1
pand   %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x14(%rbp)
movss  -0x10(%rbp),%xmm1
movss  -0x14(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jbe    1197 <pythag+0x87>
jmp    1174 <pythag+0x64>
movss  -0x10(%rbp),%xmm0
movss  %xmm0,-0x18(%rbp)
movss  -0x14(%rbp),%xmm0
movss  %xmm0,-0x10(%rbp)
movss  -0x18(%rbp),%xmm0
movss  %xmm0,-0x14(%rbp)
jmp    1197 <pythag+0x87>
movss  -0x10(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    11b8 <pythag+0xa8>
jp     11b8 <pythag+0xa8>
jmp    120d <pythag+0xfd>
movss  -0x10(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
movss  -0x14(%rbp),%xmm0
mulss  -0x14(%rbp),%xmm0
movss  -0x10(%rbp),%xmm1
mulss  -0x10(%rbp),%xmm1
divss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe26(%rip),%xmm0
addsd  %xmm1,%xmm0
call   1030 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0x20(%rbp),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    1215 <pythag+0x105>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x4(%rbp)
movss  -0x4(%rbp),%xmm0
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<pythag>:
push   %rax
movaps 0xee8(%rip),%xmm2
andps  %xmm2,%xmm0
andps  %xmm2,%xmm1
movaps %xmm1,%xmm3
maxss  %xmm0,%xmm3
xorps  %xmm2,%xmm2
ucomiss %xmm2,%xmm3
jne    112f <pythag+0x1f>
jnp    117b <pythag+0x6b>
minss  %xmm1,%xmm0
xorps  %xmm2,%xmm2
cvtss2sd %xmm3,%xmm2
mulss  %xmm0,%xmm0
mulss  %xmm3,%xmm3
divss  %xmm3,%xmm0
cvtss2sd %xmm0,%xmm0
addsd  0xebe(%rip),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     1161 <pythag+0x51>
sqrtsd %xmm0,%xmm0
jmp    1170 <pythag+0x60>
movsd  %xmm2,(%rsp)
call   1030 <sqrt@plt>
movsd  (%rsp),%xmm2
mulsd  %xmm2,%xmm0
xorps  %xmm2,%xmm2
cvtsd2ss %xmm0,%xmm2
movaps %xmm2,%xmm0
pop    %rax
ret

LLVM IR (-O0):
define float @pythag(float %0, float %1) {
%3 = alloca float
%4 = alloca float
%5 = alloca float
%6 = alloca float
%7 = alloca float
%8 = alloca float
store float %0, float* %4
store float %1, float* %5
%9 = load float, float* %4
%10 = fpext float %9 to double
%11 = call double @llvm.fabs.f64(double %10)
%12 = fptrunc double %11 to float
store float %12, float* %6
%13 = load float, float* %5
%14 = fpext float %13 to double
%15 = call double @llvm.fabs.f64(double %14)
%16 = fptrunc double %15 to float
store float %16, float* %7
%17 = load float, float* %6
%18 = load float, float* %7
%19 = fcmp olt float %17, %18
br i1 %19, label %20, label %26
20: ; preds = %2
br label %21
21: ; preds = %20
%22 = load float, float* %6
store float %22, float* %8
%23 = load float, float* %7
store float %23, float* %6
%24 = load float, float* %8
store float %24, float* %7
br label %25
25: ; preds = %21
br label %26
26: ; preds = %25, %2
%27 = load float, float* %6
%28 = fpext float %27 to double
%29 = fcmp une double %28, 0.000000e+00
br i1 %29, label %30, label %45
30: ; preds = %26
%31 = load float, float* %6
%32 = fpext float %31 to double
%33 = load float, float* %7
%34 = load float, float* %7
%35 = fmul float %33, %34
%36 = load float, float* %6
%37 = load float, float* %6
%38 = fmul float %36, %37
%39 = fdiv float %35, %38
%40 = fpext float %39 to double
%41 = fadd double 1.000000e+00, %40
%42 = call double @sqrt(double %41)
%43 = fmul double %32, %42
%44 = fptrunc double %43 to float
store float %44, float* %3
br label %46
45: ; preds = %26
store float 0.000000e+00, float* %3
br label %46
46: ; preds = %45, %30
%47 = load float, float* %3
ret float %47
}

LLVM IR (-O2):
define float @pythag(float %0, float %1) {
%3 = tail call float @llvm.fabs.f32(float %0)
%4 = tail call float @llvm.fabs.f32(float %1)
%5 = fcmp olt float %3, %4
%6 = select i1 %5, float %4, float %3
%7 = fcmp une float %6, 0.000000e+00
br i1 %7, label %8, label %19
8: ; preds = %2
%9 = select i1 %5, float %3, float %4
%10 = fpext float %6 to double
%11 = fmul float %9, %9
%12 = fmul float %6, %6
%13 = fdiv float %11, %12
%14 = fpext float %13 to double
%15 = fadd double %14, 1.000000e+00
%16 = tail call double @sqrt(double %15)
%17 = fmul double %16, %10
%18 = fptrunc double %17 to float
br label %19
19: ; preds = %2, %8
%20 = phi float [ %18, %8 ], [ 0.000000e+00, %2 ]
ret float %20
}