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

```asm
<read_matrix>:
  endbr64
  push   %r15
  mov    %ecx,%r10d
  mov    %edx,%r15d
  push   %r14
  sub    %edx,%r10d
  mov    %r8d,%r14d
  push   %r13
  add    $0x1,%r10d
  push   %r12
  mov    %r9d,%r12d
  push   %rbp
  mov    %rsi,%rbp
  push   %rbx
  mov    %ecx,%ebx
  sub    $0x88,%rsp
  mov    %rdi,(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x78(%rsp)
  xor    %eax,%eax
  mov    %r9d,%eax
  sub    %r8d,%eax
  add    $0x1,%eax
  imul   %eax,%r10d
  movslq %r10d,%rdi
  mov    %r10d,0x8(%rsp)
  shl    $0x2,%rdi
  call   <malloc@plt>
  mov    (%rsp),%rdi
  lea    D4(%rip),%rsi
  mov    %rax,%r13
  call   <fopen@plt>
  mov    0x8(%rsp),%r10d
  test   %rax,%rax
  mov    %rax,%r9
  je     L3
L11:
  mov    %r9,%rcx
  mov    %r10d,%edx
  mov    $0x4,%esi
  mov    %r13,%rdi
  mov    %r9,0x8(%rsp)
  call   <fread@plt>
  cmp    %r14d,%r12d
  mov    0x8(%rsp),%r9
  jl     L39
  cmp    %r15d,%ebx
  jl     L39
  movslq %r14d,%rax
  sub    %r14d,%r12d
  movslq %r15d,%rdx
  sub    %r15d,%ebx
  lea    0x1(%rax,%r12,1),%r11
  mov    %ebx,%r15d
  xor    %r10d,%r10d
  add    $0x1,%ebx
  lea    0x0(,%rax,8),%rcx
  mov    %r13,%rax
  shl    $0x3,%r11
  lea    0x4(%r13),%r12
  neg    %rax
  add    %rax,%rax
  lea    (%rax,%rdx,8),%r14
  nopl   0x0(%rax)
L69:
  movslq %r10d,%rdx
  lea    (%r15,%rdx,1),%rsi
  lea    0x0(%r13,%rdx,4),%rax
  neg    %rdx
  lea    (%r12,%rsi,4),%rdi
  lea    (%r14,%rdx,8),%rsi
  add    %rbp,%rsi
  nopw   0x0(%rax,%rax,1)
L42:
  mov    (%rsi,%rax,2),%rdx
  pxor   %xmm0,%xmm0
  add    $0x4,%rax
  cvtss2sd -0x4(%rax),%xmm0
  movsd  %xmm0,(%rdx,%rcx,1)
  cmp    %rax,%rdi
  jne    L42
  add    $0x8,%rcx
  add    %ebx,%r10d
  cmp    %r11,%rcx
  jne    L69
L39:
  mov    %r9,%rdi
  call   <fclose@plt>
  test   %rax,%rax
  jne    L86
L38:
  mov    0x78(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L71
  add    $0x88,%rsp
  mov    %r13,%rdi
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    10c0 <free@plt>
L86:
  mov    (%rsp),%rdx
  lea    0x20(%rsp),%r12
  lea    D48(%rip),%rsi
  mov    %r12,%rdi
  call   <sprintf@plt>
  mov    %r12,%rdx
  mov    $0x13f,%edi
  lea    D40(%rip),%rsi
  call   <errhandler@plt>
  jmp    L38
L3:
  lea    0x20(%rsp),%r11
  mov    (%rsp),%rdx
  mov    %r10d,0x14(%rsp)
  lea    D19(%rip),%rsi
  mov    %r11,%rdi
  mov    %r11,0x8(%rsp)
  mov    %rax,0x18(%rsp)
  call   <sprintf@plt>
  mov    0x8(%rsp),%rdx
  mov    $0x12f,%edi
  lea    D40(%rip),%rsi
  call   <errhandler@plt>
  mov    0x18(%rsp),%r9
  mov    0x14(%rsp),%r10d
  jmp    L11
L71:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: retrieve_data
arguments: label: D4
type: string
name: retrieve_data
arguments: label: D48
type: string
name: retrieve_data
arguments: label: D40
type: string
name: retrieve_data
arguments: label: D19
type: string

tool: D4: r

tool: D48: Inputfile %s cannot be closed

tool: D40: /scratch/repos/new/home/jordi_armengol_estape/c-scraper/outputs/2022-01-22/02-19-57/repos/grantathon/computational_fluid_dynamics/refs/heads/master/project/ns_sim_parallel/helper.c

tool: D19: Can not read file %s !!!

assistant: ```c
void read_matrix(const char *szFileName, double **m, int nrl, int nrh, int ncl, int nch) {
  int i, j;
  FILE *fh = 0;
  int nSize = (nrh - nrl + 1) * (nch - ncl + 1);
  float *tmp = (float *)malloc((size_t)(nSize * sizeof(float)));
  int k = 0;
  fh = fopen(szFileName, "r");
  if (fh == NULL) {
    char szBuff[80];
    sprintf(szBuff, "Can not read file %s !!!", szFileName);
    errhandler(303, "/scratch/repos/new/home/jordi_armengol_estape/c-scraper/outputs/2022-01-22/02-19-57/repos/grantathon/computational_fluid_dynamics/refs/heads/master/project/ns_sim_parallel/helper.c", szBuff);
  }
  fread(tmp, sizeof(float), nSize, fh);
  for (j = ncl; j <= nch; j++)
    for (i = nrl; i <= nrh; i++)
      m[i][j] = tmp[k++];
  if (fclose(fh)) {
    char szBuff[80];
    sprintf(szBuff, "Inputfile %s cannot be closed", szFileName);
    errhandler(319, "/scratch/repos/new/home/jordi_armengol_estape/c-scraper/outputs/2022-01-22/02-19-57/repos/grantathon/computational_fluid_dynamics/refs/heads/master/project/ns_sim_parallel/helper.c", szBuff);
  };
  free(tmp);
}
```