#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 */
int abs(int);
double sqrt(double);
double ynnm(int n, int m) {
  int pm, k;
  double a, ynnm;
  a = 1.0 / sqrt(4.0 * 3.14159265358979323846);
  pm = abs(m);
  if (n < pm) {
    ynnm = 0;
  } else if (n == 0) {
    ynnm = a;
  } else {
    ynnm = a;
    for (k = 1; k <= n; k++) {
      ynnm = sqrt((2.0 * k + 1) / (8.0 * k)) * ynnm;
    }
    if (n != pm) {
      for (k = n - 1; k >= pm; k--) {
        ynnm = sqrt(1.0 * (n + k + 1.0) / (n - k)) * ynnm;
      }
    }
  }
  return ynnm;
}


x86 assembly (clang, -O0):
<ynnm>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
movsd  0xed2(%rip),%xmm0
call   1040 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  0xeba(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x8(%rbp),%edi
call   1030 <abs@plt>
mov    %eax,-0xc(%rbp)
mov    -0x4(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1173 <ynnm+0x53>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
jmp    126d <ynnm+0x14d>
cmpl   $0x0,-0x4(%rbp)
jne    118c <ynnm+0x6c>
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
jmp    1268 <ynnm+0x148>
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
movl   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     11f0 <ynnm+0xd0>
cvtsi2sdl -0x10(%rbp),%xmm0
addsd  %xmm0,%xmm0
movsd  0xe46(%rip),%xmm1
addsd  %xmm1,%xmm0
cvtsi2sdl -0x10(%rbp),%xmm2
movsd  0xe45(%rip),%xmm1
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
call   1040 <sqrt@plt>
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    119d <ynnm+0x7d>
mov    -0x4(%rbp),%eax
cmp    -0xc(%rbp),%eax
je     1263 <ynnm+0x143>
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jl     125e <ynnm+0x13e>
mov    -0x4(%rbp),%eax
add    -0x10(%rbp),%eax
cvtsi2sd %eax,%xmm1
movsd  0xddd(%rip),%xmm0
addsd  %xmm0,%xmm1
movsd  0xdd1(%rip),%xmm0
mulsd  %xmm1,%xmm0
mov    -0x4(%rbp),%eax
sub    -0x10(%rbp),%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
call   1040 <sqrt@plt>
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x10(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x10(%rbp)
jmp    1205 <ynnm+0xe5>
jmp    1263 <ynnm+0x143>
jmp    1268 <ynnm+0x148>
jmp    126d <ynnm+0x14d>
movsd  -0x20(%rbp),%xmm0
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<ynnm>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x10,%rsp
mov    %esi,%r15d
neg    %r15d
cmovs  %esi,%r15d
xorpd  %xmm2,%xmm2
cmp    %edi,%r15d
jg     131d <ynnm+0x20d>
mov    %edi,%r14d
test   %edi,%edi
je     1161 <ynnm+0x51>
jle    116e <ynnm+0x5e>
cmp    $0x1,%r14d
jne    117b <ynnm+0x6b>
movsd  0xebe(%rip),%xmm1
movsd  0xeae(%rip),%xmm2
test   $0x1,%r14b
jne    127e <ynnm+0x16e>
jmp    12bd <ynnm+0x1ad>
movsd  0xe97(%rip),%xmm2
jmp    131d <ynnm+0x20d>
movsd  0xe8a(%rip),%xmm2
jmp    12bd <ynnm+0x1ad>
mov    %r14d,%r12d
and    $0xfffffffe,%r12d
mov    $0x1,%eax
movsd  0xe71(%rip),%xmm2
movsd  0xe71(%rip),%xmm3
movsd  0xe71(%rip),%xmm4
xorpd  %xmm5,%xmm5
jmp    11ee <ynnm+0xde>
cs nopw 0x0(%rax,%rax,1)
nop
movsd  %xmm6,0x8(%rsp)
call   1030 <sqrt@plt>
movsd  0x8(%rsp),%xmm6
xorpd  %xmm5,%xmm5
movsd  0xe43(%rip),%xmm4
movsd  0xe33(%rip),%xmm3
movsd  (%rsp),%xmm2
mulsd  %xmm6,%xmm2
mulsd  %xmm0,%xmm2
lea    0x1(%rbx),%eax
cmp    %r12d,%ebx
je     126e <ynnm+0x15e>
mov    %eax,%ebp
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
movapd %xmm1,%xmm0
addsd  %xmm1,%xmm0
addsd  %xmm3,%xmm0
mulsd  %xmm4,%xmm1
divsd  %xmm1,%xmm0
ucomisd %xmm5,%xmm0
movsd  %xmm2,(%rsp)
jb     1220 <ynnm+0x110>
xorps  %xmm6,%xmm6
sqrtsd %xmm0,%xmm6
jmp    123d <ynnm+0x12d>
nop
call   1030 <sqrt@plt>
xorpd  %xmm5,%xmm5
movsd  0xddf(%rip),%xmm4
movsd  0xdcf(%rip),%xmm3
movapd %xmm0,%xmm6
lea    0x1(%rbp),%ebx
xorps  %xmm1,%xmm1
cvtsi2sd %ebx,%xmm1
movapd %xmm1,%xmm0
addsd  %xmm1,%xmm0
addsd  %xmm3,%xmm0
mulsd  %xmm4,%xmm1
divsd  %xmm1,%xmm0
ucomisd %xmm5,%xmm0
jb     11b0 <ynnm+0xa0>
sqrtsd %xmm0,%xmm0
jmp    11d5 <ynnm+0xc5>
add    $0x2,%ebp
xorps  %xmm1,%xmm1
cvtsi2sd %ebp,%xmm1
test   $0x1,%r14b
je     12bd <ynnm+0x1ad>
movapd %xmm1,%xmm0
addsd  %xmm1,%xmm0
addsd  0xd7a(%rip),%xmm0
mulsd  0xd7a(%rip),%xmm1
divsd  %xmm1,%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     12aa <ynnm+0x19a>
sqrtsd %xmm0,%xmm0
jmp    12b9 <ynnm+0x1a9>
movsd  %xmm2,(%rsp)
call   1030 <sqrt@plt>
movsd  (%rsp),%xmm2
mulsd  %xmm0,%xmm2
cmp    %r14d,%r15d
jge    131d <ynnm+0x20d>
mov    %r14d,%eax
sub    %r15d,%eax
mov    %r15d,%ebp
not    %ebp
test   $0x1,%al
jne    12db <ynnm+0x1cb>
mov    %r14d,%ebx
add    %r14d,%ebp
je     131d <ynnm+0x20d>
jmp    132e <ynnm+0x21e>
lea    (%r14,%r14,1),%eax
add    $0xffffffff,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
addsd  0xd17(%rip),%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     1301 <ynnm+0x1f1>
sqrtsd %xmm0,%xmm0
jmp    1310 <ynnm+0x200>
movsd  %xmm2,(%rsp)
call   1030 <sqrt@plt>
movsd  (%rsp),%xmm2
lea    -0x1(%r14),%ebx
mulsd  %xmm0,%xmm2
add    %r14d,%ebp
jne    132e <ynnm+0x21e>
movapd %xmm2,%xmm0
add    $0x10,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret
mov    %r14d,%ebp
sub    %ebx,%ebp
add    $0x2,%ebp
movsd  0xcca(%rip),%xmm3
xorpd  %xmm4,%xmm4
jmp    1382 <ynnm+0x272>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  %xmm5,0x8(%rsp)
call   1030 <sqrt@plt>
movsd  0x8(%rsp),%xmm5
xorpd  %xmm4,%xmm4
movsd  0xc9b(%rip),%xmm3
movsd  (%rsp),%xmm2
mulsd  %xmm5,%xmm2
mulsd  %xmm0,%xmm2
add    $0x2,%ebp
cmp    %r15d,%ebx
jle    131d <ynnm+0x20d>
lea    (%r14,%rbx,1),%eax
add    $0xffffffff,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
lea    -0x1(%rbp),%eax
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
addsd  %xmm3,%xmm0
divsd  %xmm1,%xmm0
ucomisd %xmm4,%xmm0
movsd  %xmm2,(%rsp)
jb     13c0 <ynnm+0x2b0>
xorps  %xmm5,%xmm5
sqrtsd %xmm0,%xmm5
jmp    13d5 <ynnm+0x2c5>
cs nopw 0x0(%rax,%rax,1)
call   1030 <sqrt@plt>
xorpd  %xmm4,%xmm4
movsd  0xc37(%rip),%xmm3
movapd %xmm0,%xmm5
mov    %ebx,%eax
add    $0xfffffffe,%ebx
add    %r14d,%eax
add    $0xfffffffe,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
xorps  %xmm1,%xmm1
cvtsi2sd %ebp,%xmm1
addsd  %xmm3,%xmm0
divsd  %xmm1,%xmm0
ucomisd %xmm4,%xmm0
jb     1350 <ynnm+0x240>
sqrtsd %xmm0,%xmm0
jmp    136d <ynnm+0x25d>

LLVM IR (-O0):
define double @ynnm(i32 %0, i32 %1) {
%3 = alloca i32
%4 = alloca i32
%5 = alloca i32
%6 = alloca i32
%7 = alloca double
%8 = alloca double
store i32 %0, i32* %3
store i32 %1, i32* %4
%9 = call double @sqrt(double 0x402921FB54442D18)
%10 = fdiv double 1.000000e+00, %9
store double %10, double* %7
%11 = load i32, i32* %4
%12 = call i32 @abs(i32 %11)
store i32 %12, i32* %5
%13 = load i32, i32* %3
%14 = load i32, i32* %5
%15 = icmp slt i32 %13, %14
br i1 %15, label %16, label %17
16: ; preds = %2
store double 0.000000e+00, double* %8
br label %74
17: ; preds = %2
%18 = load i32, i32* %3
%19 = icmp eq i32 %18, 0
br i1 %19, label %20, label %22
20: ; preds = %17
%21 = load double, double* %7
store double %21, double* %8
br label %73
22: ; preds = %17
%23 = load double, double* %7
store double %23, double* %8
store i32 1, i32* %6
br label %24
24: ; preds = %39, %22
%25 = load i32, i32* %6
%26 = load i32, i32* %3
%27 = icmp sle i32 %25, %26
br i1 %27, label %28, label %42
28: ; preds = %24
%29 = load i32, i32* %6
%30 = sitofp i32 %29 to double
%31 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %30, double 1.000000e+00)
%32 = load i32, i32* %6
%33 = sitofp i32 %32 to double
%34 = fmul double 8.000000e+00, %33
%35 = fdiv double %31, %34
%36 = call double @sqrt(double %35)
%37 = load double, double* %8
%38 = fmul double %36, %37
store double %38, double* %8
br label %39
39: ; preds = %28
%40 = load i32, i32* %6
%41 = add nsw i32 %40, 1
store i32 %41, i32* %6
br label %24
42: ; preds = %24
%43 = load i32, i32* %3
%44 = load i32, i32* %5
%45 = icmp ne i32 %43, %44
br i1 %45, label %46, label %72
46: ; preds = %42
%47 = load i32, i32* %3
%48 = sub nsw i32 %47, 1
store i32 %48, i32* %6
br label %49
49: ; preds = %68, %46
%50 = load i32, i32* %6
%51 = load i32, i32* %5
%52 = icmp sge i32 %50, %51
br i1 %52, label %53, label %71
53: ; preds = %49
%54 = load i32, i32* %3
%55 = load i32, i32* %6
%56 = add nsw i32 %54, %55
%57 = sitofp i32 %56 to double
%58 = fadd double %57, 1.000000e+00
%59 = fmul double 1.000000e+00, %58
%60 = load i32, i32* %3
%61 = load i32, i32* %6
%62 = sub nsw i32 %60, %61
%63 = sitofp i32 %62 to double
%64 = fdiv double %59, %63
%65 = call double @sqrt(double %64)
%66 = load double, double* %8
%67 = fmul double %65, %66
store double %67, double* %8
br label %68
68: ; preds = %53
%69 = load i32, i32* %6
%70 = add nsw i32 %69, -1
store i32 %70, i32* %6
br label %49
71: ; preds = %49
br label %72
72: ; preds = %71, %42
br label %73
73: ; preds = %72, %20
br label %74
74: ; preds = %73, %16
%75 = load double, double* %8
ret double %75
}

LLVM IR (-O2):
define double @ynnm(i32 %0, i32 %1) {
%3 = tail call i32 @llvm.abs.i32(i32 %1, i1 true)
%4 = icmp sgt i32 %3, %0
br i1 %4, label %90, label %5
5: ; preds = %2
%6 = icmp eq i32 %0, 0
br i1 %6, label %90, label %7
7: ; preds = %5
%8 = icmp slt i32 %0, 1
br i1 %8, label %47, label %9
9: ; preds = %7
%10 = and i32 %0, 1
%11 = icmp eq i32 %0, 1
br i1 %11, label %36, label %12
12: ; preds = %9
%13 = and i32 %0, -2
br label %14
14: ; preds = %14, %12
%15 = phi double [ 0x3FD20DD750429B6D, %12 ], [ %30, %14 ]
%16 = phi i32 [ 1, %12 ], [ %31, %14 ]
%17 = phi i32 [ 0, %12 ], [ %32, %14 ]
%18 = sitofp i32 %16 to double
%19 = tail call double @llvm.fmuladd.f64(double %18, double 2.000000e+00, double 1.000000e+00)
%20 = fmul double %18, 8.000000e+00
%21 = fdiv double %19, %20
%22 = tail call double @sqrt(double %21)
%23 = fmul double %15, %22
%24 = add nuw i32 %16, 1
%25 = sitofp i32 %24 to double
%26 = tail call double @llvm.fmuladd.f64(double %25, double 2.000000e+00, double 1.000000e+00)
%27 = fmul double %25, 8.000000e+00
%28 = fdiv double %26, %27
%29 = tail call double @sqrt(double %28)
%30 = fmul double %23, %29
%31 = add nuw i32 %16, 2
%32 = add i32 %17, 2
%33 = icmp eq i32 %32, %13
br i1 %33, label %34, label %14
34: ; preds = %14
%35 = sitofp i32 %31 to double
br label %36
36: ; preds = %34, %9
%37 = phi double [ undef, %9 ], [ %30, %34 ]
%38 = phi double [ 0x3FD20DD750429B6D, %9 ], [ %30, %34 ]
%39 = phi double [ 1.000000e+00, %9 ], [ %35, %34 ]
%40 = icmp eq i32 %10, 0
br i1 %40, label %47, label %41
41: ; preds = %36
%42 = tail call double @llvm.fmuladd.f64(double %39, double 2.000000e+00, double 1.000000e+00)
%43 = fmul double %39, 8.000000e+00
%44 = fdiv double %42, %43
%45 = tail call double @sqrt(double %44)
%46 = fmul double %38, %45
br label %47
47: ; preds = %41, %36, %7
%48 = phi double [ 0x3FD20DD750429B6D, %7 ], [ %37, %36 ], [ %46, %41 ]
%49 = icmp slt i32 %3, %0
br i1 %49, label %50, label %90
50: ; preds = %47
%51 = sub i32 %0, %3
%52 = xor i32 %3, -1
%53 = and i32 %51, 1
%54 = icmp eq i32 %53, 0
br i1 %54, label %62, label %55
55: ; preds = %50
%56 = add nsw i32 %0, -1
%57 = add nsw i32 %56, %0
%58 = sitofp i32 %57 to double
%59 = fadd double %58, 1.000000e+00
%60 = tail call double @sqrt(double %59)
%61 = fmul double %48, %60
br label %62
62: ; preds = %55, %50
%63 = phi i32 [ %0, %50 ], [ %56, %55 ]
%64 = phi double [ %48, %50 ], [ %61, %55 ]
%65 = phi double [ undef, %50 ], [ %61, %55 ]
%66 = sub i32 0, %0
%67 = icmp eq i32 %52, %66
br i1 %67, label %90, label %68
68: ; preds = %62, %68
%69 = phi i32 [ %80, %68 ], [ %63, %62 ]
%70 = phi double [ %88, %68 ], [ %64, %62 ]
%71 = add nsw i32 %69, -1
%72 = add nsw i32 %71, %0
%73 = sitofp i32 %72 to double
%74 = fadd double %73, 1.000000e+00
%75 = sub nsw i32 %0, %71
%76 = sitofp i32 %75 to double
%77 = fdiv double %74, %76
%78 = tail call double @sqrt(double %77)
%79 = fmul double %70, %78
%80 = add nsw i32 %69, -2
%81 = add nsw i32 %80, %0
%82 = sitofp i32 %81 to double
%83 = fadd double %82, 1.000000e+00
%84 = sub nsw i32 %0, %80
%85 = sitofp i32 %84 to double
%86 = fdiv double %83, %85
%87 = tail call double @sqrt(double %86)
%88 = fmul double %79, %87
%89 = icmp sgt i32 %80, %3
br i1 %89, label %68, label %90
90: ; preds = %62, %68, %5, %2, %47
%91 = phi double [ %48, %47 ], [ 0.000000e+00, %2 ], [ 0x3FD20DD750429B6D, %5 ], [ %65, %62 ], [ %88, %68 ]
ret double %91
}