#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */
typedef size_t int64_t;
typedef size_t int32_t;
typedef scalar_t__ int16_t;

/* Variables and functions */
void acf(double *x, int32_t N, int64_t k, int16_t norm, double *rxx) {
  int64_t i, n;
  double sum = 0, mean = 0;
  for (i = 0; i < N; i++)
    sum += x[i];
  mean = sum / N;
  for (i = 0; i <= k; i++) {
    rxx[i] = 0.0;
    for (n = i; n < N; n++)
      rxx[i] += (x[n] - mean) * (x[n - i] - mean);
  }
  if (norm) {
    for (i = 1; i <= k; i++)
      rxx[i] /= rxx[0];
    rxx[0] = 1.0;
  }
}


x86 assembly (clang, -O0):
<acf>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
movq   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x10(%rbp),%rax
jae    1166 <acf+0x66>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    1130 <acf+0x30>
movsd  -0x40(%rbp),%xmm0
movaps 0xe8e(%rip),%xmm2
movq   -0x10(%rbp),%xmm1
punpckldq %xmm2,%xmm1
movapd 0xe8d(%rip),%xmm2
subpd  %xmm2,%xmm1
movaps %xmm1,%xmm2
unpckhpd %xmm1,%xmm1
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movq   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x18(%rbp),%rax
ja     1240 <acf+0x140>
mov    -0x28(%rbp),%rax
mov    -0x30(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x10(%rbp),%rax
jae    122a <acf+0x12a>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x48(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
sub    -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
subsd  -0x48(%rbp),%xmm2
mov    -0x28(%rbp),%rax
mov    -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    11c9 <acf+0xc9>
jmp    122f <acf+0x12f>
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    11a3 <acf+0xa3>
cmpq   $0x0,-0x20(%rbp)
je     12a0 <acf+0x1a0>
movq   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x18(%rbp),%rax
ja     1290 <acf+0x190>
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm1
mov    -0x28(%rbp),%rax
mov    -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    1253 <acf+0x153>
mov    -0x28(%rbp),%rax
movsd  0xd84(%rip),%xmm0
movsd  %xmm0,(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<acf>:
test   %rsi,%rsi
je     111e <acf+0x1e>
lea    -0x1(%rsi),%rax
mov    %esi,%r9d
and    $0x7,%r9d
cmp    $0x7,%rax
jae    1124 <acf+0x24>
xorpd  %xmm0,%xmm0
xor    %eax,%eax
jmp    1178 <acf+0x78>
xorpd  %xmm0,%xmm0
jmp    119f <acf+0x9f>
mov    %rsi,%r10
and    $0xfffffffffffffff8,%r10
xorpd  %xmm0,%xmm0
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
addsd  (%rdi,%rax,8),%xmm0
addsd  0x8(%rdi,%rax,8),%xmm0
addsd  0x10(%rdi,%rax,8),%xmm0
addsd  0x18(%rdi,%rax,8),%xmm0
addsd  0x20(%rdi,%rax,8),%xmm0
addsd  0x28(%rdi,%rax,8),%xmm0
addsd  0x30(%rdi,%rax,8),%xmm0
addsd  0x38(%rdi,%rax,8),%xmm0
add    $0x8,%rax
cmp    %rax,%r10
jne    1140 <acf+0x40>
test   %r9,%r9
je     119f <acf+0x9f>
lea    (%rdi,%rax,8),%r10
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
addsd  (%r10,%rax,8),%xmm0
add    $0x1,%rax
cmp    %rax,%r9
jne    1190 <acf+0x90>
movq   %rsi,%xmm1
punpckldq 0xe54(%rip),%xmm1
subpd  0xe5c(%rip),%xmm1
movapd %xmm1,%xmm2
unpckhpd %xmm1,%xmm2
addsd  %xmm1,%xmm2
divsd  %xmm2,%xmm0
xor    %r9d,%r9d
jmp    11dc <acf+0xdc>
nopl   0x0(%rax)
lea    0x1(%r9),%rax
cmp    %rdx,%r9
mov    %rax,%r9
je     1230 <acf+0x130>
movq   $0x0,(%r8,%r9,8)
cmp    %rsi,%r9
jae    11d0 <acf+0xd0>
xorpd  %xmm1,%xmm1
xor    %r10d,%r10d
mov    %r9,%rax
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%rdi,%rax,8),%xmm2
subsd  %xmm0,%xmm2
movsd  (%rdi,%r10,1),%xmm3
subsd  %xmm0,%xmm3
mulsd  %xmm2,%xmm3
addsd  %xmm3,%xmm1
movsd  %xmm1,(%r8,%r9,8)
add    $0x1,%rax
add    $0x8,%r10
cmp    %rsi,%rax
jb     1200 <acf+0x100>
jmp    11d0 <acf+0xd0>
test   %rcx,%rcx
je     1316 <acf+0x216>
test   %rdx,%rdx
je     1309 <acf+0x209>
movsd  (%r8),%xmm0
mov    $0x1,%eax
cmp    $0x1,%rdx
je     12dd <acf+0x1dd>
mov    %rdx,%rax
and    $0xfffffffffffffffe,%rax
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x2(%rax),%rsi
mov    %rsi,%rcx
shr    %rcx
add    $0x1,%rcx
test   %rsi,%rsi
je     1317 <acf+0x217>
mov    %rcx,%rdi
and    $0xfffffffffffffffe,%rdi
mov    $0x1,%esi
nopl   0x0(%rax,%rax,1)
movupd (%r8,%rsi,8),%xmm2
movupd 0x10(%r8,%rsi,8),%xmm3
divpd  %xmm1,%xmm2
movupd %xmm2,(%r8,%rsi,8)
divpd  %xmm1,%xmm3
movupd %xmm3,0x10(%r8,%rsi,8)
add    $0x4,%rsi
add    $0xfffffffffffffffe,%rdi
jne    1290 <acf+0x190>
test   $0x1,%cl
je     12d1 <acf+0x1d1>
movupd (%r8,%rsi,8),%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,(%r8,%rsi,8)
cmp    %rdx,%rax
je     1309 <acf+0x209>
mov    %rdx,%rax
or     $0x1,%rax
add    $0x1,%rdx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  (%r8,%rax,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%r8,%rax,8)
add    $0x1,%rax
cmp    %rax,%rdx
jne    12f0 <acf+0x1f0>
movabs $0x3ff0000000000000,%rax
mov    %rax,(%r8)
ret
mov    $0x1,%esi
test   $0x1,%cl
jne    12c1 <acf+0x1c1>
jmp    12d1 <acf+0x1d1>

LLVM IR (-O0):
define void @acf(double* %0, i64 %1, i64 %2, i64 %3, double* %4) {
%6 = alloca double*
%7 = alloca i64
%8 = alloca i64
%9 = alloca i64
%10 = alloca double*
%11 = alloca i64
%12 = alloca i64
%13 = alloca double
%14 = alloca double
store double* %0, double** %6
store i64 %1, i64* %7
store i64 %2, i64* %8
store i64 %3, i64* %9
store double* %4, double** %10
store double 0.000000e+00, double* %13
store double 0.000000e+00, double* %14
store i64 0, i64* %11
br label %15
15: ; preds = %26, %5
%16 = load i64, i64* %11
%17 = load i64, i64* %7
%18 = icmp ult i64 %16, %17
br i1 %18, label %19, label %29
19: ; preds = %15
%20 = load double*, double** %6
%21 = load i64, i64* %11
%22 = getelementptr inbounds double, double* %20, i64 %21
%23 = load double, double* %22
%24 = load double, double* %13
%25 = fadd double %24, %23
store double %25, double* %13
br label %26
26: ; preds = %19
%27 = load i64, i64* %11
%28 = add i64 %27, 1
store i64 %28, i64* %11
br label %15
29: ; preds = %15
%30 = load double, double* %13
%31 = load i64, i64* %7
%32 = uitofp i64 %31 to double
%33 = fdiv double %30, %32
store double %33, double* %14
store i64 0, i64* %11
br label %34
34: ; preds = %71, %29
%35 = load i64, i64* %11
%36 = load i64, i64* %8
%37 = icmp ule i64 %35, %36
br i1 %37, label %38, label %74
38: ; preds = %34
%39 = load double*, double** %10
%40 = load i64, i64* %11
%41 = getelementptr inbounds double, double* %39, i64 %40
store double 0.000000e+00, double* %41
%42 = load i64, i64* %11
store i64 %42, i64* %12
br label %43
43: ; preds = %67, %38
%44 = load i64, i64* %12
%45 = load i64, i64* %7
%46 = icmp ult i64 %44, %45
br i1 %46, label %47, label %70
47: ; preds = %43
%48 = load double*, double** %6
%49 = load i64, i64* %12
%50 = getelementptr inbounds double, double* %48, i64 %49
%51 = load double, double* %50
%52 = load double, double* %14
%53 = fsub double %51, %52
%54 = load double*, double** %6
%55 = load i64, i64* %12
%56 = load i64, i64* %11
%57 = sub i64 %55, %56
%58 = getelementptr inbounds double, double* %54, i64 %57
%59 = load double, double* %58
%60 = load double, double* %14
%61 = fsub double %59, %60
%62 = load double*, double** %10
%63 = load i64, i64* %11
%64 = getelementptr inbounds double, double* %62, i64 %63
%65 = load double, double* %64
%66 = call double @llvm.fmuladd.f64(double %53, double %61, double %65)
store double %66, double* %64
br label %67
67: ; preds = %47
%68 = load i64, i64* %12
%69 = add i64 %68, 1
store i64 %69, i64* %12
br label %43
70: ; preds = %43
br label %71
71: ; preds = %70
%72 = load i64, i64* %11
%73 = add i64 %72, 1
store i64 %73, i64* %11
br label %34
74: ; preds = %34
%75 = load i64, i64* %9
%76 = icmp ne i64 %75, 0
br i1 %76, label %77, label %97
77: ; preds = %74
store i64 1, i64* %11
br label %78
78: ; preds = %91, %77
%79 = load i64, i64* %11
%80 = load i64, i64* %8
%81 = icmp ule i64 %79, %80
br i1 %81, label %82, label %94
82: ; preds = %78
%83 = load double*, double** %10
%84 = getelementptr inbounds double, double* %83, i64 0
%85 = load double, double* %84
%86 = load double*, double** %10
%87 = load i64, i64* %11
%88 = getelementptr inbounds double, double* %86, i64 %87
%89 = load double, double* %88
%90 = fdiv double %89, %85
store double %90, double* %88
br label %91
91: ; preds = %82
%92 = load i64, i64* %11
%93 = add i64 %92, 1
store i64 %93, i64* %11
br label %78
94: ; preds = %78
%95 = load double*, double** %10
%96 = getelementptr inbounds double, double* %95, i64 0
store double 1.000000e+00, double* %96
br label %97
97: ; preds = %94, %74
ret void
}

LLVM IR (-O2):
define void @acf(double* nocapture %0, i64 %1, i64 %2, i64 %3, double* nocapture %4) {
%6 = icmp eq i64 %1, 0
br i1 %6, label %66, label %7
7: ; preds = %5
%8 = add i64 %1, -1
%9 = and i64 %1, 7
%10 = icmp ult i64 %8, 7
br i1 %10, label %51, label %11
11: ; preds = %7
%12 = and i64 %1, -8
br label %13
13: ; preds = %13, %11
%14 = phi double [ 0.000000e+00, %11 ], [ %47, %13 ]
%15 = phi i64 [ 0, %11 ], [ %48, %13 ]
%16 = phi i64 [ 0, %11 ], [ %49, %13 ]
%17 = getelementptr inbounds double, double* %0, i64 %15
%18 = load double, double* %17
%19 = fadd double %14, %18
%20 = or i64 %15, 1
%21 = getelementptr inbounds double, double* %0, i64 %20
%22 = load double, double* %21
%23 = fadd double %19, %22
%24 = or i64 %15, 2
%25 = getelementptr inbounds double, double* %0, i64 %24
%26 = load double, double* %25
%27 = fadd double %23, %26
%28 = or i64 %15, 3
%29 = getelementptr inbounds double, double* %0, i64 %28
%30 = load double, double* %29
%31 = fadd double %27, %30
%32 = or i64 %15, 4
%33 = getelementptr inbounds double, double* %0, i64 %32
%34 = load double, double* %33
%35 = fadd double %31, %34
%36 = or i64 %15, 5
%37 = getelementptr inbounds double, double* %0, i64 %36
%38 = load double, double* %37
%39 = fadd double %35, %38
%40 = or i64 %15, 6
%41 = getelementptr inbounds double, double* %0, i64 %40
%42 = load double, double* %41
%43 = fadd double %39, %42
%44 = or i64 %15, 7
%45 = getelementptr inbounds double, double* %0, i64 %44
%46 = load double, double* %45
%47 = fadd double %43, %46
%48 = add nuw i64 %15, 8
%49 = add i64 %16, 8
%50 = icmp eq i64 %49, %12
br i1 %50, label %51, label %13
51: ; preds = %13, %7
%52 = phi double [ undef, %7 ], [ %47, %13 ]
%53 = phi double [ 0.000000e+00, %7 ], [ %47, %13 ]
%54 = phi i64 [ 0, %7 ], [ %48, %13 ]
%55 = icmp eq i64 %9, 0
br i1 %55, label %66, label %56
56: ; preds = %51, %56
%57 = phi double [ %62, %56 ], [ %53, %51 ]
%58 = phi i64 [ %63, %56 ], [ %54, %51 ]
%59 = phi i64 [ %64, %56 ], [ 0, %51 ]
%60 = getelementptr inbounds double, double* %0, i64 %58
%61 = load double, double* %60
%62 = fadd double %57, %61
%63 = add nuw i64 %58, 1
%64 = add i64 %59, 1
%65 = icmp eq i64 %64, %9
br i1 %65, label %66, label %56
66: ; preds = %51, %56, %5
%67 = phi double [ 0.000000e+00, %5 ], [ %52, %51 ], [ %62, %56 ]
%68 = uitofp i64 %1 to double
%69 = fdiv double %67, %68
br label %70
70: ; preds = %66, %87
%71 = phi i64 [ 0, %66 ], [ %88, %87 ]
%72 = getelementptr inbounds double, double* %4, i64 %71
store double 0.000000e+00, double* %72
%73 = icmp ult i64 %71, %1
br i1 %73, label %74, label %87
74: ; preds = %70, %74
%75 = phi double [ %84, %74 ], [ 0.000000e+00, %70 ]
%76 = phi i64 [ %85, %74 ], [ %71, %70 ]
%77 = getelementptr inbounds double, double* %0, i64 %76
%78 = load double, double* %77
%79 = fsub double %78, %69
%80 = sub nuw i64 %76, %71
%81 = getelementptr inbounds double, double* %0, i64 %80
%82 = load double, double* %81
%83 = fsub double %82, %69
%84 = tail call double @llvm.fmuladd.f64(double %79, double %83, double %75)
store double %84, double* %72
%85 = add nuw i64 %76, 1
%86 = icmp ult i64 %85, %1
br i1 %86, label %74, label %87
87: ; preds = %74, %70
%88 = add nuw i64 %71, 1
%89 = icmp eq i64 %71, %2
br i1 %89, label %90, label %70
90: ; preds = %87
%91 = icmp eq i64 %3, 0
br i1 %91, label %150, label %92
92: ; preds = %90
%93 = icmp eq i64 %2, 0
br i1 %93, label %149, label %94
94: ; preds = %92
%95 = load double, double* %4
%96 = icmp eq i64 %2, 1
br i1 %96, label %140, label %97
97: ; preds = %94
%98 = and i64 %2, -2
%99 = or i64 %2, 1
%100 = insertelement <2 x double> poison, double %95, i64 0
%101 = shufflevector <2 x double> %100, <2 x double> poison, <2 x i32> zeroinitializer
%102 = add i64 %98, -2
%103 = lshr exact i64 %102, 1
%104 = add nuw i64 %103, 1
%105 = and i64 %104, 1
%106 = icmp eq i64 %102, 0
br i1 %106, label %129, label %107
107: ; preds = %97
%108 = and i64 %104, -2
br label %109
109: ; preds = %109, %107
%110 = phi i64 [ 0, %107 ], [ %124, %109 ]
%111 = phi i64 [ 0, %107 ], [ %125, %109 ]
%112 = or i64 %110, 1
%113 = getelementptr inbounds double, double* %4, i64 %112
%114 = bitcast double* %113 to <2 x double>*
%115 = load <2 x double>, <2 x double>* %114
%116 = fdiv <2 x double> %115, %101
%117 = bitcast double* %113 to <2 x double>*
store <2 x double> %116, <2 x double>* %117
%118 = or i64 %110, 3
%119 = getelementptr inbounds double, double* %4, i64 %118
%120 = bitcast double* %119 to <2 x double>*
%121 = load <2 x double>, <2 x double>* %120
%122 = fdiv <2 x double> %121, %101
%123 = bitcast double* %119 to <2 x double>*
store <2 x double> %122, <2 x double>* %123
%124 = add nuw i64 %110, 4
%125 = add nuw i64 %111, 2
%126 = icmp eq i64 %125, %108
br i1 %126, label %127, label %109
127: ; preds = %109
%128 = or i64 %124, 1
br label %129
129: ; preds = %127, %97
%130 = phi i64 [ 1, %97 ], [ %128, %127 ]
%131 = icmp eq i64 %105, 0
br i1 %131, label %138, label %132
132: ; preds = %129
%133 = getelementptr inbounds double, double* %4, i64 %130
%134 = bitcast double* %133 to <2 x double>*
%135 = load <2 x double>, <2 x double>* %134
%136 = fdiv <2 x double> %135, %101
%137 = bitcast double* %133 to <2 x double>*
store <2 x double> %136, <2 x double>* %137
br label %138
138: ; preds = %129, %132
%139 = icmp eq i64 %98, %2
br i1 %139, label %149, label %140
140: ; preds = %94, %138
%141 = phi i64 [ 1, %94 ], [ %99, %138 ]
br label %142
142: ; preds = %140, %142
%143 = phi i64 [ %147, %142 ], [ %141, %140 ]
%144 = getelementptr inbounds double, double* %4, i64 %143
%145 = load double, double* %144
%146 = fdiv double %145, %95
store double %146, double* %144
%147 = add nuw i64 %143, 1
%148 = icmp eq i64 %143, %2
br i1 %148, label %149, label %142
149: ; preds = %142, %138, %92
store double 1.000000e+00, double* %4
br label %150
150: ; preds = %149, %90
ret void
}