user: What is the c source code of the assembly code below?

```asm
<scriptedSphereMotionSignedDistance>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x58,%rsp
  mov    %edx,0x40(%rsp)
  mov    %rcx,0x28(%rsp)
  mov    %r9,0x48(%rsp)
  test   %edi,%edi
  jle    L137
  mov    %esi,%ebp
  mov    0x90(%rsp),%rsi
  lea    (%rdi,%rdi,2),%eax
  mov    %r8,%r15
  mov    %eax,0x44(%rsp)
  pxor   %xmm2,%xmm2
  mov    %rsi,0x30(%rsp)
  movl   $0x0,0x3c(%rsp)
  movl   $0x0,0x38(%rsp)
  nopl   (%rax)
L111:
  movslq 0x40(%rsp),%rax
  test   %eax,%eax
  jle    L168
  mov    0x48(%rsp),%rdi
  mov    %rax,0x20(%rsp)
  xor    %r13d,%r13d
  mov    %r15,%rbx
  movslq 0x38(%rsp),%rax
  lea    (%rdi,%rax,8),%r12
  movslq 0x3c(%rsp),%rax
  mov    0x98(%rsp),%rdi
  mov    %r12,%r15
  mov    %r13,%r12
  lea    (%rdi,%rax,8),%r14
  lea    -0x1(%rbp),%eax
  mov    %eax,0x14(%rsp)
  mov    %r14,%r13
  nop
L170:
  test   %ebp,%ebp
  jle    L53
  lea    (%r12,%r12,2),%eax
  mov    0x14(%rsp),%r9d
  movapd %xmm2,%xmm1
  cltq
  lea    (%rbx,%rax,8),%r8
  xor    %eax,%eax
  nopl   (%rax)
L96:
  movsd  (%r15,%rax,8),%xmm0
  subsd  (%r8,%rax,8),%xmm0
  mov    %rax,%rdx
  add    $0x1,%rax
  mulsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm1
  cmp    %rdx,%r9
  jne    L96
  ucomisd %xmm1,%xmm2
  ja     L142
  mov    0x28(%rsp),%rax
  sqrtsd %xmm1,%xmm1
  subsd  (%rax,%r12,8),%xmm1
  test   %r12,%r12
  je     L185
L126:
  movapd %xmm1,%xmm0
  movsd  %xmm1,0x18(%rsp)
  call   <fabs@plt>
  movsd  0x8(%rsp),%xmm0
  mov    %rax,%r14
  call   <fabs@plt>
  movsd  0x18(%rsp),%xmm1
  pxor   %xmm2,%xmm2
  cmp    %rax,%r14
  jl     L177
L136:
  add    $0x1,%r12
  cmp    %r12,0x20(%rsp)
  jne    L170
  mov    %rbx,%r15
L168:
  mov    0x30(%rsp),%rax
  movsd  0x8(%rsp),%xmm5
  addl   $0x3,0x38(%rsp)
  movsd  %xmm5,(%rax)
  add    $0x8,%rax
  mov    %rax,0x30(%rsp)
  mov    0x38(%rsp),%eax
  add    %ebp,0x3c(%rsp)
  cmp    0x44(%rsp),%eax
  jne    L111
L137:
  add    $0x58,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopw   0x0(%rax,%rax,1)
L185:
  movsd  %xmm1,0x8(%rsp)
  movsd  D11(%rip),%xmm1
  addsd  0x8(%rsp),%xmm1
  xor    %eax,%eax
  mov    0x14(%rsp),%r8d
  nopl   0x0(%rax,%rax,1)
L94:
  movsd  (%r15,%rax,8),%xmm0
  subsd  (%rbx,%rax,8),%xmm0
  mov    %rax,%rdx
  divsd  %xmm1,%xmm0
  movsd  %xmm0,0x0(%r13,%rax,8)
  add    $0x1,%rax
  cmp    %rdx,%r8
  jne    L94
  jmp    L136
  nop
L177:
  movsd  %xmm1,0x8(%rsp)
  test   %ebp,%ebp
  jle    L136
L169:
  movsd  D11(%rip),%xmm1
  lea    (%r12,%r12,2),%eax
  addsd  0x8(%rsp),%xmm1
  cltq
  mov    0x14(%rsp),%r9d
  lea    (%rbx,%rax,8),%r8
  xor    %eax,%eax
  nopl   (%rax)
L48:
  movsd  (%r15,%rax,8),%xmm0
  subsd  (%r8,%rax,8),%xmm0
  mov    %rax,%rdx
  divsd  %xmm1,%xmm0
  movsd  %xmm0,0x0(%r13,%rax,8)
  add    $0x1,%rax
  cmp    %rdx,%r9
  jne    L48
  jmp    L136
  nopl   0x0(%rax,%rax,1)
L53:
  mov    0x28(%rsp),%rax
  movapd %xmm2,%xmm1
  subsd  (%rax,%r12,8),%xmm1
  test   %r12,%r12
  jne    L126
  movsd  %xmm1,0x8(%rsp)
  jmp    L136
L142:
  movapd %xmm1,%xmm0
  call   <sqrt@plt>
  mov    0x28(%rsp),%rax
  test   %r12,%r12
  pxor   %xmm2,%xmm2
  movapd %xmm0,%xmm1
  subsd  (%rax,%r12,8),%xmm1
  je     L185
  movapd %xmm1,%xmm0
  movsd  %xmm1,0x18(%rsp)
  call   <fabs@plt>
  movsd  0x8(%rsp),%xmm0
  mov    %rax,%r14
  call   <fabs@plt>
  movsd  0x18(%rsp),%xmm1
  pxor   %xmm2,%xmm2
  cmp    %r14,%rax
  jle    L136
  movsd  %xmm1,0x8(%rsp)
  jmp    L169
```

tool_calls: name: read_data
arguments: data_label: D11
type_infered: byte[8]

tool: {"D11": [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;
  }
}
```