user: What is the c source code of the assembly code below:

```asm
<R2Quaternion>:
  endbr64
  movss  (%rdi),%xmm0
  movss  D25(%rip),%xmm1
  push   %r12
  xor    %r12d,%r12d
  movss  0x10(%rdi),%xmm5
  movss  0x20(%rdi),%xmm4
  push   %rbp
  mov    %rsi,%rbp
  movaps %xmm0,%xmm2
  push   %rbx
  mov    %rdi,%rbx
  addss  %xmm1,%xmm2
  subss  %xmm0,%xmm1
  movaps %xmm5,%xmm0
  movaps %xmm2,%xmm3
  subss  %xmm5,%xmm2
  addss  %xmm1,%xmm0
  addss  %xmm5,%xmm3
  subss  %xmm5,%xmm1
  subss  %xmm4,%xmm2
  subss  %xmm4,%xmm0
  addss  %xmm4,%xmm3
  addss  %xmm4,%xmm1
  comiss %xmm3,%xmm2
  jbe    L4
  movaps %xmm2,%xmm3
  mov    $0x1,%r12d
L4:
  comiss %xmm3,%xmm0
  ja     L53
  comiss %xmm3,%xmm1
  ja     L45
  movaps %xmm3,%xmm0
  call   <sqrtf@plt>
  pxor   %xmm1,%xmm1
  movss  D23(%rip),%xmm0
  add    %eax,%eax
  cvtsi2ss %eax,%xmm1
  mulss  %xmm1,%xmm0
  test   %r12d,%r12d
  jne    L17
  movss  %xmm0,0x0(%rbp)
  movss  0x14(%rbx),%xmm2
  subss  0x1c(%rbx),%xmm2
  divss  %xmm1,%xmm2
  movss  %xmm2,0x4(%rbp)
  movss  0x18(%rbx),%xmm2
  subss  0x8(%rbx),%xmm2
  divss  %xmm1,%xmm2
  movss  %xmm2,0x8(%rbp)
  movss  0x4(%rbx),%xmm2
  subss  0xc(%rbx),%xmm2
  divss  %xmm1,%xmm2
  pxor   %xmm1,%xmm1
  comiss %xmm0,%xmm1
  movss  %xmm2,0xc(%rbp)
  ja     L41
L69:
  pop    %rbx
  pop    %rbp
  pop    %r12
  ret
  nopl   0x0(%rax,%rax,1)
L53:
  comiss %xmm0,%xmm1
  jbe    L66
L45:
  movaps %xmm1,%xmm0
  call   <sqrtf@plt>
  pxor   %xmm2,%xmm2
  movss  D23(%rip),%xmm0
  add    %eax,%eax
  cvtsi2ss %eax,%xmm2
  mulss  %xmm2,%xmm0
  movss  %xmm0,0xc(%rbp)
  movss  0x4(%rbx),%xmm0
  subss  0xc(%rbx),%xmm0
  divss  %xmm2,%xmm0
  movss  %xmm0,0x0(%rbp)
  movss  0x8(%rbx),%xmm1
  addss  0x18(%rbx),%xmm1
  divss  %xmm2,%xmm1
  movss  %xmm1,0x4(%rbp)
  movss  0x14(%rbx),%xmm1
  addss  0x1c(%rbx),%xmm1
  divss  %xmm2,%xmm1
  movss  %xmm1,0x8(%rbp)
L37:
  pxor   %xmm1,%xmm1
  comiss %xmm0,%xmm1
  jbe    L69
L41:
  movss  D40(%rip),%xmm1
  pop    %rbx
  xorps  %xmm1,%xmm0
  movss  %xmm0,0x0(%rbp)
  movss  0x4(%rbp),%xmm0
  xorps  %xmm1,%xmm0
  movss  %xmm0,0x4(%rbp)
  movss  0x8(%rbp),%xmm0
  xorps  %xmm1,%xmm0
  movss  %xmm0,0x8(%rbp)
  movss  0xc(%rbp),%xmm0
  xorps  %xmm1,%xmm0
  movss  %xmm0,0xc(%rbp)
  pop    %rbp
  pop    %r12
  ret
  nopw   0x0(%rax,%rax,1)
L66:
  call   <sqrtf@plt>
  pxor   %xmm2,%xmm2
  movss  D23(%rip),%xmm0
  add    %eax,%eax
  cvtsi2ss %eax,%xmm2
  mulss  %xmm2,%xmm0
  movss  %xmm0,0x8(%rbp)
  movss  0x18(%rbx),%xmm0
  subss  0x8(%rbx),%xmm0
  divss  %xmm2,%xmm0
  movss  %xmm0,0x0(%rbp)
  movss  0x4(%rbx),%xmm1
  addss  0xc(%rbx),%xmm1
  divss  %xmm2,%xmm1
  movss  %xmm1,0x4(%rbp)
  movss  0x14(%rbx),%xmm1
  addss  0x1c(%rbx),%xmm1
  divss  %xmm2,%xmm1
  movss  %xmm1,0xc(%rbp)
  jmp    L37
L17:
  movss  %xmm0,0x4(%rbp)
  movss  0x14(%rbx),%xmm0
  subss  0x1c(%rbx),%xmm0
  divss  %xmm1,%xmm0
  movss  %xmm0,0x0(%rbp)
  movss  0x4(%rbx),%xmm2
  addss  0xc(%rbx),%xmm2
  divss  %xmm1,%xmm2
  movss  %xmm2,0x8(%rbp)
  movss  0x8(%rbx),%xmm2
  addss  0x18(%rbx),%xmm2
  divss  %xmm1,%xmm2
  movss  %xmm2,0xc(%rbp)
  jmp    L37
```

