user: Translate the below assembly code into C programming language:

```asm
<workload>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  lea    -0x60000(%rsp),%r11
L98:
  sub    $0x1000,%rsp
  orq    $0x0,(%rsp)
  cmp    %r11,%rsp
  jne    L98
  sub    $0x68,%rsp
  mov    %rsi,0x38(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x60058(%rsp)
  xor    %eax,%eax
  lea    0x40050(%rsp),%r10
  lea    0x50(%rsp),%rax
  mov    %rax,0x18(%rsp)
  mov    %rdi,0x20(%rsp)
  mov    %rdx,0x28(%rsp)
  movl   $0x0,0xc(%rsp)
  lea    0x60050(%rsp),%rax
  mov    %rax,0x10(%rsp)
  lea    0x20050(%rsp),%rax
  mov    %rax,0x30(%rsp)
  mov    %rax,0x40(%rsp)
  mov    %r10,0x48(%rsp)
  lea    0x60450(%rsp),%r8
L107:
  mov    0x38(%rsp),%rbp
  mov    0x28(%rsp),%rbx
  mov    $0x0,%r11d
L7:
  mov    %r10,%rax
L63:
  mov    %rax,%rdx
  add    $0x400,%rax
L156:
  movq   $0x0,(%rdx)
  add    $0x8,%rdx
  cmp    %rax,%rdx
  jne    L156
  cmp    %rax,0x10(%rsp)
  jne    L63
  mov    0x20(%rsp),%r14
  mov    %rbp,%r13
  mov    $0x0,%r12d
L94:
  mov    0x18(%rsp),%rdx
  mov    %r14,%rcx
L137:
  mov    $0x0,%eax
L154:
  movsd  (%rcx,%rax,1),%xmm0
  movsd  %xmm0,(%rdx,%rax,1)
  add    $0x8,%rax
  cmp    $0x400,%rax
  jne    L154
  add    $0x2000,%rcx
  add    $0x400,%rdx
  cmp    0x40(%rsp),%rdx
  jne    L137
  mov    0x30(%rsp),%rdx
  mov    %r13,%rcx
L83:
  mov    $0x0,%eax
L16:
  movsd  (%rcx,%rax,1),%xmm0
  movsd  %xmm0,(%rdx,%rax,1)
  add    $0x8,%rax
  cmp    $0x400,%rax
  jne    L16
  add    $0x2000,%rcx
  add    $0x400,%rdx
  cmp    0x48(%rsp),%rdx
  jne    L83
  mov    0x18(%rsp),%r15
  mov    $0x0,%r9d
L135:
  lea    0x400(%r10),%rcx
  mov    %r9,%rdi
  shl    $0xa,%rdi
  add    0x30(%rsp),%rdi
  mov    %r15,%rsi
L108:
  movsd  (%rsi),%xmm1
  lea    -0x400(%rcx),%rax
  mov    %rdi,%rdx
L74:
  movapd %xmm1,%xmm0
  mulsd  (%rdx),%xmm0
  addsd  (%rax),%xmm0
  movsd  %xmm0,(%rax)
  add    $0x8,%rax
  add    $0x8,%rdx
  cmp    %rcx,%rax
  jne    L74
  add    $0x400,%rsi
  add    $0x400,%rcx
  cmp    %r8,%rcx
  jne    L108
  add    $0x1,%r9
  add    $0x8,%r15
  cmp    $0x80,%r9
  jne    L135
  add    $0x100000,%r13
  sub    $0xffffffffffffff80,%r12
  add    $0x400,%r14
  cmp    $0x400,%r12
  jne    L94
  mov    %rbx,%rcx
  mov    %r10,%rdx
L110:
  mov    $0x0,%eax
L31:
  movsd  (%rdx,%rax,1),%xmm0
  movsd  %xmm0,(%rcx,%rax,1)
  add    $0x8,%rax
  cmp    $0x400,%rax
  jne    L31
  add    $0x400,%rdx
  add    $0x2000,%rcx
  cmp    0x10(%rsp),%rdx
  jne    L110
  sub    $0xffffffffffffff80,%r11
  add    $0x400,%rbx
  add    $0x400,%rbp
  cmp    $0x400,%r11
  jne    L7
  subl   $0xffffff80,0xc(%rsp)
  mov    0xc(%rsp),%eax
  addq   $0x100000,0x20(%rsp)
  addq   $0x100000,0x28(%rsp)
  cmp    $0x400,%eax
  jne    L107
  mov    0x60058(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L150
  add    $0x60068,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L150:
  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;
}
```