user: 请将这段汇编代码片段转换成等效的C代码：

```asm
<anomaly>:
  endbr64
  push   %r13
  mov    %rdi,%r13
  push   %r12
  mov    %rsi,%r12
  push   %rbp
  push   %rbx
  sub    $0x48,%rsp
  movsd  D115(%rip),%xmm6
  movsd  %xmm0,0x20(%rsp)
  comisd %xmm1,%xmm6
  movsd  %xmm1,0x8(%rsp)
  jbe    L76
  movapd %xmm0,%xmm7
  movsd  D129(%rip),%xmm0
  movapd %xmm7,%xmm1
  divsd  %xmm0,%xmm1
  cvttsd2si %xmm1,%rax
  pxor   %xmm1,%xmm1
  cvtsi2sd %rax,%xmm1
  mulsd  %xmm0,%xmm1
  subsd  %xmm1,%xmm7
  comisd D131(%rip),%xmm7
  movsd  %xmm7,0x28(%rsp)
  jbe    L28
  subsd  %xmm0,%xmm7
  movsd  %xmm7,0x28(%rsp)
L28:
  movsd  D53(%rip),%xmm1
  movsd  0x28(%rsp),%xmm6
  comisd %xmm6,%xmm1
  ja     L62
L35:
  movsd  0x28(%rsp),%xmm1
  lea    0x38(%rsp),%rbp
  lea    0x30(%rsp),%rbx
  jmp    L24
L17:
  movsd  0x18(%rsp),%xmm5
  mulsd  0x8(%rsp),%xmm5
  movsd  D115(%rip),%xmm3
  movsd  D28(%rip),%xmm6
  movsd  0x20(%rsp),%xmm0
  subsd  %xmm5,%xmm3
  maxsd  %xmm3,%xmm6
  divsd  %xmm6,%xmm0
  subsd  %xmm0,%xmm1
L24:
  movapd %xmm1,%xmm0
  mov    %rbx,%rsi
  mov    %rbp,%rdi
  movsd  %xmm1,0x10(%rsp)
  call   <sincos@plt>
  movsd  0x8(%rsp),%xmm2
  mulsd  0x38(%rsp),%xmm2
  movsd  0x10(%rsp),%xmm1
  movsd  0x30(%rsp),%xmm5
  movapd %xmm1,%xmm0
  movsd  %xmm5,0x18(%rsp)
  subsd  %xmm2,%xmm0
  subsd  0x28(%rsp),%xmm0
  movsd  %xmm0,0x20(%rsp)
  call   <fabs@plt>
  movsd  0x10(%rsp),%xmm1
  test   %eax,%eax
  jg     L17
  movsd  0x8(%rsp),%xmm6
  movsd  D115(%rip),%xmm0
  movsd  %xmm1,0x10(%rsp)
  movsd  D115(%rip),%xmm2
  addsd  %xmm6,%xmm0
  subsd  %xmm6,%xmm2
  divsd  %xmm2,%xmm0
  call   <sqrt@plt>
  movsd  0x10(%rsp),%xmm1
  mov    %eax,%ebx
  movapd %xmm1,%xmm0
  movsd  %xmm1,0x8(%rsp)
  mulsd  D1(%rip),%xmm0
  call   <tan@plt>
  imul   %eax,%ebx
  mov    %ebx,%edi
  call   <atan@plt>
  pxor   %xmm0,%xmm0
  movsd  0x8(%rsp),%xmm1
  add    %eax,%eax
  cvtsi2sd %eax,%xmm0
L66:
  movsd  %xmm0,0x0(%r13)
  movsd  %xmm1,(%r12)
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  ret
  nopl   0x0(%rax)
L62:
  addsd  %xmm6,%xmm0
  movsd  %xmm0,0x28(%rsp)
  jmp    L35
  nop
L76:
  movsd  0x20(%rsp),%xmm0
  call   <fabs@plt>
  pxor   %xmm5,%xmm5
  movsd  0x8(%rsp),%xmm6
  xor    %esi,%esi
  cvtsi2sd %eax,%xmm5
  movapd %xmm6,%xmm2
  movapd %xmm6,%xmm7
  subsd  D115(%rip),%xmm7
  mulsd  %xmm6,%xmm2
  movapd %xmm5,%xmm0
  movapd %xmm5,%xmm1
  movq   %xmm5,%rbx
  movsd  %xmm7,0x28(%rsp)
  mulsd  D58(%rip),%xmm0
  divsd  %xmm7,%xmm1
  divsd  %xmm2,%xmm0
  movsd  %xmm1,0x10(%rsp)
  cvttsd2si %xmm0,%edi
  call   <pow@plt>
  movsd  0x10(%rsp),%xmm1
  mov    D115(%rip),%rax
  minsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  jmp    L54
  nopl   (%rax)
L0:
  movapd %xmm1,%xmm0
  movsd  %xmm1,0x10(%rsp)
  call   <sinh@plt>
  movsd  0x10(%rsp),%xmm1
  movsd  %xmm0,0x18(%rsp)
  movapd %xmm1,%xmm0
  call   <cosh@plt>
  movsd  0x8(%rsp),%xmm4
  movsd  0x18(%rsp),%xmm3
  movapd %xmm0,%xmm2
  movq   %rbx,%xmm0
  movsd  0x10(%rsp),%xmm1
  mulsd  %xmm4,%xmm3
  mulsd  %xmm4,%xmm2
  subsd  D115(%rip),%xmm2
  subsd  %xmm3,%xmm0
  addsd  %xmm1,%xmm0
  divsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
L54:
  movsd  %xmm1,0x10(%rsp)
  call   <fabs@plt>
  movsd  0x10(%rsp),%xmm1
  test   %eax,%eax
  jg     L0
  pxor   %xmm0,%xmm0
  comisd 0x20(%rsp),%xmm0
  jbe    L44
  xorpd  D145(%rip),%xmm1
L44:
  movsd  D115(%rip),%xmm0
  movsd  %xmm1,0x10(%rsp)
  addsd  0x8(%rsp),%xmm0
  divsd  0x28(%rsp),%xmm0
  call   <sqrt@plt>
  movsd  0x10(%rsp),%xmm1
  mov    %eax,%ebx
  movapd %xmm1,%xmm0
  movsd  %xmm1,0x8(%rsp)
  mulsd  D1(%rip),%xmm0
  call   <tanh@plt>
  mov    %ebx,%edi
  imul   %eax,%edi
  call   <atan@plt>
  pxor   %xmm0,%xmm0
  movsd  0x8(%rsp),%xmm1
  add    %eax,%eax
  cvtsi2sd %eax,%xmm0
  jmp    L66
```