tool_calls: name: parse_data
arguments: data_tag: D25
type: byte[4]
name: parse_data
arguments: data_tag: D23
type: byte[4]
name: parse_data
arguments: data_tag: D40
type: byte[4]

tool: {"D25": [0x00,0x00,0x80,0x3f]}

tool: {"D23": [0x00,0x00,0x80,0x3e]}

tool: {"D40": [0x00,0x00,0x00,0x80]}

assistant: ```c
void R2Quaternion(float R[3][3], float q[4]) {
  float m[4], mag;
  uint8_t index, i;
  m[0] = 1 + R[0][0] + R[1][1] + R[2][2];
  m[1] = 1 + R[0][0] - R[1][1] - R[2][2];
  m[2] = 1 - R[0][0] + R[1][1] - R[2][2];
  m[3] = 1 - R[0][0] - R[1][1] + R[2][2];
  index = 0;
  mag = m[0];
  for (i = 1; i < 4; i++) {
    if (m[i] > mag) {
      mag = m[i];
      index = i;
    }
  }
  mag = 2 * sqrtf(mag);
  if (index == 0) {
    q[0] = mag / 4;
    q[1] = (R[1][2] - R[2][1]) / mag;
    q[2] = (R[2][0] - R[0][2]) / mag;
    q[3] = (R[0][1] - R[1][0]) / mag;
  } else if (index == 1) {
    q[1] = mag / 4;
    q[0] = (R[1][2] - R[2][1]) / mag;
    q[2] = (R[0][1] + R[1][0]) / mag;
    q[3] = (R[0][2] + R[2][0]) / mag;
  } else if (index == 2) {
    q[2] = mag / 4;
    q[0] = (R[2][0] - R[0][2]) / mag;
    q[1] = (R[0][1] + R[1][0]) / mag;
    q[3] = (R[1][2] + R[2][1]) / mag;
  } else {
    q[3] = mag / 4;
    q[0] = (R[0][1] - R[1][0]) / mag;
    q[1] = (R[0][2] + R[2][0]) / mag;
    q[2] = (R[1][2] + R[2][1]) / mag;
  }
  if (q[0] < 0) {
    q[0] = -q[0];
    q[1] = -q[1];
    q[2] = -q[2];
    q[3] = -q[3];
  }
}
```