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

```asm
<tpd9>:
  endbr64
  push   %rbx
  mov    %rdx,%rbx
  mov    %ecx,%edx
  lea    0x3(%rdi),%ecx
  movslq %ecx,%rcx
  sub    $0x20,%rsp
  cmpl   $0x3c,(%rsi,%rcx,4)
  jne    L93
  cmp    $0x2,%edx
  jle    L7
L93:
  add    $0x20,%rsp
  mov    $0x1,%eax
  pop    %rbx
  ret
  nopw   0x0(%rax,%rax,1)
L7:
  mov    0x14(%rsi),%ecx
  movsd  (%r8),%xmm2
  mov    %rsi,%rax
  movsd  0x8(%r8),%xmm1
  test   %ecx,%ecx
  je     L34
  movsd  0x10(%rbx),%xmm3
  movsd  0x8(%rbx),%xmm0
  add    $0x30,%rbx
  mulsd  %xmm1,%xmm3
  mulsd  -0x8(%rbx),%xmm1
  mulsd  %xmm2,%xmm0
  mulsd  -0x10(%rbx),%xmm2
  addsd  -0x30(%rbx),%xmm0
  addsd  -0x18(%rbx),%xmm2
  addsd  %xmm2,%xmm1
  movapd %xmm0,%xmm2
  addsd  %xmm3,%xmm2
L34:
  test   %edi,%edi
  jne    L29
L6:
  movsd  0x188(%rbx),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  0x140(%rbx),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  0xf8(%rbx),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  0xc0(%rbx),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  0x88(%rbx),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  0x60(%rbx),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  0x38(%rbx),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  0x20(%rbx),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  0x8(%rbx),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  (%rbx),%xmm0
  movsd  %xmm0,(%r9)
  cmp    $0x1,%edx
  je     L84
  movsd  0x190(%rbx),%xmm4
  movsd  0x198(%rbx),%xmm3
  mulsd  %xmm2,%xmm4
  mulsd  %xmm1,%xmm3
  addsd  0x148(%rbx),%xmm3
  addsd  %xmm4,%xmm3
  movsd  0x1a0(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  0x150(%rbx),%xmm4
  mulsd  %xmm2,%xmm3
  mulsd  %xmm1,%xmm4
  addsd  0x100(%rbx),%xmm4
  addsd  %xmm4,%xmm3
  movsd  0x1a8(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  0x158(%rbx),%xmm4
  mulsd  %xmm2,%xmm3
  mulsd  %xmm1,%xmm4
  addsd  0x108(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  0xc8(%rbx),%xmm4
  addsd  %xmm4,%xmm3
  movsd  0x1b0(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  0x160(%rbx),%xmm4
  mulsd  %xmm2,%xmm3
  mulsd  %xmm1,%xmm4
  addsd  0x110(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  0xd0(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  0x90(%rbx),%xmm4
  addsd  %xmm4,%xmm3
  movsd  0x1b8(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  0x168(%rbx),%xmm4
  mulsd  %xmm2,%xmm3
  mulsd  %xmm1,%xmm4
  addsd  0x118(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  0xd8(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  0x98(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  0x68(%rbx),%xmm4
  addsd  %xmm4,%xmm3
  movsd  0x1c0(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  0x170(%rbx),%xmm4
  mulsd  %xmm2,%xmm3
  mulsd  %xmm1,%xmm4
  addsd  0x120(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  0xe0(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  0xa0(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  0x70(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  0x40(%rbx),%xmm4
  addsd  %xmm4,%xmm3
  movsd  0x1c8(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  0x178(%rbx),%xmm4
  mulsd  %xmm2,%xmm3
  mulsd  %xmm1,%xmm4
  addsd  0x128(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  0xe8(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  0xa8(%rbx),%xmm4
  mov    0x18(%rax),%eax
  mulsd  %xmm1,%xmm4
  addsd  0x78(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  0x48(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  0x28(%rbx),%xmm4
  addsd  %xmm4,%xmm3
  movsd  0x1d0(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  0x180(%rbx),%xmm4
  mulsd  %xmm2,%xmm3
  mulsd  %xmm1,%xmm4
  addsd  0x130(%rbx),%xmm4
  mulsd  %xmm1,%xmm3
  mulsd  %xmm1,%xmm4
  addsd  0xf0(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  0xb0(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  0x80(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  0x50(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  0x30(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  0x10(%rbx),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  %xmm4,%xmm3
  addsd  %xmm0,%xmm3
  movsd  %xmm3,(%r9)
  test   %eax,%eax
  jne    L64
L84:
  xor    %eax,%eax
L96:
  add    $0x20,%rsp
  pop    %rbx
  ret
  nopw   0x0(%rax,%rax,1)
L29:
  movslq 0xc(%rax),%rcx
  lea    (%rbx,%rcx,8),%rbx
  jmp    L6
  nopl   (%rax)
L64:
  mulsd  %xmm2,%xmm2
  pxor   %xmm0,%xmm0
  mulsd  %xmm1,%xmm1
  addsd  %xmm1,%xmm2
  ucomisd %xmm2,%xmm0
  ja     L32
  movapd %xmm2,%xmm1
  sqrtsd %xmm1,%xmm1
L20:
  movsd  0x1d8(%rbx),%xmm0
  xor    %eax,%eax
  mulsd  %xmm2,%xmm0
  addsd  0x138(%rbx),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  0xb8(%rbx),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  0x58(%rbx),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  0x18(%rbx),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  %xmm3,%xmm0
  movsd  %xmm0,(%r9)
  jmp    L96
L32:
  movapd %xmm2,%xmm0
  mov    %r9,0x18(%rsp)
  movsd  %xmm3,0x10(%rsp)
  movsd  %xmm2,0x8(%rsp)
  call   <sqrt@plt>
  mov    0x18(%rsp),%r9
  movsd  0x10(%rsp),%xmm3
  movsd  0x8(%rsp),%xmm2
  movapd %xmm0,%xmm1
  jmp    L20
```

