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


x86 assembly (clang, -O0):
<bessi1_salahat_exponScaled>:
push   %rbp
mov    %rsp,%rbp
sub    $0xa0,%rsp
movsd  %xmm0,-0x18(%rbp)
movsd  0xed8(%rip),%xmm0
ucomisd -0x18(%rbp),%xmm0
jbe    11f8 <bessi1_salahat_exponScaled+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 <bessi1_salahat_exponScaled+0x36a>
movsd  0xe08(%rip),%xmm0
ucomisd -0x18(%rbp),%xmm0
jbe    12d0 <bessi1_salahat_exponScaled+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 <bessi1_salahat_exponScaled+0x36a>
movsd  0xd38(%rip),%xmm0
ucomisd -0x18(%rbp),%xmm0
jbe    13a8 <bessi1_salahat_exponScaled+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 <bessi1_salahat_exponScaled+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):
<bessi1_salahat_exponScaled>:
sub    $0x18,%rsp
movapd %xmm0,%xmm1
movsd  0xee0(%rip),%xmm0
ucomisd %xmm1,%xmm0
movsd  %xmm1,0x8(%rsp)
jbe    11bc <bessi1_salahat_exponScaled+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 <bessi1_salahat_exponScaled+0x26f>
movsd  0xe44(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    125a <bessi1_salahat_exponScaled+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 <bessi1_salahat_exponScaled+0x26f>
movsd  0xdae(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    12f8 <bessi1_salahat_exponScaled+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 <bessi1_salahat_exponScaled+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 @bessi1_salahat_exponScaled(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, 0xBFCF8A0902DE00D0
%9 = call double @expl(double %8)
%10 = load double, double* %3
%11 = fmul double %10, 0xBF9AB9F559B3D080
%12 = call double @expl(double %11)
%13 = fmul double 0xBF6F7917181F8687, %12
%14 = call double @llvm.fmuladd.f64(double 0xBFA53835FC66427E, double %9, double %13)
%15 = load double, double* %3
%16 = fmul double %15, 0xBFFB70A3D70A3D70
%17 = call double @expl(double %16)
%18 = call double @llvm.fmuladd.f64(double 0xBFE86BEDFA43FE5C, double %17, double %14)
%19 = load double, double* %3
%20 = fmul double %19, -7.657000e-01
%21 = call double @expl(double %20)
%22 = call double @llvm.fmuladd.f64(double 0xBFC7588C0836854D, 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, -5.290000e-01
%30 = call double @expl(double %29)
%31 = load double, double* %3
%32 = fmul double %31, -1.644000e+02
%33 = call double @expl(double %32)
%34 = fmul double 0xC054346B26BF876A, %33
%35 = call double @llvm.fmuladd.f64(double 0xBFC20F0CE33D519A, double %30, double %34)
%36 = load double, double* %3
%37 = fmul double %36, 0xBF8E4F765FD8ADC0
%38 = call double @expl(double %37)
%39 = call double @llvm.fmuladd.f64(double 0xBF5AEA6711871113, double %38, double %35)
%40 = load double, double* %3
%41 = fmul double %40, 0xBFC28240B780346C
%42 = call double @expl(double %41)
%43 = call double @llvm.fmuladd.f64(double 0xBF934EF2008DEF6A, 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, 0xBF93C36113404EA0
%51 = call double @expl(double %50)
%52 = load double, double* %3
%53 = fmul double %52, 0xBFC0FF9724745390
%54 = call double @expl(double %53)
%55 = fmul double 0xBF8CB17C0E6E57CE, %54
%56 = call double @llvm.fmuladd.f64(double 0xBF61B940530159C3, double %51, double %55)
%57 = load double, double* %3
%58 = fmul double %57, 0x3FB45A1CAC083120
%59 = call double @expl(double %58)
%60 = call double @llvm.fmuladd.f64(double 0xBEEE029F6C9757EC, double %59, double %56)
%61 = load double, double* %3
%62 = fmul double %61, 0x3FA3B645A1CAC080
%63 = call double @expl(double %62)
%64 = call double @llvm.fmuladd.f64(double 0x3F20736AB0CDE8A9, 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, 0x3FC26E978D4FDF38
%69 = call double @expl(double %68)
%70 = load double, double* %3
%71 = fmul double %70, 0xBF747AE147AE1480
%72 = call double @expl(double %71)
%73 = fmul double 0xBF361A1DB877AB37, %72
%74 = call double @llvm.fmuladd.f64(double 0x3DF7D93136F5714B, double %69, double %73)
%75 = load double, double* %3
%76 = fmul double %75, -4.314000e-01
%77 = call double @expl(double %76)
%78 = call double @llvm.fmuladd.f64(double 0xBF982B184B6A257D, double %77, double %74)
%79 = load double, double* %3
%80 = fmul double %79, 0xBFABA5E353F7CEE0
%81 = call double @expl(double %80)
%82 = call double @llvm.fmuladd.f64(double 0xBF7167AC0BB7185F, 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 @bessi1_salahat_exponScaled(double %0) {
%2 = fcmp olt double %0, 1.150000e+01
br i1 %2, label %3, label %16
3: ; preds = %1
%4 = fmul double %0, 0xBFCF8A0902DE00D0
%5 = tail call double @expl(double %4)
%6 = fmul double %0, 0xBF9AB9F559B3D080
%7 = tail call double @expl(double %6)
%8 = fmul double %7, 0xBF6F7917181F8687
%9 = tail call double @llvm.fmuladd.f64(double %5, double 0xBFA53835FC66427E, double %8)
%10 = fmul double %0, 0xBFFB70A3D70A3D70
%11 = tail call double @expl(double %10)
%12 = tail call double @llvm.fmuladd.f64(double %11, double 0xBFE86BEDFA43FE5C, double %9)
%13 = fmul double %0, -7.657000e-01
%14 = tail call double @expl(double %13)
%15 = tail call double @llvm.fmuladd.f64(double %14, double 0xBFC7588C0836854D, 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, -5.290000e-01
%20 = tail call double @expl(double %19)
%21 = fmul double %0, -1.644000e+02
%22 = tail call double @expl(double %21)
%23 = fmul double %22, 0xC054346B26BF876A
%24 = tail call double @llvm.fmuladd.f64(double %20, double 0xBFC20F0CE33D519A, double %23)
%25 = fmul double %0, 0xBF8E4F765FD8ADC0
%26 = tail call double @expl(double %25)
%27 = tail call double @llvm.fmuladd.f64(double %26, double 0xBF5AEA6711871113, double %24)
%28 = fmul double %0, 0xBFC28240B780346C
%29 = tail call double @expl(double %28)
%30 = tail call double @llvm.fmuladd.f64(double %29, double 0xBF934EF2008DEF6A, 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, 0xBF93C36113404EA0
%35 = tail call double @expl(double %34)
%36 = fmul double %0, 0xBFC0FF9724745390
%37 = tail call double @expl(double %36)
%38 = fmul double %37, 0xBF8CB17C0E6E57CE
%39 = tail call double @llvm.fmuladd.f64(double %35, double 0xBF61B940530159C3, double %38)
%40 = fmul double %0, 0x3FB45A1CAC083120
%41 = tail call double @expl(double %40)
%42 = tail call double @llvm.fmuladd.f64(double %41, double 0xBEEE029F6C9757EC, double %39)
%43 = fmul double %0, 0x3FA3B645A1CAC080
%44 = tail call double @expl(double %43)
%45 = tail call double @llvm.fmuladd.f64(double %44, double 0x3F20736AB0CDE8A9, double %42)
br label %59
46: ; preds = %31
%47 = fmul double %0, 0x3FC26E978D4FDF38
%48 = tail call double @expl(double %47)
%49 = fmul double %0, 0xBF747AE147AE1480
%50 = tail call double @expl(double %49)
%51 = fmul double %50, 0xBF361A1DB877AB37
%52 = tail call double @llvm.fmuladd.f64(double %48, double 0x3DF7D93136F5714B, double %51)
%53 = fmul double %0, -4.314000e-01
%54 = tail call double @expl(double %53)
%55 = tail call double @llvm.fmuladd.f64(double %54, double 0xBF982B184B6A257D, double %52)
%56 = fmul double %0, 0xBFABA5E353F7CEE0
%57 = tail call double @expl(double %56)
%58 = tail call double @llvm.fmuladd.f64(double %57, double 0xBF7167AC0BB7185F, 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
}