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

```asm
<neiEK>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  add    $0xffffffffffffff80,%rsp
  mov    %rdi,-0x58(%rbp)
  mov    %esi,-0x5c(%rbp)
  mov    %edx,-0x60(%rbp)
  mov    %ecx,-0x64(%rbp)
  mov    %r8d,-0x68(%rbp)
  mov    %r9d,-0x6c(%rbp)
  movsd  %xmm0,-0x78(%rbp)
  mov    0x10(%rbp),%eax
  sub    $0x1,%eax
  mov    %eax,-0x2c(%rbp)
  mov    -0x68(%rbp),%eax
  cltq
  shl    $0x3,%eax
  mov    %eax,%edi
  call   <malloc@plt>
  mov    %rax,-0x10(%rbp)
  movl   $0x0,-0x38(%rbp)
  jmp    L631
L543:
  mov    -0x38(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x5c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,%rdx
  mov    -0x38(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0x10(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x38(%rbp)
L631:
  mov    -0x38(%rbp),%eax
  cmp    -0x68(%rbp),%eax
  jl     L543
  mov    -0x60(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  divsd  -0x78(%rbp),%xmm0
  cvttsd2si %xmm0,%eax
  and    -0x2c(%rbp),%eax
  mov    %eax,-0x28(%rbp)
  mov    -0x64(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  divsd  -0x78(%rbp),%xmm0
  cvttsd2si %xmm0,%eax
  and    -0x2c(%rbp),%eax
  mov    %eax,-0x24(%rbp)
  movl   $0x0,-0x3c(%rbp)
  movl   $0x0,-0x34(%rbp)
  jmp    L441
L317:
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    0x30(%rbp),%rax
  add    %rdx,%rax
  movl   $0x0,(%rax)
  addl   $0x1,-0x34(%rbp)
L441:
  mov    -0x34(%rbp),%eax
  cmp    -0x6c(%rbp),%eax
  jl     L317
  jmp    L63
L570:
  cmpl   $0x0,-0x3c(%rbp)
  je     L627
  mov    -0x3c(%rbp),%eax
  add    %eax,%eax
  jmp    L86
L627:
  mov    $0x1,%eax
L86:
  mov    %eax,-0x20(%rbp)
  mov    -0x28(%rbp),%eax
  sub    -0x3c(%rbp),%eax
  mov    %eax,-0x44(%rbp)
  jmp    L346
L76:
  mov    -0x44(%rbp),%eax
  and    -0x2c(%rbp),%eax
  mov    %eax,-0x1c(%rbp)
  mov    -0x24(%rbp),%eax
  sub    -0x3c(%rbp),%eax
  mov    %eax,-0x40(%rbp)
  jmp    L508
L624:
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x20(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x40(%rbp),%eax
  and    -0x2c(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0x30(%rbp)
  jmp    L191
L217:
  mov    -0x10(%rbp),%rax
  movsd  (%rax),%xmm0
  mov    -0x58(%rbp),%rax
  mov    (%rax),%rdx
  mov    -0x30(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movq   D1(%rip),%xmm1
  andpd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movl   $0x1,-0x38(%rbp)
  jmp    L478
L384:
  mov    -0x38(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x38(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x30(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movq   D1(%rip),%xmm1
  andpd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x8(%rbp),%xmm0
  comisd -0x18(%rbp),%xmm0
  jbe    L377
  movsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x6c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  movsd  -0x18(%rbp),%xmm0
  comisd %xmm1,%xmm0
  ja     L20
L377:
  addl   $0x1,-0x38(%rbp)
L478:
  mov    -0x38(%rbp),%eax
  cmp    -0x68(%rbp),%eax
  jl     L384
  jmp    L67
L20:
  nop
L67:
  movl   $0x0,-0x34(%rbp)
  jmp    L311
L501:
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  comisd -0x18(%rbp),%xmm0
  jb     L566
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    0x30(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%edx
  add    $0x1,%edx
  mov    %edx,(%rax)
L566:
  addl   $0x1,-0x34(%rbp)
L311:
  mov    -0x34(%rbp),%eax
  cmp    -0x6c(%rbp),%eax
  jl     L501
  mov    -0x30(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    0x28(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0x30(%rbp)
L191:
  cmpl   $0xffffffff,-0x30(%rbp)
  jne    L217
  mov    -0x20(%rbp),%eax
  add    %eax,-0x40(%rbp)
L508:
  mov    -0x24(%rbp),%edx
  mov    -0x3c(%rbp),%eax
  add    %edx,%eax
  cmp    %eax,-0x40(%rbp)
  jle    L624
  addl   $0x1,-0x44(%rbp)
L346:
  mov    -0x28(%rbp),%edx
  mov    -0x3c(%rbp),%eax
  add    %edx,%eax
  cmp    %eax,-0x44(%rbp)
  jle    L76
  mov    -0x28(%rbp),%eax
  sub    -0x3c(%rbp),%eax
  mov    %eax,-0x44(%rbp)
  jmp    L514
L245:
  mov    -0x44(%rbp),%eax
  and    -0x2c(%rbp),%eax
  mov    %eax,-0x1c(%rbp)
  mov    -0x24(%rbp),%eax
  sub    -0x3c(%rbp),%eax
  add    $0x1,%eax
  mov    %eax,-0x40(%rbp)
  jmp    L85
L392:
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x20(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x40(%rbp),%eax
  and    -0x2c(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0x30(%rbp)
  jmp    L189
L352:
  mov    -0x10(%rbp),%rax
  movsd  (%rax),%xmm0
  mov    -0x58(%rbp),%rax
  mov    (%rax),%rdx
  mov    -0x30(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movq   D1(%rip),%xmm1
  andpd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movl   $0x1,-0x38(%rbp)
  jmp    L129
L194:
  mov    -0x38(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x38(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x30(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movq   D1(%rip),%xmm1
  andpd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x8(%rbp),%xmm0
  comisd -0x18(%rbp),%xmm0
  jbe    L137
  movsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x6c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  movsd  -0x18(%rbp),%xmm0
  comisd %xmm1,%xmm0
  ja     L79
L137:
  addl   $0x1,-0x38(%rbp)
L129:
  mov    -0x38(%rbp),%eax
  cmp    -0x68(%rbp),%eax
  jl     L194
  jmp    L15
L79:
  nop
L15:
  movl   $0x0,-0x34(%rbp)
  jmp    L258
L259:
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  comisd -0x18(%rbp),%xmm0
  jb     L166
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    0x30(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%edx
  add    $0x1,%edx
  mov    %edx,(%rax)
L166:
  addl   $0x1,-0x34(%rbp)
L258:
  mov    -0x34(%rbp),%eax
  cmp    -0x6c(%rbp),%eax
  jl     L259
  mov    -0x30(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    0x28(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0x30(%rbp)
L189:
  cmpl   $0xffffffff,-0x30(%rbp)
  jne    L352
  addl   $0x1,-0x40(%rbp)
L85:
  mov    -0x24(%rbp),%edx
  mov    -0x3c(%rbp),%eax
  add    %edx,%eax
  cmp    %eax,-0x40(%rbp)
  jl     L392
  mov    -0x20(%rbp),%eax
  add    %eax,-0x44(%rbp)
L514:
  mov    -0x28(%rbp),%edx
  mov    -0x3c(%rbp),%eax
  add    %edx,%eax
  cmp    %eax,-0x44(%rbp)
  jle    L245
  addl   $0x1,-0x3c(%rbp)
  mov    0x10(%rbp),%eax
  mov    %eax,%edx
  shr    $0x1f,%edx
  add    %edx,%eax
  sar    %eax
  cmp    %eax,-0x3c(%rbp)
  je     L609
L63:
  mov    -0x6c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x3c(%rbp),%eax
  sub    $0x1,%eax
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  mulsd  -0x78(%rbp),%xmm1
  comisd %xmm1,%xmm0
  ja     L570
  jmp    L154
L609:
  nop
L154:
  mov    0x10(%rbp),%eax
  mov    %eax,%edx
  shr    $0x1f,%edx
  add    %edx,%eax
  sar    %eax
  cmp    %eax,-0x3c(%rbp)
  jne    L3
  mov    -0x28(%rbp),%eax
  sub    -0x3c(%rbp),%eax
  mov    %eax,-0x44(%rbp)
  jmp    L24
L365:
  mov    -0x44(%rbp),%eax
  and    -0x2c(%rbp),%eax
  mov    %eax,-0x1c(%rbp)
  mov    -0x24(%rbp),%eax
  sub    -0x3c(%rbp),%eax
  mov    %eax,-0x40(%rbp)
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x20(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x40(%rbp),%eax
  and    -0x2c(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0x30(%rbp)
  jmp    L107
L23:
  mov    -0x10(%rbp),%rax
  movsd  (%rax),%xmm0
  mov    -0x58(%rbp),%rax
  mov    (%rax),%rdx
  mov    -0x30(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movq   D1(%rip),%xmm1
  andpd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movl   $0x1,-0x38(%rbp)
  jmp    L35
L614:
  mov    -0x38(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x38(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x30(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movq   D1(%rip),%xmm1
  andpd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x8(%rbp),%xmm0
  comisd -0x18(%rbp),%xmm0
  jbe    L176
  movsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x6c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  movsd  -0x18(%rbp),%xmm0
  comisd %xmm1,%xmm0
  ja     L597
L176:
  addl   $0x1,-0x38(%rbp)
L35:
  mov    -0x38(%rbp),%eax
  cmp    -0x68(%rbp),%eax
  jl     L614
  jmp    L309
L597:
  nop
L309:
  movl   $0x0,-0x34(%rbp)
  jmp    L328
L142:
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  comisd -0x18(%rbp),%xmm0
  jb     L17
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    0x30(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%edx
  add    $0x1,%edx
  mov    %edx,(%rax)
L17:
  addl   $0x1,-0x34(%rbp)
L328:
  mov    -0x34(%rbp),%eax
  cmp    -0x6c(%rbp),%eax
  jl     L142
  mov    -0x30(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    0x28(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0x30(%rbp)
L107:
  cmpl   $0xffffffff,-0x30(%rbp)
  jne    L23
  addl   $0x1,-0x44(%rbp)
L24:
  mov    -0x28(%rbp),%edx
  mov    -0x3c(%rbp),%eax
  add    %edx,%eax
  cmp    %eax,-0x44(%rbp)
  jl     L365
  mov    -0x28(%rbp),%eax
  sub    -0x3c(%rbp),%eax
  mov    %eax,-0x44(%rbp)
  mov    -0x44(%rbp),%eax
  and    -0x2c(%rbp),%eax
  mov    %eax,-0x1c(%rbp)
  mov    -0x24(%rbp),%eax
  sub    -0x3c(%rbp),%eax
  add    $0x1,%eax
  mov    %eax,-0x40(%rbp)
  jmp    L163
L512:
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x20(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x40(%rbp),%eax
  and    -0x2c(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0x30(%rbp)
  jmp    L633
L548:
  mov    -0x10(%rbp),%rax
  movsd  (%rax),%xmm0
  mov    -0x58(%rbp),%rax
  mov    (%rax),%rdx
  mov    -0x30(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movq   D1(%rip),%xmm1
  andpd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movl   $0x1,-0x38(%rbp)
  jmp    L134
L622:
  mov    -0x38(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x38(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x30(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movq   D1(%rip),%xmm1
  andpd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x8(%rbp),%xmm0
  comisd -0x18(%rbp),%xmm0
  jbe    L489
  movsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x6c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  movsd  -0x18(%rbp),%xmm0
  comisd %xmm1,%xmm0
  ja     L407
L489:
  addl   $0x1,-0x38(%rbp)
L134:
  mov    -0x38(%rbp),%eax
  cmp    -0x68(%rbp),%eax
  jl     L622
  jmp    L603
L407:
  nop
L603:
  movl   $0x0,-0x34(%rbp)
  jmp    L498
L221:
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  comisd -0x18(%rbp),%xmm0
  jb     L535
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    0x30(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%edx
  add    $0x1,%edx
  mov    %edx,(%rax)
L535:
  addl   $0x1,-0x34(%rbp)
L498:
  mov    -0x34(%rbp),%eax
  cmp    -0x6c(%rbp),%eax
  jl     L221
  mov    -0x30(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    0x28(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0x30(%rbp)
L633:
  cmpl   $0xffffffff,-0x30(%rbp)
  jne    L548
  addl   $0x1,-0x40(%rbp)
L163:
  mov    -0x24(%rbp),%edx
  mov    -0x3c(%rbp),%eax
  add    %edx,%eax
  cmp    %eax,-0x40(%rbp)
  jl     L512
L3:
  mov    -0x10(%rbp),%rax
  mov    %rax,%rdi
  call   <free@plt>
  movl   $0x0,-0x34(%rbp)
  jmp    L136
L436:
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    0x30(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%edx
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rcx
  mov    0x30(%rbp),%rax
  add    %rcx,%rax
  sub    $0x1,%edx
  mov    %edx,(%rax)
  addl   $0x1,-0x34(%rbp)
L136:
  mov    -0x34(%rbp),%eax
  cmp    -0x6c(%rbp),%eax
  jl     L436
  nop
  nop
  leave
  ret
```

