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

```asm
<quaddobl_scaler>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbx
  sub    $0x38,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,-0x38(%rbp)
  xor    %eax,%eax
  mov    $0x0,%eax
  call   <syscon_read_quaddobl_system@plt>
  mov    %eax,-0x50(%rbp)
  cmpl   $0x0,-0x50(%rbp)
  jne    L16
  lea    D4(%rip),%rax
  mov    %rax,%rdi
  call   <puts@plt>
  mov    $0x0,%eax
  call   <syscon_write_quaddobl_system@plt>
  mov    %eax,-0x50(%rbp)
  lea    -0x54(%rbp),%rax
  mov    %rax,%rdi
  call   <syscon_number_of_quaddobl_polynomials@plt>
  mov    %eax,-0x50(%rbp)
  cmpl   $0x0,-0x50(%rbp)
  jne    L16
  mov    %rsp,%rax
  mov    %rax,%rbx
  mov    -0x54(%rbp),%eax
  shl    $0x4,%eax
  add    $0x8,%eax
  movslq %eax,%rdx
  sub    $0x1,%rdx
  mov    %rdx,-0x48(%rbp)
  movslq %eax,%rdx
  mov    %rdx,%r14
  mov    $0x0,%r15d
  movslq %eax,%rdx
  mov    %rdx,%r12
  mov    $0x0,%r13d
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    $0x10,%eax
  sub    $0x1,%rax
  add    %rdx,%rax
  mov    $0x10,%esi
  mov    $0x0,%edx
  div    %rsi
  imul   $0x10,%rax,%rax
  mov    %rax,%rcx
  and    $0xfffffffffffff000,%rcx
  mov    %rsp,%rdx
  sub    %rcx,%rdx
L40:
  cmp    %rdx,%rsp
  je     L9
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L40
L9:
  mov    %rax,%rdx
  and    $0xfff,%edx
  sub    %rdx,%rsp
  mov    %rax,%rdx
  and    $0xfff,%edx
  test   %rdx,%rdx
  je     L20
  and    $0xfff,%eax
  sub    $0x8,%rax
  add    %rsp,%rax
  orq    $0x0,(%rax)
L20:
  mov    %rsp,%rax
  add    $0x7,%rax
  shr    $0x3,%rax
  shl    $0x3,%rax
  mov    %rax,-0x40(%rbp)
  mov    -0x40(%rbp),%rax
  mov    %rax,%rsi
  mov    $0x2,%edi
  call   <quaddobl_scale_system@plt>
  mov    %eax,-0x50(%rbp)
  mov    -0x54(%rbp),%eax
  shl    $0x4,%eax
  mov    %eax,%edx
  mov    -0x40(%rbp),%rax
  movslq %edx,%rdx
  mov    (%rax,%rdx,8),%rax
  movq   %rax,%xmm0
  lea    D39(%rip),%rax
  mov    %rax,%rdi
  mov    $0x1,%eax
  call   <printf@plt>
  lea    D43(%rip),%rax
  mov    %rax,%rdi
  call   <puts@plt>
  movl   $0x0,-0x4c(%rbp)
  jmp    L13
L63:
  mov    -0x4c(%rbp),%eax
  lea    0x7(%rax),%edx
  mov    -0x40(%rbp),%rax
  movslq %edx,%rdx
  movsd  (%rax,%rdx,8),%xmm6
  mov    -0x4c(%rbp),%eax
  lea    0x6(%rax),%edx
  mov    -0x40(%rbp),%rax
  movslq %edx,%rdx
  movsd  (%rax,%rdx,8),%xmm5
  mov    -0x4c(%rbp),%eax
  lea    0x5(%rax),%edx
  mov    -0x40(%rbp),%rax
  movslq %edx,%rdx
  movsd  (%rax,%rdx,8),%xmm4
  mov    -0x4c(%rbp),%eax
  lea    0x4(%rax),%edx
  mov    -0x40(%rbp),%rax
  movslq %edx,%rdx
  movsd  (%rax,%rdx,8),%xmm3
  mov    -0x4c(%rbp),%eax
  lea    0x3(%rax),%edx
  mov    -0x40(%rbp),%rax
  movslq %edx,%rdx
  movsd  (%rax,%rdx,8),%xmm2
  mov    -0x4c(%rbp),%eax
  lea    0x2(%rax),%edx
  mov    -0x40(%rbp),%rax
  movslq %edx,%rdx
  movsd  (%rax,%rdx,8),%xmm1
  mov    -0x4c(%rbp),%eax
  lea    0x1(%rax),%edx
  mov    -0x40(%rbp),%rax
  movslq %edx,%rdx
  movsd  (%rax,%rdx,8),%xmm0
  mov    -0x40(%rbp),%rax
  mov    -0x4c(%rbp),%edx
  movslq %edx,%rdx
  mov    (%rax,%rdx,8),%rdx
  mov    -0x4c(%rbp),%eax
  lea    0x7(%rax),%ecx
  test   %eax,%eax
  cmovs  %ecx,%eax
  sar    $0x3,%eax
  movapd %xmm6,%xmm7
  movapd %xmm5,%xmm6
  movapd %xmm4,%xmm5
  movapd %xmm3,%xmm4
  movapd %xmm2,%xmm3
  movapd %xmm1,%xmm2
  movapd %xmm0,%xmm1
  movq   %rdx,%xmm0
  mov    %eax,%esi
  lea    D27(%rip),%rax
  mov    %rax,%rdi
  mov    $0x8,%eax
  call   <printf@plt>
  addl   $0x4,-0x4c(%rbp)
L13:
  mov    -0x54(%rbp),%eax
  shl    $0x4,%eax
  cmp    %eax,-0x4c(%rbp)
  jl     L63
  lea    D4(%rip),%rax
  mov    %rax,%rdi
  call   <puts@plt>
  mov    $0x0,%eax
  call   <syscon_write_quaddobl_system@plt>
  mov    %eax,-0x50(%rbp)
  mov    %rbx,%rsp
L16:
  mov    -0x50(%rbp),%eax
  mov    -0x38(%rbp),%rdx
  sub    %fs:0x28,%rdx
  je     L15
  call   <__stack_chk_fail@plt>
L15:
  lea    -0x28(%rbp),%rsp
  pop    %rbx
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  pop    %rbp
  ret
```

