void EvalPtsW(int n, double *eval_pts) {
  int i;
  double twoN;
  twoN = (double)(2 * n);
  for (i = 0; i < n; i++)
    eval_pts[i] = ((2.0 * ((double)i) + 1.0) * 3.14159265358979323846) / twoN;
}


x86 assembly (clang, -O0):
<EvalPtsW>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    -0x4(%rbp),%eax
shl    $0x1,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x20(%rbp)
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    116e <EvalPtsW+0x6e>
cvtsi2sdl -0x14(%rbp),%xmm0
addsd  %xmm0,%xmm0
movsd  0xeca(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  0xeb6(%rip),%xmm1
mulsd  %xmm1,%xmm0
divsd  -0x20(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1121 <EvalPtsW+0x21>
pop    %rbp
ret

x86 assembly (clang, -O2):
<EvalPtsW>:
test   %edi,%edi
jle    1235 <EvalPtsW+0x135>
lea    (%rdi,%rdi,1),%eax
cvtsi2sd %eax,%xmm0
mov    %edi,%eax
cmp    $0x1,%edi
jne    111d <EvalPtsW+0x1d>
xor    %ecx,%ecx
jmp    11f5 <EvalPtsW+0xf5>
mov    %eax,%ecx
and    $0xfffffffe,%ecx
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x2(%rcx),%rdx
mov    %rdx,%r8
shr    %r8
add    $0x1,%r8
test   %rdx,%rdx
je     1236 <EvalPtsW+0x136>
mov    %r8,%rdx
and    $0xfffffffffffffffe,%rdx
movdqa 0xec0(%rip),%xmm2
xor    %edi,%edi
movapd 0xec6(%rip),%xmm3
movapd 0xece(%rip),%xmm4
movdqa 0xed6(%rip),%xmm5
movdqa 0xede(%rip),%xmm6
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
cvtdq2pd %xmm2,%xmm7
addpd  %xmm7,%xmm7
addpd  %xmm3,%xmm7
mulpd  %xmm4,%xmm7
divpd  %xmm1,%xmm7
movupd %xmm7,(%rsi,%rdi,8)
movdqa %xmm2,%xmm7
paddd  %xmm5,%xmm7
cvtdq2pd %xmm7,%xmm7
addpd  %xmm7,%xmm7
addpd  %xmm3,%xmm7
mulpd  %xmm4,%xmm7
divpd  %xmm1,%xmm7
movupd %xmm7,0x10(%rsi,%rdi,8)
add    $0x4,%rdi
paddd  %xmm6,%xmm2
add    $0xfffffffffffffffe,%rdx
jne    1180 <EvalPtsW+0x80>
cvtdq2pd %xmm2,%xmm2
test   $0x1,%r8b
je     11f0 <EvalPtsW+0xf0>
addpd  %xmm2,%xmm2
addpd  0xe41(%rip),%xmm2
mulpd  0xe49(%rip),%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,(%rsi,%rdi,8)
cmp    %rax,%rcx
je     1235 <EvalPtsW+0x135>
movsd  0xe03(%rip),%xmm1
movsd  0xe03(%rip),%xmm2
cs nopw 0x0(%rax,%rax,1)
nop
xorps  %xmm3,%xmm3
cvtsi2sd %ecx,%xmm3
addsd  %xmm3,%xmm3
addsd  %xmm1,%xmm3
mulsd  %xmm2,%xmm3
divsd  %xmm0,%xmm3
movsd  %xmm3,(%rsi,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    1210 <EvalPtsW+0x110>
ret
xorpd  %xmm2,%xmm2
movhpd 0xdbe(%rip),%xmm2
xor    %edi,%edi
test   $0x1,%r8b
jne    11d3 <EvalPtsW+0xd3>
jmp    11f0 <EvalPtsW+0xf0>

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

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