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

```asm
<R_Schrodinger_Fixed_E>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x58,%rsp
  mov    %rdx,%rbp
  mov    %ecx,0x3c(%rsp)
  movsd  %xmm0,0x30(%rsp)
  mov    %r8,%r12
  mov    %r9,%rbx
  lea    0x1(%rsi),%eax
  pxor   %xmm7,%xmm7
  cvtsi2sd %eax,%xmm7
  movsd  %xmm7,0x20(%rsp)
  imul   %esi,%eax
  pxor   %xmm7,%xmm7
  cvtsi2sd %eax,%xmm7
  movsd  %xmm7,0x10(%rsp)
  mov    $0x0,%eax
  call   <N_LogGrid@plt>
  mov    %eax,%r14d
  mov    %eax,0x48(%rsp)
  mov    $0x0,%eax
  call   <log_amesh_LogGrid@plt>
  pxor   %xmm6,%xmm6
  cvtsi2sd %eax,%xmm6
  movsd  %xmm6,0x18(%rsp)
  mulsd  %xmm6,%xmm6
  movsd  %xmm6,0x28(%rsp)
  mov    $0x0,%eax
  call   <r_LogGrid@plt>
  mov    %eax,%r13d
  mov    $0x0,%eax
  call   <alloc_LogGrid@plt>
  mov    %eax,%r15d
  cltq
  mov    %rax,0x40(%rsp)
  mov    $0x0,%eax
  call   <alloc_LogGrid@plt>
  mov    %eax,0x4c(%rsp)
  cltq
  mov    %rax,0x8(%rsp)
  test   %r14d,%r14d
  jle    L94
  movslq %r13d,%rsi
  lea    -0x1(%r14),%eax
  lea    0x8(%rsi,%rax,8),%rcx
  mov    %rsi,%rax
  movslq %r15d,%rdx
  sub    %rsi,%rdx
  sub    %rsi,%rbp
L0:
  movsd  (%rax),%xmm1
  movsd  0x0(%rbp,%rax,1),%xmm0
  subsd  0x30(%rsp),%xmm0
  addsd  %xmm0,%xmm0
  mulsd  %xmm1,%xmm1
  mulsd  %xmm1,%xmm0
  addsd  0x10(%rsp),%xmm0
  mulsd  0x28(%rsp),%xmm0
  movsd  %xmm0,(%rdx,%rax,1)
  add    $0x8,%rax
  cmp    %rcx,%rax
  jne    L0
L94:
  movslq %r13d,%r13
  lea    0x20(%r13),%r14
  mov    %r13,%rbp
  neg    %r13
L116:
  movsd  0x20(%rsp),%xmm1
  movsd  0x0(%rbp),%xmm0
  call   <pow@plt>
  lea    (%r12,%r13,1),%rax
  movsd  %xmm0,(%rax,%rbp,1)
  lea    (%rbx,%r13,1),%rax
  movsd  0x20(%rsp),%xmm1
  mulsd  0x18(%rsp),%xmm1
  mulsd  %xmm0,%xmm1
  movsd  %xmm1,(%rax,%rbp,1)
  add    $0x8,%rbp
  cmp    %r14,%rbp
  jne    L116
  movslq %r15d,%rdx
  lea    0x20(%rdx),%r15
  mov    %rdx,%rax
  movslq 0x4c(%rsp),%rcx
  mov    %rcx,%r8
  sub    %rdx,%r8
  mov    %r12,%rdi
  sub    %rdx,%rdi
  mov    %rbx,%rsi
  sub    %rdx,%rsi
L102:
  movsd  (%rdi,%rax,1),%xmm0
  mulsd  (%rax),%xmm0
  movsd  0x18(%rsp),%xmm1
  mulsd  (%rsi,%rax,1),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%r8,%rax,1)
  add    $0x8,%rax
  cmp    %r15,%rax
  jne    L102
  mov    0x3c(%rsp),%eax
  cmp    $0x3,%eax
  jle    L147
  add    $0x20,%rcx
  mov    %rcx,0x10(%rsp)
  mov    %eax,%eax
  mov    %rax,0x28(%rsp)
  mov    $0x3,%ebp
  jmp    L39
L72:
  add    $0x1,%rbp
  add    $0x8,%r15
  addq   $0x8,0x10(%rsp)
  cmp    %rbp,0x28(%rsp)
  je     L147
L39:
  mov    %ebp,%r14d
  mov    %rbx,%rdx
  mov    %r12,%rsi
  mov    %ebp,%edi
  mov    $0x0,%eax
  call   <Predictor_Out@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  %xmm0,0x8(%r12,%rbp,8)
  mov    0x8(%rsp),%rdx
  mov    %rbx,%rsi
  mov    %ebp,%edi
  mov    $0x0,%eax
  call   <Predictor_Out@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  %xmm0,0x8(%rbx,%rbp,8)
  mov    $0x6,%r13d
L44:
  movsd  (%r15),%xmm0
  mulsd  0x8(%r12,%rbp,8),%xmm0
  movsd  0x18(%rsp),%xmm1
  mulsd  0x8(%rbx,%rbp,8),%xmm1
  addsd  %xmm1,%xmm0
  mov    0x10(%rsp),%rax
  movsd  %xmm0,(%rax)
  mov    0x8(%rsp),%rdx
  mov    %rbx,%rsi
  mov    %r14d,%edi
  mov    $0x0,%eax
  call   <Corrector_Out@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  %xmm0,0x8(%rbx,%rbp,8)
  mov    %rbx,%rdx
  mov    %r12,%rsi
  mov    %r14d,%edi
  mov    $0x0,%eax
  call   <Corrector_Out@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  %xmm0,0x8(%r12,%rbp,8)
  sub    $0x1,%r13d
  jne    L44
  jmp    L72
L147:
  mov    %r12,%rsi
  movsd  0x20(%rsp),%xmm0
  mov    0x3c(%rsp),%edi
  mov    $0x1,%eax
  call   <Norm_LogGrid@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm0,%xmm1
  ja     L35
  sqrtsd %xmm0,%xmm0
L62:
  movsd  D9(%rip),%xmm1
  divsd  %xmm0,%xmm1
  mov    0x3c(%rsp),%eax
  test   %eax,%eax
  js     L25
  lea    0x1(%rax),%edx
  shl    $0x3,%rdx
  mov    $0x0,%eax
L50:
  movapd %xmm1,%xmm0
  mulsd  (%r12,%rax,1),%xmm0
  movsd  %xmm0,(%r12,%rax,1)
  movapd %xmm1,%xmm0
  mulsd  (%rbx,%rax,1),%xmm0
  movsd  %xmm0,(%rbx,%rax,1)
  add    $0x8,%rax
  cmp    %rdx,%rax
  jne    L50
L25:
  mov    0x3c(%rsp),%eax
  add    $0x1,%eax
  mov    0x48(%rsp),%edx
  cmp    %eax,%edx
  jle    L133
  cltq
L19:
  movq   $0x0,(%r12,%rax,8)
  movq   $0x0,(%rbx,%rax,8)
  add    $0x1,%rax
  cmp    %eax,%edx
  jg     L19
L133:
  mov    0x8(%rsp),%rdi
  mov    $0x0,%eax
  call   <dealloc_LogGrid@plt>
  mov    0x40(%rsp),%rdi
  mov    $0x0,%eax
  call   <dealloc_LogGrid@plt>
  add    $0x58,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L35:
  call   <sqrt@plt>
  jmp    L62
```

