user: Translate the below assembly code into C programming language:

```asm
<tpd5>:
  endbr64
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x10,%rsp
  mov    %rdx,%rbx
  mov    %ecx,%edx
  lea    0x4(%rdi),%ecx
  movslq %ecx,%rcx
  cmpl   $0x18,(%rsi,%rcx,4)
  jne    L68
  mov    %rsi,%rax
  mov    %r9,%rbp
  cmp    $0x2,%edx
  jle    L9
L68:
  mov    $0x1,%ecx
L44:
  mov    %ecx,%eax
  add    $0x10,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  ret
L9:
  movsd  (%r8),%xmm3
  movsd  0x8(%r8),%xmm2
  cmpl   $0x0,0x18(%rsi)
  je     L60
  movapd %xmm3,%xmm0
  mulsd  0x8(%rbx),%xmm0
  addsd  (%rbx),%xmm0
  movapd %xmm2,%xmm1
  mulsd  0x10(%rbx),%xmm1
  mulsd  0x28(%rbx),%xmm2
  mulsd  0x20(%rbx),%xmm3
  addsd  0x18(%rbx),%xmm3
  addsd  %xmm3,%xmm2
  add    $0x30,%rbx
  movapd %xmm0,%xmm3
  addsd  %xmm1,%xmm3
L60:
  test   %edi,%edi
  je     L30
  movslq 0x10(%rax),%rcx
  lea    (%rbx,%rcx,8),%rbx
L30:
  movapd %xmm3,%xmm0
  mulsd  0x88(%rbx),%xmm0
  addsd  0x60(%rbx),%xmm0
  mulsd  %xmm3,%xmm0
  addsd  0x38(%rbx),%xmm0
  mulsd  %xmm3,%xmm0
  addsd  0x20(%rbx),%xmm0
  mulsd  %xmm3,%xmm0
  addsd  0x8(%rbx),%xmm0
  mulsd  %xmm3,%xmm0
  addsd  (%rbx),%xmm0
  movsd  %xmm0,0x0(%rbp)
  mov    $0x0,%ecx
  cmp    $0x1,%edx
  je     L44
  movapd %xmm2,%xmm1
  mulsd  0x98(%rbx),%xmm1
  addsd  0x68(%rbx),%xmm1
  movapd %xmm3,%xmm4
  mulsd  0x90(%rbx),%xmm4
  addsd  %xmm4,%xmm1
  mulsd  %xmm3,%xmm1
  movapd %xmm2,%xmm4
  mulsd  0xa0(%rbx),%xmm4
  addsd  0x70(%rbx),%xmm4
  mulsd  %xmm2,%xmm4
  addsd  0x40(%rbx),%xmm4
  addsd  %xmm4,%xmm1
  mulsd  %xmm3,%xmm1
  movapd %xmm2,%xmm4
  mulsd  0xa8(%rbx),%xmm4
  addsd  0x78(%rbx),%xmm4
  mulsd  %xmm2,%xmm4
  addsd  0x48(%rbx),%xmm4
  mulsd  %xmm2,%xmm4
  addsd  0x28(%rbx),%xmm4
  addsd  %xmm4,%xmm1
  mulsd  %xmm3,%xmm1
  mulsd  %xmm2,%xmm1
  movapd %xmm2,%xmm4
  mulsd  0xb0(%rbx),%xmm4
  addsd  0x80(%rbx),%xmm4
  mulsd  %xmm2,%xmm4
  addsd  0x50(%rbx),%xmm4
  mulsd  %xmm2,%xmm4
  addsd  0x30(%rbx),%xmm4
  mulsd  %xmm2,%xmm4
  addsd  0x10(%rbx),%xmm4
  mulsd  %xmm2,%xmm4
  addsd  %xmm4,%xmm1
  addsd  %xmm0,%xmm1
  movq   %xmm1,%r12
  movsd  %xmm1,0x0(%rbp)
  mov    0x1c(%rax),%ecx
  test   %ecx,%ecx
  je     L44
  mulsd  %xmm3,%xmm3
  mulsd  %xmm2,%xmm2
  addsd  %xmm2,%xmm3
  movapd %xmm3,%xmm5
  movsd  %xmm3,0x8(%rsp)
  pxor   %xmm0,%xmm0
  ucomisd %xmm3,%xmm0
  ja     L31
  sqrtsd %xmm5,%xmm5
  movapd %xmm5,%xmm1
L75:
  movsd  0x8(%rsp),%xmm6
  movapd %xmm6,%xmm0
  mulsd  0xb8(%rbx),%xmm0
  addsd  0x58(%rbx),%xmm0
  mulsd  %xmm6,%xmm0
  addsd  0x18(%rbx),%xmm0
  mulsd  %xmm1,%xmm0
  movq   %r12,%xmm7
  addsd  %xmm7,%xmm0
  movsd  %xmm0,0x0(%rbp)
  mov    $0x0,%ecx
  jmp    L44
L31:
  movsd  0x8(%rsp),%xmm0
  call   <sqrt@plt>
  movapd %xmm0,%xmm1
  jmp    L75
```

assistant: ```c
int tpd5(int inverse, const int i[], const double p[], int Nhat, const double rawcrd[], double *discrd) {
  double r, s, u, v;
  if (i[4 + inverse] != 24 || 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])))));
  if (Nhat == 1)
    return 0;
  *discrd += v * (p[2] + v * (p[6] + v * (p[10] + v * (p[16] + v * (p[22]))))) + u * (p[5] + v * (p[9] + v * (p[15] + v * (p[21]))) + u * (p[8] + v * (p[14] + v * (p[20])) + u * (p[13] + v * (p[19]) + u * (p[18])))) * v;
  if (i[7]) {
    s = u * u + v * v;
    r = sqrt(s);
    *discrd += r * (p[3] + s * (p[11] + s * (p[23])));
  }
  return 0;
}
```