tool_calls: name: extract_symbolic_data
arguments: tag: D1
type: qword

tool: D1: 9223372036854775807

assistant: ```c
void neiEK(double **x, int i, int comp1, int comp2, int dim, int K, int bs, double epsgrid, double *eps, int **box, int *lis, int *nx) {
  int ix, iy, ix1, iy1, ix2, jj, step, d, ik;
  int el, ib = bs - 1;
  double dd, dy;
  double *xx;
  xx = (double *)malloc(dim * sizeof(double));
  for (d = 0; d < dim; d++)
    xx[d] = x[d][i];
  ix = (int)(xx[comp1] / epsgrid) & ib;
  iy = (int)(xx[comp2] / epsgrid) & ib;
  jj = 0;
  for (ik = 0; ik < K; ik++)
    nx[ik] = 0;
  while (eps[K - 1] > epsgrid * (jj - 1)) {
    step = (jj) ? 2 * jj : 1;
    for (ix1 = ix - jj; ix1 <= ix + jj; ix1++) {
      ix2 = ix1 & ib;
      for (iy1 = iy - jj; iy1 <= iy + jj; iy1 += step) {
        el = box[ix2][iy1 & ib];
        while (el != -1) {
          dd = fabs(xx[0] - x[0][el]);
          for (d = 1; d < dim; d++)
            if ((dy = fabs(xx[d] - x[d][el])) > dd)
              if ((dd = dy) > eps[K - 1])
                break;
          ;
          for (ik = 0; ik < K; ik++) {
            if (dd <= eps[ik])
              nx[ik]++;
          }
          el = lis[el];
        }
      }
    }
    for (ix1 = ix - jj; ix1 <= ix + jj; ix1 += step) {
      ix2 = ix1 & ib;
      for (iy1 = iy - jj + 1; iy1 <= iy + jj - 1; iy1++) {
        el = box[ix2][iy1 & ib];
        while (el != -1) {
          dd = fabs(xx[0] - x[0][el]);
          for (d = 1; d < dim; d++)
            if ((dy = fabs(xx[d] - x[d][el])) > dd)
              if ((dd = dy) > eps[K - 1])
                break;
          for (ik = 0; ik < K; ik++) {
            if (dd <= eps[ik])
              nx[ik]++;
          }
          el = lis[el];
        }
      }
    }
    jj++;
    if (jj == (bs / 2))
      break;
  }
  if (jj == (bs / 2)) {
    for (ix1 = ix - jj; ix1 < ix + jj; ix1++) {
      ix2 = ix1 & ib;
      iy1 = iy - jj;
      el = box[ix2][iy1 & ib];
      while (el != -1) {
        dd = fabs(xx[0] - x[0][el]);
        for (d = 1; d < dim; d++)
          if ((dy = fabs(xx[d] - x[d][el])) > dd)
            if ((dd = dy) > eps[K - 1])
              break;
        for (ik = 0; ik < K; ik++) {
          if (dd <= eps[ik])
            nx[ik]++;
        }
        el = lis[el];
      }
    }
    ix1 = ix - jj;
    ix2 = ix1 & ib;
    for (iy1 = iy - jj + 1; iy1 <= iy + jj - 1; iy1++) {
      el = box[ix2][iy1 & ib];
      while (el != -1) {
        dd = fabs(xx[0] - x[0][el]);
        for (d = 1; d < dim; d++)
          if ((dy = fabs(xx[d] - x[d][el])) > dd)
            if ((dd = dy) > eps[K - 1])
              break;
        for (ik = 0; ik < K; ik++) {
          if (dd <= eps[ik])
            nx[ik]++;
        }
        el = lis[el];
      }
    }
  }
  free(xx);
  for (ik = 0; ik < K; ik++)
    nx[ik] -= 1;
}
```