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

```asm
<paw_scattering_test>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r12
  push   %rbx
  sub    $0x310,%rsp
  movsd  %xmm0,-0x308(%rbp)
  movsd  %xmm1,-0x310(%rbp)
  mov    %edi,-0x314(%rbp)
  mov    %esi,-0x318(%rbp)
  movsd  %xmm2,-0x320(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x18(%rbp)
  xor    %eax,%eax
  mov    $0x0,%eax
  call   <paw_N_LogGrid@plt>
  mov    %eax,-0x2f0(%rbp)
  mov    $0x0,%eax
  call   <paw_r_LogGrid@plt>
  mov    %rax,-0x2e8(%rbp)
  mov    $0x0,%eax
  call   <paw_log_amesh_LogGrid@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x2e0(%rbp)
  mov    $0x0,%eax
  call   <paw_alloc_LogGrid@plt>
  mov    %rax,-0x2d8(%rbp)
  mov    $0x0,%eax
  call   <paw_alloc_LogGrid@plt>
  mov    %rax,-0x2d0(%rbp)
  mov    $0x0,%eax
  call   <paw_alloc_LogGrid@plt>
  mov    %rax,-0x2c8(%rbp)
  mov    $0x0,%eax
  call   <paw_alloc_LogGrid@plt>
  mov    %rax,-0x2c0(%rbp)
  mov    -0x314(%rbp),%eax
  mov    %eax,%edi
  call   <paw_alloc_1d_array@plt>
  mov    %rax,-0x2b8(%rbp)
  mov    -0x314(%rbp),%eax
  mov    %eax,%edi
  call   <paw_alloc_1d_array@plt>
  mov    %rax,-0x2b0(%rbp)
  mov    -0x314(%rbp),%eax
  mov    %eax,%edi
  call   <paw_alloc_1d_array@plt>
  mov    %rax,-0x2a8(%rbp)
  movsd  -0x310(%rbp),%xmm0
  subsd  -0x308(%rbp),%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x314(%rbp),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x2a0(%rbp)
  mov    $0x0,%eax
  call   <paw_get_kohn_sham_potential@plt>
  mov    %rax,-0x298(%rbp)
  mov    -0x320(%rbp),%rax
  movq   %rax,%xmm0
  call   <paw_get_grid_index@plt>
  mov    %eax,-0x2ec(%rbp)
  movl   $0x0,-0x2f8(%rbp)
  jmp    L27
L23:
  pxor   %xmm0,%xmm0
  cvtsi2sdl -0x2f8(%rbp),%xmm0
  mulsd  -0x2a0(%rbp),%xmm0
  movsd  -0x308(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x288(%rbp)
  mov    -0x2f8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x2a8(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x288(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x2c0(%rbp),%rdx
  mov    -0x2c8(%rbp),%rsi
  movsd  -0x288(%rbp),%xmm0
  mov    -0x320(%rbp),%rcx
  mov    -0x318(%rbp),%eax
  movapd %xmm0,%xmm1
  movq   %rcx,%xmm0
  mov    %eax,%edi
  call   <paw_solve_paw_scattering@plt>
  mov    -0x2d0(%rbp),%r8
  mov    -0x2d8(%rbp),%rcx
  mov    -0x288(%rbp),%rdi
  mov    -0x2ec(%rbp),%edx
  mov    -0x298(%rbp),%rsi
  mov    -0x318(%rbp),%eax
  movq   %rdi,%xmm0
  mov    %eax,%edi
  call   <paw_R_Schrodinger_Fixed_E@plt>
  mov    -0x2ec(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x2d0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x2ec(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x2d8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  mov    -0x2ec(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x2e8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mulsd  %xmm2,%xmm1
  mulsd  -0x2e0(%rbp),%xmm1
  mov    -0x2f8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x2b8(%rbp),%rax
  add    %rdx,%rax
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x2ec(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x2c0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x2ec(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x2c8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  mov    -0x2ec(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x2e8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mulsd  %xmm2,%xmm1
  mulsd  -0x2e0(%rbp),%xmm1
  mov    -0x2f8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x2b0(%rbp),%rax
  add    %rdx,%rax
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x2f8(%rbp)
L27:
  mov    -0x314(%rbp),%eax
  cmp    -0x2f8(%rbp),%eax
  jg     L23
  mov    $0x0,%eax
  call   <paw_debug@plt>
  test   %rax,%rax
  je     L54
  mov    -0x318(%rbp),%eax
  mov    %eax,%edi
  call   <paw_spd_Name@plt>
  mov    %rax,%r12
  mov    $0x0,%eax
  call   <paw_get_atom_name@plt>
  mov    %rax,%rbx
  mov    $0x0,%eax
  call   <paw_sdir@plt>
  mov    %rax,%rdx
  lea    -0x280(%rbp),%rax
  mov    %r12,%r8
  mov    %rbx,%rcx
  lea    D48(%rip),%rsi
  mov    %rax,%rdi
  call   <sprintf@plt>
  lea    -0x280(%rbp),%rax
  lea    D131(%rip),%rdx
  mov    %rdx,%rsi
  mov    %rax,%rdi
  call   <fopen@plt>
  mov    %rax,-0x290(%rbp)
  movl   $0x0,-0x2f4(%rbp)
  jmp    L3
L55:
  mov    -0x2f4(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x2b0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x2f4(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x2b8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x2f4(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x2a8(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x290(%rbp),%rax
  movapd %xmm1,%xmm2
  movapd %xmm0,%xmm1
  movq   %rdx,%xmm0
  lea    D89(%rip),%rdx
  mov    %rdx,%rsi
  mov    %rax,%rdi
  mov    $0x3,%eax
  call   <fprintf@plt>
  addl   $0x1,-0x2f4(%rbp)
L3:
  mov    -0x314(%rbp),%eax
  cmp    -0x2f4(%rbp),%eax
  jg     L55
  mov    -0x290(%rbp),%rax
  mov    %rax,%rdi
  call   <fclose@plt>
  mov    -0x318(%rbp),%eax
  mov    %eax,%edi
  call   <paw_spd_Name@plt>
  mov    %rax,%r12
  mov    $0x0,%eax
  call   <paw_get_atom_name@plt>
  mov    %rax,%rbx
  mov    $0x0,%eax
  call   <paw_sdir@plt>
  mov    %rax,%rdx
  lea    -0x150(%rbp),%rax
  mov    %r12,%r8
  mov    %rbx,%rcx
  lea    D77(%rip),%rsi
  mov    %rax,%rdi
  call   <sprintf@plt>
  lea    -0x150(%rbp),%rax
  mov    %rax,%rsi
  lea    D85(%rip),%rax
  mov    %rax,%rdi
  call   <printf@plt>
  lea    -0x150(%rbp),%rax
  lea    D131(%rip),%rdx
  mov    %rdx,%rsi
  mov    %rax,%rdi
  call   <fopen@plt>
  mov    %rax,-0x290(%rbp)
  mov    -0x290(%rbp),%rax
  mov    %rax,%rcx
  mov    $0x16,%edx
  mov    $0x1,%esi
  lea    D68(%rip),%rax
  mov    %rax,%rdi
  call   <fwrite@plt>
  mov    -0x290(%rbp),%rax
  mov    %rax,%rcx
  mov    $0xd,%edx
  mov    $0x1,%esi
  lea    D84(%rip),%rax
  mov    %rax,%rdi
  call   <fwrite@plt>
  mov    -0x290(%rbp),%rax
  mov    %rax,%rcx
  mov    $0xf,%edx
  mov    $0x1,%esi
  lea    D64(%rip),%rax
  mov    %rax,%rdi
  call   <fwrite@plt>
  movsd  -0x310(%rbp),%xmm0
  mov    -0x308(%rbp),%rdx
  mov    -0x290(%rbp),%rax
  movapd %xmm0,%xmm1
  movq   %rdx,%xmm0
  lea    D166(%rip),%rdx
  mov    %rdx,%rsi
  mov    %rax,%rdi
  mov    $0x2,%eax
  call   <fprintf@plt>
  mov    -0x290(%rbp),%rax
  mov    %rax,%rcx
  mov    $0xa,%edx
  mov    $0x1,%esi
  lea    D47(%rip),%rax
  mov    %rax,%rdi
  call   <fwrite@plt>
  mov    -0x290(%rbp),%rax
  mov    %rax,%rcx
  mov    $0xd,%edx
  mov    $0x1,%esi
  lea    D84(%rip),%rax
  mov    %rax,%rdi
  call   <fwrite@plt>
  mov    -0x290(%rbp),%rax
  mov    %rax,%rcx
  mov    $0x1a,%edx
  mov    $0x1,%esi
  lea    D90(%rip),%rax
  mov    %rax,%rdi
  call   <fwrite@plt>
  mov    -0x320(%rbp),%rdx
  mov    -0x290(%rbp),%rax
  movq   %rdx,%xmm0
  lea    D81(%rip),%rdx
  mov    %rdx,%rsi
  mov    %rax,%rdi
  mov    $0x1,%eax
  call   <fprintf@plt>
  mov    -0x318(%rbp),%eax
  mov    %eax,%edi
  call   <paw_spd_Name@plt>
  mov    %rax,%rbx
  mov    $0x0,%eax
  call   <paw_get_atom_name@plt>
  mov    %rax,%rdx
  mov    -0x290(%rbp),%rax
  mov    %rbx,%rcx
  lea    D14(%rip),%rsi
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <fprintf@plt>
  lea    -0x280(%rbp),%rdx
  mov    -0x290(%rbp),%rax
  lea    D66(%rip),%rcx
  mov    %rcx,%rsi
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <fprintf@plt>
  mov    -0x290(%rbp),%rax
  mov    %rax,%rcx
  mov    $0x1a,%edx
  mov    $0x1,%esi
  lea    D140(%rip),%rax
  mov    %rax,%rdi
  call   <fwrite@plt>
  mov    -0x290(%rbp),%rax
  mov    %rax,%rsi
  mov    $0xa,%edi
  call   <fputc@plt>
  mov    -0x290(%rbp),%rax
  mov    %rax,%rcx
  mov    $0x9,%edx
  mov    $0x1,%esi
  lea    D53(%rip),%rax
  mov    %rax,%rdi
  call   <fwrite@plt>
  mov    -0x290(%rbp),%rax
  mov    %rax,%rdi
  call   <fclose@plt>
L54:
  nop
  mov    -0x18(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L26
  call   <__stack_chk_fail@plt>
L26:
  add    $0x310,%rsp
  pop    %rbx
  pop    %r12
  pop    %rbp
  ret
```

