void rick_ab2cs(rdata, n) double *rdata;
int n;
{
  int i;
  double en;
  en = 1.0 / (double)n;
  rdata[0] *= en;
  rdata[1] = 0.0;
  en *= 2.0;
  for (i = 2; i < n; i++)
    rdata[i] *= en;
}


x86 assembly (clang, -O0):
<rick_ab2cs>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
cvtsi2sdl -0xc(%rbp),%xmm1
movsd  0xef0(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x18(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mulsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
mov    -0x8(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x8(%rax)
movsd  0xeba(%rip),%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
movl   $0x2,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1188 <rick_ab2cs+0x88>
movsd  -0x18(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1157 <rick_ab2cs+0x57>
pop    %rbp
ret

x86 assembly (clang, -O2):
<rick_ab2cs>:
cvtsi2sd %esi,%xmm1
movsd  0xef4(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  (%rdi),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rdi)
movq   $0x0,0x8(%rdi)
cmp    $0x3,%esi
jl     1217 <rick_ab2cs+0x117>
addsd  %xmm0,%xmm0
mov    %esi,%eax
lea    -0x2(%rax),%r8
mov    $0x2,%ecx
cmp    $0x4,%r8
jb     1200 <rick_ab2cs+0x100>
mov    %r8,%rcx
and    $0xfffffffffffffffc,%rcx
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x4(%rcx),%rdx
mov    %rdx,%r9
shr    $0x2,%r9
add    $0x1,%r9
test   %rdx,%rdx
je     1218 <rick_ab2cs+0x118>
mov    %r9,%rdx
and    $0xfffffffffffffffe,%rdx
mov    $0x2,%esi
nopl   0x0(%rax)
movupd (%rdi,%rsi,8),%xmm2
movupd 0x10(%rdi,%rsi,8),%xmm3
movupd 0x20(%rdi,%rsi,8),%xmm4
movupd 0x30(%rdi,%rsi,8),%xmm5
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,(%rdi,%rsi,8)
movupd %xmm3,0x10(%rdi,%rsi,8)
mulpd  %xmm1,%xmm4
mulpd  %xmm1,%xmm5
movupd %xmm4,0x20(%rdi,%rsi,8)
movupd %xmm5,0x30(%rdi,%rsi,8)
add    $0x8,%rsi
add    $0xfffffffffffffffe,%rdx
jne    1180 <rick_ab2cs+0x80>
test   $0x1,%r9b
je     11ec <rick_ab2cs+0xec>
movupd (%rdi,%rsi,8),%xmm2
movupd 0x10(%rdi,%rsi,8),%xmm3
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,(%rdi,%rsi,8)
movupd %xmm3,0x10(%rdi,%rsi,8)
cmp    %rcx,%r8
je     1217 <rick_ab2cs+0x117>
or     $0x2,%rcx
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rdi,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rdi,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    1200 <rick_ab2cs+0x100>
ret
mov    $0x2,%esi
test   $0x1,%r9b
jne    11ce <rick_ab2cs+0xce>
jmp    11ec <rick_ab2cs+0xec>

LLVM IR (-O0):
define void @rick_ab2cs(double* %0, i32 %1) {
%3 = alloca double*
%4 = alloca i32
%5 = alloca i32
%6 = alloca double
store double* %0, double** %3
store i32 %1, i32* %4
%7 = load i32, i32* %4
%8 = sitofp i32 %7 to double
%9 = fdiv double 1.000000e+00, %8
store double %9, double* %6
%10 = load double, double* %6
%11 = load double*, double** %3
%12 = getelementptr inbounds double, double* %11, i64 0
%13 = load double, double* %12
%14 = fmul double %13, %10
store double %14, double* %12
%15 = load double*, double** %3
%16 = getelementptr inbounds double, double* %15, i64 1
store double 0.000000e+00, double* %16
%17 = load double, double* %6
%18 = fmul double %17, 2.000000e+00
store double %18, double* %6
store i32 2, i32* %5
br label %19
19: ; preds = %31, %2
%20 = load i32, i32* %5
%21 = load i32, i32* %4
%22 = icmp slt i32 %20, %21
br i1 %22, label %23, label %34
23: ; preds = %19
%24 = load double, double* %6
%25 = load double*, double** %3
%26 = load i32, i32* %5
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds double, double* %25, i64 %27
%29 = load double, double* %28
%30 = fmul double %29, %24
store double %30, double* %28
br label %31
31: ; preds = %23
%32 = load i32, i32* %5
%33 = add nsw i32 %32, 1
store i32 %33, i32* %5
br label %19
34: ; preds = %19
ret void
}

LLVM IR (-O2):
define void @rick_ab2cs(double* nocapture %0, i32 %1) {
%3 = sitofp i32 %1 to double
%4 = fdiv double 1.000000e+00, %3
%5 = load double, double* %0
%6 = fmul double %4, %5
store double %6, double* %0
%7 = getelementptr inbounds double, double* %0, i64 1
store double 0.000000e+00, double* %7
%8 = fmul double %4, 2.000000e+00
%9 = icmp sgt i32 %1, 2
br i1 %9, label %10, label %83
10: ; preds = %2
%11 = zext i32 %1 to i64
%12 = add nsw i64 %11, -2
%13 = icmp ult i64 %12, 4
br i1 %13, label %74, label %14
14: ; preds = %10
%15 = and i64 %12, -4
%16 = or i64 %15, 2
%17 = insertelement <2 x double> poison, double %8, i64 0
%18 = shufflevector <2 x double> %17, <2 x double> poison, <2 x i32> zeroinitializer
%19 = insertelement <2 x double> poison, double %8, i64 0
%20 = shufflevector <2 x double> %19, <2 x double> poison, <2 x i32> zeroinitializer
%21 = add nsw i64 %15, -4
%22 = lshr exact i64 %21, 2
%23 = add nuw nsw i64 %22, 1
%24 = and i64 %23, 1
%25 = icmp eq i64 %21, 0
br i1 %25, label %58, label %26
26: ; preds = %14
%27 = and i64 %23, 9223372036854775806
br label %28
28: ; preds = %28, %26
%29 = phi i64 [ 0, %26 ], [ %53, %28 ]
%30 = phi i64 [ 0, %26 ], [ %54, %28 ]
%31 = or i64 %29, 2
%32 = getelementptr inbounds double, double* %0, i64 %31
%33 = bitcast double* %32 to <2 x double>*
%34 = load <2 x double>, <2 x double>* %33
%35 = getelementptr inbounds double, double* %32, i64 2
%36 = bitcast double* %35 to <2 x double>*
%37 = load <2 x double>, <2 x double>* %36
%38 = fmul <2 x double> %18, %34
%39 = fmul <2 x double> %20, %37
%40 = bitcast double* %32 to <2 x double>*
store <2 x double> %38, <2 x double>* %40
%41 = bitcast double* %35 to <2 x double>*
store <2 x double> %39, <2 x double>* %41
%42 = or i64 %29, 6
%43 = getelementptr inbounds double, double* %0, i64 %42
%44 = bitcast double* %43 to <2 x double>*
%45 = load <2 x double>, <2 x double>* %44
%46 = getelementptr inbounds double, double* %43, i64 2
%47 = bitcast double* %46 to <2 x double>*
%48 = load <2 x double>, <2 x double>* %47
%49 = fmul <2 x double> %18, %45
%50 = fmul <2 x double> %20, %48
%51 = bitcast double* %43 to <2 x double>*
store <2 x double> %49, <2 x double>* %51
%52 = bitcast double* %46 to <2 x double>*
store <2 x double> %50, <2 x double>* %52
%53 = add nuw i64 %29, 8
%54 = add i64 %30, 2
%55 = icmp eq i64 %54, %27
br i1 %55, label %56, label %28
56: ; preds = %28
%57 = or i64 %53, 2
br label %58
58: ; preds = %56, %14
%59 = phi i64 [ 2, %14 ], [ %57, %56 ]
%60 = icmp eq i64 %24, 0
br i1 %60, label %72, label %61
61: ; preds = %58
%62 = getelementptr inbounds double, double* %0, i64 %59
%63 = bitcast double* %62 to <2 x double>*
%64 = load <2 x double>, <2 x double>* %63
%65 = getelementptr inbounds double, double* %62, i64 2
%66 = bitcast double* %65 to <2 x double>*
%67 = load <2 x double>, <2 x double>* %66
%68 = fmul <2 x double> %18, %64
%69 = fmul <2 x double> %20, %67
%70 = bitcast double* %62 to <2 x double>*
store <2 x double> %68, <2 x double>* %70
%71 = bitcast double* %65 to <2 x double>*
store <2 x double> %69, <2 x double>* %71
br label %72
72: ; preds = %58, %61
%73 = icmp eq i64 %12, %15
br i1 %73, label %83, label %74
74: ; preds = %10, %72
%75 = phi i64 [ 2, %10 ], [ %16, %72 ]
br label %76
76: ; preds = %74, %76
%77 = phi i64 [ %81, %76 ], [ %75, %74 ]
%78 = getelementptr inbounds double, double* %0, i64 %77
%79 = load double, double* %78
%80 = fmul double %8, %79
store double %80, double* %78
%81 = add nuw nsw i64 %77, 1
%82 = icmp eq i64 %81, %11
br i1 %82, label %83, label %76
83: ; preds = %76, %72, %2
ret void
}