user: Please show me the C implementation that matches the assembly code listed below:

```asm
<nn_single_full>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0xc8,%rsp
  mov    %edi,0x48(%rsp)
  mov    %esi,%r14d
  mov    %esi,0xb0(%rsp)
  mov    %edx,%r12d
  mov    %edx,0xb4(%rsp)
  mov    %ecx,%r13d
  mov    %ecx,0x98(%rsp)
  mov    %r8,0xb8(%rsp)
  mov    %r9,%rbp
  mov    $0x400,%eax
  cmp    %eax,%esi
  mov    %eax,%r15d
  cmovle %esi,%r15d
  mov    %r15d,0x9c(%rsp)
  cmp    %eax,%edx
  cmovle %edx,%eax
  mov    %eax,%ebx
  mov    %r15d,%edi
  imul   %eax,%edi
  call   <fvec_new@plt>
  mov    %rax,0x70(%rsp)
  movslq %r14d,%rax
  mov    %rax,0x90(%rsp)
  test   %rax,%rax
  jle    L218
  movslq %r15d,%rdi
  mov    %rdi,0x80(%rsp)
  mov    %rax,0x68(%rsp)
  movslq %r13d,%rax
  shl    $0x2,%rax
  mov    %rdi,%rsi
  imul   %rax,%rsi
  mov    %rsi,0xa0(%rsp)
  mov    %rbp,0x88(%rsp)
  shl    $0x2,%rdi
  mov    %rdi,0xa8(%rsp)
  mov    0x108(%rsp),%r13
  mov    0x110(%rsp),%r15
  movslq %ebx,%rdx
  imul   %rdx,%rax
  mov    %rax,0x50(%rsp)
  lea    0x0(,%rdx,4),%rax
  mov    %rax,0x58(%rsp)
  movq   $0x0,0x78(%rsp)
  movq   $0x0,0x60(%rsp)
  movslq %r12d,%rax
  mov    %rax,0x38(%rsp)
  mov    %rdx,0x18(%rsp)
  jmp    L120
L169:
  mov    %rbx,%r9
  mov    0x40(%rsp),%r8
  mov    0x30(%rsp),%rcx
  mov    0x4c(%rsp),%edx
  mov    %eax,%esi
  mov    0x98(%rsp),%edi
  call   <compute_cross_distances@plt>
  jmp    L40
L136:
  mov    $0x0,%esi
  movslq 0xc(%rsp),%rcx
  test   %r14,%r14
  jg     L240
L146:
  mov    0x18(%rsp),%rdi
  add    %rdi,%rbp
  sub    %rdi,0x10(%rsp)
  mov    0x50(%rsp),%rdi
  add    %rdi,0x30(%rsp)
  add    %rdi,0x28(%rsp)
  add    0x58(%rsp),%r12
  cmp    0x38(%rsp),%rbp
  jge    L217
L183:
  mov    0x10(%rsp),%rdi
  mov    0x18(%rsp),%rax
  cmp    %rax,%rdi
  cmovle %rdi,%rax
  mov    %rax,0x20(%rsp)
  mov    %eax,0xc(%rsp)
  cmpl   $0x2,0x48(%rsp)
  je     L169
  sub    $0x8,%rsp
  push   %rbx
  mov    0x50(%rsp),%r9
  mov    0x38(%rsp),%r8
  mov    0x5c(%rsp),%ecx
  mov    0x1c(%rsp),%edx
  mov    0xa8(%rsp),%esi
  mov    0x58(%rsp),%edi
  call   <compute_cross_distances_alt@plt>
  add    $0x10,%rsp
L40:
  cmpq   $0x0,0x100(%rsp)
  jne    L136
  test   %r14,%r14
  jle    L146
L234:
  movslq 0x20(%rsp),%r11
  shl    $0x2,%r11
  mov    %rbx,%rdx
  mov    $0x0,%esi
  movslq 0xc(%rsp),%rcx
  jmp    L14
L2:
  mov    %rcx,%rax
  imul   %rsi,%rax
  lea    (%rbx,%rax,4),%rdx
  mov    $0x0,%eax
L30:
  movss  (%rdx,%rax,4),%xmm0
  mulss  (%r12,%rax,4),%xmm0
  movss  %xmm0,(%rdx,%rax,4)
  add    $0x1,%rax
  cmp    %rcx,%rax
  jne    L30
L211:
  add    $0x1,%rsi
  cmp    %r14,%rsi
  jge    L234
L240:
  test   %rcx,%rcx
  jg     L2
  jmp    L211
L70:
  add    $0x1,%rax
  cmp    %rax,%rcx
  je     L148
L238:
  movss  (%rdx,%rax,4),%xmm0
  comiss %xmm0,%xmm1
  jbe    L70
  lea    (%rax,%rbp,1),%edi
  movaps %xmm0,%xmm1
  jmp    L70
L148:
  mov    %edi,(%r10,%r8,4)
  movss  %xmm1,(%r9,%r8,4)
  add    $0x1,%rsi
  add    %r11,%rdx
  cmp    %r14,%rsi
  jge    L146
L14:
  mov    %r13,%r10
  mov    %rsi,%r8
  mov    0x0(%r13,%rsi,4),%edi
  mov    %r15,%r9
  movss  (%r15,%rsi,4),%xmm1
  test   %rcx,%rcx
  jle    L148
  mov    $0x0,%eax
  jmp    L238
L217:
  cmpq   $0x0,0x118(%rsp)
  je     L257
  pxor   %xmm0,%xmm0
  cvtsi2sdq 0x60(%rsp),%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2sdl 0xb0(%rsp),%xmm1
  divsd  %xmm1,%xmm0
  mov    0x120(%rsp),%rdi
  call   *0x118(%rsp)
L257:
  mov    0x80(%rsp),%rax
  add    %rax,0x60(%rsp)
  mov    0x60(%rsp),%rbx
  add    %rax,0x78(%rsp)
  movslq 0x9c(%rsp),%rax
  sub    %rax,0x68(%rsp)
  mov    0xa0(%rsp),%rdi
  add    %rdi,0x88(%rsp)
  mov    0xa8(%rsp),%rax
  add    %rax,%r13
  add    %rax,%r15
  cmp    %rbx,0x90(%rsp)
  jle    L218
L120:
  mov    0x68(%rsp),%rax
  mov    0x80(%rsp),%rdx
  cmp    %rdx,%rax
  cmovle %rax,%rdx
  mov    %edx,0x4c(%rsp)
  movslq %edx,%r14
  test   %r14,%r14
  jle    L116
  mov    0x78(%rsp),%rbx
  lea    0x0(,%rbx,4),%rax
  mov    %r14,%rdx
  add    %rbx,%rdx
  shl    $0x2,%rdx
  mov    0x108(%rsp),%rcx
  mov    0x110(%rsp),%rsi
L84:
  movl   $0xffffffff,(%rcx,%rax,1)
  movss  D1(%rip),%xmm2
  movss  %xmm2,(%rsi,%rax,1)
  add    $0x4,%rax
  cmp    %rdx,%rax
  jne    L84
L116:
  cmpq   $0x0,0x38(%rsp)
  jle    L217
  mov    0x88(%rsp),%rax
  mov    %rax,0x40(%rsp)
  movslq 0xb4(%rsp),%rax
  mov    %rax,0x10(%rsp)
  mov    0xb8(%rsp),%rax
  mov    %rax,0x30(%rsp)
  mov    0x100(%rsp),%r12
  mov    %rax,0x28(%rsp)
  mov    $0x0,%ebp
  mov    0x70(%rsp),%rbx
  jmp    L183
L218:
  mov    0x70(%rsp),%rdi
  call   <free@plt>
  cmpq   $0x0,0x118(%rsp)
  je     L49
  movsd  D3(%rip),%xmm0
  mov    0x120(%rsp),%rdi
  call   *0x118(%rsp)
L49:
  add    $0xc8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: extract_symbolic_data
arguments: data_tag: D1
data_type: f32
name: extract_symbolic_data
arguments: data_tag: D3
data_type: byte[8]

tool: D1: 1e+30

tool: {"D3": [0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0x3f]}

assistant: ```c
void nn_single_full(int distance_type, int n1, int n2, int d, const float *mat2, const float *mat1, const float *vw_weights, int *vw, float *vwdis, void (*peek_fun)(void *arg, double frac), void *peek_arg) {
  int step1 = ((n1) < (1024) ? (n1) : (1024)), step2 = ((n2) < (1024) ? (n2) : (1024)), m1, m2, imin;
  long i1, i2, j1, j2;
  float dmin;
  float *dists = fvec_new(step1 * step2), *dline;
  for (i1 = 0; i1 < n1; i1 += step1) {
    m1 = ((step1) < (n1 - i1) ? (step1) : (n1 - i1));
    for (j1 = 0; j1 < m1; j1++) {
      vw[j1 + i1] = -1;
      vwdis[j1 + i1] = 1e30;
    }
    for (i2 = 0; i2 < n2; i2 += step2) {
      m2 = ((step2) < (n2 - i2) ? (step2) : (n2 - i2));
      if (distance_type == 2) {
        compute_cross_distances(d, m2, m1, mat2 + i2 * d, mat1 + i1 * d, dists);
      } else {
        compute_cross_distances_alt(distance_type, d, m2, m1, mat2 + i2 * d, mat1 + i1 * d, dists);
      }
      if (vw_weights) {
        for (j1 = 0; j1 < m1; j1++) {
          for (j2 = 0; j2 < m2; j2++) {
            dists[j1 * m2 + j2] *= vw_weights[j2 + i2];
          }
        }
      }
      for (j1 = 0; j1 < m1; j1++) {
        dline = dists + j1 * m2;
        imin = vw[i1 + j1];
        dmin = vwdis[i1 + j1];
        for (j2 = 0; j2 < m2; j2++) {
          if (dline[j2] < dmin) {
            imin = j2 + i2;
            dmin = dline[j2];
          }
        }
        vw[i1 + j1] = imin;
        vwdis[i1 + j1] = dmin;
      }
    }
    if (peek_fun) {
      (*peek_fun)(peek_arg, i1 / (double)n1);
    }
  }
  free(dists);
  if (peek_fun) {
    (*peek_fun)(peek_arg, 1);
  }
}
```