tool_calls: name: read_data
arguments: data_label: D9
datatype: f64

tool: D9: 1

assistant: ```c
void R_Schrodinger_Fixed_E(n, l, v, match, E, u, uprime) int n, l;
double v[];
int match;
double E;
double u[], uprime[];
{
  int i, j, Ngrid;
  double log_amesh, log_amesh2, gamma, L2, r2, sum, *r, *f_upp, *upp;
  gamma = ((double)(l + 1));
  L2 = ((double)(l * (l + 1)));
  Ngrid = N_LogGrid();
  log_amesh = log_amesh_LogGrid();
  log_amesh2 = log_amesh * log_amesh;
  r = (double *)r_LogGrid();
  f_upp = (double *)alloc_LogGrid();
  upp = (double *)alloc_LogGrid();
  for (i = 0; i < Ngrid; ++i) {
    r2 = r[i];
    r2 = r2 * r2;
    f_upp[i] = log_amesh2 * (L2 + 2.0 * (v[i] - E) * r2);
  }
  for (i = 0; i < 4; ++i) {
    u[i] = pow(r[i], gamma);
    uprime[i] = log_amesh * gamma * u[i];
  }
  for (i = 0; i < 4; ++i)
    upp[i] = log_amesh * uprime[i] + f_upp[i] * u[i];
  for (i = 3; i < match; ++i) {
    u[i + 1] = Predictor_Out(i, u, uprime);
    uprime[i + 1] = Predictor_Out(i, uprime, upp);
    for (j = 0; j < 6; ++j) {
      upp[i + 1] = log_amesh * uprime[i + 1] + f_upp[i + 1] * u[i + 1];
      uprime[i + 1] = Corrector_Out(i, uprime, upp);
      u[i + 1] = Corrector_Out(i, u, uprime);
    }
  }
  sum = Norm_LogGrid(match, gamma, u);
  sum = 1.0 / sqrt(sum);
  for (i = 0; i <= match; ++i) {
    u[i] = sum * u[i];
    uprime[i] = sum * uprime[i];
  }
  for (i = match + 1; i < Ngrid; ++i) {
    u[i] = 0.0;
    uprime[i] = 0.0;
  }
  dealloc_LogGrid(upp);
  dealloc_LogGrid(f_upp);
  return;
}
```