#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 fabs(float);
float sqrt(double);
float pythag(float a, float b) {
  float absa, absb;
  absa = fabs(a);
  absb = fabs(b);
  if (absa > absb)
    return absa * sqrt(1.0 + (absb / absa) * (absb / absa));
  else
    return (absb == 0.0 ? 0.0
                        : absb * sqrt(1.0 + (absa / absb) * (absa / absb)));
}


x86 assembly (clang, -O0):
<pythag>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
movss  %xmm0,-0x8(%rbp)
movss  %xmm1,-0xc(%rbp)
movss  -0x8(%rbp),%xmm0
call   1030 <fabs@plt>
movss  %xmm0,-0x10(%rbp)
movss  -0xc(%rbp),%xmm0
call   1030 <fabs@plt>
movss  %xmm0,-0x14(%rbp)
movss  -0x10(%rbp),%xmm0
ucomiss -0x14(%rbp),%xmm0
jbe    11ac <pythag+0x8c>
movss  -0x10(%rbp),%xmm0
movss  %xmm0,-0x18(%rbp)
movss  -0x14(%rbp),%xmm0
divss  -0x10(%rbp),%xmm0
movss  -0x14(%rbp),%xmm1
divss  -0x10(%rbp),%xmm1
mulss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe73(%rip),%xmm0
addsd  %xmm1,%xmm0
call   1040 <sqrt@plt>
movaps %xmm0,%xmm1
movss  -0x18(%rbp),%xmm0
mulss  %xmm1,%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    122f <pythag+0x10f>
movss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    11d5 <pythag+0xb5>
jp     11d5 <pythag+0xb5>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
jmp    1221 <pythag+0x101>
movss  -0x14(%rbp),%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x10(%rbp),%xmm0
divss  -0x14(%rbp),%xmm0
movss  -0x10(%rbp),%xmm1
divss  -0x14(%rbp),%xmm1
mulss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xdfd(%rip),%xmm0
addsd  %xmm1,%xmm0
call   1040 <sqrt@plt>
movaps %xmm0,%xmm1
movss  -0x24(%rbp),%xmm0
mulss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x4(%rbp)
movss  -0x4(%rbp),%xmm0
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<pythag>:
sub    $0x18,%rsp
movaps %xmm0,%xmm3
movaps 0xee2(%rip),%xmm2
andps  %xmm2,%xmm3
andps  %xmm1,%xmm2
ucomiss %xmm2,%xmm3
jbe    1137 <pythag+0x27>
divss  %xmm3,%xmm2
mulss  %xmm2,%xmm2
movaps %xmm3,(%rsp)
jmp    1158 <pythag+0x48>
xorps  %xmm0,%xmm0
ucomiss %xmm0,%xmm1
jne    1146 <pythag+0x36>
jp     1146 <pythag+0x36>
add    $0x18,%rsp
ret
movaps %xmm2,%xmm0
divss  %xmm2,%xmm3
mulss  %xmm3,%xmm3
movaps %xmm3,%xmm2
movaps %xmm0,(%rsp)
xorps  %xmm0,%xmm0
cvtss2sd %xmm2,%xmm0
addsd  0xea9(%rip),%xmm0
call   1030 <sqrt@plt>
mulss  (%rsp),%xmm0
add    $0x18,%rsp
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
store float %0, float* %4
store float %1, float* %5
%8 = load float, float* %4
%9 = call float @fabs(float %8)
store float %9, float* %6
%10 = load float, float* %5
%11 = call float @fabs(float %10)
store float %11, float* %7
%12 = load float, float* %6
%13 = load float, float* %7
%14 = fcmp ogt float %12, %13
br i1 %14, label %15, label %28
15: ; preds = %2
%16 = load float, float* %6
%17 = load float, float* %7
%18 = load float, float* %6
%19 = fdiv float %17, %18
%20 = load float, float* %7
%21 = load float, float* %6
%22 = fdiv float %20, %21
%23 = fmul float %19, %22
%24 = fpext float %23 to double
%25 = fadd double 1.000000e+00, %24
%26 = call float @sqrt(double %25)
%27 = fmul float %16, %26
store float %27, float* %3
br label %50
28: ; preds = %2
%29 = load float, float* %7
%30 = fpext float %29 to double
%31 = fcmp oeq double %30, 0.000000e+00
br i1 %31, label %32, label %33
32: ; preds = %28
br label %47
33: ; preds = %28
%34 = load float, float* %7
%35 = load float, float* %6
%36 = load float, float* %7
%37 = fdiv float %35, %36
%38 = load float, float* %6
%39 = load float, float* %7
%40 = fdiv float %38, %39
%41 = fmul float %37, %40
%42 = fpext float %41 to double
%43 = fadd double 1.000000e+00, %42
%44 = call float @sqrt(double %43)
%45 = fmul float %34, %44
%46 = fpext float %45 to double
br label %47
47: ; preds = %33, %32
%48 = phi double [ 0.000000e+00, %32 ], [ %46, %33 ]
%49 = fptrunc double %48 to float
store float %49, float* %3
br label %50
50: ; preds = %47, %15
%51 = load float, float* %3
ret float %51
}

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 ogt float %3, %4
br i1 %5, label %6, label %9
6: ; preds = %2
%7 = fdiv float %4, %3
%8 = fmul float %7, %7
br label %14
9: ; preds = %2
%10 = fcmp oeq float %1, 0.000000e+00
br i1 %10, label %21, label %11
11: ; preds = %9
%12 = fdiv float %3, %4
%13 = fmul float %12, %12
br label %14
14: ; preds = %6, %11
%15 = phi float [ %13, %11 ], [ %8, %6 ]
%16 = phi float [ %4, %11 ], [ %3, %6 ]
%17 = fpext float %15 to double
%18 = fadd double %17, 1.000000e+00
%19 = tail call float @sqrt(double %18)
%20 = fmul float %16, %19
br label %21
21: ; preds = %14, %9
%22 = phi float [ 0.000000e+00, %9 ], [ %20, %14 ]
ret float %22
}