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

```asm
<workload>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  lea    -0x60000(%rsp),%r11
L3:
  sub    $0x1000,%rsp
  orq    $0x0,(%rsp)
  cmp    %r11,%rsp
  jne    L3
  sub    $0x68,%rsp
  xor    %r13d,%r13d
  mov    %rdi,0x20(%rsp)
  lea    0x20050(%rsp),%r12
  lea    0x40050(%rsp),%r15
  mov    %rsi,0x28(%rsp)
  lea    0x60450(%rsp),%rbx
  mov    %rdx,0x30(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x60058(%rsp)
  xor    %eax,%eax
  lea    0x40450(%rsp),%rax
  movl   $0x0,0xc(%rsp)
  mov    %rax,0x40(%rsp)
  lea    0x50(%rsp),%rax
  mov    %rax,0x38(%rsp)
L9:
  mov    0x30(%rsp),%rdi
  lea    0x0(,%r13,8),%rax
  mov    %r13,0x48(%rsp)
  xor    %ebp,%ebp
  add    %rax,%rdi
  add    0x20(%rsp),%rax
  mov    %rdi,0x10(%rsp)
  mov    %rax,0x18(%rsp)
L72:
  mov    $0x20000,%edx
  xor    %esi,%esi
  mov    %r15,%rdi
  xor    %r13d,%r13d
  call   <memset@plt>
  mov    0x28(%rsp),%rax
  mov    0x18(%rsp),%r9
  lea    (%rax,%rbp,8),%r10
L12:
  mov    0x38(%rsp),%rax
  mov    %r9,%rdx
  mov    %rax,%r11
L60:
  mov    (%rdx),%rcx
  lea    0x8(%rax),%rdi
  mov    %rdx,%rsi
  add    $0x2000,%rdx
  and    $0xfffffffffffffff8,%rdi
  mov    %rcx,(%rax)
  mov    -0x1c08(%rdx),%rcx
  mov    %rcx,0x3f8(%rax)
  mov    %rax,%rcx
  add    $0x400,%rax
  sub    %rdi,%rcx
  sub    %rcx,%rsi
  add    $0x400,%ecx
  shr    $0x3,%ecx
  rep movsq %ds:(%rsi),%es:(%rdi)
  cmp    %rax,%r12
  jne    L60
  mov    %r10,%rdx
  mov    %r12,%rax
L81:
  mov    (%rdx),%rcx
  lea    0x8(%rax),%rdi
  mov    %rdx,%rsi
  add    $0x2000,%rdx
  and    $0xfffffffffffffff8,%rdi
  mov    %rcx,(%rax)
  mov    -0x1c08(%rdx),%rcx
  mov    %rcx,0x3f8(%rax)
  mov    %rax,%rcx
  add    $0x400,%rax
  sub    %rdi,%rcx
  sub    %rcx,%rsi
  add    $0x400,%ecx
  shr    $0x3,%ecx
  rep movsq %ds:(%rsi),%es:(%rdi)
  cmp    %rax,%r15
  jne    L81
  xor    %r14d,%r14d
L26:
  mov    %r14,%rdi
  mov    0x40(%rsp),%rcx
  mov    %r15,%r8
  mov    %r11,%rsi
  shl    $0xa,%rdi
  add    %r12,%rdi
L85:
  movsd  (%rsi),%xmm1
  lea    -0x400(%rcx),%rax
  mov    %rdi,%rdx
L10:
  movsd  (%rdx),%xmm0
  add    $0x8,%rax
  add    $0x8,%rdx
  mulsd  %xmm1,%xmm0
  addsd  -0x8(%rax),%xmm0
  movsd  %xmm0,-0x8(%rax)
  cmp    %rcx,%rax
  jne    L10
  lea    0x400(%rax),%rcx
  add    $0x400,%rsi
  cmp    %rbx,%rcx
  jne    L85
  add    $0x1,%r14
  add    $0x8,%r11
  cmp    $0x80,%r14
  jne    L26
  sub    $0xffffffffffffff80,%r13
  add    $0x100000,%r10
  add    $0x400,%r9
  cmp    $0x400,%r13
  jne    L12
  mov    0x10(%rsp),%rdx
  lea    0x20000(%r15),%r9
L76:
  mov    (%r8),%rax
  lea    0x8(%rdx),%rdi
  mov    %r8,%rsi
  add    $0x400,%r8
  and    $0xfffffffffffffff8,%rdi
  mov    %rax,(%rdx)
  mov    -0x8(%r8),%rax
  mov    %rax,0x3f8(%rdx)
  mov    %rdx,%rax
  add    $0x2000,%rdx
  sub    %rdi,%rax
  lea    0x400(%rax),%ecx
  sub    %rax,%rsi
  shr    $0x3,%ecx
  rep movsq %ds:(%rsi),%es:(%rdi)
  cmp    %r8,%r9
  jne    L76
  addq   $0x400,0x10(%rsp)
  sub    $0xffffffffffffff80,%rbp
  cmp    $0x400,%rbp
  jne    L72
  subl   $0xffffff80,0xc(%rsp)
  mov    0x48(%rsp),%r13
  mov    0xc(%rsp),%eax
  add    $0x20000,%r13
  cmp    $0x400,%eax
  jne    L9
  mov    0x60058(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L52
  add    $0x60068,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L52:
  call   <__stack_chk_fail@plt>
```

