void taper_norm(g, dt, nt, tap_per) float *g, *dt, *tap_per;
int nt;
{
  float fac, df, arg;
  int i;
  int ntap;
  ntap = nt * (*tap_per);
  for (i = 0; i < nt - ntap; i++)
    g[i] = g[i] * (*dt);
  if (ntap > 0) {
    df = 3.14159 / (float)(ntap);
    for (i = nt - ntap; i < nt; i++) {
      arg = (i - (nt - (ntap + 1))) * df;
      fac = (*dt) * 0.5 * (1.0 + cos(arg));
      g[i] = g[i] * fac;
    }
  }
}


x86 assembly (clang, -O0):
<taper_norm>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
cvtsi2ssl -0x14(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mulss  (%rax),%xmm0
cvttss2si %xmm0,%eax
mov    %eax,-0x34(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    -0x34(%rbp),%ecx
cmp    %ecx,%eax
jge    1183 <taper_norm+0x73>
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mulss  (%rax),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1142 <taper_norm+0x32>
cmpl   $0x0,-0x34(%rbp)
jle    125c <taper_norm+0x14c>
cvtsi2ssl -0x34(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe62(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x28(%rbp)
mov    -0x14(%rbp),%eax
sub    -0x34(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1257 <taper_norm+0x147>
mov    -0x30(%rbp),%eax
mov    -0x14(%rbp),%ecx
mov    -0x34(%rbp),%edx
add    $0x1,%edx
sub    %edx,%ecx
sub    %ecx,%eax
cvtsi2ss %eax,%xmm0
mulss  -0x28(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe1e(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movss  -0x2c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
call   1030 <cos@plt>
movaps %xmm0,%xmm2
movsd  -0x40(%rbp),%xmm0
movsd  0xdef(%rip),%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x24(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mulss  -0x24(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11b4 <taper_norm+0xa4>
jmp    125c <taper_norm+0x14c>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<taper_norm>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x10,%rsp
cvtsi2ss %edx,%xmm0
mov    %rsi,%r14
mov    %rdi,%r15
mulss  (%rcx),%xmm0
cvttss2si %xmm0,%r8d
mov    %edx,%eax
sub    %r8d,%eax
test   %eax,%eax
jle    12b8 <taper_norm+0x1a8>
mov    %eax,%ebx
cmp    $0x8,%eax
jb     115d <taper_norm+0x4d>
lea    0x4(%r14),%rcx
cmp    %r15,%rcx
jbe    1208 <taper_norm+0xf8>
lea    (%r15,%rbx,4),%rcx
cmp    %r14,%rcx
jbe    1208 <taper_norm+0xf8>
xor    %edi,%edi
mov    %rdi,%rcx
not    %rcx
add    %rbx,%rcx
mov    %rbx,%rsi
and    $0x3,%rsi
je     119b <taper_norm+0x8b>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movss  (%r15,%rdi,4),%xmm0
mulss  (%r14),%xmm0
movss  %xmm0,(%r15,%rdi,4)
add    $0x1,%rdi
add    $0xffffffffffffffff,%rsi
jne    1180 <taper_norm+0x70>
cmp    $0x3,%rcx
jb     12b8 <taper_norm+0x1a8>
cs nopw 0x0(%rax,%rax,1)
nop
movss  (%r15,%rdi,4),%xmm0
mulss  (%r14),%xmm0
movss  0x4(%r15,%rdi,4),%xmm1
movss  %xmm0,(%r15,%rdi,4)
mulss  (%r14),%xmm1
movss  %xmm1,0x4(%r15,%rdi,4)
movss  0x8(%r15,%rdi,4),%xmm0
mulss  (%r14),%xmm0
movss  %xmm0,0x8(%r15,%rdi,4)
movss  0xc(%r15,%rdi,4),%xmm0
mulss  (%r14),%xmm0
movss  %xmm0,0xc(%r15,%rdi,4)
add    $0x4,%rdi
cmp    %rdi,%rbx
jne    11b0 <taper_norm+0xa0>
jmp    12b8 <taper_norm+0x1a8>
mov    %ebx,%edi
and    $0xfffffff8,%edi
lea    -0x8(%rdi),%rcx
mov    %rcx,%r9
shr    $0x3,%r9
add    $0x1,%r9
test   %rcx,%rcx
je     1354 <taper_norm+0x244>
mov    %r9,%rbp
and    $0xfffffffffffffffe,%rbp
movss  (%r14),%xmm0
shufps $0x0,%xmm0,%xmm0
xor    %ecx,%ecx
nopw   0x0(%rax,%rax,1)
movups (%r15,%rcx,4),%xmm1
movups 0x10(%r15,%rcx,4),%xmm2
movups 0x20(%r15,%rcx,4),%xmm3
movups 0x30(%r15,%rcx,4),%xmm4
mulps  %xmm0,%xmm1
mulps  %xmm0,%xmm2
movups %xmm1,(%r15,%rcx,4)
movups %xmm2,0x10(%r15,%rcx,4)
mulps  %xmm0,%xmm3
mulps  %xmm0,%xmm4
movups %xmm3,0x20(%r15,%rcx,4)
movups %xmm4,0x30(%r15,%rcx,4)
add    $0x10,%rcx
add    $0xfffffffffffffffe,%rbp
jne    1240 <taper_norm+0x130>
test   $0x1,%r9b
je     12af <taper_norm+0x19f>
movups (%r15,%rcx,4),%xmm0
movups 0x10(%r15,%rcx,4),%xmm1
movss  (%r14),%xmm2
shufps $0x0,%xmm2,%xmm2
mulps  %xmm2,%xmm0
mulps  %xmm1,%xmm2
movups %xmm0,(%r15,%rcx,4)
movups %xmm2,0x10(%r15,%rcx,4)
cmp    %rbx,%rdi
jne    115f <taper_norm+0x4f>
test   %r8d,%r8d
jle    1347 <taper_norm+0x237>
xorps  %xmm0,%xmm0
cvtsi2sd %r8d,%xmm0
movsd  0xd2f(%rip),%xmm1
divsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
movss  %xmm0,0x4(%rsp)
movslq %eax,%rbx
movslq %edx,%r12
mov    $0x1,%ebp
nopl   (%rax)
xorps  %xmm0,%xmm0
cvtsi2ss %ebp,%xmm0
mulss  0x4(%rsp),%xmm0
movss  (%r14),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  0xcfa(%rip),%xmm1
movsd  %xmm1,0x8(%rsp)
cvtss2sd %xmm0,%xmm0
call   1030 <cos@plt>
addsd  0xceb(%rip),%xmm0
mulsd  0x8(%rsp),%xmm0
cvtsd2ss %xmm0,%xmm0
mulss  (%r15,%rbx,4),%xmm0
movss  %xmm0,(%r15,%rbx,4)
add    $0x1,%rbx
add    $0x1,%ebp
cmp    %r12,%rbx
jl     12f0 <taper_norm+0x1e0>
add    $0x10,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %ecx,%ecx
test   $0x1,%r9b
jne    128a <taper_norm+0x17a>
jmp    12af <taper_norm+0x19f>

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

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