#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 expl(double);
long double bessi0_salahat(double x) {
  if (x < 11.5)
    return 0.1682 * expl(x * 0.7536) + 0.1472 * expl(x * 0.9739) +
           0.4450 * expl(x * -0.715) + 0.2382 * expl(x * 0.2343);
  else if (x < 20)
    return 0.2667 * expl(x * 0.4710) + 0.4916 * expl(x * -163.4) +
           0.1110 * expl(x * 0.9852) + 0.1304 * expl(x * 0.8554);
  else if (x < 37.25)
    return 0.1121 * expl(x * 0.9807) + 0.1055 * expl(x * 0.8672) +
           -0.00018 * expl(x * 1.0795) + 0.00326 * expl(x * 1.0385);
  else
    return 2.41E-9 * expl(x * 1.144) + 0.06745 * expl(x * 0.995) +
           0.05471 * expl(x * 0.5686) + 0.07869 * expl(x * 0.946);
}


x86 assembly (clang, -O0):
<bessi0_salahat>:
push   %rbp
mov    %rsp,%rbp
sub    $0xa0,%rsp
movsd  %xmm0,-0x18(%rbp)
movsd  0xed8(%rip),%xmm0
ucomisd -0x18(%rbp),%xmm0
jbe    11f8 <bessi0_salahat+0xe8>
movsd  -0x18(%rbp),%xmm0
movsd  0xf98(%rip),%xmm1
mulsd  %xmm1,%xmm0
call   1030 <expl@plt>
movsd  %xmm0,-0x50(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  0xf85(%rip),%xmm1
mulsd  %xmm1,%xmm0
call   1030 <expl@plt>
movaps %xmm0,%xmm1
movsd  -0x50(%rbp),%xmm0
movsd  0xf74(%rip),%xmm2
mulsd  %xmm2,%xmm1
movsd  0xf70(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  0xf5e(%rip),%xmm1
mulsd  %xmm1,%xmm0
call   1030 <expl@plt>
movsd  -0x48(%rbp),%xmm1
movsd  0xf50(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  0xf3e(%rip),%xmm1
mulsd  %xmm1,%xmm0
call   1030 <expl@plt>
movsd  -0x40(%rbp),%xmm1
movsd  0xf30(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
fldl   -0x38(%rbp)
fstpt  -0x10(%rbp)
jmp    147a <bessi0_salahat+0x36a>
movsd  0xe08(%rip),%xmm0
ucomisd -0x18(%rbp),%xmm0
jbe    12d0 <bessi0_salahat+0x1c0>
movsd  -0x18(%rbp),%xmm0
movsd  0xe80(%rip),%xmm1
mulsd  %xmm1,%xmm0
call   1030 <expl@plt>
movsd  %xmm0,-0x68(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  0xe6d(%rip),%xmm1
mulsd  %xmm1,%xmm0
call   1030 <expl@plt>
movaps %xmm0,%xmm1
movsd  -0x68(%rbp),%xmm0
movsd  0xe5c(%rip),%xmm2
mulsd  %xmm2,%xmm1
movsd  0xe58(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  0xe46(%rip),%xmm1
mulsd  %xmm1,%xmm0
call   1030 <expl@plt>
movsd  -0x60(%rbp),%xmm1
movsd  0xe38(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  0xe26(%rip),%xmm1
mulsd  %xmm1,%xmm0
call   1030 <expl@plt>
movsd  -0x58(%rbp),%xmm1
movsd  0xe18(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
fldl   -0x30(%rbp)
fstpt  -0x10(%rbp)
jmp    147a <bessi0_salahat+0x36a>
movsd  0xd38(%rip),%xmm0
ucomisd -0x18(%rbp),%xmm0
jbe    13a8 <bessi0_salahat+0x298>
movsd  -0x18(%rbp),%xmm0
movsd  0xd68(%rip),%xmm1
mulsd  %xmm1,%xmm0
call   1030 <expl@plt>
movsd  %xmm0,-0x80(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  0xd55(%rip),%xmm1
mulsd  %xmm1,%xmm0
call   1030 <expl@plt>
movaps %xmm0,%xmm1
movsd  -0x80(%rbp),%xmm0
movsd  0xd44(%rip),%xmm2
mulsd  %xmm2,%xmm1
movsd  0xd40(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  0xd2e(%rip),%xmm1
mulsd  %xmm1,%xmm0
call   1030 <expl@plt>
movsd  -0x78(%rbp),%xmm1
movsd  0xd20(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  0xd0e(%rip),%xmm1
mulsd  %xmm1,%xmm0
call   1030 <expl@plt>
movsd  -0x70(%rbp),%xmm1
movsd  0xd00(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
fldl   -0x28(%rbp)
fstpt  -0x10(%rbp)
jmp    147a <bessi0_salahat+0x36a>
movsd  -0x18(%rbp),%xmm0
movsd  0xc63(%rip),%xmm1
mulsd  %xmm1,%xmm0
call   1030 <expl@plt>
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x18(%rbp),%xmm0
movsd  0xc4d(%rip),%xmm1
mulsd  %xmm1,%xmm0
call   1030 <expl@plt>
movaps %xmm0,%xmm1
movsd  -0x98(%rbp),%xmm0
ff
movsd  0xc39(%rip),%xmm2
mulsd  %xmm2,%xmm1
movsd  0xc35(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x18(%rbp),%xmm0
movsd  0xc20(%rip),%xmm1
mulsd  %xmm1,%xmm0
call   1030 <expl@plt>
movsd  -0x90(%rbp),%xmm1
ff
movsd  0xc0f(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x18(%rbp),%xmm0
movsd  0xbfa(%rip),%xmm1
mulsd  %xmm1,%xmm0
call   1030 <expl@plt>
movsd  -0x88(%rbp),%xmm1
ff
movsd  0xbe9(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
fldl   -0x20(%rbp)
fstpt  -0x10(%rbp)
fldt   -0x10(%rbp)
add    $0xa0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<bessi0_salahat>:
sub    $0x18,%rsp
movapd %xmm0,%xmm1
movsd  0xee0(%rip),%xmm0
ucomisd %xmm1,%xmm0
movsd  %xmm1,0x8(%rsp)
jbe    11bc <bessi0_salahat+0xac>
movsd  0xfa0(%rip),%xmm0
mulsd  %xmm1,%xmm0
call   1030 <expl@plt>
movsd  %xmm0,(%rsp)
movsd  0xf92(%rip),%xmm0
mulsd  0x8(%rsp),%xmm0
call   1030 <expl@plt>
mulsd  0xf87(%rip),%xmm0
movsd  (%rsp),%xmm1
mulsd  0xf82(%rip),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rsp)
movsd  0xf79(%rip),%xmm0
mulsd  0x8(%rsp),%xmm0
call   1030 <expl@plt>
mulsd  0xf6e(%rip),%xmm0
addsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
movsd  0x8(%rsp),%xmm0
mulsd  0xf5e(%rip),%xmm0
call   1030 <expl@plt>
mulsd  0xf59(%rip),%xmm0
jmp    137f <bessi0_salahat+0x26f>
movsd  0xe44(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    125a <bessi0_salahat+0x14a>
movsd  0xec2(%rip),%xmm0
mulsd  %xmm1,%xmm0
call   1030 <expl@plt>
movsd  %xmm0,(%rsp)
movsd  0xeb4(%rip),%xmm0
mulsd  0x8(%rsp),%xmm0
call   1030 <expl@plt>
mulsd  0xea9(%rip),%xmm0
movsd  (%rsp),%xmm1
mulsd  0xea4(%rip),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rsp)
movsd  0xe9b(%rip),%xmm0
mulsd  0x8(%rsp),%xmm0
call   1030 <expl@plt>
mulsd  0xe90(%rip),%xmm0
addsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
movsd  0x8(%rsp),%xmm0
mulsd  0xe80(%rip),%xmm0
call   1030 <expl@plt>
mulsd  0xe7b(%rip),%xmm0
jmp    137f <bessi0_salahat+0x26f>
movsd  0xdae(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    12f8 <bessi0_salahat+0x1e8>
movsd  0xde4(%rip),%xmm0
mulsd  %xmm1,%xmm0
call   1030 <expl@plt>
movsd  %xmm0,(%rsp)
movsd  0xdd6(%rip),%xmm0
mulsd  0x8(%rsp),%xmm0
call   1030 <expl@plt>
mulsd  0xdcb(%rip),%xmm0
movsd  (%rsp),%xmm1
mulsd  0xdc6(%rip),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rsp)
movsd  0xdbd(%rip),%xmm0
mulsd  0x8(%rsp),%xmm0
call   1030 <expl@plt>
mulsd  0xdb2(%rip),%xmm0
addsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
movsd  0x8(%rsp),%xmm0
mulsd  0xda2(%rip),%xmm0
call   1030 <expl@plt>
mulsd  0xd9d(%rip),%xmm0
jmp    137f <bessi0_salahat+0x26f>
movsd  0xd18(%rip),%xmm0
mulsd  %xmm1,%xmm0
call   1030 <expl@plt>
movsd  %xmm0,(%rsp)
movsd  0xd0a(%rip),%xmm0
mulsd  0x8(%rsp),%xmm0
call   1030 <expl@plt>
mulsd  0xcff(%rip),%xmm0
movsd  (%rsp),%xmm1
mulsd  0xcfa(%rip),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rsp)
movsd  0xcf1(%rip),%xmm0
mulsd  0x8(%rsp),%xmm0
call   1030 <expl@plt>
mulsd  0xce6(%rip),%xmm0
addsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
movsd  0x8(%rsp),%xmm0
mulsd  0xcd6(%rip),%xmm0
call   1030 <expl@plt>
mulsd  0xcd1(%rip),%xmm0
addsd  (%rsp),%xmm0
movsd  %xmm0,0x10(%rsp)
fldl   0x10(%rsp)
add    $0x18,%rsp
ret

LLVM IR (-O0):
define x86_fp80 @bessi0_salahat(double %0) {
%2 = alloca x86_fp80
%3 = alloca double
store double %0, double* %3
%4 = load double, double* %3
%5 = fcmp olt double %4, 1.150000e+01
br i1 %5, label %6, label %24
6: ; preds = %1
%7 = load double, double* %3
%8 = fmul double %7, 7.536000e-01
%9 = call double @expl(double %8)
%10 = load double, double* %3
%11 = fmul double %10, 0x3FEF2A305532617C
%12 = call double @expl(double %11)
%13 = fmul double 1.472000e-01, %12
%14 = call double @llvm.fmuladd.f64(double 1.682000e-01, double %9, double %13)
%15 = load double, double* %3
%16 = fmul double %15, -7.150000e-01
%17 = call double @expl(double %16)
%18 = call double @llvm.fmuladd.f64(double 4.450000e-01, double %17, double %14)
%19 = load double, double* %3
%20 = fmul double %19, 2.343000e-01
%21 = call double @expl(double %20)
%22 = call double @llvm.fmuladd.f64(double 2.382000e-01, double %21, double %18)
%23 = fpext double %22 to x86_fp80
store x86_fp80 %23, x86_fp80* %2
br label %84
24: ; preds = %1
%25 = load double, double* %3
%26 = fcmp olt double %25, 2.000000e+01
br i1 %26, label %27, label %45
27: ; preds = %24
%28 = load double, double* %3
%29 = fmul double %28, 4.710000e-01
%30 = call double @expl(double %29)
%31 = load double, double* %3
%32 = fmul double %31, -1.634000e+02
%33 = call double @expl(double %32)
%34 = fmul double 4.916000e-01, %33
%35 = call double @llvm.fmuladd.f64(double 2.667000e-01, double %30, double %34)
%36 = load double, double* %3
%37 = fmul double %36, 0x3FEF86C226809D49
%38 = call double @expl(double %37)
%39 = call double @llvm.fmuladd.f64(double 1.110000e-01, double %38, double %35)
%40 = load double, double* %3
%41 = fmul double %40, 8.554000e-01
%42 = call double @expl(double %41)
%43 = call double @llvm.fmuladd.f64(double 1.304000e-01, double %42, double %39)
%44 = fpext double %43 to x86_fp80
store x86_fp80 %44, x86_fp80* %2
br label %84
45: ; preds = %24
%46 = load double, double* %3
%47 = fcmp olt double %46, 3.725000e+01
br i1 %47, label %48, label %66
48: ; preds = %45
%49 = load double, double* %3
%50 = fmul double %49, 9.807000e-01
%51 = call double @expl(double %50)
%52 = load double, double* %3
%53 = fmul double %52, 0x3FEBC01A36E2EB1C
%54 = call double @expl(double %53)
%55 = fmul double 1.055000e-01, %54
%56 = call double @llvm.fmuladd.f64(double 1.121000e-01, double %51, double %55)
%57 = load double, double* %3
%58 = fmul double %57, 1.079500e+00
%59 = call double @expl(double %58)
%60 = call double @llvm.fmuladd.f64(double -1.800000e-04, double %59, double %56)
%61 = load double, double* %3
%62 = fmul double %61, 1.038500e+00
%63 = call double @expl(double %62)
%64 = call double @llvm.fmuladd.f64(double 3.260000e-03, double %63, double %60)
%65 = fpext double %64 to x86_fp80
store x86_fp80 %65, x86_fp80* %2
br label %84
66: ; preds = %45
%67 = load double, double* %3
%68 = fmul double %67, 1.144000e+00
%69 = call double @expl(double %68)
%70 = load double, double* %3
%71 = fmul double %70, 0x3FEFD70A3D70A3D7
%72 = call double @expl(double %71)
%73 = fmul double 0x3FB14467381D7DBF, %72
%74 = call double @llvm.fmuladd.f64(double 2.410000e-09, double %69, double %73)
%75 = load double, double* %3
%76 = fmul double %75, 5.686000e-01
%77 = call double @expl(double %76)
%78 = call double @llvm.fmuladd.f64(double 5.471000e-02, double %77, double %74)
%79 = load double, double* %3
%80 = fmul double %79, 0x3FEE45A1CAC08312
%81 = call double @expl(double %80)
%82 = call double @llvm.fmuladd.f64(double 0x3FB425072085B185, double %81, double %78)
%83 = fpext double %82 to x86_fp80
store x86_fp80 %83, x86_fp80* %2
br label %84
84: ; preds = %66, %48, %27, %6
%85 = load x86_fp80, x86_fp80* %2
ret x86_fp80 %85
}

LLVM IR (-O2):
define x86_fp80 @bessi0_salahat(double %0) {
%2 = fcmp olt double %0, 1.150000e+01
br i1 %2, label %3, label %16
3: ; preds = %1
%4 = fmul double %0, 7.536000e-01
%5 = tail call double @expl(double %4)
%6 = fmul double %0, 0x3FEF2A305532617C
%7 = tail call double @expl(double %6)
%8 = fmul double %7, 1.472000e-01
%9 = tail call double @llvm.fmuladd.f64(double %5, double 1.682000e-01, double %8)
%10 = fmul double %0, -7.150000e-01
%11 = tail call double @expl(double %10)
%12 = tail call double @llvm.fmuladd.f64(double %11, double 4.450000e-01, double %9)
%13 = fmul double %0, 2.343000e-01
%14 = tail call double @expl(double %13)
%15 = tail call double @llvm.fmuladd.f64(double %14, double 2.382000e-01, double %12)
br label %59
16: ; preds = %1
%17 = fcmp olt double %0, 2.000000e+01
br i1 %17, label %18, label %31
18: ; preds = %16
%19 = fmul double %0, 4.710000e-01
%20 = tail call double @expl(double %19)
%21 = fmul double %0, -1.634000e+02
%22 = tail call double @expl(double %21)
%23 = fmul double %22, 4.916000e-01
%24 = tail call double @llvm.fmuladd.f64(double %20, double 2.667000e-01, double %23)
%25 = fmul double %0, 0x3FEF86C226809D49
%26 = tail call double @expl(double %25)
%27 = tail call double @llvm.fmuladd.f64(double %26, double 1.110000e-01, double %24)
%28 = fmul double %0, 8.554000e-01
%29 = tail call double @expl(double %28)
%30 = tail call double @llvm.fmuladd.f64(double %29, double 1.304000e-01, double %27)
br label %59
31: ; preds = %16
%32 = fcmp olt double %0, 3.725000e+01
br i1 %32, label %33, label %46
33: ; preds = %31
%34 = fmul double %0, 9.807000e-01
%35 = tail call double @expl(double %34)
%36 = fmul double %0, 0x3FEBC01A36E2EB1C
%37 = tail call double @expl(double %36)
%38 = fmul double %37, 1.055000e-01
%39 = tail call double @llvm.fmuladd.f64(double %35, double 1.121000e-01, double %38)
%40 = fmul double %0, 1.079500e+00
%41 = tail call double @expl(double %40)
%42 = tail call double @llvm.fmuladd.f64(double %41, double -1.800000e-04, double %39)
%43 = fmul double %0, 1.038500e+00
%44 = tail call double @expl(double %43)
%45 = tail call double @llvm.fmuladd.f64(double %44, double 3.260000e-03, double %42)
br label %59
46: ; preds = %31
%47 = fmul double %0, 1.144000e+00
%48 = tail call double @expl(double %47)
%49 = fmul double %0, 0x3FEFD70A3D70A3D7
%50 = tail call double @expl(double %49)
%51 = fmul double %50, 0x3FB14467381D7DBF
%52 = tail call double @llvm.fmuladd.f64(double %48, double 2.410000e-09, double %51)
%53 = fmul double %0, 5.686000e-01
%54 = tail call double @expl(double %53)
%55 = tail call double @llvm.fmuladd.f64(double %54, double 5.471000e-02, double %52)
%56 = fmul double %0, 0x3FEE45A1CAC08312
%57 = tail call double @expl(double %56)
%58 = tail call double @llvm.fmuladd.f64(double %57, double 0x3FB425072085B185, double %55)
br label %59
59: ; preds = %46, %33, %18, %3
%60 = phi double [ %15, %3 ], [ %30, %18 ], [ %45, %33 ], [ %58, %46 ]
%61 = fpext double %60 to x86_fp80
ret x86_fp80 %61
}