double radiation_logarithm(double Z) {
  const int i = (int)Z;
  if ((i >= 1) && (i <= 92)) {
    static double Lz[92] = {
        202.4, 151.9, 159.9, 172.3, 177.9, 178.3, 176.6, 173.4, 170.0, 165.8,
        165.8, 167.1, 169.1, 170.8, 172.2, 173.4, 174.3, 174.8, 175.1, 175.6,
        176.2, 176.8, 0.0,   0.0,   0.0,   175.8, 0.0,   0.0,   173.1, 0.0,
        0.0,   173.0, 0.0,   0.0,   173.5, 0.0,   0.0,   0.0,   0.0,   0.0,
        0.0,   175.9, 0.0,   0.0,   0.0,   0.0,   0.0,   0.0,   0.0,   177.4,
        0.0,   0.0,   178.6, 0.0,   0.0,   0.0,   0.0,   0.0,   0.0,   0.0,
        0.0,   0.0,   0.0,   0.0,   0.0,   0.0,   0.0,   0.0,   0.0,   0.0,
        0.0,   0.0,   0.0,   177.6, 0.0,   0.0,   0.0,   0.0,   0.0,   0.0,
        0.0,   178.0, 0.0,   0.0,   0.0,   0.0,   0.0,   0.0,   0.0,   0.0,
        0.0,   179.8};
    const double l = Lz[i - 1];
    if (l > 0)
      return l;
  }
  return 182.7;
}


x86 assembly (clang, -O0):
<radiation_logarithm>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x10(%rbp)
cvttsd2si -0x10(%rbp),%eax
mov    %eax,-0x14(%rbp)
cmpl   $0x1,-0x14(%rbp)
jl     1165 <radiation_logarithm+0x65>
cmpl   $0x5c,-0x14(%rbp)
jg     1165 <radiation_logarithm+0x65>
mov    -0x14(%rbp),%eax
sub    $0x1,%eax
movslq %eax,%rcx
lea    0x2efb(%rip),%rax
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1160 <radiation_logarithm+0x60>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1172 <radiation_logarithm+0x72>
jmp    1165 <radiation_logarithm+0x65>
movsd  0xe93(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<radiation_logarithm>:
cvttsd2si %xmm0,%eax
add    $0xffffffff,%eax
cmp    $0x5b,%eax
ja     1124 <radiation_logarithm+0x24>
mov    %eax,%eax
lea    0xefb(%rip),%rcx
movsd  (%rcx,%rax,8),%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
ja     112c <radiation_logarithm+0x2c>
movsd  0xed4(%rip),%xmm0
ret

LLVM IR (-O0):
@radiation_logarithm.Lz = internal global [92 x double] [double 2.024000e+02, double 1.519000e+02, double 1.599000e+02, double 1.723000e+02, double 1.779000e+02, double 1.783000e+02, double 1.766000e+02, double 1.734000e+02, double 1.700000e+02, double 1.658000e+02, double 1.658000e+02, double 1.671000e+02, double 1.691000e+02, double 1.708000e+02, double 1.722000e+02, double 1.734000e+02, double 1.743000e+02, double 1.748000e+02, double 1.751000e+02, double 1.756000e+02, double 1.762000e+02, double 1.768000e+02, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 1.758000e+02, double 0.000000e+00, double 0.000000e+00, double 1.731000e+02, double 0.000000e+00, double 0.000000e+00, double 1.730000e+02, double 0.000000e+00, double 0.000000e+00, double 1.735000e+02, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 1.759000e+02, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 1.774000e+02, double 0.000000e+00, double 0.000000e+00, double 1.786000e+02, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 1.776000e+02, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 1.780000e+02, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 1.798000e+02]
define double @radiation_logarithm(double %0) {
%2 = alloca double
%3 = alloca double
%4 = alloca i32
%5 = alloca double
store double %0, double* %3
%6 = load double, double* %3
%7 = fptosi double %6 to i32
store i32 %7, i32* %4
%8 = load i32, i32* %4
%9 = icmp sge i32 %8, 1
br i1 %9, label %10, label %24
10: ; preds = %1
%11 = load i32, i32* %4
%12 = icmp sle i32 %11, 92
br i1 %12, label %13, label %24
13: ; preds = %10
%14 = load i32, i32* %4
%15 = sub nsw i32 %14, 1
%16 = sext i32 %15 to i64
%17 = getelementptr inbounds [92 x double], [92 x double]* @radiation_logarithm.Lz, i64 0, i64 %16
%18 = load double, double* %17
store double %18, double* %5
%19 = load double, double* %5
%20 = fcmp ogt double %19, 0.000000e+00
br i1 %20, label %21, label %23
21: ; preds = %13
%22 = load double, double* %5
store double %22, double* %2
br label %25
23: ; preds = %13
br label %24
24: ; preds = %23, %10, %1
store double 1.827000e+02, double* %2
br label %25
25: ; preds = %24, %21
%26 = load double, double* %2
ret double %26
}

LLVM IR (-O2):
@radiation_logarithm.Lz = internal unnamed_addr constant [92 x double] [double 2.024000e+02, double 1.519000e+02, double 1.599000e+02, double 1.723000e+02, double 1.779000e+02, double 1.783000e+02, double 1.766000e+02, double 1.734000e+02, double 1.700000e+02, double 1.658000e+02, double 1.658000e+02, double 1.671000e+02, double 1.691000e+02, double 1.708000e+02, double 1.722000e+02, double 1.734000e+02, double 1.743000e+02, double 1.748000e+02, double 1.751000e+02, double 1.756000e+02, double 1.762000e+02, double 1.768000e+02, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 1.758000e+02, double 0.000000e+00, double 0.000000e+00, double 1.731000e+02, double 0.000000e+00, double 0.000000e+00, double 1.730000e+02, double 0.000000e+00, double 0.000000e+00, double 1.735000e+02, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 1.759000e+02, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 1.774000e+02, double 0.000000e+00, double 0.000000e+00, double 1.786000e+02, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 1.776000e+02, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 1.780000e+02, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 1.798000e+02]
define double @radiation_logarithm(double %0) {
%2 = fptosi double %0 to i32
%3 = add i32 %2, -1
%4 = icmp ult i32 %3, 92
br i1 %4, label %5, label %10
5: ; preds = %1
%6 = zext i32 %3 to i64
%7 = getelementptr inbounds [92 x double], [92 x double]* @radiation_logarithm.Lz, i64 0, i64 %6
%8 = load double, double* %7
%9 = fcmp ogt double %8, 0.000000e+00
br i1 %9, label %11, label %10
10: ; preds = %5, %1
br label %11
11: ; preds = %5, %10
%12 = phi double [ 1.827000e+02, %10 ], [ %8, %5 ]
ret double %12
}