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

```asm
<mse_pfdr_deltaapprox>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x38,%rsp
  movsd  %xmm0,0x28(%rsp)
  movsd  %xmm1,0x8(%rsp)
  mov    %edi,%r12d
  mov    %esi,%ebp
  mov    %edx,%ebx
  mov    %ecx,%r13d
  mov    %r8d,%r14d
  mov    $0x0,%edi
  comisd D67(%rip),%xmm0
  jb     L1
  movapd %xmm0,%xmm4
  movsd  D52(%rip),%xmm0
  comisd %xmm4,%xmm0
  setae  %dil
  movzbl %dil,%edi
L1:
  call   <assert@plt>
  movsd  0x8(%rsp),%xmm7
  mov    $0x0,%edi
  comisd D67(%rip),%xmm7
  jb     L12
  movsd  D52(%rip),%xmm0
  comisd %xmm7,%xmm0
  setae  %dil
  movzbl %dil,%edi
L12:
  call   <assert@plt>
  mov    $0x1,%edi
  call   <assert@plt>
  mov    $0x1,%edi
  call   <assert@plt>
  mov    $0x1,%edi
  call   <assert@plt>
  mov    $0x1,%edi
  call   <assert@plt>
  lea    (%r12,%rbp,1),%r15d
  lea    (%rbx,%r13,1),%eax
  add    %r15d,%eax
  cmp    %r14d,%eax
  sete   %dil
  movzbl %dil,%edi
  call   <assert@plt>
  mov    %r15d,%r15d
  pxor   %xmm0,%xmm0
  cvtsi2sd %r15,%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,0x10(%rsp)
  mov    %ebp,%eax
  sub    %ebx,%eax
  pxor   %xmm5,%xmm5
  cvtsi2sd %rax,%xmm5
  movsd  %xmm5,0x18(%rsp)
  lea    (%r12,%rbx,1),%eax
  pxor   %xmm6,%xmm6
  cvtsi2sd %rax,%xmm6
  movsd  %xmm6,0x20(%rsp)
  mov    %r14d,%r14d
  movsd  0x8(%rsp),%xmm7
  xorpd  D12(%rip),%xmm7
  movapd %xmm7,%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2sd %r14,%xmm1
  call   <pow1pm1@plt>
  movsd  D52(%rip),%xmm0
  subsd  0x28(%rsp),%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  mulsd  %xmm1,%xmm0
  movsd  D62(%rip),%xmm1
  divsd  %xmm0,%xmm1
  mov    %r12d,%r12d
  pxor   %xmm0,%xmm0
  cvtsi2sd %r12,%xmm0
  movsd  0x18(%rsp),%xmm5
  mulsd  %xmm5,%xmm0
  mulsd  %xmm5,%xmm0
  mov    %ebp,%ebp
  pxor   %xmm2,%xmm2
  cvtsi2sd %rbp,%xmm2
  movsd  0x20(%rsp),%xmm6
  mulsd  %xmm6,%xmm2
  mulsd  %xmm6,%xmm2
  addsd  %xmm2,%xmm0
  mov    %ebx,%ebx
  pxor   %xmm2,%xmm2
  cvtsi2sd %rbx,%xmm2
  movsd  0x10(%rsp),%xmm3
  mulsd  %xmm3,%xmm2
  addsd  %xmm2,%xmm0
  mulsd  %xmm1,%xmm1
  movapd %xmm3,%xmm2
  mulsd  %xmm3,%xmm2
  divsd  %xmm2,%xmm1
  mulsd  %xmm1,%xmm0
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: read_data
arguments: data_tag: D67
data_type: f64
name: read_data
arguments: data_tag: D52
data_type: f64
name: read_data
arguments: data_tag: D12
data_type: byte[16]
name: read_data
arguments: data_tag: D62
data_type: byte[8]

tool: D67: 0

tool: D52: 1

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

tool: {"D62": [0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0xbf]}

assistant: ```c
double mse_pfdr_deltaapprox(const double lambda, const double gamma, const unsigned int a, const unsigned int b, const unsigned int c, const unsigned int d, const unsigned int m) {
  assert(lambda >= 0. && lambda <= 1.);
  assert(gamma >= 0. && gamma <= 1.);
  assert(a >= 0);
  assert(b >= 0);
  assert(c >= 0);
  assert(d >= 0);
  assert(a + b + c + d == m);
  const double apb = a + b;
  const double apb_sqr = apb * apb;
  const double bmc = b - c;
  const double apc = a + c;
  const double Klambda = -1. / ((1. - lambda) * pow1pm1(-gamma, (double)m));
  return Klambda * Klambda / (apb_sqr * apb_sqr) * (a * bmc * bmc + b * apc * apc + c * apb_sqr);
}
```