assistant: ```c
int tpd9(int inverse, const int i[], const double p[], int Nhat, const double rawcrd[], double *discrd) {
  double r, s, u, v;
  if (i[3 + inverse] != 60 || 2 < Nhat) {
    return 1;
  }
  u = rawcrd[0];
  v = rawcrd[1];
  if (i[5]) {
    r = p[0] + p[1] * u + p[2] * v;
    v = p[3] + p[4] * u + p[5] * v;
    u = r;
    p += 6;
  }
  if (inverse)
    p += i[3];
  *discrd = p[0] + u * (p[1] + u * (p[4] + u * (p[7] + u * (p[12] + u * (p[17] + u * (p[24] + u * (p[31] + u * (p[40] + u * (p[49])))))))));
  if (Nhat == 1)
    return 0;
  *discrd += v * (p[2] + v * (p[6] + v * (p[10] + v * (p[16] + v * (p[22] + v * (p[30] + v * (p[38] + v * (p[48] + v * (p[58]))))))))) + u * (p[5] + v * (p[9] + v * (p[15] + v * (p[21] + v * (p[29] + v * (p[37] + v * (p[47] + v * (p[57]))))))) + u * (p[8] + v * (p[14] + v * (p[20] + v * (p[28] + v * (p[36] + v * (p[46] + v * (p[56])))))) + u * (p[13] + v * (p[19] + v * (p[27] + v * (p[35] + v * (p[45] + v * (p[55]))))) + u * (p[18] + v * (p[26] + v * (p[34] + v * (p[44] + v * (p[54])))) + u * (p[25] + v * (p[33] + v * (p[43] + v * (p[53]))) + u * (p[32] + v * (p[42] + v * (p[52])) + u * (p[41] + v * (p[51]) + u * (p[50])))))))) * v;
  if (i[6]) {
    s = u * u + v * v;
    r = sqrt(s);
    *discrd += r * (p[3] + s * (p[11] + s * (p[23] + s * (p[39] + s * (p[59])))));
  }
  return 0;
}
```