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

```asm
<Thomas_Fermi>:
  endbr64
  push   %r15
  xor    %eax,%eax
  push   %r14
  movq   %xmm0,%r14
  push   %r12
  push   %rbp
  push   %rbx
  mov    %rdi,%rbx
  sub    $0x10,%rsp
  call   <N_LogGrid@plt>
  mov    %eax,%ebp
  xor    %eax,%eax
  call   <r_LogGrid@plt>
  test   %ebp,%ebp
  jle    L18
  movslq %eax,%r12
  lea    -0x1(%rbp),%eax
  movq   %r14,%xmm7
  divsd  D22(%rip),%xmm7
  lea    0x8(%r12,%rax,8),%rbp
  movq   %xmm7,%r15
  sub    %r12,%rbx
L36:
  mov    D68(%rip),%rax
  movsd  (%r12),%xmm3
  movq   %r15,%xmm0
  movq   %rax,%xmm1
  movsd  %xmm3,(%rsp)
  call   <pow@plt>
  movsd  (%rsp),%xmm3
  pxor   %xmm5,%xmm5
  movapd %xmm0,%xmm1
  mulsd  %xmm3,%xmm1
  ucomisd %xmm1,%xmm5
  ja     L33
  movapd %xmm1,%xmm4
  sqrtsd %xmm4,%xmm4
L32:
  movsd  D91(%rip),%xmm0
  movsd  D6(%rip),%xmm2
  movsd  D113(%rip),%xmm6
  mov    D102(%rip),%rax
  mulsd  %xmm1,%xmm0
  subsd  %xmm0,%xmm2
  movsd  D54(%rip),%xmm0
  mulsd  %xmm1,%xmm2
  subsd  %xmm2,%xmm0
  movsd  D113(%rip),%xmm2
  mulsd  %xmm4,%xmm0
  addsd  %xmm0,%xmm2
  movsd  D43(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  D37(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  D55(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movq   %r14,%xmm1
  addsd  %xmm2,%xmm0
  divsd  %xmm0,%xmm1
  movq   %rax,%xmm0
  comisd %xmm1,%xmm6
  ja     L3
  xorpd  D26(%rip),%xmm1
  movapd %xmm1,%xmm0
L3:
  divsd  %xmm3,%xmm0
  movsd  %xmm0,(%rbx,%r12,1)
  add    $0x8,%r12
  cmp    %r12,%rbp
  jne    L36
L18:
  add    $0x10,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r14
  pop    %r15
  ret
L33:
  movapd %xmm1,%xmm0
  movsd  %xmm3,0x8(%rsp)
  movsd  %xmm1,(%rsp)
  call   <sqrt@plt>
  movsd  0x8(%rsp),%xmm3
  movsd  (%rsp),%xmm1
  movapd %xmm0,%xmm4
  jmp    L32
```

tool_calls: name: retrieve_data
arguments: label: D22
type_infered: f64
name: retrieve_data
arguments: label: D68
type_infered: byte[8]
name: retrieve_data
arguments: label: D91
type_infered: f64
name: retrieve_data
arguments: label: D6
type_infered: f64
name: retrieve_data
arguments: label: D113
type_infered: f64
name: retrieve_data
arguments: label: D102
type_infered: byte[8]
name: retrieve_data
arguments: label: D54
type_infered: f64
name: retrieve_data
arguments: label: D43
type_infered: f64
name: retrieve_data
arguments: label: D37
type_infered: f64
name: retrieve_data
arguments: label: D55
type_infered: f64
name: retrieve_data
arguments: label: D26
type_infered: byte[16]

tool: D22: 0.69395656

tool: {"D68": [0x55,0x55,0x55,0x55,0x55,0x55,0xd5,0x3f]}

tool: D91: 0.007298

tool: D6: 0.1486

tool: D113: 1

tool: {"D102": [0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0xbf]}

tool: D54: 0.02747

tool: D43: 0.006944

tool: D37: 0.2302

tool: D55: 1.243

tool: {"D26": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00]}

assistant: ```c
void Thomas_Fermi(Z, Vtmp) double Z;
double Vtmp[];
{
  int i, Ngrid;
  double *rgrid, x, t;
  Ngrid = N_LogGrid();
  rgrid = r_LogGrid();
  for (i = 0; i < Ngrid; ++i) {
    x = rgrid[i] * pow((Z / 0.69395656), (1.0 / 3.0));
    t = Z / (1.0 + sqrt(x) * (0.02747 - x * (0.1486 - 0.007298 * x)) + x * (1.243 + x * (0.2302 + 0.006944 * x)));
    if (t < 1.0)
      t = 1.0;
    Vtmp[i] = -t / rgrid[i];
  }
}
```