double plgndr_a(l, m, t)
int l, m;
double t;
{
  int i, ll;
  double x, fact, pll, pmm, pmmp1, somx2, plgndr;
  const double two = 2.0;
  const double one = 1.0;
  x = cos(t);
  pmm = one;
  if (m > 0) {
    somx2 = sqrt((one - x) * (one + x));
    fact = one;
    for (i = 1; i <= m; i++) {
      pmm = -pmm * fact * somx2;
      fact = fact + two;
    }
  }
  if (l == m)
    plgndr = pmm;
  else {
    pmmp1 = x * (2 * m + 1) * pmm;
    if (l == m + 1)
      plgndr = pmmp1;
    else {
      for (ll = m + 2; ll <= l; ll++) {
        pll = (x * (2 * ll - 1) * pmmp1 - (ll + m - 1) * pmm) / (ll - m);
        pmm = pmmp1;
        pmmp1 = pll;
      }
      plgndr = pll;
    }
  }
  return plgndr;
}


x86 assembly (clang, -O0):
<plgndr_a>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
movsd  %xmm0,-0x10(%rbp)
movsd  0xecd(%rip),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  0xeb8(%rip),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x10(%rbp),%xmm0
call   1030 <cos@plt>
movsd  %xmm0,-0x20(%rbp)
movsd  0xe9c(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
cmpl   $0x0,-0x8(%rbp)
jle    120b <plgndr_a+0xeb>
movsd  0xe85(%rip),%xmm0
subsd  -0x20(%rbp),%xmm0
movsd  0xe78(%rip),%xmm1
addsd  -0x20(%rbp),%xmm1
mulsd  %xmm1,%xmm0
call   1040 <sqrt@plt>
movsd  %xmm0,-0x48(%rbp)
movsd  0xe5d(%rip),%xmm0
movsd  %xmm0,-0x28(%rbp)
movl   $0x1,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0x8(%rbp),%eax
jg     1206 <plgndr_a+0xe6>
movsd  -0x38(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  -0x28(%rbp),%xmm0
mulsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  0xe1a(%rip),%xmm0
addsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    11af <plgndr_a+0x8f>
jmp    120b <plgndr_a+0xeb>
mov    -0x4(%rbp),%eax
cmp    -0x8(%rbp),%eax
jne    1226 <plgndr_a+0x106>
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
jmp    130b <plgndr_a+0x1eb>
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%eax
shl    $0x1,%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm1,%xmm0
mulsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x4(%rbp),%eax
mov    -0x8(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jne    1266 <plgndr_a+0x146>
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
jmp    1306 <plgndr_a+0x1e6>
mov    -0x8(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     12fc <plgndr_a+0x1dc>
movsd  -0x20(%rbp),%xmm0
mov    -0x18(%rbp),%eax
shl    $0x1,%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x40(%rbp),%xmm2
mov    -0x18(%rbp),%eax
add    -0x8(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
mulsd  -0x38(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%eax
sub    -0x8(%rbp),%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    126f <plgndr_a+0x14f>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
jmp    130b <plgndr_a+0x1eb>
movsd  -0x50(%rbp),%xmm0
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<plgndr_a>:
push   %rbp
push   %rbx
push   %rax
mov    %esi,%ebx
mov    %edi,%ebp
call   1030 <cos@plt>
movapd %xmm0,%xmm5
test   %ebx,%ebx
jle    1160 <plgndr_a+0x40>
movsd  0xec4(%rip),%xmm6
movapd %xmm6,%xmm1
subsd  %xmm5,%xmm1
movapd %xmm5,%xmm0
addsd  %xmm6,%xmm0
mulsd  %xmm1,%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     116d <plgndr_a+0x4d>
sqrtsd %xmm0,%xmm0
jmp    1184 <plgndr_a+0x64>
movsd  0xe98(%rip),%xmm1
jmp    1235 <plgndr_a+0x115>
movsd  %xmm5,(%rsp)
call   1040 <sqrt@plt>
movsd  0xe81(%rip),%xmm6
movsd  (%rsp),%xmm5
lea    -0x1(%rbx),%ecx
mov    %ebx,%eax
and    $0x3,%eax
movapd %xmm6,%xmm1
cmp    $0x3,%ecx
jb     1209 <plgndr_a+0xe9>
mov    %ebx,%ecx
and    $0xfffffffc,%ecx
movapd 0xe6e(%rip),%xmm2
xorpd  %xmm0,%xmm2
movsd  0xe52(%rip),%xmm6
movapd 0xe5a(%rip),%xmm3
movsd  0xe4a(%rip),%xmm4
movapd %xmm6,%xmm1
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
xorpd  %xmm3,%xmm1
mulsd  %xmm6,%xmm1
addsd  %xmm4,%xmm6
mulsd  %xmm2,%xmm1
mulsd  %xmm6,%xmm1
addsd  %xmm4,%xmm6
mulsd  %xmm2,%xmm1
mulsd  %xmm6,%xmm1
addsd  %xmm4,%xmm6
mulsd  %xmm2,%xmm1
mulsd  %xmm6,%xmm1
mulsd  %xmm0,%xmm1
addsd  %xmm4,%xmm6
add    $0xfffffffc,%ecx
jne    11d0 <plgndr_a+0xb0>
test   %eax,%eax
je     1235 <plgndr_a+0x115>
movapd 0xdfb(%rip),%xmm2
movsd  0xdeb(%rip),%xmm3
nopl   (%rax)
xorpd  %xmm2,%xmm1
mulsd  %xmm6,%xmm1
mulsd  %xmm0,%xmm1
addsd  %xmm3,%xmm6
add    $0xffffffff,%eax
jne    1220 <plgndr_a+0x100>
mov    %ebp,%eax
sub    %ebx,%eax
jne    1241 <plgndr_a+0x121>
movapd %xmm1,%xmm0
jmp    12ba <plgndr_a+0x19a>
lea    (%rbx,%rbx,1),%ecx
add    $0x1,%ecx
xorps  %xmm0,%xmm0
cvtsi2sd %ecx,%xmm0
mulsd  %xmm5,%xmm0
mulsd  %xmm1,%xmm0
lea    0x1(%rbx),%ecx
cmp    %ebp,%ecx
je     12ba <plgndr_a+0x19a>
lea    0x2(%rbx),%ecx
cmp    %ebp,%ecx
jg     12ba <plgndr_a+0x19a>
mov    %ebx,%ecx
add    %ebx,%ebx
add    %ecx,%ecx
add    $0x3,%ecx
mov    $0x1,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
lea    (%rbx,%rdx,1),%esi
xorps  %xmm2,%xmm2
cvtsi2sd %esi,%xmm2
mulsd  %xmm1,%xmm2
movapd %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsi2sd %ecx,%xmm0
mulsd  %xmm5,%xmm0
mulsd  %xmm1,%xmm0
add    $0x1,%edx
xorps  %xmm3,%xmm3
cvtsi2sd %edx,%xmm3
subsd  %xmm2,%xmm0
divsd  %xmm3,%xmm0
add    $0x2,%ecx
cmp    %edx,%eax
jne    1280 <plgndr_a+0x160>
add    $0x8,%rsp
pop    %rbx
pop    %rbp
ret

LLVM IR (-O0):
define double @plgndr_a(i32 %0, i32 %1, double %2) {
%4 = alloca i32
%5 = alloca i32
%6 = alloca double
%7 = alloca i32
%8 = alloca i32
%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 i32 %0, i32* %4
store i32 %1, i32* %5
store double %2, double* %6
store double 2.000000e+00, double* %16
store double 1.000000e+00, double* %17
%18 = load double, double* %6
%19 = call double @cos(double %18)
store double %19, double* %9
store double 1.000000e+00, double* %12
%20 = load i32, i32* %5
%21 = icmp sgt i32 %20, 0
br i1 %21, label %22, label %46
22: ; preds = %3
%23 = load double, double* %9
%24 = fsub double 1.000000e+00, %23
%25 = load double, double* %9
%26 = fadd double 1.000000e+00, %25
%27 = fmul double %24, %26
%28 = call double @sqrt(double %27)
store double %28, double* %14
store double 1.000000e+00, double* %10
store i32 1, i32* %7
br label %29
29: ; preds = %42, %22
%30 = load i32, i32* %7
%31 = load i32, i32* %5
%32 = icmp sle i32 %30, %31
br i1 %32, label %33, label %45
33: ; preds = %29
%34 = load double, double* %12
%35 = fneg double %34
%36 = load double, double* %10
%37 = fmul double %35, %36
%38 = load double, double* %14
%39 = fmul double %37, %38
store double %39, double* %12
%40 = load double, double* %10
%41 = fadd double %40, 2.000000e+00
store double %41, double* %10
br label %42
42: ; preds = %33
%43 = load i32, i32* %7
%44 = add nsw i32 %43, 1
store i32 %44, i32* %7
br label %29
45: ; preds = %29
br label %46
46: ; preds = %45, %3
%47 = load i32, i32* %4
%48 = load i32, i32* %5
%49 = icmp eq i32 %47, %48
br i1 %49, label %50, label %52
50: ; preds = %46
%51 = load double, double* %12
store double %51, double* %15
br label %104
52: ; preds = %46
%53 = load double, double* %9
%54 = load i32, i32* %5
%55 = mul nsw i32 2, %54
%56 = add nsw i32 %55, 1
%57 = sitofp i32 %56 to double
%58 = fmul double %53, %57
%59 = load double, double* %12
%60 = fmul double %58, %59
store double %60, double* %13
%61 = load i32, i32* %4
%62 = load i32, i32* %5
%63 = add nsw i32 %62, 1
%64 = icmp eq i32 %61, %63
br i1 %64, label %65, label %67
65: ; preds = %52
%66 = load double, double* %13
store double %66, double* %15
br label %103
67: ; preds = %52
%68 = load i32, i32* %5
%69 = add nsw i32 %68, 2
store i32 %69, i32* %8
br label %70
70: ; preds = %98, %67
%71 = load i32, i32* %8
%72 = load i32, i32* %4
%73 = icmp sle i32 %71, %72
br i1 %73, label %74, label %101
74: ; preds = %70
%75 = load double, double* %9
%76 = load i32, i32* %8
%77 = mul nsw i32 2, %76
%78 = sub nsw i32 %77, 1
%79 = sitofp i32 %78 to double
%80 = fmul double %75, %79
%81 = load double, double* %13
%82 = load i32, i32* %8
%83 = load i32, i32* %5
%84 = add nsw i32 %82, %83
%85 = sub nsw i32 %84, 1
%86 = sitofp i32 %85 to double
%87 = load double, double* %12
%88 = fmul double %86, %87
%89 = fneg double %88
%90 = call double @llvm.fmuladd.f64(double %80, double %81, double %89)
%91 = load i32, i32* %8
%92 = load i32, i32* %5
%93 = sub nsw i32 %91, %92
%94 = sitofp i32 %93 to double
%95 = fdiv double %90, %94
store double %95, double* %11
%96 = load double, double* %13
store double %96, double* %12
%97 = load double, double* %11
store double %97, double* %13
br label %98
98: ; preds = %74
%99 = load i32, i32* %8
%100 = add nsw i32 %99, 1
store i32 %100, i32* %8
br label %70
101: ; preds = %70
%102 = load double, double* %11
store double %102, double* %15
br label %103
103: ; preds = %101, %65
br label %104
104: ; preds = %103, %50
%105 = load double, double* %15
ret double %105
}

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