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

```asm
<scriptedSphereMotionSignedDistance>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x88,%rsp
  mov    %edx,0x78(%rsp)
  mov    %rcx,0x58(%rsp)
  test   %edi,%edi
  jle    L237
  mov    %esi,%r12d
  lea    (%rdi,%rdi,2),%eax
  xor    %r15d,%r15d
  mov    %r8,%rbx
  mov    %eax,0x7c(%rsp)
  mov    %r12d,%eax
  mov    %r12d,%r13d
  pxor   %xmm3,%xmm3
  shr    %eax
  and    $0xfffffffe,%r13d
  mov    %r15d,%edx
  mov    %r9,%rcx
  sub    $0x1,%eax
  mov    0xc0(%rsp),%rsi
  movl   $0x0,0x70(%rsp)
  add    $0x1,%rax
  shl    $0x4,%rax
  mov    %rsi,0x68(%rsp)
  mov    %r12d,%esi
  mov    %rax,%r14
  movslq %r13d,%rax
  lea    (%r8,%rax,8),%rax
  mov    %r13d,%r8d
  mov    %rax,0x48(%rsp)
  nopl   0x0(%rax)
L25:
  movslq 0x78(%rsp),%rax
  test   %eax,%eax
  jle    L258
  movslq 0x70(%rsp),%rdi
  mov    0xc8(%rsp),%r9
  mov    %rax,0x50(%rsp)
  movslq %edx,%rax
  shl    $0x3,%rax
  cmp    $0x1,%esi
  mov    %rbx,%rbp
  mov    %rdi,%r15
  lea    (%rcx,%rax,1),%r12
  lea    (%r9,%rdi,8),%rdi
  setne  %r11b
  lea    0x8(%rcx,%rax,1),%rax
  cmp    %rax,%rdi
  setne  %r10b
  mov    %r10d,%eax
  and    %r11d,%eax
  mov    %al,0x67(%rsp)
  lea    (%rdx,%r8,1),%eax
  cltq
  lea    (%rcx,%rax,8),%r9
  lea    (%r15,%r8,1),%eax
  mov    0xc8(%rsp),%r15
  cltq
  lea    (%r15,%rax,8),%r15
  lea    0x8(%rbx),%rax
  cmp    %rax,%rdi
  setne  %al
  xor    %r13d,%r13d
  and    %r11d,%eax
  and    %r10d,%eax
  mov    %al,0x66(%rsp)
  lea    -0x1(%rsi),%eax
  mov    %eax,0x74(%rsp)
  nopl   0x0(%rax,%rax,1)
L182:
  test   %esi,%esi
  jle    L70
  lea    0x0(%r13,%r13,2),%eax
  cmp    $0x1,%esi
  je     L3
  xor    %r10d,%r10d
  movapd %xmm3,%xmm1
  nopl   (%rax)
L152:
  movupd (%r12,%r10,1),%xmm0
  movupd 0x0(%rbp,%r10,1),%xmm5
  add    $0x10,%r10
  subpd  %xmm5,%xmm0
  mulpd  %xmm0,%xmm0
  addsd  %xmm0,%xmm1
  unpckhpd %xmm0,%xmm0
  addsd  %xmm0,%xmm1
  cmp    %r14,%r10
  jne    L152
  mov    %r8d,%r11d
  cmp    %r8d,%esi
  je     L0
L210:
  lea    (%rdx,%r11,1),%r10d
  add    %r11d,%eax
  movslq %r10d,%r10
  cltq
  movsd  (%rcx,%r10,8),%xmm0
  subsd  (%rbx,%rax,8),%xmm0
  mulsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm1
L0:
  ucomisd %xmm1,%xmm3
  ja     L107
  mov    0x58(%rsp),%rax
  sqrtsd %xmm1,%xmm1
  subsd  (%rax,%r13,8),%xmm1
  test   %r13,%r13
  je     L21
L282:
  movapd %xmm1,%xmm0
  mov    %rcx,0x38(%rsp)
  mov    %esi,0x60(%rsp)
  mov    %r8d,0x30(%rsp)
  mov    %r9,0x28(%rsp)
  mov    %rdi,0x20(%rsp)
  mov    %edx,0x18(%rsp)
  movsd  %xmm1,0x10(%rsp)
  call   <fabs@plt>
  movsd  0x40(%rsp),%xmm0
  mov    %rax,0x8(%rsp)
  call   <fabs@plt>
  cmp    %rax,0x8(%rsp)
  movsd  0x10(%rsp),%xmm1
  pxor   %xmm3,%xmm3
  mov    0x18(%rsp),%edx
  mov    0x20(%rsp),%rdi
  mov    0x28(%rsp),%r9
  mov    0x30(%rsp),%r8d
  mov    0x60(%rsp),%esi
  mov    0x38(%rsp),%rcx
  jl     L174
L11:
  add    $0x1,%r13
  add    $0x18,%rbp
  cmp    0x50(%rsp),%r13
  jne    L182
L258:
  mov    0x68(%rsp),%rax
  movsd  0x40(%rsp),%xmm7
  add    $0x3,%edx
  add    %esi,0x70(%rsp)
  movsd  %xmm7,(%rax)
  add    $0x8,%rax
  mov    %rax,0x68(%rsp)
  cmp    0x7c(%rsp),%edx
  jne    L25
L237:
  add    $0x88,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   (%rax)
L21:
  cmpb   $0x0,0x66(%rsp)
  je     L142
  movsd  D15(%rip),%xmm2
  xor    %eax,%eax
  addsd  %xmm1,%xmm2
  movapd %xmm2,%xmm4
  unpcklpd %xmm4,%xmm4
  nopl   0x0(%rax)
L167:
  movupd (%r12,%rax,1),%xmm0
  movupd (%rbx,%rax,1),%xmm6
  subpd  %xmm6,%xmm0
  divpd  %xmm4,%xmm0
  movups %xmm0,(%rdi,%rax,1)
  add    $0x10,%rax
  cmp    %rax,%r14
  jne    L167
  cmp    %r8d,%esi
  je     L15
  mov    0x48(%rsp),%rax
  movsd  (%r9),%xmm0
  subsd  (%rax),%xmm0
  divsd  %xmm2,%xmm0
  movsd  %xmm0,(%r15)
L15:
  movsd  %xmm1,0x40(%rsp)
  jmp    L11
  nopw   0x0(%rax,%rax,1)
L174:
  test   %esi,%esi
  jle    L15
L295:
  movsd  D15(%rip),%xmm2
  lea    0x8(%rbp),%r10
  lea    0x0(%r13,%r13,2),%eax
  addsd  %xmm1,%xmm2
  cmp    %r10,%rdi
  je     L296
  cmpb   $0x0,0x67(%rsp)
  je     L296
  movapd %xmm2,%xmm4
  xor    %r10d,%r10d
  unpcklpd %xmm4,%xmm4
  nopl   0x0(%rax,%rax,1)
L24:
  movupd (%r12,%r10,1),%xmm0
  movupd 0x0(%rbp,%r10,1),%xmm7
  subpd  %xmm7,%xmm0
  divpd  %xmm4,%xmm0
  movups %xmm0,(%rdi,%r10,1)
  add    $0x10,%r10
  cmp    %r14,%r10
  jne    L24
  cmp    %r8d,%esi
  je     L15
  add    %r8d,%eax
  movsd  (%r9),%xmm0
  cltq
  subsd  (%rbx,%rax,8),%xmm0
  divsd  %xmm2,%xmm0
  movsd  %xmm0,(%r15)
  jmp    L15
  nopl   (%rax)
L70:
  mov    0x58(%rsp),%rax
  movapd %xmm3,%xmm1
  subsd  (%rax,%r13,8),%xmm1
  test   %r13,%r13
  jne    L282
  jmp    L15
L3:
  movapd %xmm3,%xmm1
  xor    %r11d,%r11d
  jmp    L210
L142:
  movsd  D15(%rip),%xmm2
  mov    0x74(%rsp),%r11d
  xor    %eax,%eax
  addsd  %xmm1,%xmm2
  nopl   0x0(%rax)
L154:
  movsd  (%r12,%rax,8),%xmm0
  subsd  (%rbx,%rax,8),%xmm0
  mov    %rax,%r10
  divsd  %xmm2,%xmm0
  movsd  %xmm0,(%rdi,%rax,8)
  add    $0x1,%rax
  cmp    %r10,%r11
  jne    L154
  jmp    L15
L296:
  mov    0x74(%rsp),%r10d
  cltq
  mov    %r14,0x8(%rsp)
  lea    (%rbx,%rax,8),%r11
  xor    %eax,%eax
  mov    %r10,%r14
  nopw   0x0(%rax,%rax,1)
L130:
  movsd  (%r12,%rax,8),%xmm0
  subsd  (%r11,%rax,8),%xmm0
  mov    %rax,%r10
  divsd  %xmm2,%xmm0
  movsd  %xmm0,(%rdi,%rax,8)
  add    $0x1,%rax
  cmp    %r10,%r14
  jne    L130
  mov    0x8(%rsp),%r14
  jmp    L15
L107:
  movapd %xmm1,%xmm0
  mov    %rcx,0x30(%rsp)
  mov    %esi,0x28(%rsp)
  mov    %r8d,0x20(%rsp)
  mov    %r9,0x18(%rsp)
  mov    %rdi,0x10(%rsp)
  mov    %edx,0x8(%rsp)
  call   <sqrt@plt>
  test   %r13,%r13
  mov    0x8(%rsp),%edx
  mov    0x28(%rsp),%esi
  mov    0x58(%rsp),%rax
  movapd %xmm0,%xmm1
  mov    0x10(%rsp),%rdi
  pxor   %xmm3,%xmm3
  mov    0x18(%rsp),%r9
  mov    0x20(%rsp),%r8d
  subsd  (%rax,%r13,8),%xmm1
  mov    0x30(%rsp),%rcx
  je     L21
  movapd %xmm1,%xmm0
  mov    %rcx,0x38(%rsp)
  mov    %esi,0x60(%rsp)
  mov    %r8d,0x30(%rsp)
  mov    %r9,0x28(%rsp)
  mov    %rdi,0x20(%rsp)
  mov    %edx,0x18(%rsp)
  movsd  %xmm1,0x10(%rsp)
  call   <fabs@plt>
  movsd  0x40(%rsp),%xmm0
  mov    %rax,0x8(%rsp)
  call   <fabs@plt>
  cmp    0x8(%rsp),%rax
  movsd  0x10(%rsp),%xmm1
  pxor   %xmm3,%xmm3
  mov    0x18(%rsp),%edx
  mov    0x20(%rsp),%rdi
  mov    0x28(%rsp),%r9
  mov    0x30(%rsp),%r8d
  mov    0x60(%rsp),%esi
  mov    0x38(%rsp),%rcx
  jle    L11
  jmp    L295
```

