void HC_puckette_lock(long len, const double *y, double *z) {
  int k;
  z[0] = y[0];
  for (k = 1; k < (len + 1) / 2; k++) {
    z[k] = y[k];
    z[len - k] = y[len - k];
    if (k > 1) {
      z[k] += y[k - 1];
      z[len - k] += y[len - (k - 1)];
    }
    if (k < (len + 1) / 2 - 1) {
      z[k] += y[k + 1];
      z[len - k] += y[len - (k + 1)];
    }
  }
  if (len % 2 == 0) {
    z[len / 2] = y[len / 2];
  }
}


x86 assembly (clang, -O0):
<HC_puckette_lock>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x1,-0x1c(%rbp)
movslq -0x1c(%rbp),%rax
mov    %rax,-0x28(%rbp)
mov    -0x8(%rbp),%rax
add    $0x1,%rax
mov    $0x2,%ecx
cqto
idiv   %rcx
mov    %rax,%rcx
mov    -0x28(%rbp),%rax
cmp    %rcx,%rax
jge    128a <HC_puckette_lock+0x18a>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
mov    -0x8(%rbp),%rcx
movslq -0x1c(%rbp),%rdx
sub    %rdx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x8(%rbp),%rcx
movslq -0x1c(%rbp),%rdx
sub    %rdx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
cmpl   $0x1,-0x1c(%rbp)
jle    11f3 <HC_puckette_lock+0xf3>
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
mov    -0x8(%rbp),%rcx
mov    -0x1c(%rbp),%edx
sub    $0x1,%edx
movslq %edx,%rdx
sub    %rdx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x8(%rbp),%rcx
movslq -0x1c(%rbp),%rdx
sub    %rdx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movslq -0x1c(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x8(%rbp),%rax
add    $0x1,%rax
mov    $0x2,%ecx
cqto
idiv   %rcx
mov    %rax,%rcx
mov    -0x30(%rbp),%rax
sub    $0x1,%rcx
cmp    %rcx,%rax
jge    1277 <HC_puckette_lock+0x177>
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
mov    -0x8(%rbp),%rcx
mov    -0x1c(%rbp),%edx
add    $0x1,%edx
movslq %edx,%rdx
sub    %rdx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x8(%rbp),%rcx
movslq -0x1c(%rbp),%rdx
sub    %rdx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    127c <HC_puckette_lock+0x17c>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1127 <HC_puckette_lock+0x27>
mov    -0x8(%rbp),%rax
mov    $0x2,%ecx
cqto
idiv   %rcx
cmp    $0x0,%rdx
jne    12e6 <HC_puckette_lock+0x1e6>
mov    -0x10(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x8(%rbp),%rax
mov    $0x2,%ecx
cqto
idiv   %rcx
mov    %rax,%rcx
mov    -0x40(%rbp),%rax
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x8(%rbp),%rax
mov    $0x2,%ecx
cqto
idiv   %rcx
mov    %rax,%rcx
mov    -0x38(%rbp),%rax
movsd  %xmm0,(%rax,%rcx,8)
pop    %rbp
ret

x86 assembly (clang, -O2):
<HC_puckette_lock>:
push   %rbx
movsd  (%rsi),%xmm0
movsd  %xmm0,(%rdx)
cmp    $0x3,%rdi
jl     121a <HC_puckette_lock+0x11a>
lea    0x1(%rdi),%rax
shr    $0x3f,%rax
lea    (%rdi,%rax,1),%r9
add    $0x1,%r9
sar    %r9
cmp    $0x3,%r9
mov    $0x2,%r8d
cmovge %r9,%r8
movsd  0x8(%rsi),%xmm0
movsd  %xmm0,0x8(%rdx)
movsd  -0x8(%rsi,%rdi,8),%xmm0
movsd  %xmm0,-0x8(%rdx,%rdi,8)
cmp    $0x5,%rdi
jb     121a <HC_puckette_lock+0x11a>
add    $0xffffffffffffffff,%r9
movsd  0x10(%rsi),%xmm0
addsd  0x8(%rdx),%xmm0
movsd  %xmm0,0x8(%rdx)
movsd  -0x10(%rsi,%rdi,8),%xmm0
addsd  -0x8(%rdx,%rdi,8),%xmm0
movsd  %xmm0,-0x8(%rdx,%rdi,8)
lea    (%rdx,%rdi,8),%r11
add    $0xfffffffffffffff0,%r11
lea    (%rsi,%rdi,8),%r10
add    $0xfffffffffffffff8,%r10
add    $0xfffffffffffffffe,%r8
xor    %eax,%eax
xor    %ecx,%ecx
jmp    11ad <HC_puckette_lock+0xad>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%rcx
add    $0xfffffffffffffff8,%rax
cmp    %rcx,%r8
je     121a <HC_puckette_lock+0x11a>
lea    0x2(%rcx),%rbx
movsd  0x10(%rsi,%rcx,8),%xmm0
movsd  %xmm0,0x10(%rdx,%rcx,8)
movsd  -0x8(%r10,%rax,1),%xmm0
movsd  %xmm0,(%r11,%rax,1)
movsd  0x8(%rsi,%rcx,8),%xmm0
addsd  0x10(%rdx,%rcx,8),%xmm0
movsd  %xmm0,0x10(%rdx,%rcx,8)
movsd  (%r10,%rax,1),%xmm0
addsd  (%r11,%rax,1),%xmm0
movsd  %xmm0,(%r11,%rax,1)
cmp    %rbx,%r9
jle    11a0 <HC_puckette_lock+0xa0>
movsd  0x18(%rsi,%rcx,8),%xmm0
addsd  0x10(%rdx,%rcx,8),%xmm0
movsd  %xmm0,0x10(%rdx,%rcx,8)
movsd  -0x10(%r10,%rax,1),%xmm0
addsd  (%r11,%rax,1),%xmm0
movsd  %xmm0,(%r11,%rax,1)
jmp    11a0 <HC_puckette_lock+0xa0>
test   $0x1,%dil
jne    123c <HC_puckette_lock+0x13c>
mov    %rdi,%rax
shr    $0x3f,%rax
add    %rdi,%rax
shl    $0x2,%rax
and    $0xfffffffffffffff8,%rax
movsd  (%rsi,%rax,1),%xmm0
movsd  %xmm0,(%rdx,%rax,1)
pop    %rbx
ret

LLVM IR (-O0):
define void @HC_puckette_lock(i64 %0, double* %1, double* %2) {
%4 = alloca i64
%5 = alloca double*
%6 = alloca double*
%7 = alloca i32
store i64 %0, i64* %4
store double* %1, double** %5
store double* %2, double** %6
%8 = load double*, double** %5
%9 = getelementptr inbounds double, double* %8, i64 0
%10 = load double, double* %9
%11 = load double*, double** %6
%12 = getelementptr inbounds double, double* %11, i64 0
store double %10, double* %12
store i32 1, i32* %7
br label %13
13: ; preds = %112, %3
%14 = load i32, i32* %7
%15 = sext i32 %14 to i64
%16 = load i64, i64* %4
%17 = add nsw i64 %16, 1
%18 = sdiv i64 %17, 2
%19 = icmp slt i64 %15, %18
br i1 %19, label %20, label %115
20: ; preds = %13
%21 = load double*, double** %5
%22 = load i32, i32* %7
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds double, double* %21, i64 %23
%25 = load double, double* %24
%26 = load double*, double** %6
%27 = load i32, i32* %7
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds double, double* %26, i64 %28
store double %25, double* %29
%30 = load double*, double** %5
%31 = load i64, i64* %4
%32 = load i32, i32* %7
%33 = sext i32 %32 to i64
%34 = sub nsw i64 %31, %33
%35 = getelementptr inbounds double, double* %30, i64 %34
%36 = load double, double* %35
%37 = load double*, double** %6
%38 = load i64, i64* %4
%39 = load i32, i32* %7
%40 = sext i32 %39 to i64
%41 = sub nsw i64 %38, %40
%42 = getelementptr inbounds double, double* %37, i64 %41
store double %36, double* %42
%43 = load i32, i32* %7
%44 = icmp sgt i32 %43, 1
br i1 %44, label %45, label %74
45: ; preds = %20
%46 = load double*, double** %5
%47 = load i32, i32* %7
%48 = sub nsw i32 %47, 1
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds double, double* %46, i64 %49
%51 = load double, double* %50
%52 = load double*, double** %6
%53 = load i32, i32* %7
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds double, double* %52, i64 %54
%56 = load double, double* %55
%57 = fadd double %56, %51
store double %57, double* %55
%58 = load double*, double** %5
%59 = load i64, i64* %4
%60 = load i32, i32* %7
%61 = sub nsw i32 %60, 1
%62 = sext i32 %61 to i64
%63 = sub nsw i64 %59, %62
%64 = getelementptr inbounds double, double* %58, i64 %63
%65 = load double, double* %64
%66 = load double*, double** %6
%67 = load i64, i64* %4
%68 = load i32, i32* %7
%69 = sext i32 %68 to i64
%70 = sub nsw i64 %67, %69
%71 = getelementptr inbounds double, double* %66, i64 %70
%72 = load double, double* %71
%73 = fadd double %72, %65
store double %73, double* %71
br label %74
74: ; preds = %45, %20
%75 = load i32, i32* %7
%76 = sext i32 %75 to i64
%77 = load i64, i64* %4
%78 = add nsw i64 %77, 1
%79 = sdiv i64 %78, 2
%80 = sub nsw i64 %79, 1
%81 = icmp slt i64 %76, %80
br i1 %81, label %82, label %111
82: ; preds = %74
%83 = load double*, double** %5
%84 = load i32, i32* %7
%85 = add nsw i32 %84, 1
%86 = sext i32 %85 to i64
%87 = getelementptr inbounds double, double* %83, i64 %86
%88 = load double, double* %87
%89 = load double*, double** %6
%90 = load i32, i32* %7
%91 = sext i32 %90 to i64
%92 = getelementptr inbounds double, double* %89, i64 %91
%93 = load double, double* %92
%94 = fadd double %93, %88
store double %94, double* %92
%95 = load double*, double** %5
%96 = load i64, i64* %4
%97 = load i32, i32* %7
%98 = add nsw i32 %97, 1
%99 = sext i32 %98 to i64
%100 = sub nsw i64 %96, %99
%101 = getelementptr inbounds double, double* %95, i64 %100
%102 = load double, double* %101
%103 = load double*, double** %6
%104 = load i64, i64* %4
%105 = load i32, i32* %7
%106 = sext i32 %105 to i64
%107 = sub nsw i64 %104, %106
%108 = getelementptr inbounds double, double* %103, i64 %107
%109 = load double, double* %108
%110 = fadd double %109, %102
store double %110, double* %108
br label %111
111: ; preds = %82, %74
br label %112
112: ; preds = %111
%113 = load i32, i32* %7
%114 = add nsw i32 %113, 1
store i32 %114, i32* %7
br label %13
115: ; preds = %13
%116 = load i64, i64* %4
%117 = srem i64 %116, 2
%118 = icmp eq i64 %117, 0
br i1 %118, label %119, label %129
119: ; preds = %115
%120 = load double*, double** %5
%121 = load i64, i64* %4
%122 = sdiv i64 %121, 2
%123 = getelementptr inbounds double, double* %120, i64 %122
%124 = load double, double* %123
%125 = load double*, double** %6
%126 = load i64, i64* %4
%127 = sdiv i64 %126, 2
%128 = getelementptr inbounds double, double* %125, i64 %127
store double %124, double* %128
br label %129
129: ; preds = %119, %115
ret void
}

LLVM IR (-O2):
define void @HC_puckette_lock(i64 %0, double* nocapture %1, double* nocapture %2) {
%4 = load double, double* %1
store double %4, double* %2
%5 = add nsw i64 %0, 1
%6 = sdiv i64 %5, 2
%7 = icmp sgt i64 %0, 2
br i1 %7, label %8, label %62
8: ; preds = %3
%9 = add nsw i64 %6, -1
%10 = call i64 @llvm.smax.i64(i64 %6, i64 2)
%11 = getelementptr inbounds double, double* %1, i64 1
%12 = load double, double* %11
%13 = getelementptr inbounds double, double* %2, i64 1
store double %12, double* %13
%14 = add nsw i64 %0, -1
%15 = getelementptr inbounds double, double* %1, i64 %14
%16 = load double, double* %15
%17 = getelementptr inbounds double, double* %2, i64 %14
store double %16, double* %17
%18 = icmp ugt i64 %0, 4
br i1 %18, label %19, label %62
19: ; preds = %8
%20 = getelementptr inbounds double, double* %1, i64 2
%21 = load double, double* %20
%22 = load double, double* %13
%23 = fadd double %21, %22
store double %23, double* %13
%24 = add nsw i64 %0, -2
%25 = getelementptr inbounds double, double* %1, i64 %24
%26 = load double, double* %25
%27 = load double, double* %17
%28 = fadd double %26, %27
store double %28, double* %17
br label %29
29: ; preds = %19, %60
%30 = phi i64 [ 2, %19 ], [ %49, %60 ]
%31 = getelementptr inbounds double, double* %1, i64 %30
%32 = load double, double* %31
%33 = getelementptr inbounds double, double* %2, i64 %30
store double %32, double* %33
%34 = sub nsw i64 %0, %30
%35 = getelementptr inbounds double, double* %1, i64 %34
%36 = load double, double* %35
%37 = getelementptr inbounds double, double* %2, i64 %34
store double %36, double* %37
%38 = add nsw i64 %30, -1
%39 = getelementptr inbounds double, double* %1, i64 %38
%40 = load double, double* %39
%41 = load double, double* %33
%42 = fadd double %40, %41
store double %42, double* %33
%43 = sub nsw i64 %0, %38
%44 = getelementptr inbounds double, double* %1, i64 %43
%45 = load double, double* %44
%46 = load double, double* %37
%47 = fadd double %45, %46
store double %47, double* %37
%48 = icmp sgt i64 %9, %30
%49 = add nuw nsw i64 %30, 1
br i1 %48, label %50, label %60
50: ; preds = %29
%51 = getelementptr inbounds double, double* %1, i64 %49
%52 = load double, double* %51
%53 = load double, double* %33
%54 = fadd double %52, %53
store double %54, double* %33
%55 = sub nsw i64 %0, %49
%56 = getelementptr inbounds double, double* %1, i64 %55
%57 = load double, double* %56
%58 = load double, double* %37
%59 = fadd double %57, %58
store double %59, double* %37
br label %60
60: ; preds = %29, %50
%61 = icmp eq i64 %49, %10
br i1 %61, label %62, label %29
62: ; preds = %60, %8, %3
%63 = and i64 %0, 1
%64 = icmp eq i64 %63, 0
br i1 %64, label %65, label %70
65: ; preds = %62
%66 = sdiv i64 %0, 2
%67 = getelementptr inbounds double, double* %1, i64 %66
%68 = load double, double* %67
%69 = getelementptr inbounds double, double* %2, i64 %66
store double %68, double* %69
br label %70
70: ; preds = %65, %62
ret void
}