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

```asm
<calc_thermal_matrix_LU>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %rbx
  sub    $0x58,%rsp
  mov    %rdi,-0x58(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x18(%rbp)
  xor    %eax,%eax
  lea    D4(%rip),%rax
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <printf@plt>
  lea    -0x24(%rbp),%rdx
  mov    -0x58(%rbp),%rax
  mov    %rdx,%rsi
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <set_up_tfield@plt>
  mov    -0x20(%rbp),%edx
  mov    -0x1c(%rbp),%eax
  imul   %edx,%eax
  mov    %eax,-0x44(%rbp)
  mov    -0x44(%rbp),%eax
  add    %eax,%eax
  add    $0x1,%eax
  mov    %eax,-0x40(%rbp)
  mov    -0x24(%rbp),%eax
  mov    -0x44(%rbp),%edx
  imul   %edx,%eax
  mov    %eax,-0x3c(%rbp)
  mov    -0x1c(%rbp),%edx
  mov    -0x44(%rbp),%ecx
  mov    -0x58(%rbp),%rax
  mov    %ecx,%esi
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <find_contained_nodes@plt>
  mov    -0x3c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  mov    %rax,%rdi
  call   <malloc@plt>
  mov    %rax,-0x38(%rbp)
  movl   $0x0,-0x4c(%rbp)
  jmp    L20
L33:
  mov    -0x4c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x4c(%rbp)
L20:
  mov    -0x4c(%rbp),%eax
  cmp    -0x3c(%rbp),%eax
  jl     L33
  mov    -0x3c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  mov    %rax,%rdi
  call   <malloc@plt>
  mov    %rax,-0x30(%rbp)
  movl   $0x0,-0x4c(%rbp)
  jmp    L63
L42:
  mov    -0x40(%rbp),%eax
  cltq
  shl    $0x3,%rax
  mov    -0x4c(%rbp),%edx
  movslq %edx,%rdx
  lea    0x0(,%rdx,8),%rcx
  mov    -0x30(%rbp),%rdx
  lea    (%rcx,%rdx,1),%rbx
  mov    %rax,%rdi
  call   <malloc@plt>
  mov    %rax,(%rbx)
  movl   $0x0,-0x48(%rbp)
  jmp    L11
L87:
  mov    -0x4c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x48(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x48(%rbp)
L11:
  mov    -0x48(%rbp),%eax
  cmp    -0x40(%rbp),%eax
  jl     L87
  addl   $0x1,-0x4c(%rbp)
L63:
  mov    -0x4c(%rbp),%eax
  cmp    -0x3c(%rbp),%eax
  jl     L42
  mov    -0x3c(%rbp),%edx
  mov    -0x38(%rbp),%rcx
  mov    -0x58(%rbp),%rax
  mov    %rcx,%rsi
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <set_up_flux@plt>
  mov    -0x1c(%rbp),%r8d
  mov    -0x44(%rbp),%edi
  mov    -0x3c(%rbp),%ecx
  mov    -0x30(%rbp),%rdx
  mov    -0x38(%rbp),%rsi
  mov    -0x58(%rbp),%rax
  mov    %r8d,%r9d
  mov    %edi,%r8d
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <set_up_coef_matrix@plt>
  mov    -0x44(%rbp),%edi
  mov    -0x3c(%rbp),%ecx
  mov    -0x38(%rbp),%rdx
  mov    -0x30(%rbp),%rsi
  mov    -0x58(%rbp),%rax
  mov    %edi,%r8d
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <LU_Decomp@plt>
  mov    -0x58(%rbp),%rax
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <find_comp_temp@plt>
  movl   $0x0,-0x4c(%rbp)
  jmp    L5
L22:
  mov    -0x4c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  mov    %rax,%rdi
  call   <free@plt>
  addl   $0x1,-0x4c(%rbp)
L5:
  mov    -0x4c(%rbp),%eax
  cmp    -0x3c(%rbp),%eax
  jl     L22
  mov    -0x30(%rbp),%rax
  mov    %rax,%rdi
  call   <free@plt>
  mov    -0x38(%rbp),%rax
  mov    %rax,%rdi
  call   <free@plt>
  nop
  mov    -0x18(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L69
  call   <__stack_chk_fail@plt>
L69:
  mov    -0x8(%rbp),%rbx
  leave
  ret
```

tool_calls: name: retrieve_data
arguments: tag: D4
type_infered: string

tool: D4: LU   

assistant: ```c
void calc_thermal_matrix_LU(design) struct Design *design;
{
  double *flux, **R;
  void set_up_tfield(), set_up_flux(), set_up_coef_matrix();
  void LU_Decomp(), find_comp_temp(), find_contained_nodes();
  int i, j, node_dim[3], tot_nodes, hbw, width;
  printf("LU   ");
  set_up_tfield(design, node_dim);
  hbw = node_dim[1] * node_dim[2];
  width = 2 * hbw + 1;
  tot_nodes = node_dim[0] * hbw;
  find_contained_nodes(design, hbw, node_dim[2]);
  flux = (double *)malloc(sizeof(double) * tot_nodes);
  for (i = 0; i < tot_nodes; i++)
    flux[i] = 0.0;
  R = (double **)malloc(sizeof(double *) * tot_nodes);
  for (i = 0; i < tot_nodes; i++) {
    R[i] = (double *)malloc(sizeof(double) * width);
    for (j = 0; j < width; ++j)
      R[i][j] = 0.0;
  }
  set_up_flux(design, flux, tot_nodes);
  set_up_coef_matrix(design, flux, R, tot_nodes, hbw, node_dim[2]);
  LU_Decomp(design, R, flux, tot_nodes, hbw);
  find_comp_temp(design);
  for (i = 0; i < tot_nodes; i++)
    free(R[i]);
  free(R);
  free(flux);
}
```