user: 请显示与下面列出的汇编代码相匹配的C语言实现：

```asm
<nullspace>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %rdi,%rbp
  mov    %rcx,-0x10(%rsp)
  movslq %edx,%r13
  lea    (%rcx,%r13,4),%rdi
  mov    %rdi,-0x18(%rsp)
  test   %edx,%edx
  jle    L113
  mov    %edx,%eax
  mov    %rcx,%rbx
  lea    -0x1(%rdx),%r15d
  lea    0x0(,%r13,8),%rcx
  mov    %rbx,-0x38(%rsp)
  mov    $0x0,%r14d
  mov    $0x0,%r9d
  mov    $0x0,%r10d
  movl   $0x0,-0x24(%rsp)
  movq   D25(%rip),%xmm6
  movapd %xmm6,%xmm3
  movsd  D29(%rip),%xmm5
  mov    %r13,-0x20(%rsp)
  lea    0x8(%rbp),%rbx
  mov    %rbx,-0x48(%rsp)
  movq   D12(%rip),%xmm4
  mov    %rsi,-0x8(%rsp)
  jmp    L32
L60:
  mov    %r10d,%esi
  jmp    L208
L215:
  lea    0x4(%rdi),%rdx
  mov    %rdx,-0x30(%rsp)
  mov    %r9d,(%rdi)
  cmp    %r9d,%eax
  jle    L269
  imul   %eax,%esi
  movslq %r8d,%r11
  lea    (%r11,%r14,1),%rdx
  lea    0x0(%rbp,%rdx,8),%rdi
  mov    %r15d,%edx
  mov    %r11,-0x40(%rsp)
  add    %r11,%rdx
  add    %r14,%rdx
  mov    -0x48(%rsp),%r11
  lea    (%r11,%rdx,8),%r11
  mov    %rdi,%rdx
  movslq %esi,%rsi
  sub    -0x40(%rsp),%rsi
L259:
  movsd  (%rdx),%xmm0
  movsd  (%rdx,%rsi,8),%xmm1
  movsd  %xmm1,(%rdx)
  movsd  %xmm0,(%rdx,%rsi,8)
  add    $0x8,%rdx
  cmp    %r11,%rdx
  jne    L259
  movsd  (%r12),%xmm1
L182:
  movsd  (%rdi),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rdi)
  add    $0x8,%rdi
  cmp    %r11,%rdi
  jne    L182
L269:
  test   %r10d,%r10d
  jle    L158
  lea    0x0(%rbp,%r14,8),%r11
  movslq %r8d,%rdi
  mov    %r15d,%edx
  add    %r14,%rdx
  mov    -0x48(%rsp),%rsi
  lea    (%rsi,%rdx,8),%rsi
  mov    $0x0,%r12d
L296:
  movsd  (%r11),%xmm1
  xorpd  %xmm4,%xmm1
  cmp    %r9d,%eax
  jle    L160
  mov    %r11,%rdx
L221:
  movapd %xmm1,%xmm0
  mulsd  (%rdx,%rdi,8),%xmm0
  addsd  (%rdx),%xmm0
  movsd  %xmm0,(%rdx)
  add    $0x8,%rdx
  cmp    %rsi,%rdx
  jne    L221
L160:
  add    $0x1,%r12d
  add    %rcx,%r11
  sub    %r13,%rdi
  add    %rcx,%rsi
  cmp    %r10d,%r12d
  jne    L296
L158:
  cmp    %ebx,%eax
  jle    L83
  lea    (%r8,%rax,1),%esi
  movslq %esi,%rsi
  lea    (%rsi,%r14,1),%rdx
  lea    0x0(%rbp,%rdx,8),%r8
  mov    -0x20(%rsp),%rdi
  neg    %rdi
  mov    %r15d,%edx
  add    %r14,%rdx
  add    %rsi,%rdx
  mov    -0x48(%rsp),%rsi
  lea    (%rsi,%rdx,8),%rsi
  mov    %ebx,%r10d
L105:
  movsd  (%r8),%xmm2
  cmp    %r9d,%eax
  jle    L21
  mov    %r8,%rdx
L312:
  movapd %xmm2,%xmm1
  mulsd  (%rdx,%rdi,8),%xmm1
  movsd  (%rdx),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rdx)
  add    $0x8,%rdx
  cmp    %rsi,%rdx
  jne    L312
L21:
  add    $0x1,%r10d
  add    %rcx,%r8
  sub    %r13,%rdi
  add    %rcx,%rsi
  cmp    %r10d,%eax
  jne    L105
  mov    %ebx,%r10d
  mov    -0x30(%rsp),%rdi
L275:
  add    $0x1,%r9d
  add    $0x1,%r14
  sub    $0x1,%r15d
  cmp    %r9d,%eax
  je     L313
L32:
  mov    %eax,%r8d
  imul   %r10d,%r8d
  lea    (%r8,%r9,1),%edx
  movslq %edx,%rdx
  lea    0x0(%rbp,%rdx,8),%r12
  movsd  (%r12),%xmm1
  andpd  %xmm6,%xmm1
  lea    0x1(%r10),%ebx
  cmp    %ebx,%eax
  jle    L60
  lea    (%r8,%rax,1),%edx
  movslq %edx,%rdx
  add    %r14,%rdx
  lea    0x0(%rbp,%rdx,8),%r11
  mov    %r10d,%esi
  mov    %ebx,%edx
L121:
  movsd  (%r11),%xmm0
  andpd  %xmm3,%xmm0
  comisd %xmm1,%xmm0
  maxsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  cmova  %edx,%esi
  add    $0x1,%edx
  add    %rcx,%r11
  cmp    %edx,%eax
  jne    L121
L208:
  comisd %xmm1,%xmm5
  jbe    L215
  mov    -0x38(%rsp),%rbx
  lea    0x4(%rbx),%rdx
  mov    %r9d,(%rbx)
  addl   $0x1,-0x24(%rsp)
  cmp    %r10d,%eax
  jle    L190
  movslq %r8d,%rsi
  add    %r14,%rsi
  lea    0x0(%rbp,%rsi,8),%r8
  mov    %r10d,%esi
L12:
  movq   $0x0,(%r8)
  add    $0x1,%esi
  add    %rcx,%r8
  cmp    %esi,%eax
  jne    L12
  mov    %rdx,-0x38(%rsp)
  jmp    L275
L190:
  mov    %rdx,-0x38(%rsp)
  jmp    L275
L83:
  mov    %ebx,%r10d
  mov    -0x30(%rsp),%rdi
  jmp    L275
L313:
  mov    -0x8(%rsp),%r8
  mov    -0x24(%rsp),%ebx
  test   %ebx,%ebx
  jle    L219
  mov    %r9d,%r14d
  mov    -0x10(%rsp),%r15
  lea    0x4(%r15),%rdx
  lea    -0x1(%rbx),%eax
  lea    (%rdx,%rax,4),%rdi
  mov    %r9d,%r12d
  sub    %ebx,%r12d
  movslq %r9d,%r9
  lea    0x0(,%r9,8),%r11
  lea    -0x1(%r12),%eax
  add    %r9,%rax
  lea    (%rdx,%rax,4),%r10
  mov    %r15,%rbx
  mov    $0x0,%r13d
  movq   D12(%rip),%xmm2
  movsd  D47(%rip),%xmm1
  pxor   %xmm3,%xmm3
  mov    %r14d,-0x48(%rsp)
  mov    -0x18(%rsp),%r14
  jmp    L229
L290:
  add    %ecx,%eax
  cltq
  movsd  %xmm0,(%r8,%rax,8)
  add    $0x4,%rdx
  cmp    %rdi,%rdx
  je     L166
L25:
  mov    (%rdx),%eax
  movapd %xmm1,%xmm0
  cmp    %esi,%eax
  je     L290
  movapd %xmm3,%xmm0
  jmp    L290
L166:
  add    $0x4,%rbx
  add    -0x48(%rsp),%r13d
  cmp    %rbx,%rdi
  je     L219
L229:
  mov    (%rbx),%esi
  test   %r12d,%r12d
  jle    L278
  mov    %r13d,%r9d
  movslq %esi,%rax
  lea    0x0(%rbp,%rax,8),%rcx
  mov    %r14,%rax
L85:
  mov    %r9d,%edx
  add    (%rax),%edx
  movslq %edx,%rdx
  movsd  (%rcx),%xmm0
  xorpd  %xmm2,%xmm0
  movsd  %xmm0,(%r8,%rdx,8)
  add    %r11,%rcx
  add    $0x4,%rax
  cmp    %r10,%rax
  jne    L85
L278:
  mov    %r13d,%ecx
  mov    %r15,%rdx
  jmp    L25
L113:
  movl   $0x0,-0x24(%rsp)
L219:
  mov    -0x24(%rsp),%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: read_data
arguments: data_label: D25
type_guessed: qword
name: read_data
arguments: data_label: D29
type_guessed: byte[8]
name: read_data
arguments: data_label: D12
type_guessed: byte[8]
name: read_data
arguments: data_label: D47
type_guessed: byte[8]

tool: D25: 9223372036854775807

tool: {"D29": [0x11,0xea,0x2d,0x81,0x99,0x97,0x71,0x3d]}

tool: {"D12": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80]}

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

assistant: ```c
int nullspace(double *matrix, double *nullspace, int n, int *buffer) {
  int *pnopivot = buffer, nonpivot = 0;
  int *ppivot = buffer + n;
  int i, j, k, l, ptr, max;
  double pivot, t;
  double tol = 1e-12;
  ptr = 0;
  i = 0;
  for (j = 0; j < n; j++) {
    pivot = fabs(matrix[n * i + j]);
    max = i;
    for (k = i + 1; k < n; k++) {
      t = fabs(matrix[n * k + j]);
      if (pivot < t) {
        pivot = t;
        max = k;
      }
    }
    if (pivot < tol) {
      *(pnopivot++) = j;
      nonpivot++;
      for (k = i; k < n; k++)
        matrix[n * k + j] = 0;
    } else {
      *(ppivot++) = j;
      for (k = j; k < n; k++) {
        t = matrix[i * n + k];
        matrix[i * n + k] = matrix[max * n + k];
        matrix[max * n + k] = t;
      }
      pivot = matrix[i * n + j];
      for (k = j; k < n; k++)
        matrix[i * n + k] /= pivot;
      for (k = 0; k < i; k++) {
        pivot = -matrix[k * n + j];
        for (l = j; l < n; l++)
          matrix[k * n + l] += pivot * matrix[i * n + l];
      }
      for (k = i + 1; k < n; k++) {
        pivot = matrix[k * n + j];
        for (l = j; l < n; l++)
          matrix[k * n + l] -= pivot * matrix[i * n + l];
      }
      i++;
    }
  }
  for (k = 0; k < nonpivot; k++) {
    j = buffer[k];
    for (l = 0; l < n - nonpivot; l++)
      nullspace[k * n + buffer[n + l]] = -matrix[l * n + j];
    for (l = 0; l < nonpivot; l++)
      nullspace[k * n + buffer[l]] = (j == buffer[l]) ? 1 : 0;
  }
  return nonpivot;
}
```