tool_calls: name: retrieve_data
arguments: label: D4
type_infered: string
name: retrieve_data
arguments: label: D39
type_infered: string
name: retrieve_data
arguments: label: D43
type_infered: string
name: retrieve_data
arguments: label: D27
type_infered: string

tool: D4: 
The system in the container : 

tool: D39: The estimated inverse condition number : %.3e


tool: D43: The scaling coefficients : 

tool: D27: %2d : %.15e %.15e %.15e %.15e  %.15e %.15e %.15e %.15e


assistant: ```c
int quaddobl_scaler(void) {
  int fail, dim, i;
  fail = syscon_read_quaddobl_system();
  if (fail == 0) {
    printf("\nThe system in the container : \n");
    fail = syscon_write_quaddobl_system();
    fail = syscon_number_of_quaddobl_polynomials(&dim);
    if (fail == 0) {
      double c[16 * dim + 8];
      fail = quaddobl_scale_system(2, c);
      printf("The estimated inverse condition number : %.3e\n", c[16 * dim]);
      printf("The scaling coefficients : \n");
      for (i = 0; i < 16 * dim; i = i + 4)
        printf("%2d : %.15e %.15e %.15e %.15e  %.15e %.15e %.15e %.15e\n", i / 8, c[i], c[i + 1], c[i + 2], c[i + 3], c[i + 4], c[i + 5], c[i + 6], c[i + 7]);
      printf("\nThe system in the container : \n");
      fail = syscon_write_quaddobl_system();
    }
  }
  return fail;
}
```