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

```asm
<fpsimple_9j>:
  endbr64
  push   %r15
  mov    %esi,%r15d
  push   %r14
  push   %r13
  mov    %edx,%r13d
  push   %r12
  mov    %r8d,%r12d
  push   %rbp
  push   %rbx
  sub    $0x88,%rsp
  mov    %rdi,0x40(%rsp)
  mov    0xc0(%rsp),%r14d
  mov    %edx,0x10(%rsp)
  mov    0xd0(%rsp),%ebx
  mov    %ecx,(%rsp)
  mov    %r9d,0x14(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x78(%rsp)
  mov    %esi,%eax
  sub    0xd8(%rsp),%eax
  mov    %eax,%edx
  neg    %edx
  cmovs  %eax,%edx
  mov    %r13d,%eax
  sub    %r14d,%eax
  mov    %eax,%esi
  neg    %esi
  cmovs  %eax,%esi
  mov    %ebx,%eax
  sub    %r8d,%eax
  mov    %eax,%edi
  neg    %edi
  cmovs  %eax,%edi
  call   <max3@plt>
  mov    0xd8(%rsp),%edx
  lea    0x0(%r13,%r14,1),%esi
  lea    (%rbx,%r12,1),%edi
  mov    %eax,%ebp
  add    %r15d,%edx
  call   <min3@plt>
  mov    (%rsp),%r11d
  cmp    %eax,%ebp
  mov    %eax,0x24(%rsp)
  jg     L9
  lea    0x58(%rsp),%rax
  movsd  D5(%rip),%xmm7
  mov    %r11d,0x20(%rsp)
  pxor   %xmm5,%xmm5
  mov    %rax,0x28(%rsp)
  lea    0x60(%rsp),%rax
  lea    0x70(%rsp),%r13
  mov    %rax,0x30(%rsp)
  lea    0x68(%rsp),%rax
  mov    %rax,0x38(%rsp)
  movsd  %xmm5,0x48(%rsp)
  movsd  %xmm5,(%rsp)
  movsd  %xmm7,0x18(%rsp)
  nopl   0x0(%rax)
L44:
  sub    $0x8,%rsp
  mov    %ebp,%r9d
  mov    %r14d,%ecx
  mov    %r15d,%edi
  push   0x30(%rsp)
  mov    0xe8(%rsp),%r8d
  mov    0x30(%rsp),%edx
  mov    0x20(%rsp),%esi
  call   <fpsimple_factor_6j@plt>
  pop    %rax
  push   0x38(%rsp)
  mov    %ebp,%r9d
  mov    0x20(%rsp),%r8d
  mov    %ebx,%ecx
  mov    %r12d,%esi
  mov    %r14d,%edi
  mov    0x24(%rsp),%edx
  call   <fpsimple_factor_6j@plt>
  pop    %rdx
  push   0x40(%rsp)
  mov    %ebp,%r9d
  mov    0xd8(%rsp),%edx
  mov    %r12d,%r8d
  mov    %r15d,%ecx
  mov    %ebx,%edi
  mov    0xe8(%rsp),%esi
  call   <fpsimple_factor_6j@plt>
  movsd  0x68(%rsp),%xmm0
  mov    %ebp,%edx
  mov    %r13,%rcx
  mulsd  0x70(%rsp),%xmm0
  movsd  0x28(%rsp),%xmm3
  mov    %r15d,%edi
  mulsd  0x78(%rsp),%xmm0
  mov    0xe8(%rsp),%esi
  movsd  %xmm3,0x80(%rsp)
  movsd  %xmm0,0x18(%rsp)
  call   <fpsimple_delta_coeff@plt>
  mov    0x20(%rsp),%esi
  mov    %ebp,%edx
  mov    %r13,%rcx
  mov    %r14d,%edi
  call   <fpsimple_delta_coeff@plt>
  mov    %ebp,%edx
  mov    %r13,%rcx
  mov    %r12d,%esi
  mov    %ebx,%edi
  call   <fpsimple_delta_coeff@plt>
  lea    0x1(%rbp),%eax
  pxor   %xmm0,%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm0
  mov    %ebp,%eax
  mov    $0x1,%edx
  add    $0x2,%ebp
  mulsd  0x80(%rsp),%xmm0
  and    $0x1,%eax
  mulsd  0x18(%rsp),%xmm0
  add    %eax,%eax
  sub    %eax,%edx
  cvtsi2sd %edx,%xmm1
  mulsd  %xmm1,%xmm0
  addsd  0x10(%rsp),%xmm0
  movsd  %xmm0,0x10(%rsp)
  pop    %rcx
  pop    %rsi
  cmp    %ebp,0x24(%rsp)
  jge    L44
  mov    0x20(%rsp),%r11d
L53:
  movsd  0x18(%rsp),%xmm6
  mov    0x10(%rsp),%esi
  mov    %r11d,%edx
  mov    %r13,%rcx
  mov    %r15d,%edi
  mov    %r11d,0x8(%rsp)
  movsd  %xmm6,0x70(%rsp)
  call   <fpsimple_delta_coeff@plt>
  mov    0x14(%rsp),%ebp
  mov    %r13,%rcx
  mov    %r14d,%edx
  mov    %r12d,%edi
  mov    %ebp,%esi
  call   <fpsimple_delta_coeff@plt>
  mov    0xd8(%rsp),%edx
  mov    0xc8(%rsp),%edi
  mov    %ebx,%esi
  mov    %r13,%rcx
  call   <fpsimple_delta_coeff@plt>
  mov    %r13,%rcx
  mov    %r12d,%esi
  mov    %r15d,%edi
  mov    0xc8(%rsp),%edx
  call   <fpsimple_delta_coeff@plt>
  mov    0x10(%rsp),%edi
  mov    %r13,%rcx
  mov    %ebx,%edx
  mov    %ebp,%esi
  call   <fpsimple_delta_coeff@plt>
  mov    0x8(%rsp),%edi
  mov    %r13,%rcx
  mov    %r14d,%esi
  mov    0xd8(%rsp),%edx
  call   <fpsimple_delta_coeff@plt>
  movsd  0x70(%rsp),%xmm0
  movsd  0x48(%rsp),%xmm7
  ucomisd %xmm0,%xmm7
  ja     L55
  sqrtsd %xmm0,%xmm0
L58:
  mulsd  (%rsp),%xmm0
  mov    0x40(%rsp),%rax
  movsd  %xmm0,(%rax)
  mov    0x78(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L50
  add    $0x88,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  cs nopw 0x0(%rax,%rax,1)
L9:
  movsd  D5(%rip),%xmm6
  pxor   %xmm7,%xmm7
  lea    0x70(%rsp),%r13
  movsd  %xmm7,0x48(%rsp)
  movsd  %xmm7,(%rsp)
  movsd  %xmm6,0x18(%rsp)
  jmp    L53
L55:
  call   <sqrt@plt>
  jmp    L58
L50:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: retrieve_data
arguments: data_tag: D5
data_type: f64

tool: D5: 1

assistant: ```c
void fpsimple_9j(double *result, int two_a, int two_b, int two_c, int two_d, int two_e, int two_f, int two_g, int two_h, int two_i) {
  int two_kmin = max3(abs(two_h - two_d), abs(two_b - two_f), abs(two_a - two_i));
  int two_kmax = min3(two_h + two_d, two_b + two_f, two_a + two_i);
  int two_k;
  double sum_prod = 0;
  for (two_k = two_kmin; two_k <= two_kmax; two_k += 2) {
    double f1, f2, f3;
    fpsimple_factor_6j(two_a, two_b, two_c, two_f, two_i, two_k, &f1);
    fpsimple_factor_6j(two_f, two_d, two_e, two_h, two_b, two_k, &f2);
    fpsimple_factor_6j(two_h, two_i, two_g, two_a, two_d, two_k, &f3);
    double factor = f1 * f2 * f3;
    double inner_delta = 1;
    fpsimple_delta_coeff(two_a, two_i, two_k, &inner_delta);
    fpsimple_delta_coeff(two_f, two_b, two_k, &inner_delta);
    fpsimple_delta_coeff(two_h, two_d, two_k, &inner_delta);
    double prod = (two_k + 1) * inner_delta * factor;
    sum_prod += (1 - 2 * ((two_k)&1)) * prod;
  }
  double prefact = 1;
  fpsimple_delta_coeff(two_a, two_b, two_c, &prefact);
  fpsimple_delta_coeff(two_d, two_e, two_f, &prefact);
  fpsimple_delta_coeff(two_g, two_h, two_i, &prefact);
  fpsimple_delta_coeff(two_a, two_d, two_g, &prefact);
  fpsimple_delta_coeff(two_b, two_e, two_h, &prefact);
  fpsimple_delta_coeff(two_c, two_f, two_i, &prefact);
  *result = sqrt(prefact) * sum_prod;
}
```