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

```asm
<pool_backwards_cpu>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x70,%rsp
  mov    %rdi,-0x58(%rbp)
  mov    %rsi,-0x60(%rbp)
  mov    %edx,-0x64(%rbp)
  mov    %ecx,-0x68(%rbp)
  mov    %r8d,-0x6c(%rbp)
  mov    %r9d,-0x70(%rbp)
  mov    -0x68(%rbp),%eax
  imul   -0x6c(%rbp),%eax
  imul   -0x64(%rbp),%eax
  cltq
  shl    $0x3,%eax
  mov    %eax,%edx
  mov    0x28(%rbp),%rax
  mov    $0x0,%esi
  mov    %rax,%rdi
  call   <memset@plt>
  mov    -0x68(%rbp),%eax
  imul   -0x6c(%rbp),%eax
  mov    %eax,-0x34(%rbp)
  movl   $0x0,-0x50(%rbp)
  jmp    L81
L55:
  movl   $0x0,-0x4c(%rbp)
  jmp    L56
L25:
  movl   $0x0,-0x48(%rbp)
  jmp    L14
L34:
  mov    -0x4c(%rbp),%eax
  imul   -0x6c(%rbp),%eax
  mov    %eax,%edx
  mov    -0x48(%rbp),%eax
  add    %eax,%edx
  mov    -0x34(%rbp),%eax
  imul   -0x50(%rbp),%eax
  add    %edx,%eax
  mov    %eax,-0x44(%rbp)
  mov    -0x44(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movl   $0x1,-0x40(%rbp)
  jmp    L44
L123:
  mov    -0x4c(%rbp),%eax
  mov    %eax,-0x24(%rbp)
  mov    -0x48(%rbp),%edx
  mov    -0x40(%rbp),%eax
  add    %edx,%eax
  mov    %eax,-0x20(%rbp)
  mov    -0x24(%rbp),%eax
  imul   -0x6c(%rbp),%eax
  mov    %eax,%edx
  mov    -0x20(%rbp),%eax
  add    %eax,%edx
  mov    -0x34(%rbp),%eax
  imul   -0x50(%rbp),%eax
  add    %edx,%eax
  mov    %eax,-0x1c(%rbp)
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x8(%rbp),%xmm0
  comisd -0x18(%rbp),%xmm0
  jbe    L134
  mov    -0x1c(%rbp),%eax
  mov    %eax,-0x44(%rbp)
  movsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,-0x18(%rbp)
L134:
  addl   $0x1,-0x40(%rbp)
L44:
  mov    -0x40(%rbp),%eax
  cmp    0x10(%rbp),%eax
  jl     L123
  movl   $0x1,-0x3c(%rbp)
  jmp    L107
L105:
  movl   $0x0,-0x38(%rbp)
  jmp    L5
L108:
  mov    -0x4c(%rbp),%edx
  mov    -0x3c(%rbp),%eax
  add    %edx,%eax
  mov    %eax,-0x30(%rbp)
  mov    -0x48(%rbp),%edx
  mov    -0x38(%rbp),%eax
  add    %edx,%eax
  mov    %eax,-0x2c(%rbp)
  mov    -0x30(%rbp),%eax
  imul   -0x6c(%rbp),%eax
  mov    %eax,%edx
  mov    -0x2c(%rbp),%eax
  add    %eax,%edx
  mov    -0x34(%rbp),%eax
  imul   -0x50(%rbp),%eax
  add    %edx,%eax
  mov    %eax,-0x28(%rbp)
  mov    -0x28(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x10(%rbp),%xmm0
  comisd -0x18(%rbp),%xmm0
  jbe    L28
  mov    -0x28(%rbp),%eax
  mov    %eax,-0x44(%rbp)
  movsd  -0x10(%rbp),%xmm0
  movsd  %xmm0,-0x18(%rbp)
L28:
  addl   $0x1,-0x38(%rbp)
L5:
  mov    -0x38(%rbp),%eax
  cmp    0x10(%rbp),%eax
  jl     L108
  addl   $0x1,-0x3c(%rbp)
L107:
  mov    -0x3c(%rbp),%eax
  cmp    -0x70(%rbp),%eax
  jl     L105
  mov    -0x60(%rbp),%rax
  lea    0x8(%rax),%rdx
  mov    %rdx,-0x60(%rbp)
  movsd  (%rax),%xmm1
  mov    -0x44(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x44(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x28(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    0x20(%rbp),%eax
  add    %eax,-0x48(%rbp)
L14:
  mov    -0x6c(%rbp),%eax
  sub    0x10(%rbp),%eax
  cmp    %eax,-0x48(%rbp)
  jle    L34
  mov    0x18(%rbp),%eax
  add    %eax,-0x4c(%rbp)
L56:
  mov    -0x68(%rbp),%eax
  sub    -0x70(%rbp),%eax
  cmp    %eax,-0x4c(%rbp)
  jle    L25
  addl   $0x1,-0x50(%rbp)
L81:
  mov    -0x50(%rbp),%eax
  cmp    -0x64(%rbp),%eax
  jl     L55
  nop
  nop
  leave
  ret
```

assistant: ```c
void pool_backwards_cpu(const double *data_im, const double *data_pooled, const int channels, const int height, const int width, const int kernel_h, const int kernel_w, const int stride_h, const int stride_w, double *data_backgrad) {
  memset(data_backgrad, 0, height * width * channels * sizeof(double));
  const int channel_size = height * width;
  for (int channel = 0; channel < channels; channel++) {
    for (int fitting_height = 0; fitting_height <= (height - kernel_h); fitting_height += stride_h) {
      for (int fitting_width = 0; fitting_width <= (width - kernel_w); fitting_width += stride_w) {
        int max_index = fitting_height * width + fitting_width + channel_size * channel;
        double max_value = data_im[max_index];
        for (int kernel_col = 1; kernel_col < kernel_w; kernel_col++) {
          int input_row = fitting_height;
          int input_col = fitting_width + kernel_col;
          int data_index = input_row * width + input_col + channel_size * channel;
          double data_value = data_im[data_index];
          if (data_value > max_value) {
            max_index = data_index;
            max_value = data_value;
          }
        }
        for (int kernel_row = 1; kernel_row < kernel_h; kernel_row++) {
          for (int kernel_col = 0; kernel_col < kernel_w; kernel_col++) {
            int input_row = fitting_height + kernel_row;
            int input_col = fitting_width + kernel_col;
            int data_index = input_row * width + input_col + channel_size * channel;
            double data_value = data_im[data_index];
            if (data_value > max_value) {
              max_index = data_index;
              max_value = data_value;
            }
          }
        }
        data_backgrad[max_index] += *(data_pooled++);
      }
    }
  }
}
```