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

```asm
<jacobi>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %rbx
  sub    $0x128,%rsp
  mov    %rdi,-0x108(%rbp)
  mov    %rsi,-0x110(%rbp)
  mov    %rdx,-0x118(%rbp)
  mov    %ecx,-0x11c(%rbp)
  movsd  %xmm0,-0x128(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x18(%rbp)
  xor    %eax,%eax
  mov    -0x108(%rbp),%rax
  mov    %rax,%rdi
  call   <len@plt>
  movq   %xmm0,%rax
  mov    %rax,-0xf8(%rbp)
  mov    -0x108(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm1
  mov    -0x110(%rbp),%rax
  add    $0x10,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  mov    -0x108(%rbp),%rax
  add    $0x10,%rax
  movsd  (%rax),%xmm2
  mov    -0x110(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm1
  mulsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0xd0(%rbp)
  mov    -0x108(%rbp),%rax
  add    $0x10,%rax
  movsd  (%rax),%xmm1
  mov    -0x110(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  mov    -0x108(%rbp),%rax
  movsd  (%rax),%xmm2
  mov    -0x110(%rbp),%rax
  add    $0x10,%rax
  movsd  (%rax),%xmm1
  mulsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0xc8(%rbp)
  mov    -0x108(%rbp),%rax
  movsd  (%rax),%xmm1
  mov    -0x110(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  mov    -0x108(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm2
  mov    -0x110(%rbp),%rax
  movsd  (%rax),%xmm1
  mulsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0xc0(%rbp)
  lea    -0xd0(%rbp),%rax
  mov    %rax,%rdi
  call   <len@plt>
  movq   %xmm0,%rax
  movsd  -0xf8(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  %xmm0,%xmm1
  movq   %rax,%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0xf0(%rbp)
  movsd  D3(%rip),%xmm0
  movsd  %xmm0,-0xe8(%rbp)
  movl   $0x0,-0xfc(%rbp)
  jmp    L46
L36:
  mov    -0xfc(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x108(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movapd %xmm0,%xmm1
  divsd  -0xf8(%rbp),%xmm1
  movsd  -0xf8(%rbp),%xmm0
  subsd  -0xe8(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  mov    -0xfc(%rbp),%eax
  cltq
  movsd  %xmm0,-0xb0(%rbp,%rax,8)
  mov    -0xfc(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x108(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movapd %xmm0,%xmm1
  divsd  -0xf8(%rbp),%xmm1
  movsd  -0xf8(%rbp),%xmm0
  addsd  -0xe8(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  mov    -0xfc(%rbp),%eax
  cltq
  movsd  %xmm0,-0x90(%rbp,%rax,8)
  addl   $0x1,-0xfc(%rbp)
L46:
  cmpl   $0x2,-0xfc(%rbp)
  jle    L36
  mov    -0x11c(%rbp),%ecx
  mov    -0x118(%rbp),%rdx
  lea    -0x70(%rbp),%rsi
  lea    -0xb0(%rbp),%rax
  mov    %rax,%rdi
  call   <force@plt>
  mov    -0x11c(%rbp),%ecx
  mov    -0x118(%rbp),%rdx
  lea    -0x50(%rbp),%rsi
  lea    -0x90(%rbp),%rax
  mov    %rax,%rdi
  call   <force@plt>
  movl   $0x0,-0xfc(%rbp)
  jmp    L45
L27:
  mov    -0xfc(%rbp),%eax
  cltq
  movsd  -0xb0(%rbp,%rax,8),%xmm0
  mov    -0xfc(%rbp),%eax
  cltq
  movsd  -0x90(%rbp,%rax,8),%xmm1
  subsd  %xmm1,%xmm0
  mov    -0xfc(%rbp),%eax
  cltq
  movsd  %xmm0,-0x30(%rbp,%rax,8)
  addl   $0x1,-0xfc(%rbp)
L45:
  cmpl   $0x2,-0xfc(%rbp)
  jle    L27
  lea    -0x70(%rbp),%rax
  mov    %rax,%rdi
  call   <len@plt>
  movq   %xmm0,%rbx
  lea    -0x50(%rbp),%rax
  mov    %rax,%rdi
  call   <len@plt>
  movq   %rbx,%xmm3
  subsd  %xmm0,%xmm3
  movsd  %xmm3,-0x130(%rbp)
  lea    -0x30(%rbp),%rax
  mov    %rax,%rdi
  call   <len@plt>
  movapd %xmm0,%xmm1
  movsd  -0x130(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0xe0(%rbp)
  movsd  -0x128(%rbp),%xmm1
  movsd  D29(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0xf0(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  addsd  -0xe0(%rbp),%xmm0
  movq   D22(%rip),%xmm2
  andpd  %xmm2,%xmm0
  divsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  mov    $0x0,%edi
  movq   %rax,%xmm0
  call   <pow@plt>
  movq   %xmm0,%rax
  mov    %rax,-0xd8(%rbp)
  movsd  -0xd8(%rbp),%xmm0
  movq   %xmm0,%rax
  mov    -0x18(%rbp),%rdx
  sub    %fs:0x28,%rdx
  je     L23
  call   <__stack_chk_fail@plt>
L23:
  movq   %rax,%xmm0
  mov    -0x8(%rbp),%rbx
  leave
  ret
```

tool_calls: name: read_data
arguments: tag: D3
type_infered: byte[8]
name: read_data
arguments: tag: D29
type_infered: byte[8]
name: read_data
arguments: tag: D22
type_infered: qword

tool: {"D3": [0x30,0xea,0xd2,0xf6,0x02,0x21,0xa1,0x43]}

tool: {"D29": [0xf1,0xd9,0x6d,0x61,0x98,0x55,0xd2,0x3d]}

tool: D22: 9223372036854775807

assistant: ```c
double jacobi(double *x, double *v, double *par, int potential, double Mcl) {
  int i;
  double R, om, dpot, delta, r;
  double omega[3], x1[3], x2[3], a1[3], a2[3], dx[3];
  r = len(x);
  omega[0] = x[1] * v[2] - x[2] * v[1];
  omega[1] = x[2] * v[0] - x[0] * v[2];
  omega[2] = x[0] * v[1] - x[1] * v[0];
  om = len(omega) / (r * r);
  delta = 0.02 * 3.08567758e19;
  for (i = 0; i < 3; i++) {
    x1[i] = x[i] / r * (r - delta);
    x2[i] = x[i] / r * (r + delta);
  }
  force(x1, a1, par, potential);
  force(x2, a2, par, potential);
  for (i = 0; i < 3; i++) {
    dx[i] = x1[i] - x2[i];
  }
  dpot = (len(a1) - len(a2)) / len(dx);
  R = pow(6.67e-11 * Mcl / fabs(om * om + dpot), 1. / 3.);
  return R;
}
```