tool_calls: name: retrieve_data
arguments: data_tag: D15
type_guessed: byte[8]

tool: {"D15": [0x11,0xea,0x2d,0x81,0x99,0x97,0x71,0x3d]}

assistant: ```c
void scriptedSphereMotionSignedDistance(const int nPoints, const double t, const int nSpace, const int nSpheres, const double *radii, const double *centers, const double *x, double *phi, double *n) {
  int i, j, k;
  double minDistance, distance, tmp;
  for (k = 0; k < nPoints; k++) {
    for (j = 0; j < nSpheres; j++) {
      tmp = 0.0;
      for (i = 0; i < nSpace; i++) {
        tmp += (x[3 * k + i] - centers[3 * j + i]) * (x[3 * k + i] - centers[3 * j + i]);
      }
      distance = sqrt(tmp) - radii[j];
      if (j == 0) {
        minDistance = distance;
        for (i = 0; i < nSpace; i++) {
          n[k * nSpace + i] = (x[3 * k + i] - centers[3 * j + i]) / (distance + 1.0e-12);
        }
      } else if (fabs(distance) < fabs(minDistance)) {
        minDistance = distance;
        for (i = 0; i < nSpace; i++) {
          n[k * nSpace + i] = (x[3 * k + i] - centers[3 * j + i]) / (distance + 1.0e-12);
        }
      }
    }
    phi[k] = minDistance;
  }
}
```