user: 请从给定的汇编指令推导出C源代码：

```asm
<tpd9>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x50,%rsp
  mov    %edi,-0x24(%rbp)
  mov    %rsi,-0x30(%rbp)
  mov    %rdx,-0x38(%rbp)
  mov    %ecx,-0x28(%rbp)
  mov    %r8,-0x40(%rbp)
  mov    %r9,-0x48(%rbp)
  mov    -0x24(%rbp),%eax
  add    $0x4,%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  cmp    $0x3c,%eax
  jne    L40
  cmpl   $0x2,-0x28(%rbp)
  jle    L7
L40:
  mov    $0x1,%eax
  jmp    L35
L7:
  mov    -0x40(%rbp),%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  mov    -0x40(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x30(%rbp),%rax
  add    $0x18,%rax
  mov    (%rax),%eax
  test   %eax,%eax
  je     L26
  mov    -0x38(%rbp),%rax
  movsd  (%rax),%xmm1
  mov    -0x38(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x20(%rbp),%xmm0
  addsd  %xmm0,%xmm1
  mov    -0x38(%rbp),%rax
  add    $0x10,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  mov    -0x38(%rbp),%rax
  add    $0x18,%rax
  movsd  (%rax),%xmm1
  mov    -0x38(%rbp),%rax
  add    $0x20,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x20(%rbp),%xmm0
  addsd  %xmm0,%xmm1
  mov    -0x38(%rbp),%rax
  add    $0x28,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x10(%rbp),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  addq   $0x30,-0x38(%rbp)
L26:
  cmpl   $0x0,-0x24(%rbp)
  je     L24
  mov    -0x30(%rbp),%rax
  add    $0x10,%rax
  mov    (%rax),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0x38(%rbp)
L24:
  mov    -0x38(%rbp),%rax
  movsd  (%rax),%xmm1
  mov    -0x38(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm2
  mov    -0x38(%rbp),%rax
  add    $0x20,%rax
  movsd  (%rax),%xmm3
  mov    -0x38(%rbp),%rax
  add    $0x38,%rax
  movsd  (%rax),%xmm4
  mov    -0x38(%rbp),%rax
  add    $0x60,%rax
  movsd  (%rax),%xmm5
  mov    -0x38(%rbp),%rax
  add    $0x88,%rax
  movsd  (%rax),%xmm6
  mov    -0x38(%rbp),%rax
  add    $0xc0,%rax
  movsd  (%rax),%xmm7
  mov    -0x38(%rbp),%rax
  add    $0xf8,%rax
  movsd  (%rax),%xmm8
  mov    -0x38(%rbp),%rax
  add    $0x140,%rax
  movsd  (%rax),%xmm9
  mov    -0x38(%rbp),%rax
  add    $0x188,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x20(%rbp),%xmm0
  addsd  %xmm9,%xmm0
  mulsd  -0x20(%rbp),%xmm0
  addsd  %xmm8,%xmm0
  mulsd  -0x20(%rbp),%xmm0
  addsd  %xmm7,%xmm0
  mulsd  -0x20(%rbp),%xmm0
  addsd  %xmm6,%xmm0
  mulsd  -0x20(%rbp),%xmm0
  addsd  %xmm5,%xmm0
  mulsd  -0x20(%rbp),%xmm0
  addsd  %xmm4,%xmm0
  mulsd  -0x20(%rbp),%xmm0
  addsd  %xmm3,%xmm0
  mulsd  -0x20(%rbp),%xmm0
  addsd  %xmm2,%xmm0
  mulsd  -0x20(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  mov    -0x48(%rbp),%rax
  movsd  %xmm0,(%rax)
  cmpl   $0x1,-0x28(%rbp)
  jne    L6
  mov    $0x0,%eax
  jmp    L35
L6:
  mov    -0x48(%rbp),%rax
  movsd  (%rax),%xmm1
  mov    -0x38(%rbp),%rax
  add    $0x10,%rax
  movsd  (%rax),%xmm2
  mov    -0x38(%rbp),%rax
  add    $0x30,%rax
  movsd  (%rax),%xmm3
  mov    -0x38(%rbp),%rax
  add    $0x50,%rax
  movsd  (%rax),%xmm4
  mov    -0x38(%rbp),%rax
  sub    $0xffffffffffffff80,%rax
  movsd  (%rax),%xmm5
  mov    -0x38(%rbp),%rax
  add    $0xb0,%rax
  movsd  (%rax),%xmm6
  mov    -0x38(%rbp),%rax
  add    $0xf0,%rax
  movsd  (%rax),%xmm7
  mov    -0x38(%rbp),%rax
  add    $0x130,%rax
  movsd  (%rax),%xmm8
  mov    -0x38(%rbp),%rax
  add    $0x180,%rax
  movsd  (%rax),%xmm9
  mov    -0x38(%rbp),%rax
  add    $0x1d0,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm9,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm8,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm7,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm6,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm5,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm4,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm3,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm2,%xmm0
  movapd %xmm0,%xmm2
  mulsd  -0x18(%rbp),%xmm2
  mov    -0x38(%rbp),%rax
  add    $0x28,%rax
  movsd  (%rax),%xmm3
  mov    -0x38(%rbp),%rax
  add    $0x48,%rax
  movsd  (%rax),%xmm4
  mov    -0x38(%rbp),%rax
  add    $0x78,%rax
  movsd  (%rax),%xmm5
  mov    -0x38(%rbp),%rax
  add    $0xa8,%rax
  movsd  (%rax),%xmm6
  mov    -0x38(%rbp),%rax
  add    $0xe8,%rax
  movsd  (%rax),%xmm7
  mov    -0x38(%rbp),%rax
  add    $0x128,%rax
  movsd  (%rax),%xmm8
  mov    -0x38(%rbp),%rax
  add    $0x178,%rax
  movsd  (%rax),%xmm9
  mov    -0x38(%rbp),%rax
  add    $0x1c8,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm9,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm8,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm7,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm6,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm5,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm4,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm0,%xmm3
  mov    -0x38(%rbp),%rax
  add    $0x40,%rax
  movsd  (%rax),%xmm4
  mov    -0x38(%rbp),%rax
  add    $0x70,%rax
  movsd  (%rax),%xmm5
  mov    -0x38(%rbp),%rax
  add    $0xa0,%rax
  movsd  (%rax),%xmm6
  mov    -0x38(%rbp),%rax
  add    $0xe0,%rax
  movsd  (%rax),%xmm7
  mov    -0x38(%rbp),%rax
  add    $0x120,%rax
  movsd  (%rax),%xmm8
  mov    -0x38(%rbp),%rax
  add    $0x170,%rax
  movsd  (%rax),%xmm9
  mov    -0x38(%rbp),%rax
  add    $0x1c0,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm9,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm8,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm7,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm6,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm5,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm0,%xmm4
  mov    -0x38(%rbp),%rax
  add    $0x68,%rax
  movsd  (%rax),%xmm5
  mov    -0x38(%rbp),%rax
  add    $0x98,%rax
  movsd  (%rax),%xmm6
  mov    -0x38(%rbp),%rax
  add    $0xd8,%rax
  movsd  (%rax),%xmm7
  mov    -0x38(%rbp),%rax
  add    $0x118,%rax
  movsd  (%rax),%xmm8
  mov    -0x38(%rbp),%rax
  add    $0x168,%rax
  movsd  (%rax),%xmm9
  mov    -0x38(%rbp),%rax
  add    $0x1b8,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm9,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm8,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm7,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm6,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm0,%xmm5
  mov    -0x38(%rbp),%rax
  add    $0x90,%rax
  movsd  (%rax),%xmm6
  mov    -0x38(%rbp),%rax
  add    $0xd0,%rax
  movsd  (%rax),%xmm7
  mov    -0x38(%rbp),%rax
  add    $0x110,%rax
  movsd  (%rax),%xmm8
  mov    -0x38(%rbp),%rax
  add    $0x160,%rax
  movsd  (%rax),%xmm9
  mov    -0x38(%rbp),%rax
  add    $0x1b0,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm9,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm8,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm7,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm0,%xmm6
  mov    -0x38(%rbp),%rax
  add    $0xc8,%rax
  movsd  (%rax),%xmm7
  mov    -0x38(%rbp),%rax
  add    $0x108,%rax
  movsd  (%rax),%xmm8
  mov    -0x38(%rbp),%rax
  add    $0x158,%rax
  movsd  (%rax),%xmm9
  mov    -0x38(%rbp),%rax
  add    $0x1a8,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm9,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm8,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm0,%xmm7
  mov    -0x38(%rbp),%rax
  add    $0x100,%rax
  movsd  (%rax),%xmm8
  mov    -0x38(%rbp),%rax
  add    $0x150,%rax
  movsd  (%rax),%xmm9
  mov    -0x38(%rbp),%rax
  add    $0x1a0,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm9,%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm0,%xmm8
  mov    -0x38(%rbp),%rax
  add    $0x148,%rax
  movsd  (%rax),%xmm9
  mov    -0x38(%rbp),%rax
  add    $0x198,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm0,%xmm9
  mov    -0x38(%rbp),%rax
  add    $0x190,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x20(%rbp),%xmm0
  addsd  %xmm9,%xmm0
  mulsd  -0x20(%rbp),%xmm0
  addsd  %xmm8,%xmm0
  mulsd  -0x20(%rbp),%xmm0
  addsd  %xmm7,%xmm0
  mulsd  -0x20(%rbp),%xmm0
  addsd  %xmm6,%xmm0
  mulsd  -0x20(%rbp),%xmm0
  addsd  %xmm5,%xmm0
  mulsd  -0x20(%rbp),%xmm0
  addsd  %xmm4,%xmm0
  mulsd  -0x20(%rbp),%xmm0
  addsd  %xmm3,%xmm0
  mulsd  -0x20(%rbp),%xmm0
  mulsd  -0x18(%rbp),%xmm0
  addsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  mov    -0x48(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x30(%rbp),%rax
  add    $0x1c,%rax
  mov    (%rax),%eax
  test   %eax,%eax
  je     L54
  movsd  -0x20(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  %xmm0,%xmm1
  movsd  -0x18(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  mov    -0x8(%rbp),%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x10(%rbp)
  mov    -0x48(%rbp),%rax
  movsd  (%rax),%xmm1
  mov    -0x38(%rbp),%rax
  add    $0x18,%rax
  movsd  (%rax),%xmm2
  mov    -0x38(%rbp),%rax
  add    $0x58,%rax
  movsd  (%rax),%xmm3
  mov    -0x38(%rbp),%rax
  add    $0xb8,%rax
  movsd  (%rax),%xmm4
  mov    -0x38(%rbp),%rax
  add    $0x138,%rax
  movsd  (%rax),%xmm5
  mov    -0x38(%rbp),%rax
  add    $0x1d8,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x8(%rbp),%xmm0
  addsd  %xmm5,%xmm0
  mulsd  -0x8(%rbp),%xmm0
  addsd  %xmm4,%xmm0
  mulsd  -0x8(%rbp),%xmm0
  addsd  %xmm3,%xmm0
  mulsd  -0x8(%rbp),%xmm0
  addsd  %xmm2,%xmm0
  mulsd  -0x10(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  mov    -0x48(%rbp),%rax
  movsd  %xmm0,(%rax)
L54:
  mov    $0x0,%eax
L35:
  leave
  ret
```

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[4 + inverse] != 60 || 2 < Nhat) {
    return 1;
  }
  u = rawcrd[0];
  v = rawcrd[1];
  if (i[6]) {
    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[4];
  *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[7]) {
    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;
}
```