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

```asm
<DOsvdPlane>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x50,%rsp
  mov    %rdi,-0x38(%rbp)
  mov    %esi,-0x3c(%rbp)
  mov    %rdx,-0x48(%rbp)
  mov    %rcx,-0x50(%rbp)
  mov    -0x3c(%rbp),%eax
  mov    $0x3,%esi
  mov    %eax,%edi
  call   <allocMatrix@plt>
  mov    %rax,-0x18(%rbp)
  mov    $0x3,%esi
  mov    $0x3,%edi
  call   <allocMatrix@plt>
  mov    %rax,-0x10(%rbp)
  mov    $0x3,%edi
  call   <allocVect@plt>
  mov    %rax,-0x8(%rbp)
  movss  D5(%rip),%xmm0
  movss  %xmm0,-0x20(%rbp)
  cmpl   $0x2,-0x3c(%rbp)
  setg   %al
  movzbl %al,%eax
  mov    %eax,%edi
  call   <assert@plt>
  mov    -0x38(%rbp),%rax
  movss  (%rax),%xmm0
  mov    -0x50(%rbp),%rax
  movss  %xmm0,(%rax)
  mov    -0x50(%rbp),%rax
  lea    0x4(%rax),%rdx
  mov    -0x38(%rbp),%rax
  movss  0x4(%rax),%xmm0
  movss  %xmm0,(%rdx)
  mov    -0x50(%rbp),%rax
  lea    0x8(%rax),%rdx
  mov    -0x38(%rbp),%rax
  movss  0x8(%rax),%xmm0
  movss  %xmm0,(%rdx)
  movl   $0x1,-0x28(%rbp)
  jmp    L94
L37:
  movl   $0x0,-0x24(%rbp)
  jmp    L43
L63:
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm1
  mov    -0x28(%rbp),%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %eax,%edx
  mov    -0x24(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  addss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  addl   $0x1,-0x24(%rbp)
L43:
  cmpl   $0x2,-0x24(%rbp)
  jle    L63
  addl   $0x1,-0x28(%rbp)
L94:
  mov    -0x28(%rbp),%eax
  cmp    -0x3c(%rbp),%eax
  jl     L37
  mov    -0x50(%rbp),%rax
  movss  (%rax),%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2ssl -0x3c(%rbp),%xmm1
  divss  %xmm1,%xmm0
  mov    -0x50(%rbp),%rax
  movss  %xmm0,(%rax)
  mov    -0x50(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2ssl -0x3c(%rbp),%xmm1
  mov    -0x50(%rbp),%rax
  add    $0x4,%rax
  divss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x50(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2ssl -0x3c(%rbp),%xmm1
  mov    -0x50(%rbp),%rax
  add    $0x8,%rax
  divss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  movl   $0x0,-0x28(%rbp)
  jmp    L106
L102:
  movl   $0x0,-0x24(%rbp)
  jmp    L84
L1:
  mov    -0x24(%rbp),%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %eax,%edx
  mov    -0x28(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mov    -0x28(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm1
  subss  %xmm1,%xmm0
  mov    -0x24(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x28(%rbp),%eax
  cltq
  add    $0x1,%rax
  shl    $0x3,%rax
  add    %rdx,%rax
  cvtss2sd %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x24(%rbp)
L84:
  mov    -0x24(%rbp),%eax
  cmp    -0x3c(%rbp),%eax
  jl     L1
  addl   $0x1,-0x28(%rbp)
L106:
  cmpl   $0x2,-0x28(%rbp)
  jle    L102
  mov    -0x10(%rbp),%rcx
  mov    -0x8(%rbp),%rdx
  mov    -0x3c(%rbp),%esi
  mov    -0x18(%rbp),%rax
  mov    %rcx,%r8
  mov    %rdx,%rcx
  mov    $0x3,%edx
  mov    %rax,%rdi
  call   <svdcmp@plt>
  movl   $0x0,-0x28(%rbp)
  jmp    L75
L25:
  mov    -0x28(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x8(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  movq   %rax,%xmm0
  call   <fabs@plt>
  movss  -0x20(%rbp),%xmm1
  comiss %xmm0,%xmm1
  jbe    L16
  mov    -0x28(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x8(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  movq   %rax,%xmm0
  call   <fabs@plt>
  movd   %xmm0,%eax
  mov    %eax,-0x20(%rbp)
  mov    -0x28(%rbp),%eax
  mov    %eax,-0x1c(%rbp)
L16:
  addl   $0x1,-0x28(%rbp)
L75:
  cmpl   $0x2,-0x28(%rbp)
  jle    L25
  mov    -0x10(%rbp),%rax
  add    $0x8,%rax
  mov    (%rax),%rdx
  mov    -0x1c(%rbp),%eax
  cltq
  add    $0x1,%rax
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  cvtsd2ss %xmm0,%xmm0
  mov    -0x48(%rbp),%rax
  movss  %xmm0,(%rax)
  mov    -0x10(%rbp),%rax
  add    $0x10,%rax
  mov    (%rax),%rdx
  mov    -0x1c(%rbp),%eax
  cltq
  add    $0x1,%rax
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x48(%rbp),%rax
  add    $0x4,%rax
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x10(%rbp),%rax
  add    $0x18,%rax
  mov    (%rax),%rdx
  mov    -0x1c(%rbp),%eax
  cltq
  add    $0x1,%rax
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x48(%rbp),%rax
  add    $0x8,%rax
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x8(%rbp),%rax
  mov    %rax,%rdi
  call   <freeVect@plt>
  mov    -0x10(%rbp),%rax
  mov    $0x3,%esi
  mov    %rax,%rdi
  call   <freeMatrix@plt>
  mov    -0x3c(%rbp),%edx
  mov    -0x18(%rbp),%rax
  mov    %edx,%esi
  mov    %rax,%rdi
  call   <freeMatrix@plt>
  nop
  leave
  ret
```

tool_calls: name: parse_data
arguments: tag: D5
type: byte[4]

tool: {"D5": [0xca,0xf2,0x49,0x71]}

assistant: ```c
void DOsvdPlane(float *pnts, int npnts, float *n, float *base) {
  int i, j;
  double **da = allocMatrix(npnts, 3);
  double **v = allocMatrix(3, 3);
  double *w = allocVect(3);
  float mn = 1E30f;
  int bestat;
  assert(npnts >= 3);
  base[0] = pnts[0];
  base[1] = pnts[1];
  base[2] = pnts[2];
  for (i = 1; i < npnts; i++) {
    for (j = 0; j < 3; j++)
      base[j] += pnts[i * 3 + j];
  }
  base[0] /= (float)(npnts);
  base[1] /= (float)(npnts);
  base[2] /= (float)(npnts);
  for (i = 0; i < 3; i++) {
    for (j = 0; j < npnts; j++)
      da[j + 1][i + 1] = pnts[j * 3 + i] - base[i];
  }
  svdcmp(da, npnts, 3, w, v);
  for (i = 0; i < 3; i++) {
    if (fabs(w[i + 1]) < mn) {
      mn = (float)fabs(w[i + 1]);
      bestat = i;
    }
  }
  n[0] = (float)v[1][bestat + 1];
  n[1] = (float)v[2][bestat + 1];
  n[2] = (float)v[3][bestat + 1];
  freeVect(w);
  freeMatrix(v, 3);
  freeMatrix(da, npnts);
}
```