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

```asm
<Accel>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x48,%rsp
  mov    %rdi,0x18(%rsp)
  mov    %rsi,0x8(%rsp)
  mov    %rdx,%r14
  mov    %rcx,0x10(%rsp)
  mov    %r8d,0x24(%rsp)
  mov    %r9d,%r12d
  mov    0x80(%rsp),%r13
  mov    %r8d,%ebp
  imul   $0x2710,%r9d,%eax
  add    %r8d,%eax
  cltq
  shl    $0x3,%rax
  lea    (%rsi,%rax,1),%rbx
  mov    %rbx,0x28(%rsp)
  lea    (%rdx,%rax,1),%rsi
  mov    %rsi,0x30(%rsp)
  add    %rcx,%rax
  mov    %rax,0x38(%rsp)
  lea    0x186a0(%r8),%r15d
  mov    $0x0,%ebx
  pxor   %xmm6,%xmm6
  movsd  %xmm6,(%rsp)
  jmp    L86
L48:
  mov    0x20(%rsp),%r9d
  mov    %r12d,%r8d
  mov    0x24(%rsp),%ecx
  mov    0x10(%rsp),%rdx
  mov    %r14,%rsi
  mov    0x8(%rsp),%rdi
  call   <Distancia@plt>
  mov    %eax,%edi
  mov    $0x3,%esi
  call   <pow@plt>
  movapd %xmm0,%xmm2
  movslq %ebp,%rax
  mov    0x8(%rsp),%rdx
  movsd  (%rdx,%rax,8),%xmm0
  mov    0x28(%rsp),%rax
  subsd  (%rax),%xmm0
  movsd  D38(%rip),%xmm1
  mov    0x18(%rsp),%rax
  mulsd  (%rax,%rbx,8),%xmm1
  mulsd  %xmm1,%xmm0
  divsd  %xmm2,%xmm0
  addsd  (%rsp),%xmm0
  movsd  %xmm0,(%rsp)
  jmp    L79
L22:
  mov    0x20(%rsp),%r9d
  mov    %r12d,%r8d
  mov    0x24(%rsp),%ecx
  mov    0x10(%rsp),%rdx
  mov    %r14,%rsi
  mov    0x8(%rsp),%rdi
  call   <Distancia@plt>
  mov    %eax,%edi
  mov    $0x3,%esi
  call   <pow@plt>
  movapd %xmm0,%xmm2
  movslq %ebp,%rax
  movsd  (%r14,%rax,8),%xmm0
  mov    0x30(%rsp),%rax
  subsd  (%rax),%xmm0
  movsd  D38(%rip),%xmm1
  mov    0x18(%rsp),%rax
  mulsd  (%rax,%rbx,8),%xmm1
  mulsd  %xmm1,%xmm0
  divsd  %xmm2,%xmm0
  addsd  (%rsp),%xmm0
  movsd  %xmm0,(%rsp)
  jmp    L28
L75:
  add    $0x1,%rbx
  add    $0x2710,%ebp
  cmp    %r15d,%ebp
  je     L56
L86:
  mov    %ebx,0x20(%rsp)
  cmp    %ebx,%r12d
  je     L75
  lea    D26(%rip),%rsi
  mov    %r13,%rdi
  call   <strcmp@plt>
  test   %rax,%rax
  je     L48
L79:
  lea    D39(%rip),%rsi
  mov    %r13,%rdi
  call   <strcmp@plt>
  test   %rax,%rax
  je     L22
L28:
  lea    D11(%rip),%rsi
  mov    %r13,%rdi
  call   <strcmp@plt>
  test   %rax,%rax
  jne    L75
  mov    0x20(%rsp),%r9d
  mov    %r12d,%r8d
  mov    0x24(%rsp),%ecx
  mov    0x10(%rsp),%rdx
  mov    %r14,%rsi
  mov    0x8(%rsp),%rdi
  call   <Distancia@plt>
  mov    %eax,%edi
  mov    $0x3,%esi
  call   <pow@plt>
  movapd %xmm0,%xmm1
  movslq %ebp,%rax
  mov    0x10(%rsp),%rcx
  movsd  (%rcx,%rax,8),%xmm0
  mov    0x38(%rsp),%rax
  subsd  (%rax),%xmm0
  mov    0x18(%rsp),%rax
  movsd  (%rax,%rbx,8),%xmm2
  mulsd  D38(%rip),%xmm2
  mulsd  %xmm2,%xmm0
  divsd  %xmm1,%xmm0
  addsd  (%rsp),%xmm0
  movsd  %xmm0,(%rsp)
  jmp    L75
L56:
  movsd  (%rsp),%xmm0
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: read_data
arguments: data_label: D38
datatype: f64
name: read_data
arguments: data_label: D26
datatype: string
name: read_data
arguments: data_label: D39
datatype: string
name: read_data
arguments: data_label: D11
datatype: string

tool: D38: 39.47841760435743

tool: D26: x

tool: D39: y

tool: D11: z

assistant: ```c
double Accel(double *masas, double *posX, double *posY, double *posZ, int t, int p, char comp[]) {
  double a = 0;
  double r3 = 0;
  for (int i = 0; i < 10; i++) {
    if (i != p) {
      if (strcmp(comp, "x") == 0) {
        r3 = pow(Distancia(posX, posY, posZ, t, p, i), 3);
        a += 39.4784176043574344753379 * masas[i] * (posX[i * 10000 + t] - posX[p * 10000 + t]) / r3;
      }
      if (strcmp(comp, "y") == 0) {
        r3 = pow(Distancia(posX, posY, posZ, t, p, i), 3);
        a += 39.4784176043574344753379 * masas[i] * (posY[i * 10000 + t] - posY[p * 10000 + t]) / r3;
      }
      if (strcmp(comp, "z") == 0) {
        r3 = pow(Distancia(posX, posY, posZ, t, p, i), 3);
        a += 39.4784176043574344753379 * masas[i] * (posZ[i * 10000 + t] - posZ[p * 10000 + t]) / r3;
      }
    }
  }
  return a;
}
```