void downhipbs(x, n, hpf, m, y) double x[], hpf[], y[];
int n, m;
{
  int n2, i, h, j, k;
  double s;
  n2 = n / 2;
  k = (m + 1) / 2;
  for (i = 0; i < n2; i++) {
    s = 0;
    j = k;
    for (h = 0; h < m; h++) {
      while (j >= n)
        j -= n;
      while (j < 0)
        j += n;
      s += hpf[h] * x[j];
      j--;
    }
    k += 2;
    y[i] = s;
  }
}


x86 assembly (clang, -O0):
<downhipbs>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0xc(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x2c(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x3c(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    1213 <downhipbs+0x113>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x38(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    11ea <downhipbs+0xea>
jmp    116e <downhipbs+0x6e>
mov    -0x38(%rbp),%eax
cmp    -0xc(%rbp),%eax
jl     118a <downhipbs+0x8a>
mov    -0xc(%rbp),%ecx
mov    -0x38(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x38(%rbp)
jmp    116e <downhipbs+0x6e>
jmp    118f <downhipbs+0x8f>
cmpl   $0x0,-0x38(%rbp)
jge    11a7 <downhipbs+0xa7>
mov    -0xc(%rbp),%eax
add    -0x38(%rbp),%eax
mov    %eax,-0x38(%rbp)
jmp    118f <downhipbs+0x8f>
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x38(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x38(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    115d <downhipbs+0x5d>
mov    -0x3c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x3c(%rbp)
movsd  -0x48(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    113c <downhipbs+0x3c>
pop    %rbp
ret

x86 assembly (clang, -O2):
<downhipbs>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
cmp    $0x2,%esi
jl     11b5 <downhipbs+0xb5>
mov    %esi,%r9d
shr    %r9d
lea    0x1(%rcx),%eax
shr    $0x1f,%eax
lea    (%rcx,%rax,1),%r10d
add    $0x1,%r10d
sar    %r10d
mov    %ecx,%r11d
mov    %esi,%r15d
neg    %r15d
xor    %r14d,%r14d
jmp    1153 <downhipbs+0x53>
cs nopw 0x0(%rax,%rax,1)
add    $0x2,%r10d
movsd  %xmm0,(%r8,%r14,8)
add    $0x1,%r14
cmp    %r9,%r14
je     11b5 <downhipbs+0xb5>
xorpd  %xmm0,%xmm0
test   %ecx,%ecx
jle    1140 <downhipbs+0x40>
xor    %r12d,%r12d
mov    %r10d,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
lea    (%rsi,%rax,1),%ebx
add    %r15d,%eax
cs nopw 0x0(%rax,%rax,1)
mov    %eax,%ebp
sub    %esi,%ebx
sub    %esi,%eax
cmp    %esi,%ebx
jge    1180 <downhipbs+0x80>
nopw   0x0(%rax,%rax,1)
add    %esi,%ebp
js     1190 <downhipbs+0x90>
movsd  (%rdx,%r12,8),%xmm1
mov    %ebp,%eax
mulsd  (%rdi,%rax,8),%xmm1
addsd  %xmm1,%xmm0
add    $0xffffffff,%ebp
add    $0x1,%r12
mov    %ebp,%eax
cmp    %r11,%r12
jne    1170 <downhipbs+0x70>
jmp    1140 <downhipbs+0x40>
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @downhipbs(double* %0, i32 %1, double* %2, i32 %3, double* %4) {
%6 = alloca double*
%7 = alloca i32
%8 = alloca double*
%9 = alloca i32
%10 = alloca double*
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca double
store double* %0, double** %6
store i32 %1, i32* %7
store double* %2, double** %8
store i32 %3, i32* %9
store double* %4, double** %10
%17 = load i32, i32* %7
%18 = sdiv i32 %17, 2
store i32 %18, i32* %11
%19 = load i32, i32* %9
%20 = add nsw i32 %19, 1
%21 = sdiv i32 %20, 2
store i32 %21, i32* %15
store i32 0, i32* %12
br label %22
22: ; preds = %75, %5
%23 = load i32, i32* %12
%24 = load i32, i32* %11
%25 = icmp slt i32 %23, %24
br i1 %25, label %26, label %78
26: ; preds = %22
store double 0.000000e+00, double* %16
%27 = load i32, i32* %15
store i32 %27, i32* %14
store i32 0, i32* %13
br label %28
28: ; preds = %64, %26
%29 = load i32, i32* %13
%30 = load i32, i32* %9
%31 = icmp slt i32 %29, %30
br i1 %31, label %32, label %67
32: ; preds = %28
br label %33
33: ; preds = %37, %32
%34 = load i32, i32* %14
%35 = load i32, i32* %7
%36 = icmp sge i32 %34, %35
br i1 %36, label %37, label %41
37: ; preds = %33
%38 = load i32, i32* %7
%39 = load i32, i32* %14
%40 = sub nsw i32 %39, %38
store i32 %40, i32* %14
br label %33
41: ; preds = %33
br label %42
42: ; preds = %45, %41
%43 = load i32, i32* %14
%44 = icmp slt i32 %43, 0
br i1 %44, label %45, label %49
45: ; preds = %42
%46 = load i32, i32* %7
%47 = load i32, i32* %14
%48 = add nsw i32 %47, %46
store i32 %48, i32* %14
br label %42
49: ; preds = %42
%50 = load double*, double** %8
%51 = load i32, i32* %13
%52 = sext i32 %51 to i64
%53 = getelementptr inbounds double, double* %50, i64 %52
%54 = load double, double* %53
%55 = load double*, double** %6
%56 = load i32, i32* %14
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds double, double* %55, i64 %57
%59 = load double, double* %58
%60 = load double, double* %16
%61 = call double @llvm.fmuladd.f64(double %54, double %59, double %60)
store double %61, double* %16
%62 = load i32, i32* %14
%63 = add nsw i32 %62, -1
store i32 %63, i32* %14
br label %64
64: ; preds = %49
%65 = load i32, i32* %13
%66 = add nsw i32 %65, 1
store i32 %66, i32* %13
br label %28
67: ; preds = %28
%68 = load i32, i32* %15
%69 = add nsw i32 %68, 2
store i32 %69, i32* %15
%70 = load double, double* %16
%71 = load double*, double** %10
%72 = load i32, i32* %12
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds double, double* %71, i64 %73
store double %70, double* %74
br label %75
75: ; preds = %67
%76 = load i32, i32* %12
%77 = add nsw i32 %76, 1
store i32 %77, i32* %12
br label %22
78: ; preds = %22
ret void
}

LLVM IR (-O2):
define void @downhipbs(double* nocapture %0, i32 %1, double* nocapture %2, i32 %3, double* nocapture writeonly %4) {
%6 = icmp sgt i32 %1, 1
br i1 %6, label %7, label %45
7: ; preds = %5
%8 = lshr i32 %1, 1
%9 = add nsw i32 %3, 1
%10 = sdiv i32 %9, 2
%11 = icmp sgt i32 %3, 0
%12 = zext i32 %8 to i64
%13 = zext i32 %3 to i64
br label %14
14: ; preds = %7, %39
%15 = phi i64 [ 0, %7 ], [ %43, %39 ]
%16 = phi i32 [ %10, %7 ], [ %41, %39 ]
br i1 %11, label %17, label %39
17: ; preds = %14, %29
%18 = phi i64 [ %37, %29 ], [ 0, %14 ]
%19 = phi double [ %35, %29 ], [ 0.000000e+00, %14 ]
%20 = phi i32 [ %36, %29 ], [ %16, %14 ]
br label %21
21: ; preds = %21, %17
%22 = phi i32 [ %24, %21 ], [ %20, %17 ]
%23 = icmp slt i32 %22, %1
%24 = sub nsw i32 %22, %1
br i1 %23, label %25, label %21
25: ; preds = %21, %25
%26 = phi i32 [ %28, %25 ], [ %22, %21 ]
%27 = icmp slt i32 %26, 0
%28 = add nsw i32 %26, %1
br i1 %27, label %25, label %29
29: ; preds = %25
%30 = getelementptr inbounds double, double* %2, i64 %18
%31 = load double, double* %30
%32 = zext i32 %26 to i64
%33 = getelementptr inbounds double, double* %0, i64 %32
%34 = load double, double* %33
%35 = tail call double @llvm.fmuladd.f64(double %31, double %34, double %19)
%36 = add nsw i32 %26, -1
%37 = add nuw nsw i64 %18, 1
%38 = icmp eq i64 %37, %13
br i1 %38, label %39, label %17
39: ; preds = %29, %14
%40 = phi double [ 0.000000e+00, %14 ], [ %35, %29 ]
%41 = add nsw i32 %16, 2
%42 = getelementptr inbounds double, double* %4, i64 %15
store double %40, double* %42
%43 = add nuw nsw i64 %15, 1
%44 = icmp eq i64 %43, %12
br i1 %44, label %45, label %14
45: ; preds = %39, %5
ret void
}