tool_calls: name: extract_symbolic_data
arguments: label: D48
type_infered: string
name: extract_symbolic_data
arguments: label: D131
type_infered: string
name: extract_symbolic_data
arguments: label: D89
type_infered: string
name: extract_symbolic_data
arguments: label: D77
type_infered: string
name: extract_symbolic_data
arguments: label: D85
type_infered: string
name: extract_symbolic_data
arguments: label: D68
type_infered: string
name: extract_symbolic_data
arguments: label: D84
type_infered: string
name: extract_symbolic_data
arguments: label: D64
type_infered: string
name: extract_symbolic_data
arguments: label: D166
type_infered: string
name: extract_symbolic_data
arguments: label: D47
type_infered: string
name: extract_symbolic_data
arguments: label: D90
type_infered: string
name: extract_symbolic_data
arguments: label: D81
type_infered: string
name: extract_symbolic_data
arguments: label: D14
type_infered: string
name: extract_symbolic_data
arguments: label: D66
type_infered: string
name: extract_symbolic_data
arguments: label: D140
type_infered: string
name: extract_symbolic_data
arguments: label: D53
type_infered: string

tool: D48: %s%s_%s_scat_test.dat

tool: D131: w+

tool: D89: %le	%le	%le


tool: D77: %s%s_%s_scat_test.plt

