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

```asm
<main>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x78,%rsp
  cmp    $0x1,%edi
  jle    L135
  mov    0x8(%rsi),%rdi
  call   <atoi@plt>
  mov    %eax,%ebx
  movslq %eax,%rbp
  mov    %eax,%r12d
  imul   %eax,%r12d
  shl    $0x3,%r12d
  mov    %r12d,%edi
  call   <malloc@plt>
  mov    %rax,%r15
  mov    %rax,0x28(%rsp)
  mov    %r12d,%edi
  call   <malloc@plt>
  mov    %rax,%r13
  mov    %rax,0x50(%rsp)
  mov    %r12d,%edi
  call   <malloc@plt>
  mov    %rax,%r14
  mov    %rax,0x40(%rsp)
  mov    %r12d,%edi
  call   <malloc@plt>
  mov    %rax,%rdx
  mov    %rax,0x58(%rsp)
  test   %ebx,%ebx
  jle    L132
  lea    0x0(,%rbp,8),%r10
  mov    %r10,0x18(%rsp)
  lea    0x8(%r15),%rbp
  lea    -0x1(%rbx),%eax
  lea    0x0(%rbp,%rax,8),%rcx
  mov    %r13,%r8
  mov    $0x0,%r11d
  mov    %ebx,%r9d
  neg    %r9
  shl    $0x3,%r9
  movsd  D43(%rip),%xmm0
  mov    %r14,%rdi
  sub    %r15,%rdi
  sub    %r15,%rdx
  mov    %rdx,%rsi
L136:
  lea    (%r9,%rcx,1),%rax
  mov    %r8,%rdx
L94:
  movsd  %xmm0,(%rax)
  movsd  %xmm0,(%rdx)
  movq   $0x0,(%rdi,%rax,1)
  movsd  %xmm0,(%rsi,%rax,1)
  add    $0x8,%rax
  add    %r10,%rdx
  cmp    %rcx,%rax
  jne    L94
  lea    0x1(%r11),%eax
  add    $0x8,%r8
  add    %r10,%rcx
  cmp    %eax,%ebx
  je     L13
  mov    %eax,%r11d
  jmp    L136
L135:
  lea    D55(%rip),%rdi
  call   <puts@plt>
  jmp    L189
L89:
  add    $0x1,%eax
  add    $0x8,%rcx
  cmp    %ebx,%eax
  jge    L157
L74:
  test   %edx,%edx
  je     L89
  ucomisd (%rcx),%xmm0
  setnp  %dl
  movzbl %dl,%edx
  cmovne %esi,%edx
  jmp    L89
L157:
  test   %edx,%edx
  jne    L72
  lea    D13(%rip),%rdi
  call   <puts@plt>
  jmp    L100
L78:
  mov    $0x0,%eax
  call   <dwalltime@plt>
  subsd  0x60(%rsp),%xmm0
  lea    D28(%rip),%rdi
  mov    $0x1,%eax
  call   <printf@plt>
  mov    0x6c(%rsp),%ebx
  imul   %ebx,%ebx
L186:
  mov    0x40(%rsp),%rcx
  mov    $0x1,%edx
  mov    $0x0,%eax
  pxor   %xmm0,%xmm0
  mov    $0x0,%esi
  jmp    L74
L34:
  divsd  %xmm2,%xmm0
  movsd  %xmm0,0x8(%rsp)
  mov    0x40(%rsp),%rax
  lea    0x8(%rax,%r14,1),%rax
  mov    %rax,0x20(%rsp)
  add    %rbp,%r14
  imul   $0xfffffffffffffff8,%r8,%rax
  mov    %rax,0x48(%rsp)
  mov    0x50(%rsp),%rax
  sub    0x28(%rsp),%rax
  mov    %rax,0x30(%rsp)
  movq   $0x0,0x10(%rsp)
  movl   $0x0,0x3c(%rsp)
L183:
  mov    0x10(%rsp),%r15
  add    0x40(%rsp),%r15
  mov    0x30(%rsp),%r13
L156:
  mov    %r15,%r12
  mov    0x10(%rsp),%rbp
  add    0x28(%rsp),%rbp
L111:
  movsd  0x0(%rbp),%xmm0
  subsd  0x8(%rsp),%xmm0
  mov    $0x2,%edi
  call   <pow@plt>
  mov    %eax,%ebx
  movsd  0x0(%rbp,%r13,1),%xmm0
  subsd  0x8(%rsp),%xmm0
  mov    $0x2,%edi
  call   <pow@plt>
  imul   %eax,%ebx
  mov    %ebx,%edi
  call   <sqrt@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %rax,%xmm0
  addsd  (%r12),%xmm0
  movsd  %xmm0,(%r12)
  add    $0x8,%rbp
  cmp    %r14,%rbp
  jne    L111
  add    0x18(%rsp),%r13
  add    $0x8,%r15
  cmp    0x20(%rsp),%r15
  jne    L156
  mov    0x3c(%rsp),%ecx
  lea    0x1(%rcx),%eax
  mov    0x18(%rsp),%rsi
  add    %rsi,0x20(%rsp)
  add    %rsi,0x10(%rsp)
  add    %rsi,%r14
  mov    0x48(%rsp),%rdx
  add    %rdx,0x30(%rsp)
  cmp    %ecx,0x68(%rsp)
  je     L78
  mov    %eax,0x3c(%rsp)
  jmp    L183
L13:
  mov    %eax,%r15d
  mov    %eax,0x6c(%rsp)
  mov    %r11d,%ebx
  mov    %r11d,0x68(%rsp)
  mov    $0x0,%eax
  call   <dwalltime@plt>
  movsd  %xmm0,0x60(%rsp)
  movslq %r15d,%r8
  mov    %ebx,%r9d
  mov    %ebx,%eax
  lea    0x0(,%rax,8),%r14
  mov    0x58(%rsp),%rbx
  lea    0x8(%rbx,%r14,1),%rdx
  pxor   %xmm2,%xmm2
  movapd %xmm2,%xmm0
  mov    $0x0,%esi
  not    %rax
  lea    0x0(,%rax,8),%rdi
  mov    0x28(%rsp),%rcx
  sub    %rbx,%rcx
  mov    0x18(%rsp),%r10
L24:
  lea    (%rdi,%rdx,1),%rax
L57:
  movsd  (%rax),%xmm1
  addsd  %xmm1,%xmm2
  mulsd  (%rcx,%rax,1),%xmm1
  addsd  %xmm1,%xmm0
  add    $0x8,%rax
  cmp    %rdx,%rax
  jne    L57
  lea    0x1(%rsi),%eax
  add    %r10,%rdx
  cmp    %esi,%r9d
  je     L34
  mov    %eax,%esi
  jmp    L24
L132:
  mov    $0x0,%eax
  call   <dwalltime@plt>
  movq   %xmm0,%rbp
  mov    $0x0,%eax
  call   <dwalltime@plt>
  movq   %rbp,%xmm5
  subsd  %xmm5,%xmm0
  lea    D28(%rip),%rdi
  mov    $0x1,%eax
  call   <printf@plt>
  imul   %ebx,%ebx
  test   %ebx,%ebx
  jg     L186
L72:
  lea    D48(%rip),%rdi
  call   <puts@plt>
L100:
  mov    0x28(%rsp),%rdi
  call   <free@plt>
  mov    0x50(%rsp),%rdi
  call   <free@plt>
  mov    0x40(%rsp),%rdi
  call   <free@plt>
  mov    0x58(%rsp),%rdi
  call   <free@plt>
L189:
  mov    $0x0,%eax
  add    $0x78,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: extract_symbolic_data
arguments: data_label: D43
data_type: f64
name: extract_symbolic_data
arguments: data_label: D55
data_type: string
name: extract_symbolic_data
arguments: data_label: D13
data_type: string
name: extract_symbolic_data
arguments: data_label: D28
data_type: string
name: extract_symbolic_data
arguments: data_label: D48
data_type: string

tool: D43: 1

tool: D55: 
 Faltan argumentos:: N dimension de la matriz 

tool: D13: Resultado erroneo

tool: D28: Tiempo en segundos: %f 


tool: D48: Resultado correcto

assistant: ```c
main(int argc, char *argv[]) {
  double *A, *B, *C, *W;
  int i, j, k, N;
  int check = 1;
  double timetick;
  if (argc < 2) {
    printf("\n Faltan argumentos:: N dimension de la matriz \n");
    return 0;
  }
  N = atoi(argv[1]);
  A = (double *)malloc(sizeof(double) * N * N);
  B = (double *)malloc(sizeof(double) * N * N);
  C = (double *)malloc(sizeof(double) * N * N);
  W = (double *)malloc(sizeof(double) * N * N);
  for (i = 0; i < N; i++) {
    for (j = 0; j < N; j++) {
      A[i * N + j] = 1.0;
      B[i + N * j] = 1.0;
      C[i * N + j] = 0.0;
      W[i * N + j] = 1.0;
    }
  }
  timetick = dwalltime();
  double total = 0;
  double temp = 0;
  for (i = 0; i < N; i++) {
    for (j = 0; j < N; j++) {
      temp += W[i * N + j];
      total += A[i * N + j] * W[i * N + j];
    }
  }
  total /= temp;
  for (i = 0; i < N; i++) {
    for (j = 0; j < N; j++) {
      for (k = 0; k < N; k++) {
        C[i * N + j] += sqrt((pow(A[i * N + k] - total, 2)) * (pow(B[j * N + k] - total, 2)));
      }
    }
  }
  printf("Tiempo en segundos: %f \n", dwalltime() - timetick);
  for (i = 0; i < N * N; i++) {
    check = check && (C[i] == 0.0);
  }
  if (check) {
    printf("Resultado correcto\n");
  } else {
    printf("Resultado erroneo\n");
  }
  free(A);
  free(B);
  free(C);
  free(W);
  return (0);
}
```