user: 请从给定的汇编指令推导出C源代码：

```asm
<workload>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  lea    -0x60000(%rsp),%r11
L31:
  sub    $0x1000,%rsp
  orq    $0x0,(%rsp)
  cmp    %r11,%rsp
  jne    L31
  sub    $0x48,%rsp
  xor    %r13d,%r13d
  xor    %r15d,%r15d
  mov    %rdi,0x10(%rsp)
  lea    0x20030(%rsp),%rbx
  lea    0x40030(%rsp),%r8
  mov    %rsi,0x18(%rsp)
  mov    %rdx,0x20(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x60038(%rsp)
  xor    %eax,%eax
  lea    0x30(%rsp),%rax
  mov    %rax,(%rsp)
L47:
  mov    0x20(%rsp),%rdx
  lea    0x0(,%r13,8),%rax
  xor    %ebp,%ebp
  lea    (%rdx,%rax,1),%r12
  add    0x10(%rsp),%rax
  mov    %rax,0x8(%rsp)
  lea    0x20000(%r8),%rax
  mov    %rax,0x28(%rsp)
L7:
  mov    $0x20000,%edx
  mov    %r8,%rdi
  xor    %esi,%esi
  call   <memset@plt>
  mov    0x8(%rsp),%rdx
  xor    %r11d,%r11d
  mov    %rax,%r8
  mov    0x18(%rsp),%rax
  lea    (%rax,%rbp,8),%r9
L93:
  mov    (%rsp),%rax
  mov    %rdx,%r10
L42:
  mov    (%r10),%rcx
  lea    0x8(%rax),%rdi
  mov    %r10,%rsi
  add    $0x2000,%r10
  and    $0xfffffffffffffff8,%rdi
  mov    %rcx,(%rax)
  mov    -0x1c08(%r10),%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,%rbx
  jne    L42
  mov    %r9,%r10
  mov    %rbx,%rax
L34:
  mov    (%r10),%rcx
  lea    0x8(%rax),%rdi
  mov    %r10,%rsi
  add    $0x2000,%r10
  and    $0xfffffffffffffff8,%rdi
  mov    %rcx,(%rax)
  mov    -0x1c08(%r10),%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,%r8
  jne    L34
  xor    %r14d,%r14d
L23:
  mov    (%rsp),%rax
  mov    %rbx,%rsi
  lea    (%r8,%r14,1),%rcx
  lea    (%rax,%r14,1),%rdi
L15:
  movsd  (%rdi),%xmm1
  mov    %r8,%r10
  xor    %eax,%eax
  unpcklpd %xmm1,%xmm1
L11:
  movapd (%rsi,%rax,1),%xmm0
  mulpd  %xmm1,%xmm0
  addpd  (%rcx,%rax,1),%xmm0
  movaps %xmm0,(%rcx,%rax,1)
  add    $0x10,%rax
  cmp    $0x400,%rax
  jne    L11
  add    $0x400,%rsi
  add    $0x8,%rdi
  cmp    %rsi,%r8
  jne    L15
  add    $0x400,%r14
  cmp    $0x20000,%r14
  jne    L23
  sub    $0xffffffffffffff80,%r11
  add    $0x100000,%r9
  add    $0x400,%rdx
  cmp    $0x400,%r11
  jne    L93
  mov    %r12,%rdx
L57:
  mov    (%r10),%rax
  lea    0x8(%rdx),%rdi
  mov    %r10,%rsi
  add    $0x400,%r10
  and    $0xfffffffffffffff8,%rdi
  mov    %rax,(%rdx)
  mov    -0x8(%r10),%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    %r10,0x28(%rsp)
  jne    L57
  sub    $0xffffffffffffff80,%rbp
  add    $0x400,%r12
  cmp    $0x400,%rbp
  jne    L7
  sub    $0xffffff80,%r15d
  add    $0x20000,%r13
  cmp    $0x400,%r15d
  jne    L47
  mov    0x60038(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L60
  add    $0x60048,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L60:
  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;
}
```