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

```asm
<sncndn>:
  endbr64
  push   %r15
  pxor   %xmm7,%xmm7
  push   %r14
  mov    %rdx,%r14
  push   %r13
  mov    %rdi,%r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0xb8,%rsp
  mov    %rsi,0x18(%rsp)
  movss  %xmm0,0x20(%rsp)
  movss  %xmm1,0x14(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0xa8(%rsp)
  xor    %eax,%eax
  ucomiss %xmm7,%xmm1
  jp     L4
  je     L38
L4:
  pxor   %xmm7,%xmm7
  comiss 0x14(%rsp),%xmm7
  ja     L28
  movss  D52(%rip),%xmm7
  movss  0x14(%rsp),%xmm3
  movss  %xmm7,0x24(%rsp)
L10:
  movss  0x24(%rsp),%xmm6
  mov    $0x1,%r15d
  lea    0x30(%rsp),%rbp
  lea    0x70(%rsp),%r12
  movss  %xmm6,(%r14)
  movss  %xmm6,0x8(%rsp)
  nopw   0x0(%rax,%rax,1)
L105:
  movss  0x8(%rsp),%xmm5
  movaps %xmm3,%xmm0
  mov    %r15d,%ebx
  movss  %xmm5,0x0(%rbp,%r15,4)
  call   <sqrt@plt>
  movss  0x8(%rsp),%xmm2
  pxor   %xmm1,%xmm1
  cvtsd2ss %xmm0,%xmm1
  movaps %xmm1,%xmm0
  movss  %xmm1,(%r12,%r15,4)
  addss  %xmm2,%xmm0
  mulss  D27(%rip),%xmm0
  movss  %xmm2,0x10(%rsp)
  movss  %xmm1,0xc(%rsp)
  movss  %xmm0,0x8(%rsp)
  movaps %xmm2,%xmm0
  subss  %xmm1,%xmm0
  call   <fabs@plt>
  movss  0x10(%rsp),%xmm2
  movss  0xc(%rsp),%xmm1
  movapd %xmm0,%xmm3
  pxor   %xmm0,%xmm0
  cvtss2sd %xmm2,%xmm0
  mulsd  D71(%rip),%xmm0
  comisd %xmm3,%xmm0
  jae    L129
  movaps %xmm1,%xmm3
  add    $0x1,%r15
  mulss  %xmm2,%xmm3
  cmp    $0xe,%r15
  jne    L105
L129:
  movss  0x20(%rsp),%xmm3
  mulss  0x8(%rsp),%xmm3
  movss  %xmm1,0x28(%rsp)
  movss  %xmm2,0x10(%rsp)
  movaps %xmm3,%xmm0
  movss  %xmm3,0xc(%rsp)
  call   <sin@plt>
  movss  0xc(%rsp),%xmm3
  movss  %xmm0,0x0(%r13)
  movaps %xmm3,%xmm0
  call   <cos@plt>
  mov    0x18(%rsp),%rax
  pxor   %xmm7,%xmm7
  movss  0x10(%rsp),%xmm2
  movss  0x28(%rsp),%xmm1
  movss  %xmm0,(%rax)
  movss  0x0(%r13),%xmm3
  ucomiss %xmm7,%xmm3
  jp     L104
  je     L76
L104:
  divss  %xmm3,%xmm0
  movss  0x8(%rsp),%xmm3
  sub    $0x1,%ebx
  movss  (%r14),%xmm4
  movslq %ebx,%rbx
  mulss  %xmm0,%xmm3
  jmp    L89
  nopl   0x0(%rax,%rax,1)
L50:
  movss  (%r12,%rbx,4),%xmm1
  movss  0x0(%rbp,%rbx,4),%xmm2
  sub    $0x1,%rbx
L89:
  mulss  %xmm3,%xmm0
  mulss  %xmm4,%xmm3
  movaps %xmm1,%xmm4
  addss  %xmm0,%xmm4
  addss  %xmm2,%xmm0
  divss  %xmm0,%xmm4
  movaps %xmm3,%xmm0
  divss  %xmm2,%xmm0
  test   %ebx,%ebx
  jne    L50
  movaps %xmm3,%xmm0
  movss  %xmm4,(%r14)
  mulss  %xmm3,%xmm0
  movss  %xmm3,0x8(%rsp)
  addss  0x24(%rsp),%xmm0
  call   <sqrt@plt>
  pxor   %xmm7,%xmm7
  movss  0x8(%rsp),%xmm3
  movapd %xmm0,%xmm1
  movsd  D19(%rip),%xmm0
  divsd  %xmm1,%xmm0
  movss  0x0(%r13),%xmm1
  comiss %xmm7,%xmm1
  cvtsd2ss %xmm0,%xmm0
  jb     L82
L135:
  movss  %xmm0,0x0(%r13)
  mulss  %xmm3,%xmm0
  mov    0x18(%rsp),%rax
  movss  %xmm0,(%rax)
L76:
  pxor   %xmm7,%xmm7
  comiss 0x14(%rsp),%xmm7
  jbe    L11
  mov    0x18(%rsp),%rax
  movss  (%r14),%xmm0
  movss  (%rax),%xmm1
  movss  %xmm1,(%r14)
  movss  %xmm0,(%rax)
  movss  0x0(%r13),%xmm0
  divss  0x2c(%rsp),%xmm0
  movss  %xmm0,0x0(%r13)
L11:
  mov    0xa8(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L123
  add    $0xb8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax)
L38:
  movss  0x20(%rsp),%xmm0
  call   <cosh@plt>
  mov    0x18(%rsp),%rax
  movapd %xmm0,%xmm1
  movsd  D19(%rip),%xmm0
  divsd  %xmm1,%xmm0
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%rax)
  movss  %xmm0,(%r14)
  movss  0x20(%rsp),%xmm0
  call   <tanh@plt>
  movss  %xmm0,0x0(%r13)
  jmp    L11
L28:
  movss  D52(%rip),%xmm6
  movss  0x14(%rsp),%xmm7
  pxor   %xmm2,%xmm2
  pxor   %xmm3,%xmm3
  movsd  D66(%rip),%xmm0
  movaps %xmm6,%xmm1
  cvtss2sd %xmm7,%xmm3
  movss  %xmm6,0x24(%rsp)
  subss  %xmm7,%xmm1
  cvtss2sd %xmm1,%xmm2
  divsd  %xmm2,%xmm0
  divsd  %xmm0,%xmm3
  movaps %xmm1,%xmm0
  cvtsd2ss %xmm3,%xmm3
  movss  %xmm3,0x8(%rsp)
  call   <sqrt@plt>
  movss  0x20(%rsp),%xmm7
  pxor   %xmm6,%xmm6
  movss  0x8(%rsp),%xmm3
  cvtsd2ss %xmm0,%xmm6
  movss  %xmm6,0x2c(%rsp)
  mulss  %xmm6,%xmm7
  movss  %xmm7,0x20(%rsp)
  jmp    L10
  cs nopw 0x0(%rax,%rax,1)
L82:
  xorps  D29(%rip),%xmm0
  jmp    L135
L123:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: extract_symbolic_data
arguments: data_tag: D52
type: f32
name: extract_symbolic_data
arguments: data_tag: D27
type: f32
name: extract_symbolic_data
arguments: data_tag: D71
type: f64
name: extract_symbolic_data
arguments: data_tag: D19
type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D66
type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D29
type: byte[16]

tool: D52: 1

tool: D27: 0.5

tool: D71: 0.0003

tool: {"D19": [0.0,1.875]}

tool: {"D66": [0.0,-1.875]}

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

assistant: ```c
void sncndn(float uu, float emmc, float *sn, float *cn, float *dn) {
  float a, b, c, d, emc, u;
  float em[14], en[14];
  int i, ii, l, bo;
  emc = emmc;
  u = uu;
  if (emc) {
    bo = (emc < 0.0);
    if (bo) {
      d = 1.0 - emc;
      emc /= -1.0 / d;
      u *= (d = sqrt(d));
    }
    a = 1.0;
    *dn = 1.0;
    for (i = 1; i <= 13; i++) {
      l = i;
      em[i] = a;
      en[i] = (emc = sqrt(emc));
      c = 0.5 * (a + emc);
      if (fabs(a - emc) <= 0.0003 * a)
        break;
      emc *= a;
      a = c;
    }
    u *= c;
    *sn = sin(u);
    *cn = cos(u);
    if (*sn) {
      a = (*cn) / (*sn);
      c *= a;
      for (ii = l; ii >= 1; ii--) {
        b = em[ii];
        a *= c;
        c *= (*dn);
        *dn = (en[ii] + a) / (b + a);
        a = c / b;
      }
      a = 1.0 / sqrt(c * c + 1.0);
      *sn = (*sn >= 0.0 ? a : -a);
      *cn = c * (*sn);
    }
    if (bo) {
      a = (*dn);
      *dn = (*cn);
      *cn = a;
      *sn /= d;
    }
  } else {
    *cn = 1.0 / cosh(u);
    *dn = (*cn);
    *sn = tanh(u);
  }
}
```