double smallj0(x)
double x;
{
  static double p[10] = {
      0.1641556014884554385346147435e+25, -0.3943559664767296636012616471e+24,
      0.2172018385924539313982287997e+23, -0.4814859952069817648285245941e+21,
      0.5345457598841972345381674607e+19, -0.3301538925689637686465426220e+17,
      0.1187390681211042949874031474e+15, -0.2479851167896144439689877514e+12,
      0.2803148940831953934479400118e+09, -0.1336625500481224741885945416e+06};
  static double q[10] = {
      0.1641556014884554385346137617e+25, 0.1603303724440893273539045602e+23,
      0.7913043777646405204323616203e+20, 0.2613165313325153278086066185e+18,
      0.6429607918826017759289213100e+15, 0.1237672982083407903483177730e+13,
      0.1893012093677918995179541438e+10, 0.2263381356781110003609399116e+07,
      0.1974019272727281783930443513e+04, 0.1000000000000000000000000000e+01};
  double xsq = x * x;
  return POLYNOM9(xsq, p) / POLYNOM9(xsq, q);
}


x86 assembly (clang, -O0):
<smallj0>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
mulsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x10(%rbp)
movsd  -0x10(%rbp),%xmm0
lea    0x2ef8(%rip),%rdi
mov    $0x1,%al
call   1030 <POLYNOM9@plt>
mov    %eax,-0x14(%rbp)
movsd  -0x10(%rbp),%xmm0
lea    0x2f32(%rip),%rdi
mov    $0x1,%al
call   1030 <POLYNOM9@plt>
mov    %eax,%ecx
mov    -0x14(%rbp),%eax
cltd
idiv   %ecx
cvtsi2sd %eax,%xmm0
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<smallj0>:
push   %rbx
sub    $0x10,%rsp
mulsd  %xmm0,%xmm0
movsd  %xmm0,0x8(%rsp)
lea    0x2f0a(%rip),%rdi
mov    $0x1,%al
call   1030 <POLYNOM9@plt>
mov    %eax,%ebx
lea    0x2f4a(%rip),%rdi
movsd  0x8(%rsp),%xmm0
mov    $0x1,%al
call   1030 <POLYNOM9@plt>
mov    %eax,%ecx
mov    %ebx,%eax
cltd
idiv   %ecx
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
add    $0x10,%rsp
pop    %rbx
ret

LLVM IR (-O0):
@smallj0.p = internal global [10 x double] [double 0x44F5B9CEE0C98FC0, double 0xC4D4E0851B1AC833, double 0x449265D0693611FB, double 0xC43A19F5E19E5D95, double 0x43D28BB8575738DA, double 0xC35D52D49531D3DE, double 0x42DAFF85FE7C0413, double 0xC24CDE884B3ACEA6, double 0x41B0B5440E154C4B, double 0xC10050F4667FA14A]
@smallj0.q = internal global [10 x double] [double 0x44F5B9CEE0C98FC0, double 0x448B2938B5D8F192, double 0x4411289F27A06D2A, double 0x438D031001801830, double 0x43024627CCD65B4E, double 0x427202B0EF643687, double 0x41DC35449F6B6306, double 0x414144AAADAB00E0, double 0x409ED813BC3AD585, double 1.000000e+00]
define double @smallj0(double %0) {
%2 = alloca double
%3 = alloca double
store double %0, double* %2
%4 = load double, double* %2
%5 = load double, double* %2
%6 = fmul double %4, %5
store double %6, double* %3
%7 = load double, double* %3
%8 = call i32 (double, double*, ...) bitcast (i32 (...)* @POLYNOM9 to i32 (double, double*, ...)*)(double %7, double* getelementptr inbounds ([10 x double], [10 x double]* @smallj0.p, i64 0, i64 0))
%9 = load double, double* %3
%10 = call i32 (double, double*, ...) bitcast (i32 (...)* @POLYNOM9 to i32 (double, double*, ...)*)(double %9, double* getelementptr inbounds ([10 x double], [10 x double]* @smallj0.q, i64 0, i64 0))
%11 = sdiv i32 %8, %10
%12 = sitofp i32 %11 to double
ret double %12
}

LLVM IR (-O2):
@smallj0.p = internal global [10 x double] [double 0x44F5B9CEE0C98FC0, double 0xC4D4E0851B1AC833, double 0x449265D0693611FB, double 0xC43A19F5E19E5D95, double 0x43D28BB8575738DA, double 0xC35D52D49531D3DE, double 0x42DAFF85FE7C0413, double 0xC24CDE884B3ACEA6, double 0x41B0B5440E154C4B, double 0xC10050F4667FA14A]
@smallj0.q = internal global [10 x double] [double 0x44F5B9CEE0C98FC0, double 0x448B2938B5D8F192, double 0x4411289F27A06D2A, double 0x438D031001801830, double 0x43024627CCD65B4E, double 0x427202B0EF643687, double 0x41DC35449F6B6306, double 0x414144AAADAB00E0, double 0x409ED813BC3AD585, double 1.000000e+00]
define double @smallj0(double %0) {
%2 = fmul double %0, %0
%3 = tail call i32 (double, double*, ...) bitcast (i32 (...)* @POLYNOM9 to i32 (double, double*, ...)*)(double %2, double* getelementptr inbounds ([10 x double], [10 x double]* @smallj0.p, i64 0, i64 0))
%4 = tail call i32 (double, double*, ...) bitcast (i32 (...)* @POLYNOM9 to i32 (double, double*, ...)*)(double %2, double* getelementptr inbounds ([10 x double], [10 x double]* @smallj0.q, i64 0, i64 0))
%5 = sdiv i32 %3, %4
%6 = sitofp i32 %5 to double
ret double %6
}