tool_calls: name: retrieve_data
arguments: label: D115
data_type: f64
name: retrieve_data
arguments: label: D129
data_type: byte[8]
name: retrieve_data
arguments: label: D131
data_type: f64[2]
name: retrieve_data
arguments: label: D53
data_type: byte[8]
name: retrieve_data
arguments: label: D28
data_type: byte[8]
name: retrieve_data
arguments: label: D1
data_type: f32[2]
name: retrieve_data
arguments: label: D58
data_type: f32[2]
name: retrieve_data
arguments: label: D145
data_type: f32[4]

tool: D115: 1

tool: {"D129": [0x18,0x2d,0x44,0x54,0xfb,0x21,0x19,0x40]}

tool: {"D131": [3.141592653589793,-3.141592653589793]}

tool: {"D53": [0x18,0x2d,0x44,0x54,0xfb,0x21,0x09,0xc0]}

tool: {"D28": [0x9a,0x99,0x99,0x99,0x99,0x99,0xb9,0x3f]}

tool: {"D1": [0.0,1.75]}

tool: {"D58": [0.0,2.375]}

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

assistant: ```c
void anomaly(double ma, double s, double *nu, double *ea) {
  double m, fea, corr;
  if (s < 1.0) {
    double dla;
    m = ma - (2 * 3.141592653589793) * (long)(ma / (2 * 3.141592653589793));
    if (m > 3.141592653589793)
      m -= (2 * 3.141592653589793);
    if (m < -3.141592653589793)
      m += (2 * 3.141592653589793);
    fea = m;
    for (;;) {
      dla = fea - (s * sin(fea)) - m;
      if (fabs(dla) < (1e-8))
        break;
      corr = 1 - (s * cos(fea));
      if (corr < .1)
        corr = .1;
      dla /= corr;
      fea -= dla;
    }
    *nu = 2 * atan(sqrt((1 + s) / (1 - s)) * tan(fea / 2));
  } else {
    double fea1;
    m = fabs(ma);
    fea = m / (s - 1.);
    fea1 = pow(6 * m / (s * s), 1. / 3.);
    if (fea1 < fea)
      fea = fea1;
    corr = 1;
    while (fabs(corr) > (1e-8)) {
      corr = (m - s * sinh(fea) + fea) / (s * cosh(fea) - 1);
      fea += corr;
    }
    if (ma < 0.)
      fea = -fea;
    *nu = 2 * atan(sqrt((s + 1) / (s - 1)) * tanh(fea / 2));
  }
  *ea = fea;
}
```