tool: D85: script_filename: %s


tool: D68: set style data lines 


tool: D84: set nolabel 


tool: D64: set autoscale 


tool: D166: set xr[%f:%f] 


tool: D47: set grid 


tool: D90: set xlabel "e (Hartree)" 


tool: D81: set ylabel "logarithmic derivative at r=%f" 


tool: D14: set title " %s %s channel scattering test


tool: D66: plot "%s" using 1:2 title "all electron",

tool: D140: "" using 1:3 title "paw" 


tool: D53: pause -1


assistant: ```c
void paw_scattering_test(double e1, double e2, int number_points, int l, double r) {
  int i;
  int k;
  int i_end_point;
  int Ngrid;
  double *rgrid;
  FILE *fp;
  double *V_ks;
  double *psi1;
  double *psi1_prime;
  double *psi;
  double *psi_prime;
  double *log_grid_ae;
  double *log_grid_paw;
  double de;
  double *e3;
  double e_test;
  double log_amesh;
  char data_filename[300];
  char script_filename[300];
  Ngrid = paw_N_LogGrid();
  rgrid = paw_r_LogGrid();
  log_amesh = paw_log_amesh_LogGrid();
  psi1 = paw_alloc_LogGrid();
  psi1_prime = paw_alloc_LogGrid();
  psi = paw_alloc_LogGrid();
  psi_prime = paw_alloc_LogGrid();
  log_grid_ae = paw_alloc_1d_array(number_points);
  log_grid_paw = paw_alloc_1d_array(number_points);
  e3 = paw_alloc_1d_array(number_points);
  de = (e2 - e1) / number_points;
  V_ks = paw_get_kohn_sham_potential();
  i_end_point = paw_get_grid_index(r);
  for (i = 0; i <= number_points - 1; i++) {
    e_test = e1 + de * i;
    e3[i] = e_test;
    paw_solve_paw_scattering(l, r, e_test, psi, psi_prime);
    paw_R_Schrodinger_Fixed_E(l, V_ks, i_end_point, e_test, psi1, psi1_prime);
    log_grid_ae[i] = psi1_prime[i_end_point - 1] / (psi1[i_end_point - 1] * rgrid[i_end_point - 1] * log_amesh);
    log_grid_paw[i] = psi_prime[i_end_point - 1] / (psi[i_end_point - 1] * rgrid[i_end_point - 1] * log_amesh);
  }
  if (paw_debug()) {
    sprintf(data_filename, "%s%s_%s_scat_test.dat", paw_sdir(), paw_get_atom_name(), paw_spd_Name(l));
    fp = fopen(data_filename, "w+");
    for (k = 0; k <= number_points - 1; k++) {
      fprintf(fp, "%le\t%le\t%le\n", e3[k], log_grid_ae[k], log_grid_paw[k]);
    }
    fclose(fp);
    sprintf(script_filename, "%s%s_%s_scat_test.plt", paw_sdir(), paw_get_atom_name(), paw_spd_Name(l));
    printf("script_filename: %s\n", script_filename);
    fp = fopen(script_filename, "w+");
    fprintf(fp, "set style data lines \n");
    fprintf(fp, "set nolabel \n");
    fprintf(fp, "set autoscale \n");
    fprintf(fp, "set xr[%f:%f] \n", e1, e2);
    fprintf(fp, "set grid \n");
    fprintf(fp, "set nolabel \n");
    fprintf(fp, "set xlabel \"e (Hartree)\" \n");
    fprintf(fp, "set ylabel \"logarithmic derivative at r=%f\" \n", r);
    fprintf(fp, "set title \" %s %s channel scattering test\n", paw_get_atom_name(), paw_spd_Name(l));
    fprintf(fp, "plot \"%s\" using 1:2 title \"all electron\",", data_filename);
    fprintf(fp, "\"\" using 1:3 title \"paw\" \n");
    fprintf(fp, "\n");
    fprintf(fp, "pause -1\n");
    fclose(fp);
  }
}
```