#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 j0(double);
double j1(double);
double jn(int n, double x) {
  int i;
  double a, b, temp;
  double xsq, t;
  if (n < 0) {
    n = -n;
    x = -x;
  }
  if (n == 0)
    return j0(x);
  if (n == 1)
    return j1(x);
  if (x == 0)
    return 0;
  if (n > x)
    goto recurs;
  a = j0(x);
  b = j1(x);
  for (i = 1; i < n; i++) {
    temp = b;
    b = (2 * i / x) * b - a;
    a = temp;
  }
  return b;
recurs:
  xsq = x * x;
  for (t = 0, i = n + 16; i > n; i--)
    t = xsq / (2 * i - t);
  t = x / (2 * n - t);
  a = t;
  b = 1;
  for (i = n - 1; i > 0; i--) {
    temp = b;
    b = (2 * i / x) * b - a;
    a = temp;
  }
  return t * j0(x) / b;
}


x86 assembly (clang, -O0):
<jn>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %edi,-0xc(%rbp)
movsd  %xmm0,-0x18(%rbp)
cmpl   $0x0,-0xc(%rbp)
jge    1163 <jn+0x43>
xor    %eax,%eax
sub    -0xc(%rbp),%eax
mov    %eax,-0xc(%rbp)
movsd  -0x18(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x18(%rbp)
cmpl   $0x0,-0xc(%rbp)
jne    1181 <jn+0x61>
movsd  -0x18(%rbp),%xmm0
call   1040 <j0@plt>
movsd  %xmm0,-0x8(%rbp)
jmp    139f <jn+0x27f>
cmpl   $0x1,-0xc(%rbp)
jne    119f <jn+0x7f>
movsd  -0x18(%rbp),%xmm0
call   1030 <j1@plt>
movsd  %xmm0,-0x8(%rbp)
jmp    139f <jn+0x27f>
movsd  -0x18(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    11c4 <jn+0xa4>
jp     11c4 <jn+0xa4>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    139f <jn+0x27f>
cvtsi2sdl -0xc(%rbp),%xmm0
ucomisd -0x18(%rbp),%xmm0
jbe    11d9 <jn+0xb9>
jmp    1278 <jn+0x158>
movsd  -0x18(%rbp),%xmm0
call   1040 <j0@plt>
movsd  %xmm0,-0x28(%rbp)
movsd  -0x18(%rbp),%xmm0
call   1030 <j1@plt>
movsd  %xmm0,-0x30(%rbp)
movl   $0x1,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1269 <jn+0x149>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x1c(%rbp),%eax
shl    $0x1,%eax
cvtsi2sd %eax,%xmm0
divsd  -0x18(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
movsd  -0x28(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    11fe <jn+0xde>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    139f <jn+0x27f>
movsd  -0x18(%rbp),%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0xc(%rbp),%eax
add    $0x10,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jle    12cf <jn+0x1af>
movsd  -0x40(%rbp),%xmm0
mov    -0x1c(%rbp),%eax
shl    $0x1,%eax
cvtsi2sd %eax,%xmm1
subsd  -0x48(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x1c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x1c(%rbp)
jmp    1298 <jn+0x178>
movsd  -0x18(%rbp),%xmm0
mov    -0xc(%rbp),%eax
shl    $0x1,%eax
cvtsi2sd %eax,%xmm1
subsd  -0x48(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  0xd02(%rip),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x1c(%rbp)
cmpl   $0x0,-0x1c(%rbp)
jle    1375 <jn+0x255>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x1c(%rbp),%eax
shl    $0x1,%eax
cvtsi2sd %eax,%xmm0
divsd  -0x18(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
movsd  -0x28(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x1c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x1c(%rbp)
jmp    130c <jn+0x1ec>
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x18(%rbp),%xmm0
call   1040 <j0@plt>
movaps %xmm0,%xmm1
movsd  -0x50(%rbp),%xmm0
mulsd  %xmm1,%xmm0
divsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<jn>:
push   %rbx
sub    $0x20,%rsp
mov    %edi,%ebx
neg    %ebx
cmovs  %edi,%ebx
test   %edi,%edi
jns    1138 <jn+0x18>
xorpd  0xec8(%rip),%xmm0
cmp    $0x1,%ebx
je     114b <jn+0x2b>
test   %ebx,%ebx
jne    1155 <jn+0x35>
add    $0x20,%rsp
pop    %rbx
jmp    1040 <j0@plt>
add    $0x20,%rsp
pop    %rbx
jmp    1030 <j1@plt>
xorpd  %xmm3,%xmm3
ucomisd %xmm3,%xmm0
jne    1165 <jn+0x45>
jnp    12c7 <jn+0x1a7>
cvtsi2sd %ebx,%xmm1
ucomisd %xmm0,%xmm1
jbe    11f9 <jn+0xd9>
movapd %xmm0,%xmm3
mulsd  %xmm0,%xmm3
lea    0x10(%rbx),%eax
lea    (%rbx,%rbx,1),%ecx
add    $0x20,%ecx
xorpd  %xmm1,%xmm1
nopl   0x0(%rax,%rax,1)
xorps  %xmm2,%xmm2
cvtsi2sd %ecx,%xmm2
subsd  %xmm1,%xmm2
movapd %xmm3,%xmm1
divsd  %xmm2,%xmm1
add    $0xffffffff,%eax
add    $0xfffffffe,%ecx
cmp    %ebx,%eax
jg     1190 <jn+0x70>
lea    (%rbx,%rbx,1),%eax
xorps  %xmm2,%xmm2
cvtsi2sd %eax,%xmm2
subsd  %xmm1,%xmm2
movapd %xmm0,%xmm1
divsd  %xmm2,%xmm1
cmp    $0x2,%ebx
movsd  %xmm1,0x10(%rsp)
jl     1247 <jn+0x127>
test   $0x1,%bl
jne    1251 <jn+0x131>
lea    (%rbx,%rbx,1),%eax
add    $0xfffffffe,%eax
xorps  %xmm3,%xmm3
cvtsi2sd %eax,%xmm3
lea    -0x1(%rbx),%eax
divsd  %xmm0,%xmm3
subsd  %xmm1,%xmm3
movsd  0xe25(%rip),%xmm1
movapd %xmm3,%xmm4
jmp    125b <jn+0x13b>
movapd %xmm0,0x10(%rsp)
call   1040 <j0@plt>
movsd  %xmm0,0x8(%rsp)
movapd 0x10(%rsp),%xmm0
call   1030 <j1@plt>
movapd %xmm0,%xmm3
cmp    $0x2,%ebx
jl     12c7 <jn+0x1a7>
movsd  0x8(%rsp),%xmm4
movapd 0x10(%rsp),%xmm2
lea    -0x1(%rbx),%eax
cmp    $0x2,%ebx
jne    12d1 <jn+0x1b1>
movsd  0xdce(%rip),%xmm0
jmp    1323 <jn+0x203>
movsd  0xdc9(%rip),%xmm4
jmp    12ac <jn+0x18c>
movsd  0xdbf(%rip),%xmm3
mov    %ebx,%eax
cmp    $0x2,%ebx
je     12ac <jn+0x18c>
lea    (%rax,%rax,1),%ecx
add    $0xfffffffc,%ecx
add    $0x1,%eax
nopl   0x0(%rax)
lea    0x2(%rcx),%edx
xorps  %xmm2,%xmm2
cvtsi2sd %edx,%xmm2
divsd  %xmm0,%xmm2
mulsd  %xmm3,%xmm2
subsd  %xmm1,%xmm2
xorps  %xmm4,%xmm4
cvtsi2sd %ecx,%xmm4
divsd  %xmm0,%xmm4
mulsd  %xmm2,%xmm4
subsd  %xmm3,%xmm4
add    $0xfffffffc,%ecx
add    $0xfffffffe,%eax
movapd %xmm4,%xmm3
movapd %xmm2,%xmm1
cmp    $0x2,%eax
ja     1270 <jn+0x150>
movsd  %xmm4,0x8(%rsp)
call   1040 <j0@plt>
movsd  0x10(%rsp),%xmm3
mulsd  %xmm0,%xmm3
divsd  0x8(%rsp),%xmm3
movapd %xmm3,%xmm0
add    $0x20,%rsp
pop    %rbx
ret
mov    %eax,%edx
and    $0xfffffffe,%edx
mov    $0x4,%ecx
nopl   0x0(%rax,%rax,1)
movapd %xmm4,%xmm0
movapd %xmm3,%xmm1
lea    -0x2(%rcx),%esi
xorps  %xmm4,%xmm4
cvtsi2sd %esi,%xmm4
divsd  %xmm2,%xmm4
mulsd  %xmm3,%xmm4
subsd  %xmm0,%xmm4
xorps  %xmm3,%xmm3
cvtsi2sd %ecx,%xmm3
divsd  %xmm2,%xmm3
mulsd  %xmm4,%xmm3
subsd  %xmm1,%xmm3
add    $0x4,%ecx
add    $0xfffffffe,%edx
jne    12e0 <jn+0x1c0>
add    $0xfffffffe,%ecx
xorps  %xmm0,%xmm0
cvtsi2sd %ecx,%xmm0
test   $0x1,%al
je     12c7 <jn+0x1a7>
divsd  %xmm2,%xmm0
mulsd  %xmm3,%xmm0
subsd  %xmm4,%xmm0
jmp    12cb <jn+0x1ab>

LLVM IR (-O0):
define double @jn(i32 %0, double %1) {
%3 = alloca double
%4 = alloca i32
%5 = alloca double
%6 = alloca i32
%7 = alloca double
%8 = alloca double
%9 = alloca double
%10 = alloca double
%11 = alloca double
store i32 %0, i32* %4
store double %1, double* %5
%12 = load i32, i32* %4
%13 = icmp slt i32 %12, 0
br i1 %13, label %14, label %19
14: ; preds = %2
%15 = load i32, i32* %4
%16 = sub nsw i32 0, %15
store i32 %16, i32* %4
%17 = load double, double* %5
%18 = fneg double %17
store double %18, double* %5
br label %19
19: ; preds = %14, %2
%20 = load i32, i32* %4
%21 = icmp eq i32 %20, 0
br i1 %21, label %22, label %25
22: ; preds = %19
%23 = load double, double* %5
%24 = call double @j0(double %23)
store double %24, double* %3
br label %124
25: ; preds = %19
%26 = load i32, i32* %4
%27 = icmp eq i32 %26, 1
br i1 %27, label %28, label %31
28: ; preds = %25
%29 = load double, double* %5
%30 = call double @j1(double %29)
store double %30, double* %3
br label %124
31: ; preds = %25
%32 = load double, double* %5
%33 = fcmp oeq double %32, 0.000000e+00
br i1 %33, label %34, label %35
34: ; preds = %31
store double 0.000000e+00, double* %3
br label %124
35: ; preds = %31
%36 = load i32, i32* %4
%37 = sitofp i32 %36 to double
%38 = load double, double* %5
%39 = fcmp ogt double %37, %38
br i1 %39, label %40, label %41
40: ; preds = %35
br label %67
41: ; preds = %35
%42 = load double, double* %5
%43 = call double @j0(double %42)
store double %43, double* %7
%44 = load double, double* %5
%45 = call double @j1(double %44)
store double %45, double* %8
store i32 1, i32* %6
br label %46
46: ; preds = %62, %41
%47 = load i32, i32* %6
%48 = load i32, i32* %4
%49 = icmp slt i32 %47, %48
br i1 %49, label %50, label %65
50: ; preds = %46
%51 = load double, double* %8
store double %51, double* %9
%52 = load i32, i32* %6
%53 = mul nsw i32 2, %52
%54 = sitofp i32 %53 to double
%55 = load double, double* %5
%56 = fdiv double %54, %55
%57 = load double, double* %8
%58 = load double, double* %7
%59 = fneg double %58
%60 = call double @llvm.fmuladd.f64(double %56, double %57, double %59)
store double %60, double* %8
%61 = load double, double* %9
store double %61, double* %7
br label %62
62: ; preds = %50
%63 = load i32, i32* %6
%64 = add nsw i32 %63, 1
store i32 %64, i32* %6
br label %46
65: ; preds = %46
%66 = load double, double* %8
store double %66, double* %3
br label %124
67: ; preds = %40
%68 = load double, double* %5
%69 = load double, double* %5
%70 = fmul double %68, %69
store double %70, double* %10
store double 0.000000e+00, double* %11
%71 = load i32, i32* %4
%72 = add nsw i32 %71, 16
store i32 %72, i32* %6
br label %73
73: ; preds = %85, %67
%74 = load i32, i32* %6
%75 = load i32, i32* %4
%76 = icmp sgt i32 %74, %75
br i1 %76, label %77, label %88
77: ; preds = %73
%78 = load double, double* %10
%79 = load i32, i32* %6
%80 = mul nsw i32 2, %79
%81 = sitofp i32 %80 to double
%82 = load double, double* %11
%83 = fsub double %81, %82
%84 = fdiv double %78, %83
store double %84, double* %11
br label %85
85: ; preds = %77
%86 = load i32, i32* %6
%87 = add nsw i32 %86, -1
store i32 %87, i32* %6
br label %73
88: ; preds = %73
%89 = load double, double* %5
%90 = load i32, i32* %4
%91 = mul nsw i32 2, %90
%92 = sitofp i32 %91 to double
%93 = load double, double* %11
%94 = fsub double %92, %93
%95 = fdiv double %89, %94
store double %95, double* %11
%96 = load double, double* %11
store double %96, double* %7
store double 1.000000e+00, double* %8
%97 = load i32, i32* %4
%98 = sub nsw i32 %97, 1
store i32 %98, i32* %6
br label %99
99: ; preds = %114, %88
%100 = load i32, i32* %6
%101 = icmp sgt i32 %100, 0
br i1 %101, label %102, label %117
102: ; preds = %99
%103 = load double, double* %8
store double %103, double* %9
%104 = load i32, i32* %6
%105 = mul nsw i32 2, %104
%106 = sitofp i32 %105 to double
%107 = load double, double* %5
%108 = fdiv double %106, %107
%109 = load double, double* %8
%110 = load double, double* %7
%111 = fneg double %110
%112 = call double @llvm.fmuladd.f64(double %108, double %109, double %111)
store double %112, double* %8
%113 = load double, double* %9
store double %113, double* %7
br label %114
114: ; preds = %102
%115 = load i32, i32* %6
%116 = add nsw i32 %115, -1
store i32 %116, i32* %6
br label %99
117: ; preds = %99
%118 = load double, double* %11
%119 = load double, double* %5
%120 = call double @j0(double %119)
%121 = fmul double %118, %120
%122 = load double, double* %8
%123 = fdiv double %121, %122
store double %123, double* %3
br label %124
124: ; preds = %117, %65, %34, %28, %22
%125 = load double, double* %3
ret double %125
}

LLVM IR (-O2):
define double @jn(i32 %0, double %1) {
%3 = icmp slt i32 %0, 0
%4 = sub nsw i32 0, %0
%5 = fneg double %1
%6 = select i1 %3, i32 %4, i32 %0
%7 = select i1 %3, double %5, double %1
switch i32 %6, label %12 [
i32 0, label %8
i32 1, label %10
]
8: ; preds = %2
%9 = tail call double @j0(double %7)
br label %114
10: ; preds = %2
%11 = tail call double @j1(double %7)
br label %114
12: ; preds = %2
%13 = fcmp oeq double %7, 0.000000e+00
br i1 %13, label %114, label %14
14: ; preds = %12
%15 = sitofp i32 %6 to double
%16 = fcmp olt double %7, %15
br i1 %16, label %46, label %17
17: ; preds = %14
%18 = tail call double @j0(double %7)
%19 = tail call double @j1(double %7)
%20 = icmp sgt i32 %6, 1
br i1 %20, label %21, label %114
21: ; preds = %17
%22 = add nsw i32 %6, -1
%23 = and i32 %22, 1
%24 = icmp eq i32 %6, 2
br i1 %24, label %104, label %25
25: ; preds = %21
%26 = and i32 %22, -2
br label %27
27: ; preds = %27, %25
%28 = phi double [ %19, %25 ], [ %42, %27 ]
%29 = phi double [ %18, %25 ], [ %36, %27 ]
%30 = phi i32 [ 1, %25 ], [ %43, %27 ]
%31 = phi i32 [ 0, %25 ], [ %44, %27 ]
%32 = shl nuw nsw i32 %30, 1
%33 = sitofp i32 %32 to double
%34 = fdiv double %33, %7
%35 = fneg double %29
%36 = tail call double @llvm.fmuladd.f64(double %34, double %28, double %35)
%37 = shl nuw i32 %30, 1
%38 = add i32 %37, 2
%39 = sitofp i32 %38 to double
%40 = fdiv double %39, %7
%41 = fneg double %28
%42 = tail call double @llvm.fmuladd.f64(double %40, double %36, double %41)
%43 = add nuw nsw i32 %30, 2
%44 = add i32 %31, 2
%45 = icmp eq i32 %44, %26
br i1 %45, label %101, label %27
46: ; preds = %14
%47 = fmul double %7, %7
%48 = add nuw nsw i32 %6, 16
br label %49
49: ; preds = %46, %49
%50 = phi double [ 0.000000e+00, %46 ], [ %55, %49 ]
%51 = phi i32 [ %48, %46 ], [ %56, %49 ]
%52 = shl nsw i32 %51, 1
%53 = sitofp i32 %52 to double
%54 = fsub double %53, %50
%55 = fdiv double %47, %54
%56 = add nsw i32 %51, -1
%57 = icmp sgt i32 %56, %6
br i1 %57, label %49, label %58
58: ; preds = %49
%59 = shl nuw nsw i32 %6, 1
%60 = sitofp i32 %59 to double
%61 = fsub double %60, %55
%62 = fdiv double %7, %61
%63 = icmp sgt i32 %6, 1
br i1 %63, label %64, label %96
64: ; preds = %58
%65 = and i32 %6, 1
%66 = icmp eq i32 %65, 0
br i1 %66, label %67, label %73
67: ; preds = %64
%68 = add nsw i32 %6, -1
%69 = shl nuw nsw i32 %68, 1
%70 = sitofp i32 %69 to double
%71 = fdiv double %70, %7
%72 = fsub double %71, %62
br label %73
73: ; preds = %67, %64
%74 = phi i32 [ %6, %64 ], [ %68, %67 ]
%75 = phi double [ 1.000000e+00, %64 ], [ %72, %67 ]
%76 = phi double [ %62, %64 ], [ 1.000000e+00, %67 ]
%77 = phi double [ undef, %64 ], [ %72, %67 ]
%78 = icmp eq i32 %6, 2
br i1 %78, label %96, label %79
79: ; preds = %73, %79
%80 = phi i32 [ %89, %79 ], [ %74, %73 ]
%81 = phi double [ %94, %79 ], [ %75, %73 ]
%82 = phi double [ %88, %79 ], [ %76, %73 ]
%83 = add nsw i32 %80, -1
%84 = shl nuw nsw i32 %83, 1
%85 = sitofp i32 %84 to double
%86 = fdiv double %85, %7
%87 = fneg double %82
%88 = tail call double @llvm.fmuladd.f64(double %86, double %81, double %87)
%89 = add nsw i32 %80, -2
%90 = shl nuw nsw i32 %89, 1
%91 = sitofp i32 %90 to double
%92 = fdiv double %91, %7
%93 = fneg double %81
%94 = tail call double @llvm.fmuladd.f64(double %92, double %88, double %93)
%95 = icmp ugt i32 %83, 2
br i1 %95, label %79, label %96
96: ; preds = %73, %79, %58
%97 = phi double [ 1.000000e+00, %58 ], [ %77, %73 ], [ %94, %79 ]
%98 = tail call double @j0(double %7)
%99 = fmul double %62, %98
%100 = fdiv double %99, %97
br label %114
101: ; preds = %27
%102 = shl nuw i32 %43, 1
%103 = sitofp i32 %102 to double
br label %104
104: ; preds = %101, %21
%105 = phi double [ undef, %21 ], [ %42, %101 ]
%106 = phi double [ %19, %21 ], [ %42, %101 ]
%107 = phi double [ %18, %21 ], [ %36, %101 ]
%108 = phi double [ 2.000000e+00, %21 ], [ %103, %101 ]
%109 = icmp eq i32 %23, 0
br i1 %109, label %114, label %110
110: ; preds = %104
%111 = fdiv double %108, %7
%112 = fneg double %107
%113 = tail call double @llvm.fmuladd.f64(double %111, double %106, double %112)
br label %114
114: ; preds = %110, %104, %17, %12, %96, %10, %8
%115 = phi double [ %9, %8 ], [ %11, %10 ], [ %100, %96 ], [ 0.000000e+00, %12 ], [ %19, %17 ], [ %105, %104 ], [ %113, %110 ]
ret double %115
}