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

```asm
<build_SPD_Phi>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x50,%rsp
  mov    %edi,-0x34(%rbp)
  mov    %rsi,-0x40(%rbp)
  mov    %rdx,-0x48(%rbp)
  mov    %rcx,-0x50(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  movsd  D2(%rip),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x34(%rbp),%edx
  mov    -0x34(%rbp),%eax
  imul   %edx,%eax
  cltq
  shl    $0x3,%eax
  mov    %eax,%edx
  mov    $0x1ab,%esi
  lea    D21(%rip),%rax
  mov    %rax,%rdi
  call   <fgls_malloc_impl@plt>
  mov    %rax,-0x10(%rbp)
  movl   $0x0,-0x24(%rbp)
  jmp    L21
L47:
  movl   $0x0,-0x28(%rbp)
  jmp    L25
L37:
  mov    -0x34(%rbp),%eax
  imul   -0x24(%rbp),%eax
  mov    %eax,%edx
  mov    -0x28(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x40(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x34(%rbp),%eax
  imul   -0x24(%rbp),%eax
  mov    %eax,%edx
  mov    -0x28(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x10(%rbp),%rax
  add    %rdx,%rax
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x28(%rbp)
L25:
  mov    -0x34(%rbp),%eax
  cmp    %eax,-0x28(%rbp)
  jl     L37
  addl   $0x1,-0x24(%rbp)
L21:
  mov    -0x34(%rbp),%eax
  cmp    %eax,-0x24(%rbp)
  jl     L47
  lea    -0x20(%rbp),%rdi
  lea    -0x34(%rbp),%rsi
  lea    -0x34(%rbp),%rdx
  lea    -0x34(%rbp),%rax
  sub    $0x8,%rsp
  lea    -0x34(%rbp),%rcx
  push   %rcx
  push   -0x50(%rbp)
  lea    -0x18(%rbp),%rcx
  push   %rcx
  lea    -0x34(%rbp),%rcx
  push   %rcx
  push   -0x40(%rbp)
  lea    -0x34(%rbp),%rcx
  push   %rcx
  push   -0x10(%rbp)
  mov    %rdi,%r9
  mov    %rsi,%r8
  mov    %rdx,%rcx
  mov    %rax,%rdx
  lea    D31(%rip),%rax
  mov    %rax,%rsi
  lea    D36(%rip),%rax
  mov    %rax,%rdi
  call   <dgemm_@plt>
  add    $0x40,%rsp
  mov    -0x10(%rbp),%rax
  mov    %rax,%rdi
  call   <free@plt>
  nop
  mov    -0x8(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L6
  call   <__stack_chk_fail@plt>
L6:
  leave
  ret
```

tool_calls: name: parse_data
arguments: label: D2
type_guessed: f64
name: parse_data
arguments: label: D21
type_guessed: string
name: parse_data
arguments: label: D31
type_guessed: string
name: parse_data
arguments: label: D36
type_guessed: string

tool: D2: 1

tool: D21: /scratch/repos/new/home/jordi_armengol_estape/c-scraper/outputs/2022-01-22/02-19-57/repos/lucasb-eyer/OmicABEL/refs/heads/master/src/fgls_chol.c

tool: D31: T

tool: D36: N

assistant: ```c
void build_SPD_Phi(int n, double *eigVecs, double *eigVals, double *Phi) {
  double ONE = 1.0, ZERO = 0.0, *tmp;
  int i, j;
  tmp = fgls_malloc_impl("/scratch/repos/new/home/jordi_armengol_estape/c-scraper/outputs/2022-01-22/02-19-57/repos/lucasb-eyer/OmicABEL/refs/heads/master/src/fgls_chol.c", 427, n * n * sizeof(double));
  for (j = 0; j < n; j++)
    for (i = 0; i < n; i++)
      tmp[j * n + i] = eigVecs[j * n + i] * eigVals[j];
  dgemm_("N", "T", &n, &n, &n, &ONE, tmp, &n, eigVecs, &n, &ZERO, Phi, &n);
  free(tmp);
}
```