void plgndr0(int n, double x, double *pln) {
  int i;
  double p0, p1, p2;
  pln[0] = p0 = 1.0;
  pln[1] = p1 = x;
  if (n < 3) {
    return;
  }
  for (i = 2; i < n; i++) {
    p2 = (x * (2 * i - 1) * p1 - (i - 1) * p0) / (double)i;
    p0 = p1;
    p1 = p2;
    pln[i] = p2;
  }
}


x86 assembly (clang, -O0):
<plgndr0>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
movsd  0xee8(%rip),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x18(%rbp),%rax
movsd  0xed7(%rip),%xmm0
movsd  %xmm0,(%rax)
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x18(%rbp),%rax
movsd  %xmm0,0x8(%rax)
cmpl   $0x3,-0x4(%rbp)
jge    114f <plgndr0+0x4f>
jmp    11ed <plgndr0+0xed>
movl   $0x2,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11ed <plgndr0+0xed>
movsd  -0x10(%rbp),%xmm0
mov    -0x1c(%rbp),%eax
shl    $0x1,%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x30(%rbp),%xmm2
mov    -0x1c(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
mulsd  -0x28(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsi2sdl -0x1c(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x38(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1156 <plgndr0+0x56>
pop    %rbp
ret

x86 assembly (clang, -O2):
<plgndr0>:
movabs $0x3ff0000000000000,%rax
mov    %rax,(%rsi)
movsd  %xmm0,0x8(%rsi)
cmp    $0x3,%edi
jl     1181 <plgndr0+0x81>
mov    %edi,%eax
add    %rax,%rax
movsd  0xedc(%rip),%xmm2
mov    $0x1,%ecx
mov    $0x4,%edx
movapd %xmm0,%xmm1
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
xorps  %xmm3,%xmm3
cvtsi2sd %ecx,%xmm3
mulsd  %xmm2,%xmm3
movapd %xmm1,%xmm2
lea    -0x1(%rdx),%edi
xorps  %xmm1,%xmm1
cvtsi2sd %edi,%xmm1
mulsd  %xmm0,%xmm1
mulsd  %xmm2,%xmm1
subsd  %xmm3,%xmm1
add    $0x1,%ecx
xorps  %xmm3,%xmm3
cvtsi2sd %ecx,%xmm3
divsd  %xmm3,%xmm1
movsd  %xmm1,(%rsi,%rdx,4)
add    $0x2,%rdx
cmp    %rdx,%rax
jne    1140 <plgndr0+0x40>
ret

LLVM IR (-O0):
define void @plgndr0(i32 %0, double %1, double* %2) {
%4 = alloca i32
%5 = alloca double
%6 = alloca double*
%7 = alloca i32
%8 = alloca double
%9 = alloca double
%10 = alloca double
store i32 %0, i32* %4
store double %1, double* %5
store double* %2, double** %6
store double 1.000000e+00, double* %8
%11 = load double*, double** %6
%12 = getelementptr inbounds double, double* %11, i64 0
store double 1.000000e+00, double* %12
%13 = load double, double* %5
store double %13, double* %9
%14 = load double*, double** %6
%15 = getelementptr inbounds double, double* %14, i64 1
store double %13, double* %15
%16 = load i32, i32* %4
%17 = icmp slt i32 %16, 3
br i1 %17, label %18, label %19
18: ; preds = %3
br label %52
19: ; preds = %3
store i32 2, i32* %7
br label %20
20: ; preds = %49, %19
%21 = load i32, i32* %7
%22 = load i32, i32* %4
%23 = icmp slt i32 %21, %22
br i1 %23, label %24, label %52
24: ; preds = %20
%25 = load double, double* %5
%26 = load i32, i32* %7
%27 = mul nsw i32 2, %26
%28 = sub nsw i32 %27, 1
%29 = sitofp i32 %28 to double
%30 = fmul double %25, %29
%31 = load double, double* %9
%32 = load i32, i32* %7
%33 = sub nsw i32 %32, 1
%34 = sitofp i32 %33 to double
%35 = load double, double* %8
%36 = fmul double %34, %35
%37 = fneg double %36
%38 = call double @llvm.fmuladd.f64(double %30, double %31, double %37)
%39 = load i32, i32* %7
%40 = sitofp i32 %39 to double
%41 = fdiv double %38, %40
store double %41, double* %10
%42 = load double, double* %9
store double %42, double* %8
%43 = load double, double* %10
store double %43, double* %9
%44 = load double, double* %10
%45 = load double*, double** %6
%46 = load i32, i32* %7
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds double, double* %45, i64 %47
store double %44, double* %48
br label %49
49: ; preds = %24
%50 = load i32, i32* %7
%51 = add nsw i32 %50, 1
store i32 %51, i32* %7
br label %20
52: ; preds = %18, %20
ret void
}

LLVM IR (-O2):
define void @plgndr0(i32 %0, double %1, double* nocapture writeonly %2) {
store double 1.000000e+00, double* %2
%4 = getelementptr inbounds double, double* %2, i64 1
store double %1, double* %4
%5 = icmp sgt i32 %0, 2
br i1 %5, label %6, label %29
6: ; preds = %3
%7 = zext i32 %0 to i64
br label %8
8: ; preds = %6, %8
%9 = phi i64 [ 2, %6 ], [ %27, %8 ]
%10 = phi double [ %1, %6 ], [ %25, %8 ]
%11 = phi double [ 1.000000e+00, %6 ], [ %10, %8 ]
%12 = trunc i64 %9 to i32
%13 = shl i32 %12, 1
%14 = add i32 %13, -1
%15 = sitofp i32 %14 to double
%16 = fmul double %15, %1
%17 = trunc i64 %9 to i32
%18 = add i32 %17, -1
%19 = sitofp i32 %18 to double
%20 = fneg double %11
%21 = fmul double %20, %19
%22 = tail call double @llvm.fmuladd.f64(double %16, double %10, double %21)
%23 = trunc i64 %9 to i32
%24 = sitofp i32 %23 to double
%25 = fdiv double %22, %24
%26 = getelementptr inbounds double, double* %2, i64 %9
store double %25, double* %26
%27 = add nuw nsw i64 %9, 1
%28 = icmp eq i64 %27, %7
br i1 %28, label %29, label %8
29: ; preds = %8, %3
ret void
}