user: Demonstrate the C code that generates the assembly output below:

```asm
<pool_forwards_cpu>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x70,%rsp
  mov    %rdi,-0x58(%rbp)
  mov    %esi,-0x5c(%rbp)
  mov    %edx,-0x60(%rbp)
  mov    %ecx,-0x64(%rbp)
  mov    %r8d,-0x68(%rbp)
  mov    %r9d,-0x6c(%rbp)
  mov    -0x60(%rbp),%eax
  imul   -0x64(%rbp),%eax
  mov    %eax,-0x38(%rbp)
  movl   $0x0,-0x50(%rbp)
  jmp    L94
L49:
  movl   $0x0,-0x4c(%rbp)
  jmp    L55
L46:
  movl   $0x0,-0x48(%rbp)
  jmp    L119
L98:
  mov    -0x4c(%rbp),%eax
  imul   -0x64(%rbp),%eax
  mov    %eax,%edx
  mov    -0x48(%rbp),%eax
  add    %eax,%edx
  mov    -0x38(%rbp),%eax
  imul   -0x50(%rbp),%eax
  add    %edx,%eax
  mov    %eax,-0x34(%rbp)
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movl   $0x1,-0x44(%rbp)
  jmp    L91
L70:
  mov    -0x4c(%rbp),%eax
  mov    %eax,-0x24(%rbp)
  mov    -0x48(%rbp),%edx
  mov    -0x44(%rbp),%eax
  add    %edx,%eax
  mov    %eax,-0x20(%rbp)
  mov    -0x24(%rbp),%eax
  imul   -0x64(%rbp),%eax
  mov    %eax,%edx
  mov    -0x20(%rbp),%eax
  add    %eax,%edx
  mov    -0x38(%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
  mov    -0x18(%rbp),%rax
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  call   <max@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x18(%rbp)
  addl   $0x1,-0x44(%rbp)
L91:
  mov    -0x44(%rbp),%eax
  cmp    -0x6c(%rbp),%eax
  jl     L70
  movl   $0x1,-0x40(%rbp)
  jmp    L12
L44:
  movl   $0x0,-0x3c(%rbp)
  jmp    L25
L72:
  mov    -0x4c(%rbp),%edx
  mov    -0x40(%rbp),%eax
  add    %edx,%eax
  mov    %eax,-0x30(%rbp)
  mov    -0x48(%rbp),%edx
  mov    -0x3c(%rbp),%eax
  add    %edx,%eax
  mov    %eax,-0x2c(%rbp)
  mov    -0x30(%rbp),%eax
  imul   -0x64(%rbp),%eax
  mov    %eax,%edx
  mov    -0x2c(%rbp),%eax
  add    %eax,%edx
  mov    -0x38(%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
  mov    -0x18(%rbp),%rax
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  call   <max@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x18(%rbp)
  addl   $0x1,-0x3c(%rbp)
L25:
  mov    -0x3c(%rbp),%eax
  cmp    -0x6c(%rbp),%eax
  jl     L72
  addl   $0x1,-0x40(%rbp)
L12:
  mov    -0x40(%rbp),%eax
  cmp    -0x68(%rbp),%eax
  jl     L44
  mov    0x20(%rbp),%rax
  lea    0x8(%rax),%rdx
  mov    %rdx,0x20(%rbp)
  movsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    0x18(%rbp),%eax
  add    %eax,-0x48(%rbp)
L119:
  mov    -0x64(%rbp),%eax
  sub    -0x6c(%rbp),%eax
  cmp    %eax,-0x48(%rbp)
  jle    L98
  mov    0x10(%rbp),%eax
  add    %eax,-0x4c(%rbp)
L55:
  mov    -0x60(%rbp),%eax
  sub    -0x68(%rbp),%eax
  cmp    %eax,-0x4c(%rbp)
  jle    L46
  addl   $0x1,-0x50(%rbp)
L94:
  mov    -0x50(%rbp),%eax
  cmp    -0x5c(%rbp),%eax
  jl     L49
  nop
  nop
  leave
  ret
```

assistant: ```c
void pool_forwards_cpu(const double *data_im, 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_pooled) {
  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];
          max_value = max(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];
            max_value = max(max_value, data_value);
          }
        }
        *(data_pooled++) = max_value;
      }
    }
  }
}
```