int leg_pol_berechnen(double t, int maxdeg, double *lp) {
  short n;
  lp[0] = 1.0;
  lp[1] = t;
  for (n = 2; n <= maxdeg; n++) {
    lp[n] = (-lp[n - 2] * (n - 1) + t * lp[n - 1] * (2 * n - 1)) / n;
  }
  return 0;
}


x86 assembly (clang, -O0):
<leg_pol_berechnen>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x8(%rbp)
mov    %edi,-0xc(%rbp)
mov    %rsi,-0x18(%rbp)
mov    -0x18(%rbp),%rax
movsd  0xee4(%rip),%xmm0
movsd  %xmm0,(%rax)
movsd  -0x8(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,0x8(%rax)
movw   $0x2,-0x1a(%rbp)
movswl -0x1a(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     11d3 <leg_pol_berechnen+0xd3>
mov    -0x18(%rbp),%rax
movswl -0x1a(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movswl -0x1a(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm2
movsd  -0x8(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movswl -0x1a(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
movswl -0x1a(%rbp),%eax
shl    $0x1,%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm3
mulsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movswl -0x1a(%rbp),%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movswq -0x1a(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1a(%rbp),%ax
add    $0x1,%ax
mov    %ax,-0x1a(%rbp)
jmp    1134 <leg_pol_berechnen+0x34>
xor    %eax,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<leg_pol_berechnen>:
movabs $0x3ff0000000000000,%rax
mov    %rax,(%rsi)
movsd  %xmm0,0x8(%rsi)
cmp    $0x2,%edi
jl     116a <leg_pol_berechnen+0x6a>
mov    $0x2,%eax
nopl   0x0(%rax)
lea    -0x1(%rax),%ecx
xorps  %xmm1,%xmm1
cvtsi2sd %ecx,%xmm1
cltq
movsd  -0x8(%rsi,%rax,8),%xmm2
mulsd  %xmm0,%xmm2
lea    -0x1(%rax,%rax,1),%ecx
xorps  %xmm3,%xmm3
cvtsi2sd %ecx,%xmm3
mulsd  %xmm2,%xmm3
mulsd  -0x10(%rsi,%rax,8),%xmm1
subsd  %xmm1,%xmm3
movswl %ax,%ecx
xorps  %xmm1,%xmm1
cvtsi2sd %ecx,%xmm1
divsd  %xmm1,%xmm3
movsd  %xmm3,(%rsi,%rax,8)
add    $0x1,%eax
cwtl
cmp    %edi,%eax
jle    1120 <leg_pol_berechnen+0x20>
xor    %eax,%eax
ret

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

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