user: Demonstrate the C code that generates the assembly output below:

```asm
<feat_postings>:
  endbr64
  push   %r15
  mov    %r8,%r15
  push   %r14
  movslq %edx,%r14
  push   %r13
  lea    0x0(,%r14,4),%edx
  mov    %r14,%r13
  push   %r12
  mov    %edi,%r12d
  mov    %rcx,%rdi
  push   %rbp
  mov    %rcx,%rbp
  push   %rbx
  mov    %rsi,%rbx
  xor    %esi,%esi
  sub    $0x18,%rsp
  movss  %xmm0,0xc(%rsp)
  call   <memset@plt>
  test   %r12d,%r12d
  jle    L91
  test   %r14d,%r14d
  jle    L25
  movss  0xc(%rsp),%xmm0
  lea    0x0(,%r14,4),%r9
  mov    %rbx,%rcx
  xor    %r8d,%r8d
  nopw   0x0(%rax,%rax,1)
L40:
  xor    %eax,%eax
  nopw   0x0(%rax,%rax,1)
L21:
  movss  (%rcx,%rax,4),%xmm1
  comiss %xmm0,%xmm1
  jbe    L71
  movslq 0x0(%rbp,%rax,4),%rdx
  mov    (%r15,%rax,8),%rsi
  lea    0x1(%rdx),%edi
  mov    %edi,0x0(%rbp,%rax,4)
  mov    %r8d,(%rsi,%rdx,4)
L71:
  add    $0x1,%rax
  cmp    %r14,%rax
  jne    L21
  add    $0x1,%r8d
  add    %r9,%rcx
  cmp    %r8d,%r12d
  jne    L40
L70:
  lea    -0x1(%r13),%eax
  cmp    $0x2,%eax
  jbe    L34
  mov    %r13d,%edx
  pxor   %xmm3,%xmm3
  pcmpeqd %xmm5,%xmm5
  mov    %rbp,%rax
  shr    $0x2,%edx
  movdqa %xmm3,%xmm4
  sub    $0x1,%edx
  shl    $0x4,%rdx
  lea    0x10(%rbp,%rdx,1),%rdx
  nopl   0x0(%rax)
L46:
  movdqu (%rax),%xmm2
  add    $0x10,%rax
  movdqa %xmm2,%xmm1
  paddd  %xmm5,%xmm1
  movdqa %xmm1,%xmm0
  psrlq  $0x20,%xmm1
  pmuludq %xmm2,%xmm0
  psrlq  $0x20,%xmm2
  pmuludq %xmm2,%xmm1
  pshufd $0x8,%xmm0,%xmm0
  pshufd $0x8,%xmm1,%xmm1
  punpckldq %xmm1,%xmm0
  movdqa %xmm4,%xmm1
  pcmpgtd %xmm0,%xmm1
  psubd  %xmm1,%xmm0
  psrad  $0x1,%xmm0
  paddd  %xmm0,%xmm3
  cmp    %rax,%rdx
  jne    L46
  movdqa %xmm3,%xmm0
  mov    %r13d,%ecx
  psrldq $0x8,%xmm0
  and    $0xfffffffc,%ecx
  paddd  %xmm0,%xmm3
  movdqa %xmm3,%xmm0
  psrldq $0x4,%xmm0
  paddd  %xmm0,%xmm3
  movd   %xmm3,%eax
  test   $0x3,%r13b
  je     L56
L48:
  movslq %ecx,%rdx
  lea    0x0(,%rdx,4),%rdi
  mov    0x0(%rbp,%rdx,4),%edx
  lea    -0x1(%rdx),%esi
  imul   %edx,%esi
  mov    %esi,%edx
  shr    $0x1f,%edx
  add    %esi,%edx
  sar    %edx
  add    %edx,%eax
  lea    0x1(%rcx),%edx
  cmp    %edx,%r13d
  jle    L56
  mov    0x4(%rbp,%rdi,1),%edx
  add    $0x2,%ecx
  lea    -0x1(%rdx),%esi
  imul   %edx,%esi
  mov    %esi,%edx
  shr    $0x1f,%edx
  add    %esi,%edx
  sar    %edx
  add    %edx,%eax
  cmp    %ecx,%r13d
  jle    L56
  mov    0x8(%rbp,%rdi,1),%edx
  lea    -0x1(%rdx),%ecx
  imul   %edx,%ecx
  mov    %ecx,%edx
  shr    $0x1f,%edx
  add    %ecx,%edx
  sar    %edx
  add    %edx,%eax
L56:
  add    $0x18,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L91:
  test   %r14d,%r14d
  jg     L70
L25:
  add    $0x18,%rsp
  xor    %eax,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L34:
  xor    %ecx,%ecx
  xor    %eax,%eax
  jmp    L48
```

assistant: ```c
int feat_postings(int N, float *feats, int D, float pthr, int *cntlist, int **occlist) {
  memset((void *)cntlist, 0, sizeof(*cntlist) * D);
  for (int n = 0; n < N; n++)
    for (int d = 0; d < D; d++)
      if (feats[n * D + d] > pthr)
        occlist[d][cntlist[d]++] = n;
  int maxdots = 0;
  for (int d = 0; d < D; d++)
    maxdots += cntlist[d] * (cntlist[d] - 1) / 2;
  return maxdots;
}
```