assistant: ```c
void workload(double m1[1024 * 1024], double m2[1024 * 1024], double prod[1024 * 1024]) {
#pragma HLS INTERFACE m_axi port = m1 offset = slave bundle = gmem
#pragma HLS INTERFACE m_axi port = m2 offset = slave bundle = gmem
#pragma HLS INTERFACE m_axi port = prod offset = slave bundle = gmem
#pragma HLS INTERFACE s_axilite port = m1 bundle = control
#pragma HLS INTERFACE s_axilite port = m2 bundle = control
#pragma HLS INTERFACE s_axilite port = prod bundle = control
#pragma HLS INTERFACE s_axilite port = return bundle = control
  double local_m1[128][128];
  double local_m2[128][128];
#pragma HLS ARRAY_PARTITION variable = local_m2 complete dim = 2
  double local_prod[128][128];
#pragma HLS ARRAY_PARTITION variable = local_prod complete dim = 2
  double mult, sum;
  int i, j, k, ii, jj, kk;
  int uu;
  for (i = 0; i < 1024; i += 128)
    for (j = 0; j < 1024; j += 128) {
      for (ii = 0; ii < 128; ii++)
        for (jj = 0; jj < 128; jj++)
#pragma HLS PIPELINE
          local_prod[ii][jj] = 0;
      for (k = 0; k < 1024; k += 128) {
        for (ii = 0; ii < 128; ii++)
          for (kk = 0; kk < 128; kk++)
#pragma HLS PIPELINE
            local_m1[ii][kk] = m1[(i + ii) * 1024 + (k + kk)];
        for (kk = 0; kk < 128; kk++)
          for (jj = 0; jj < 128; jj++)
#pragma HLS PIPELINE
            local_m2[kk][jj] = m2[(k + kk) * 1024 + (j + jj)];
        for (kk = 0; kk < 128; kk++)
          for (ii = 0; ii < 128; ii++)
            for (jj = 0; jj < 128; jj += 128)
#pragma HLS PIPELINE II = 20
              for (uu = 0; uu < 128; uu++)
#pragma HLS UNROLL
#pragma HLS DEPENDENCE variable = "local_prod" inter false
                local_prod[ii][jj + uu] += local_m1[ii][kk] * local_m2[kk][jj + uu];
      }
      for (ii = 0; ii < 128; ii++)
        for (jj = 0; jj < 128; jj++)
#pragma HLS PIPELINE
          prod[(i + ii) * 1024 + (j + jj)] = local_prod[ii][jj];
    }
  return;
}
```