void wf_gsb1pol(x, order, basis) double x;
int order;
double *basis;
{
  int i;
  basis[0] = 1.0;
  if (order == 1)
    return;
  basis[1] = x;
  if (order == 2)
    return;
  for (i = 2; i < order; i++)
    basis[i] = x * basis[i - 1];
  return;
}


x86 assembly (clang, -O0):
<wf_gsb1pol>:
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)
cmpl   $0x1,-0xc(%rbp)
jne    112f <wf_gsb1pol+0x2f>
jmp    1196 <wf_gsb1pol+0x96>
movsd  -0x8(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,0x8(%rax)
cmpl   $0x2,-0xc(%rbp)
jne    114c <wf_gsb1pol+0x4c>
jmp    1196 <wf_gsb1pol+0x96>
movl   $0x2,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1191 <wf_gsb1pol+0x91>
movsd  -0x8(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%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    1153 <wf_gsb1pol+0x53>
jmp    1196 <wf_gsb1pol+0x96>
pop    %rbp
ret

x86 assembly (clang, -O2):
<wf_gsb1pol>:
movabs $0x3ff0000000000000,%rax
mov    %rax,(%rsi)
cmp    $0x1,%edi
je     11c2 <wf_gsb1pol+0xc2>
movsd  %xmm0,0x8(%rsi)
cmp    $0x3,%edi
jl     11c2 <wf_gsb1pol+0xc2>
mov    %edi,%edi
lea    -0x2(%rdi),%rcx
add    $0xfffffffffffffffd,%rdi
mov    %ecx,%eax
and    $0x3,%eax
mov    $0x2,%edx
movapd %xmm0,%xmm1
cmp    $0x3,%rdi
jb     1199 <wf_gsb1pol+0x99>
and    $0xfffffffffffffffc,%rcx
neg    %rcx
mov    $0x2,%edx
movapd %xmm0,%xmm1
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rsi,%rdx,8)
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rsi,%rdx,8)
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x10(%rsi,%rdx,8)
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x18(%rsi,%rdx,8)
lea    (%rcx,%rdx,1),%rdi
add    $0x4,%rdi
add    $0x4,%rdx
cmp    $0x2,%rdi
jne    1160 <wf_gsb1pol+0x60>
test   %rax,%rax
je     11c2 <wf_gsb1pol+0xc2>
lea    (%rsi,%rdx,8),%rcx
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rcx,%rdx,8)
add    $0x1,%rdx
cmp    %rdx,%rax
jne    11b0 <wf_gsb1pol+0xb0>
ret

LLVM IR (-O0):
define void @wf_gsb1pol(double %0, i32 %1, double* %2) {
%4 = alloca double
%5 = alloca i32
%6 = alloca double*
%7 = alloca i32
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 i32, i32* %5
%11 = icmp eq i32 %10, 1
br i1 %11, label %12, label %13
12: ; preds = %3
br label %42
13: ; preds = %3
%14 = load double, double* %4
%15 = load double*, double** %6
%16 = getelementptr inbounds double, double* %15, i64 1
store double %14, double* %16
%17 = load i32, i32* %5
%18 = icmp eq i32 %17, 2
br i1 %18, label %19, label %20
19: ; preds = %13
br label %42
20: ; preds = %13
store i32 2, i32* %7
br label %21
21: ; preds = %38, %20
%22 = load i32, i32* %7
%23 = load i32, i32* %5
%24 = icmp slt i32 %22, %23
br i1 %24, label %25, label %41
25: ; preds = %21
%26 = load double, double* %4
%27 = load double*, double** %6
%28 = load i32, i32* %7
%29 = sub nsw i32 %28, 1
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds double, double* %27, i64 %30
%32 = load double, double* %31
%33 = fmul double %26, %32
%34 = load double*, double** %6
%35 = load i32, i32* %7
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds double, double* %34, i64 %36
store double %33, double* %37
br label %38
38: ; preds = %25
%39 = load i32, i32* %7
%40 = add nsw i32 %39, 1
store i32 %40, i32* %7
br label %21
41: ; preds = %21
br label %42
42: ; preds = %41, %19, %12
ret void
}

LLVM IR (-O2):
define void @wf_gsb1pol(double %0, i32 %1, double* nocapture %2) {
store double 1.000000e+00, double* %2
%4 = icmp eq i32 %1, 1
br i1 %4, label %47, label %5
5: ; preds = %3
%6 = getelementptr inbounds double, double* %2, i64 1
store double %0, double* %6
%7 = icmp sgt i32 %1, 2
br i1 %7, label %8, label %47
8: ; preds = %5
%9 = zext i32 %1 to i64
%10 = add nsw i64 %9, -2
%11 = add nsw i64 %9, -3
%12 = and i64 %10, 3
%13 = icmp ult i64 %11, 3
br i1 %13, label %34, label %14
14: ; preds = %8
%15 = and i64 %10, -4
br label %16
16: ; preds = %16, %14
%17 = phi double [ %0, %14 ], [ %29, %16 ]
%18 = phi i64 [ 2, %14 ], [ %31, %16 ]
%19 = phi i64 [ 0, %14 ], [ %32, %16 ]
%20 = fmul double %17, %0
%21 = getelementptr inbounds double, double* %2, i64 %18
store double %20, double* %21
%22 = or i64 %18, 1
%23 = fmul double %20, %0
%24 = getelementptr inbounds double, double* %2, i64 %22
store double %23, double* %24
%25 = add nuw nsw i64 %18, 2
%26 = fmul double %23, %0
%27 = getelementptr inbounds double, double* %2, i64 %25
store double %26, double* %27
%28 = add nuw nsw i64 %18, 3
%29 = fmul double %26, %0
%30 = getelementptr inbounds double, double* %2, i64 %28
store double %29, double* %30
%31 = add nuw nsw i64 %18, 4
%32 = add i64 %19, 4
%33 = icmp eq i64 %32, %15
br i1 %33, label %34, label %16
34: ; preds = %16, %8
%35 = phi double [ %0, %8 ], [ %29, %16 ]
%36 = phi i64 [ 2, %8 ], [ %31, %16 ]
%37 = icmp eq i64 %12, 0
br i1 %37, label %47, label %38
38: ; preds = %34, %38
%39 = phi double [ %42, %38 ], [ %35, %34 ]
%40 = phi i64 [ %44, %38 ], [ %36, %34 ]
%41 = phi i64 [ %45, %38 ], [ 0, %34 ]
%42 = fmul double %39, %0
%43 = getelementptr inbounds double, double* %2, i64 %40
store double %42, double* %43
%44 = add nuw nsw i64 %40, 1
%45 = add i64 %41, 1
%46 = icmp eq i64 %45, %12
br i1 %46, label %47, label %38
47: ; preds = %34, %38, %5, %3
ret void
}