void alphabgrnd(int n, double y[]) {
  int i, j, m;
  double s = 65536.0, w = 2053.0, v = 13849.0, r = 0.0, t = 0.0;
  for (i = 0; i < n; i++) {
    t = 0.0;
    for (j = 0; j < 12; j++) {
      r = w * r + v;
      m = (int)(r / s);
      r = r - m * s;
      t = t + r / s;
    }
    y[i] = 0.5 * (t - 6.0);
  }
  return;
}


x86 assembly (clang, -O0):
<alphabgrnd>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
movsd  0xefd(%rip),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  0xee8(%rip),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  0xed3(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    122c <alphabgrnd+0x12c>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x18(%rbp)
cmpl   $0xc,-0x18(%rbp)
jge    11f4 <alphabgrnd+0xf4>
movsd  -0x30(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm2
movsd  -0x38(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
divsd  -0x28(%rbp),%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x1c(%rbp)
movsd  -0x40(%rbp),%xmm1
cvtsi2sdl -0x1c(%rbp),%xmm0
movsd  -0x28(%rbp),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm1
divsd  -0x28(%rbp),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1164 <alphabgrnd+0x64>
movsd  -0x48(%rbp),%xmm1
movsd  0xe1f(%rip),%xmm0
subsd  %xmm0,%xmm1
movsd  0xe0b(%rip),%xmm0
mulsd  %xmm1,%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    1149 <alphabgrnd+0x49>
pop    %rbp
ret

x86 assembly (clang, -O2):
<alphabgrnd>:
test   %edi,%edi
jle    137d <alphabgrnd+0x27d>
mov    %edi,%eax
xorpd  %xmm8,%xmm8
xor    %ecx,%ecx
movsd  0xee7(%rip),%xmm1
movsd  0xee7(%rip),%xmm2
movsd  0xee7(%rip),%xmm3
movsd  0xee7(%rip),%xmm4
movsd  0xee6(%rip),%xmm9
movsd  0xee5(%rip),%xmm10
xorpd  %xmm6,%xmm6
nopw   0x0(%rax,%rax,1)
mulsd  %xmm1,%xmm6
addsd  %xmm2,%xmm6
movapd %xmm6,%xmm7
mulsd  %xmm3,%xmm7
cvttpd2dq %xmm7,%xmm7
cvtdq2pd %xmm7,%xmm7
mulsd  %xmm4,%xmm7
subsd  %xmm7,%xmm6
movapd %xmm6,%xmm7
mulsd  %xmm3,%xmm7
mulsd  %xmm1,%xmm6
addsd  %xmm2,%xmm6
movapd %xmm6,%xmm0
mulsd  %xmm3,%xmm0
cvttpd2dq %xmm0,%xmm0
cvtdq2pd %xmm0,%xmm0
addsd  %xmm8,%xmm7
mulsd  %xmm4,%xmm0
subsd  %xmm0,%xmm6
movapd %xmm6,%xmm0
mulsd  %xmm3,%xmm0
addsd  %xmm7,%xmm0
mulsd  %xmm1,%xmm6
addsd  %xmm2,%xmm6
movapd %xmm6,%xmm7
mulsd  %xmm3,%xmm7
cvttpd2dq %xmm7,%xmm7
cvtdq2pd %xmm7,%xmm7
mulsd  %xmm4,%xmm7
subsd  %xmm7,%xmm6
movapd %xmm6,%xmm7
mulsd  %xmm3,%xmm7
mulsd  %xmm1,%xmm6
addsd  %xmm2,%xmm6
movapd %xmm6,%xmm5
mulsd  %xmm3,%xmm5
cvttpd2dq %xmm5,%xmm5
addsd  %xmm0,%xmm7
cvtdq2pd %xmm5,%xmm0
mulsd  %xmm4,%xmm0
subsd  %xmm0,%xmm6
movapd %xmm6,%xmm0
mulsd  %xmm3,%xmm0
addsd  %xmm7,%xmm0
mulsd  %xmm1,%xmm6
addsd  %xmm2,%xmm6
movapd %xmm6,%xmm5
mulsd  %xmm3,%xmm5
cvttpd2dq %xmm5,%xmm5
cvtdq2pd %xmm5,%xmm5
mulsd  %xmm4,%xmm5
subsd  %xmm5,%xmm6
movapd %xmm6,%xmm5
mulsd  %xmm3,%xmm5
mulsd  %xmm1,%xmm6
addsd  %xmm2,%xmm6
movapd %xmm6,%xmm7
mulsd  %xmm3,%xmm7
cvttpd2dq %xmm7,%xmm7
cvtdq2pd %xmm7,%xmm7
addsd  %xmm0,%xmm5
mulsd  %xmm4,%xmm7
subsd  %xmm7,%xmm6
movapd %xmm6,%xmm0
mulsd  %xmm3,%xmm0
addsd  %xmm5,%xmm0
mulsd  %xmm1,%xmm6
addsd  %xmm2,%xmm6
movapd %xmm6,%xmm5
mulsd  %xmm3,%xmm5
cvttpd2dq %xmm5,%xmm5
cvtdq2pd %xmm5,%xmm5
mulsd  %xmm4,%xmm5
subsd  %xmm5,%xmm6
movapd %xmm6,%xmm5
mulsd  %xmm3,%xmm5
mulsd  %xmm1,%xmm6
addsd  %xmm2,%xmm6
movapd %xmm6,%xmm7
mulsd  %xmm3,%xmm7
cvttpd2dq %xmm7,%xmm7
addsd  %xmm0,%xmm5
cvtdq2pd %xmm7,%xmm0
mulsd  %xmm4,%xmm0
subsd  %xmm0,%xmm6
movapd %xmm6,%xmm0
mulsd  %xmm3,%xmm0
addsd  %xmm5,%xmm0
mulsd  %xmm1,%xmm6
addsd  %xmm2,%xmm6
movapd %xmm6,%xmm5
mulsd  %xmm3,%xmm5
cvttpd2dq %xmm5,%xmm5
cvtdq2pd %xmm5,%xmm5
mulsd  %xmm4,%xmm5
subsd  %xmm5,%xmm6
movapd %xmm6,%xmm5
mulsd  %xmm3,%xmm5
mulsd  %xmm1,%xmm6
addsd  %xmm2,%xmm6
movapd %xmm6,%xmm7
mulsd  %xmm3,%xmm7
cvttpd2dq %xmm7,%xmm7
cvtdq2pd %xmm7,%xmm7
addsd  %xmm0,%xmm5
mulsd  %xmm4,%xmm7
subsd  %xmm7,%xmm6
movapd %xmm6,%xmm0
mulsd  %xmm3,%xmm0
addsd  %xmm5,%xmm0
mulsd  %xmm1,%xmm6
addsd  %xmm2,%xmm6
movapd %xmm6,%xmm5
mulsd  %xmm3,%xmm5
cvttpd2dq %xmm5,%xmm5
cvtdq2pd %xmm5,%xmm5
mulsd  %xmm4,%xmm5
subsd  %xmm5,%xmm6
movapd %xmm6,%xmm5
mulsd  %xmm3,%xmm5
mulsd  %xmm1,%xmm6
addsd  %xmm2,%xmm6
movapd %xmm6,%xmm7
mulsd  %xmm3,%xmm7
cvttpd2dq %xmm7,%xmm7
addsd  %xmm0,%xmm5
cvtdq2pd %xmm7,%xmm0
mulsd  %xmm4,%xmm0
subsd  %xmm0,%xmm6
movapd %xmm6,%xmm0
mulsd  %xmm3,%xmm0
addsd  %xmm5,%xmm0
addsd  %xmm9,%xmm0
mulsd  %xmm10,%xmm0
movsd  %xmm0,(%rsi,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    1150 <alphabgrnd+0x50>
ret

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

LLVM IR (-O2):
define void @alphabgrnd(i32 %0, double* nocapture writeonly %1) {
%3 = icmp sgt i32 %0, 0
br i1 %3, label %4, label %110
4: ; preds = %2
%5 = zext i32 %0 to i64
br label %6
6: ; preds = %4, %6
%7 = phi i64 [ 0, %4 ], [ %108, %6 ]
%8 = phi double [ 0.000000e+00, %4 ], [ %102, %6 ]
%9 = tail call double @llvm.fmuladd.f64(double %8, double 2.053000e+03, double 1.384900e+04)
%10 = fmul double %9, 0x3EF0000000000000
%11 = fptosi double %10 to i32
%12 = sitofp i32 %11 to double
%13 = fneg double %12
%14 = tail call double @llvm.fmuladd.f64(double %13, double 6.553600e+04, double %9)
%15 = fmul double %14, 0x3EF0000000000000
%16 = fadd double %15, 0.000000e+00
%17 = tail call double @llvm.fmuladd.f64(double %14, double 2.053000e+03, double 1.384900e+04)
%18 = fmul double %17, 0x3EF0000000000000
%19 = fptosi double %18 to i32
%20 = sitofp i32 %19 to double
%21 = fneg double %20
%22 = tail call double @llvm.fmuladd.f64(double %21, double 6.553600e+04, double %17)
%23 = fmul double %22, 0x3EF0000000000000
%24 = fadd double %16, %23
%25 = tail call double @llvm.fmuladd.f64(double %22, double 2.053000e+03, double 1.384900e+04)
%26 = fmul double %25, 0x3EF0000000000000
%27 = fptosi double %26 to i32
%28 = sitofp i32 %27 to double
%29 = fneg double %28
%30 = tail call double @llvm.fmuladd.f64(double %29, double 6.553600e+04, double %25)
%31 = fmul double %30, 0x3EF0000000000000
%32 = fadd double %24, %31
%33 = tail call double @llvm.fmuladd.f64(double %30, double 2.053000e+03, double 1.384900e+04)
%34 = fmul double %33, 0x3EF0000000000000
%35 = fptosi double %34 to i32
%36 = sitofp i32 %35 to double
%37 = fneg double %36
%38 = tail call double @llvm.fmuladd.f64(double %37, double 6.553600e+04, double %33)
%39 = fmul double %38, 0x3EF0000000000000
%40 = fadd double %32, %39
%41 = tail call double @llvm.fmuladd.f64(double %38, double 2.053000e+03, double 1.384900e+04)
%42 = fmul double %41, 0x3EF0000000000000
%43 = fptosi double %42 to i32
%44 = sitofp i32 %43 to double
%45 = fneg double %44
%46 = tail call double @llvm.fmuladd.f64(double %45, double 6.553600e+04, double %41)
%47 = fmul double %46, 0x3EF0000000000000
%48 = fadd double %40, %47
%49 = tail call double @llvm.fmuladd.f64(double %46, double 2.053000e+03, double 1.384900e+04)
%50 = fmul double %49, 0x3EF0000000000000
%51 = fptosi double %50 to i32
%52 = sitofp i32 %51 to double
%53 = fneg double %52
%54 = tail call double @llvm.fmuladd.f64(double %53, double 6.553600e+04, double %49)
%55 = fmul double %54, 0x3EF0000000000000
%56 = fadd double %48, %55
%57 = tail call double @llvm.fmuladd.f64(double %54, double 2.053000e+03, double 1.384900e+04)
%58 = fmul double %57, 0x3EF0000000000000
%59 = fptosi double %58 to i32
%60 = sitofp i32 %59 to double
%61 = fneg double %60
%62 = tail call double @llvm.fmuladd.f64(double %61, double 6.553600e+04, double %57)
%63 = fmul double %62, 0x3EF0000000000000
%64 = fadd double %56, %63
%65 = tail call double @llvm.fmuladd.f64(double %62, double 2.053000e+03, double 1.384900e+04)
%66 = fmul double %65, 0x3EF0000000000000
%67 = fptosi double %66 to i32
%68 = sitofp i32 %67 to double
%69 = fneg double %68
%70 = tail call double @llvm.fmuladd.f64(double %69, double 6.553600e+04, double %65)
%71 = fmul double %70, 0x3EF0000000000000
%72 = fadd double %64, %71
%73 = tail call double @llvm.fmuladd.f64(double %70, double 2.053000e+03, double 1.384900e+04)
%74 = fmul double %73, 0x3EF0000000000000
%75 = fptosi double %74 to i32
%76 = sitofp i32 %75 to double
%77 = fneg double %76
%78 = tail call double @llvm.fmuladd.f64(double %77, double 6.553600e+04, double %73)
%79 = fmul double %78, 0x3EF0000000000000
%80 = fadd double %72, %79
%81 = tail call double @llvm.fmuladd.f64(double %78, double 2.053000e+03, double 1.384900e+04)
%82 = fmul double %81, 0x3EF0000000000000
%83 = fptosi double %82 to i32
%84 = sitofp i32 %83 to double
%85 = fneg double %84
%86 = tail call double @llvm.fmuladd.f64(double %85, double 6.553600e+04, double %81)
%87 = fmul double %86, 0x3EF0000000000000
%88 = fadd double %80, %87
%89 = tail call double @llvm.fmuladd.f64(double %86, double 2.053000e+03, double 1.384900e+04)
%90 = fmul double %89, 0x3EF0000000000000
%91 = fptosi double %90 to i32
%92 = sitofp i32 %91 to double
%93 = fneg double %92
%94 = tail call double @llvm.fmuladd.f64(double %93, double 6.553600e+04, double %89)
%95 = fmul double %94, 0x3EF0000000000000
%96 = fadd double %88, %95
%97 = tail call double @llvm.fmuladd.f64(double %94, double 2.053000e+03, double 1.384900e+04)
%98 = fmul double %97, 0x3EF0000000000000
%99 = fptosi double %98 to i32
%100 = sitofp i32 %99 to double
%101 = fneg double %100
%102 = tail call double @llvm.fmuladd.f64(double %101, double 6.553600e+04, double %97)
%103 = fmul double %102, 0x3EF0000000000000
%104 = fadd double %96, %103
%105 = fadd double %104, -6.000000e+00
%106 = fmul double %105, 5.000000e-01
%107 = getelementptr inbounds double, double* %1, i64 %7
store double %106, double* %107
%108 = add nuw nsw i64 %7, 1
%109 = icmp eq i64 %108, %5
br i1 %109, label %110, label %6
110: ; preds = %6, %2
ret void
}