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

```asm
<DQ2QuatTrans>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x30,%rsp
  mov    %rdi,-0x18(%rbp)
  mov    %rsi,-0x20(%rbp)
  mov    %rdx,-0x28(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movl   $0x0,-0x10(%rbp)
  jmp    L11
L29:
  mov    -0x18(%rbp),%rax
  mov    -0x10(%rbp),%edx
  movslq %edx,%rdx
  movsd  (%rax,%rdx,8),%xmm1
  mov    -0x18(%rbp),%rax
  mov    -0x10(%rbp),%edx
  movslq %edx,%rdx
  movsd  (%rax,%rdx,8),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x8(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  addl   $0x1,-0x10(%rbp)
L11:
  cmpl   $0x3,-0x10(%rbp)
  jle    L29
  mov    -0x8(%rbp),%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x8(%rbp)
  movl   $0x0,-0xc(%rbp)
  jmp    L10
L3:
  mov    -0x18(%rbp),%rax
  mov    -0xc(%rbp),%edx
  movslq %edx,%rdx
  movsd  (%rax,%rdx,8),%xmm0
  mov    -0xc(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  divsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0xc(%rbp)
L10:
  cmpl   $0x3,-0xc(%rbp)
  jle    L3
  mov    -0x18(%rbp),%rax
  add    $0x20,%rax
  movsd  (%rax),%xmm0
  movq   D22(%rip),%xmm1
  xorpd  %xmm0,%xmm1
  mov    -0x18(%rbp),%rax
  movsd  0x8(%rax),%xmm0
  mulsd  %xmm0,%xmm1
  mov    -0x18(%rbp),%rax
  add    $0x20,%rax
  movsd  0x8(%rax),%xmm2
  mov    -0x18(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  mov    -0x18(%rbp),%rax
  add    $0x20,%rax
  movsd  0x10(%rax),%xmm2
  mov    -0x18(%rbp),%rax
  movsd  0x18(%rax),%xmm1
  mulsd  %xmm1,%xmm2
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  mov    -0x18(%rbp),%rax
  add    $0x20,%rax
  movsd  0x18(%rax),%xmm2
  mov    -0x18(%rbp),%rax
  movsd  0x10(%rax),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  addsd  %xmm0,%xmm0
  divsd  -0x8(%rbp),%xmm0
  mov    -0x28(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x18(%rbp),%rax
  add    $0x20,%rax
  movsd  (%rax),%xmm0
  movq   D22(%rip),%xmm1
  xorpd  %xmm0,%xmm1
  mov    -0x18(%rbp),%rax
  movsd  0x10(%rax),%xmm0
  mulsd  %xmm0,%xmm1
  mov    -0x18(%rbp),%rax
  add    $0x20,%rax
  movsd  0x8(%rax),%xmm2
  mov    -0x18(%rbp),%rax
  movsd  0x18(%rax),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  mov    -0x18(%rbp),%rax
  add    $0x20,%rax
  movsd  0x10(%rax),%xmm2
  mov    -0x18(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  mov    -0x18(%rbp),%rax
  add    $0x20,%rax
  movsd  0x18(%rax),%xmm2
  mov    -0x18(%rbp),%rax
  movsd  0x8(%rax),%xmm1
  mulsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  addsd  %xmm0,%xmm0
  mov    -0x28(%rbp),%rax
  add    $0x8,%rax
  divsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x18(%rbp),%rax
  add    $0x20,%rax
  movsd  (%rax),%xmm0
  movq   D22(%rip),%xmm1
  xorpd  %xmm0,%xmm1
  mov    -0x18(%rbp),%rax
  movsd  0x18(%rax),%xmm0
  mulsd  %xmm1,%xmm0
  mov    -0x18(%rbp),%rax
  add    $0x20,%rax
  movsd  0x8(%rax),%xmm2
  mov    -0x18(%rbp),%rax
  movsd  0x10(%rax),%xmm1
  mulsd  %xmm1,%xmm2
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  mov    -0x18(%rbp),%rax
  add    $0x20,%rax
  movsd  0x10(%rax),%xmm2
  mov    -0x18(%rbp),%rax
  movsd  0x8(%rax),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  mov    -0x18(%rbp),%rax
  add    $0x20,%rax
  movsd  0x18(%rax),%xmm2
  mov    -0x18(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  addsd  %xmm0,%xmm0
  mov    -0x28(%rbp),%rax
  add    $0x10,%rax
  divsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  nop
  leave
  ret
```

tool_calls: name: unpack_data
arguments: data_label: D22
type: f32[2]

tool: {"D22": [0.0,-0.0]}

assistant: ```c
void DQ2QuatTrans(const double dq[2][4], double q0[4], double t[3]) {
  double len = 0.0;
  for (int i = 0; i < 4; i++)
    len += dq[0][i] * dq[0][i];
  len = sqrt(len);
  for (int i = 0; i < 4; i++)
    q0[i] = dq[0][i] / len;
  t[0] = 2.0 * (-dq[1][0] * dq[0][1] + dq[1][1] * dq[0][0] - dq[1][2] * dq[0][3] + dq[1][3] * dq[0][2]) / len;
  t[1] = 2.0 * (-dq[1][0] * dq[0][2] + dq[1][1] * dq[0][3] + dq[1][2] * dq[0][0] - dq[1][3] * dq[0][1]) / len;
  t[2] = 2.0 * (-dq[1][0] * dq[0][3] - dq[1][1] * dq[0][2] + dq[1][2] * dq[0][1] + dq[1][3] * dq[0][0]) / len;
}
```