#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 alphaB_PPB(double);
double exp(double);
double pow(int, double);
double sqrt(double);
double rec_HRecFast_dxedlna(double xe, double nH, double H, double TM,
                            double TR, double energy_rate) {
  double RLya, alphaB, four_betaB, C;
  double chi_ion_H;
  RLya = 4.662899067555897e15 * H / nH / (1. - xe);
  alphaB = 1.14 * alphaB_PPB(TM);
  four_betaB = 3.016103031869581e21 * TR * sqrt(TR) *
               exp(-0.25 * 13.598286071938324 / TR) * alphaB;
  C = (3. * RLya + 8.2206) / (3. * RLya + 8.2206 + four_betaB);
  if (xe < 1.)
    chi_ion_H = 0.369202 * pow(1. - pow(xe, 0.463929), 1.70237);
  else
    chi_ion_H = 0.;
  return (-nH * xe * xe * alphaB +
          four_betaB * (1. - xe) * exp(-10.198714553953742 / TR)) *
             C / H +
         chi_ion_H / nH * energy_rate *
             (1. / 13.598286071938324 + (1. - C) / 10.198714553953742) / H;
}


x86 assembly (clang, -O0):
<rec_HRecFast_dxedlna>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
movsd  %xmm2,-0x18(%rbp)
movsd  %xmm3,-0x20(%rbp)
movsd  %xmm4,-0x28(%rbp)
movsd  %xmm5,-0x30(%rbp)
movsd  0xebf(%rip),%xmm0
mulsd  -0x18(%rbp),%xmm0
divsd  -0x10(%rbp),%xmm0
movsd  0xe7d(%rip),%xmm1
subsd  -0x8(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x20(%rbp),%xmm0
call   1060 <alphaB_PPB@plt>
movaps %xmm0,%xmm1
movsd  0xe82(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  0xe69(%rip),%xmm0
mulsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x28(%rbp),%xmm0
call   1050 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0x68(%rbp),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  0xe34(%rip),%xmm0
divsd  -0x28(%rbp),%xmm0
call   1030 <exp@plt>
movaps %xmm0,%xmm1
movsd  -0x60(%rbp),%xmm0
mulsd  %xmm1,%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  0xdf7(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xdf3(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  -0x38(%rbp),%xmm1
movsd  0xdda(%rip),%xmm2
mulsd  %xmm2,%xmm1
movsd  0xdd6(%rip),%xmm2
addsd  %xmm2,%xmm1
addsd  -0x48(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  0xdac(%rip),%xmm0
ucomisd -0x8(%rbp),%xmm0
jbe    12aa <rec_HRecFast_dxedlna+0x16a>
cvttsd2si -0x8(%rbp),%edi
movsd  0xddc(%rip),%xmm0
call   1040 <pow@plt>
movaps %xmm0,%xmm1
movsd  0xd84(%rip),%xmm0
subsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
movsd  0xdb4(%rip),%xmm0
call   1040 <pow@plt>
movaps %xmm0,%xmm1
movsd  0xd9c(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
jmp    12b2 <rec_HRecFast_dxedlna+0x172>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x10(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  -0x8(%rbp),%xmm0
mulsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  0xd0c(%rip),%xmm1
subsd  -0x8(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  0xd56(%rip),%xmm0
divsd  -0x28(%rbp),%xmm0
call   1030 <exp@plt>
movsd  -0x80(%rbp),%xmm1
movsd  -0x78(%rbp),%xmm2
movaps %xmm0,%xmm3
movsd  -0x70(%rbp),%xmm0
mulsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mulsd  -0x50(%rbp),%xmm0
divsd  -0x18(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm1
divsd  -0x10(%rbp),%xmm1
mulsd  -0x30(%rbp),%xmm1
movsd  0xcad(%rip),%xmm3
subsd  -0x50(%rbp),%xmm3
movsd  0xcf8(%rip),%xmm2
divsd  %xmm2,%xmm3
movsd  0xce4(%rip),%xmm2
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm1
divsd  -0x18(%rbp),%xmm1
addsd  %xmm1,%xmm0
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<rec_HRecFast_dxedlna>:
sub    $0x78,%rsp
movsd  %xmm5,0x48(%rsp)
movsd  %xmm4,0x8(%rsp)
movapd %xmm0,%xmm5
movsd  0xea4(%rip),%xmm4
movapd %xmm2,0x60(%rsp)
mulsd  %xmm2,%xmm4
movsd  %xmm1,0x28(%rsp)
divsd  %xmm1,%xmm4
movsd  0xe90(%rip),%xmm0
movsd  %xmm5,0x18(%rsp)
subsd  %xmm5,%xmm0
movsd  %xmm0,0x10(%rsp)
divsd  %xmm0,%xmm4
movapd %xmm4,0x30(%rsp)
movapd %xmm3,%xmm0
call   1060 <alphaB_PPB@plt>
movsd  0x8(%rsp),%xmm1
mulsd  0xe67(%rip),%xmm0
movsd  %xmm0,0x20(%rsp)
movsd  0xe61(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,(%rsp)
xorpd  %xmm2,%xmm2
xorpd  %xmm0,%xmm0
movapd %xmm0,0x50(%rsp)
ucomisd %xmm2,%xmm1
jb     11dd <rec_HRecFast_dxedlna+0x9d>
xorps  %xmm0,%xmm0
sqrtsd %xmm1,%xmm0
jmp    11ec <rec_HRecFast_dxedlna+0xac>
movapd %xmm1,%xmm0
call   1050 <sqrt@plt>
movsd  0x8(%rsp),%xmm1
movsd  (%rsp),%xmm2
mulsd  %xmm0,%xmm2
movsd  %xmm2,(%rsp)
movsd  0xe1e(%rip),%xmm0
divsd  %xmm1,%xmm0
call   1030 <exp@plt>
movapd %xmm0,%xmm2
mulsd  (%rsp),%xmm2
mulsd  0x20(%rsp),%xmm2
movapd 0x30(%rsp),%xmm1
mulsd  0xe00(%rip),%xmm1
addsd  0xe00(%rip),%xmm1
movapd %xmm1,%xmm0
addsd  %xmm2,%xmm0
divsd  %xmm0,%xmm1
movapd %xmm1,0x30(%rsp)
movsd  0xdbe(%rip),%xmm0
movsd  0x18(%rsp),%xmm1
ucomisd %xmm1,%xmm0
jbe    12a2 <rec_HRecFast_dxedlna+0x162>
cvttsd2si %xmm1,%edi
movsd  0xdd6(%rip),%xmm0
movsd  %xmm2,(%rsp)
call   1040 <pow@plt>
movsd  0xd94(%rip),%xmm1
subsd  %xmm0,%xmm1
cvttsd2si %xmm1,%edi
movsd  0xdbc(%rip),%xmm0
call   1040 <pow@plt>
movsd  (%rsp),%xmm2
movsd  0x18(%rsp),%xmm1
mulsd  0xdac(%rip),%xmm0
movapd %xmm0,0x50(%rsp)
movsd  0x28(%rsp),%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,(%rsp)
movsd  0x10(%rsp),%xmm0
mulsd  %xmm2,%xmm0
movsd  %xmm0,0x10(%rsp)
movsd  0xd83(%rip),%xmm0
divsd  0x8(%rsp),%xmm0
call   1030 <exp@plt>
mulsd  0x10(%rsp),%xmm0
movsd  (%rsp),%xmm1
mulsd  0x20(%rsp),%xmm1
subsd  %xmm1,%xmm0
movapd 0x50(%rsp),%xmm1
divsd  0x28(%rsp),%xmm1
mulsd  0x48(%rsp),%xmm1
unpcklpd %xmm0,%xmm1
movapd 0x30(%rsp),%xmm0
movsd  0xcf7(%rip),%xmm2
subsd  %xmm0,%xmm2
divsd  0xd3b(%rip),%xmm2
addsd  0xd3b(%rip),%xmm2
unpcklpd %xmm0,%xmm2
mulpd  %xmm1,%xmm2
movapd 0x60(%rsp),%xmm0
unpcklpd %xmm0,%xmm0
divpd  %xmm0,%xmm2
movapd %xmm2,%xmm0
unpckhpd %xmm2,%xmm0
addsd  %xmm2,%xmm0
add    $0x78,%rsp
ret

LLVM IR (-O0):
define double @rec_HRecFast_dxedlna(double %0, double %1, double %2, double %3, double %4, double %5) {
%7 = alloca double
%8 = alloca double
%9 = alloca double
%10 = alloca double
%11 = alloca double
%12 = alloca double
%13 = alloca double
%14 = alloca double
%15 = alloca double
%16 = alloca double
%17 = alloca double
store double %0, double* %7
store double %1, double* %8
store double %2, double* %9
store double %3, double* %10
store double %4, double* %11
store double %5, double* %12
%18 = load double, double* %9
%19 = fmul double 0x433090E1CA8B9039, %18
%20 = load double, double* %8
%21 = fdiv double %19, %20
%22 = load double, double* %7
%23 = fsub double 1.000000e+00, %22
%24 = fdiv double %21, %23
store double %24, double* %13
%25 = load double, double* %10
%26 = call double @alphaB_PPB(double %25)
%27 = fmul double 1.140000e+00, %26
store double %27, double* %14
%28 = load double, double* %11
%29 = fmul double 0x4464701ACFEF92A1, %28
%30 = load double, double* %11
%31 = call double @sqrt(double %30)
%32 = fmul double %29, %31
%33 = load double, double* %11
%34 = fdiv double 0xC00B32528D51413B, %33
%35 = call double @exp(double %34)
%36 = fmul double %32, %35
%37 = load double, double* %14
%38 = fmul double %36, %37
store double %38, double* %15
%39 = load double, double* %13
%40 = call double @llvm.fmuladd.f64(double 3.000000e+00, double %39, double 8.220600e+00)
%41 = load double, double* %13
%42 = call double @llvm.fmuladd.f64(double 3.000000e+00, double %41, double 8.220600e+00)
%43 = load double, double* %15
%44 = fadd double %42, %43
%45 = fdiv double %40, %44
store double %45, double* %16
%46 = load double, double* %7
%47 = fcmp olt double %46, 1.000000e+00
br i1 %47, label %48, label %56
48: ; preds = %6
%49 = load double, double* %7
%50 = fptosi double %49 to i32
%51 = call double @pow(i32 %50, double 4.639290e-01)
%52 = fsub double 1.000000e+00, %51
%53 = fptosi double %52 to i32
%54 = call double @pow(i32 %53, double 1.702370e+00)
%55 = fmul double 3.692020e-01, %54
store double %55, double* %17
br label %57
56: ; preds = %6
store double 0.000000e+00, double* %17
br label %57
57: ; preds = %56, %48
%58 = load double, double* %8
%59 = fneg double %58
%60 = load double, double* %7
%61 = fmul double %59, %60
%62 = load double, double* %7
%63 = fmul double %61, %62
%64 = load double, double* %14
%65 = load double, double* %15
%66 = load double, double* %7
%67 = fsub double 1.000000e+00, %66
%68 = fmul double %65, %67
%69 = load double, double* %11
%70 = fdiv double 0xC02465BDE9FCF0EC, %69
%71 = call double @exp(double %70)
%72 = fmul double %68, %71
%73 = call double @llvm.fmuladd.f64(double %63, double %64, double %72)
%74 = load double, double* %16
%75 = fmul double %73, %74
%76 = load double, double* %9
%77 = fdiv double %75, %76
%78 = load double, double* %17
%79 = load double, double* %8
%80 = fdiv double %78, %79
%81 = load double, double* %12
%82 = fmul double %80, %81
%83 = load double, double* %16
%84 = fsub double 1.000000e+00, %83
%85 = fdiv double %84, 0x402465BDE9FCF0EC
%86 = fadd double 0x3FB2D36E4F11087F, %85
%87 = fmul double %82, %86
%88 = load double, double* %9
%89 = fdiv double %87, %88
%90 = fadd double %77, %89
ret double %90
}

LLVM IR (-O2):
define double @rec_HRecFast_dxedlna(double %0, double %1, double %2, double %3, double %4, double %5) {
%7 = fmul double %2, 0x433090E1CA8B9039
%8 = fdiv double %7, %1
%9 = fsub double 1.000000e+00, %0
%10 = fdiv double %8, %9
%11 = tail call double @alphaB_PPB(double %3)
%12 = fmul double %11, 1.140000e+00
%13 = fmul double %4, 0x4464701ACFEF92A1
%14 = tail call double @sqrt(double %4)
%15 = fmul double %13, %14
%16 = fdiv double 0xC00B32528D51413B, %4
%17 = tail call double @exp(double %16)
%18 = fmul double %15, %17
%19 = fmul double %12, %18
%20 = tail call double @llvm.fmuladd.f64(double %10, double 3.000000e+00, double 8.220600e+00)
%21 = fadd double %20, %19
%22 = fdiv double %20, %21
%23 = fcmp olt double %0, 1.000000e+00
br i1 %23, label %24, label %31
24: ; preds = %6
%25 = fptosi double %0 to i32
%26 = tail call double @pow(i32 %25, double 4.639290e-01)
%27 = fsub double 1.000000e+00, %26
%28 = fptosi double %27 to i32
%29 = tail call double @pow(i32 %28, double 1.702370e+00)
%30 = fmul double %29, 3.692020e-01
br label %31
31: ; preds = %6, %24
%32 = phi double [ %30, %24 ], [ 0.000000e+00, %6 ]
%33 = fneg double %1
%34 = fmul double %33, %0
%35 = fmul double %34, %0
%36 = fmul double %9, %19
%37 = fdiv double 0xC02465BDE9FCF0EC, %4
%38 = tail call double @exp(double %37)
%39 = fmul double %36, %38
%40 = tail call double @llvm.fmuladd.f64(double %35, double %12, double %39)
%41 = fdiv double %32, %1
%42 = fmul double %41, %5
%43 = fsub double 1.000000e+00, %22
%44 = fdiv double %43, 0x402465BDE9FCF0EC
%45 = fadd double %44, 0x3FB2D36E4F11087F
%46 = insertelement <2 x double> poison, double %22, i64 1
%47 = insertelement <2 x double> %46, double %45, i64 0
%48 = insertelement <2 x double> poison, double %42, i64 0
%49 = insertelement <2 x double> %48, double %40, i64 1
%50 = fmul <2 x double> %47, %49
%51 = insertelement <2 x double> poison, double %2, i64 0
%52 = shufflevector <2 x double> %51, <2 x double> poison, <2 x i32> zeroinitializer
%53 = fdiv <2 x double> %50, %52
%54 = shufflevector <2 x double> %53, <2 x double> poison, <2 x i32> <i32 1, i32 undef>
%55 = fadd <2 x double> %53, %54
%56 = extractelement <2 x double> %55, i64 0
ret double %56
}