// File: /workspace/github_c_makefile_star2_O0/YunyiQian[P]multitel3/fft.c
// Function: filter
void	filter(complex *d, int n, float f1, float f2, float dt, int sgn) {
     int i, if1, if2;
     float a;
     dt = 0.5/dt/n;
     if1 = rint(f1/dt);
     if2 = rint(f2/dt); if (if2>=n) if2=n-1;
     if (if2<=if1) {
	fprintf(stderr, "filter freq. wrong f2<f1\n");
	return;
     }
     dt = PI/(if2-if1);
     for(a=0.,i=if1;i<if2;i++,a+=dt) d[i] = dmltp(0.5*(1-sgn*cos(a)), d[i]);
     if (sgn<0) {
	for(i=if2;i<n;i++) d[i] = Zero;
	d[0].y = 0.;
     } else  {
        d[0].x = 0.;
        for(i=1;i<if1;i++) d[i] = Zero;
     }
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0x80, %rsp
movq %rdi, -0x8(%rbp)
movl %esi, -0xc(%rbp)
movss %xmm0, -0x10(%rbp)
movss %xmm1, -0x14(%rbp)
movss %xmm2, -0x18(%rbp)
movl %edx, -0x1c(%rbp)
movss -0x18(%rbp), %xmm0
cvtss2sd %xmm0, %xmm1
movsd 0x4ee3(%rip), %xmm0     # 0x89c8
divsd %xmm1, %xmm0
cvtsi2sdl -0xc(%rbp), %xmm1
divsd %xmm1, %xmm0
cvtsd2ss %xmm0, %xmm0
movss %xmm0, -0x18(%rbp)
movss -0x10(%rbp), %xmm0
divss -0x18(%rbp), %xmm0
cvtss2sd %xmm0, %xmm0
callq 0x1120
cvttsd2si %xmm0, %eax
movl %eax, -0x24(%rbp)
movss -0x14(%rbp), %xmm0
divss -0x18(%rbp), %xmm0
cvtss2sd %xmm0, %xmm0
callq 0x1120
cvttsd2si %xmm0, %eax
movl %eax, -0x28(%rbp)
movl -0x28(%rbp), %eax
cmpl -0xc(%rbp), %eax
jl 0x3b40
movl -0xc(%rbp), %eax
subl $0x1, %eax
movl %eax, -0x28(%rbp)
movl -0x28(%rbp), %eax
cmpl -0x24(%rbp), %eax
jg 0x3b65
movq 0x74a9(%rip), %rax      # 0xaff8
movq (%rax), %rdi
leaq 0x4ed1(%rip), %rsi      # 0x8a2a
movb $0x0, %al
callq 0x1100
jmp 0x3ce3
movl -0x28(%rbp), %eax
subl -0x24(%rbp), %eax
cvtsi2sd %eax, %xmm1
movsd 0x4e39(%rip), %xmm0     # 0x89b0
divsd %xmm1, %xmm0
cvtsd2ss %xmm0, %xmm0
movss %xmm0, -0x18(%rbp)
xorps %xmm0, %xmm0
movss %xmm0, -0x2c(%rbp)
movl -0x24(%rbp), %eax
movl %eax, -0x20(%rbp)
movl -0x20(%rbp), %eax
cmpl -0x28(%rbp), %eax
jge 0x3c2f
movq -0x8(%rbp), %rax
movq %rax, -0x50(%rbp)
movslq -0x20(%rbp), %rax
movq %rax, -0x58(%rbp)
cvtsi2sdl -0x1c(%rbp), %xmm0
movsd %xmm0, -0x60(%rbp)
movss -0x2c(%rbp), %xmm0
cvtss2sd %xmm0, %xmm0
callq 0x10d0
movsd -0x60(%rbp), %xmm1
mulsd %xmm0, %xmm1
movsd 0x4471(%rip), %xmm0     # 0x8048
subsd %xmm1, %xmm0
movsd 0x4de5(%rip), %xmm1     # 0x89c8
mulsd %xmm1, %xmm0
cvtsd2ss %xmm0, %xmm0
movq -0x8(%rbp), %rax
movslq -0x20(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm1
callq 0x4750
movq -0x58(%rbp), %rcx
movq -0x50(%rbp), %rax
movlpd %xmm0, -0x34(%rbp)
movq -0x34(%rbp), %rdx
movq %rdx, (%rax,%rcx,8)
movl -0x20(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x20(%rbp)
movss -0x18(%rbp), %xmm0
addss -0x2c(%rbp), %xmm0
movss %xmm0, -0x2c(%rbp)
jmp 0x3b92
cmpl $0x0, -0x1c(%rbp)
jge 0x3c8c
movl -0x28(%rbp), %eax
movl %eax, -0x20(%rbp)
movl -0x20(%rbp), %eax
cmpl -0xc(%rbp), %eax
jge 0x3c7e
movq -0x8(%rbp), %rax
movq %rax, -0x68(%rbp)
movslq -0x20(%rbp), %rax
movq %rax, -0x70(%rbp)
xorps %xmm1, %xmm1
movaps %xmm1, %xmm0
callq 0x48c0
movq -0x70(%rbp), %rcx
movq -0x68(%rbp), %rax
movlpd %xmm0, -0x3c(%rbp)
movq -0x3c(%rbp), %rdx
movq %rdx, (%rax,%rcx,8)
movl -0x20(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x20(%rbp)
jmp 0x3c3b
movq -0x8(%rbp), %rax
xorps %xmm0, %xmm0
movss %xmm0, 0x4(%rax)
jmp 0x3ce3
movq -0x8(%rbp), %rax
xorps %xmm0, %xmm0
movss %xmm0, (%rax)
movl $0x1, -0x20(%rbp)
movl -0x20(%rbp), %eax
cmpl -0x24(%rbp), %eax
jge 0x3ce1
movq -0x8(%rbp), %rax
movq %rax, -0x78(%rbp)
movslq -0x20(%rbp), %rax
movq %rax, -0x80(%rbp)
xorps %xmm1, %xmm1
movaps %xmm1, %xmm0
callq 0x48c0
movq -0x80(%rbp), %rcx
movq -0x78(%rbp), %rax
movlpd %xmm0, -0x44(%rbp)
movq -0x44(%rbp), %rdx
movq %rdx, (%rax,%rcx,8)
movl -0x20(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x20(%rbp)
jmp 0x3c9e
jmp 0x3ce3
addq $0x80, %rsp
popq %rbp
retq
nopl (%rax)