user: What would the C source look like for the assembly code provided here?

```asm
<separable_convolution>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %ecx,-0x4(%rsp)
  mov    %r9d,-0x1c(%rsp)
  test   %ecx,%ecx
  jle    L21
  mov    %esi,%r14d
  mov    %edx,%esi
  mov    %r8,%rbp
  mov    %r14d,%eax
  imul   %edx,%eax
  mov    %eax,-0x8(%rsp)
  movslq %edx,%rax
  shl    $0x2,%rax
  mov    %rax,-0x28(%rsp)
  movl   $0x0,-0xc(%rsp)
  movl   $0x0,-0x10(%rsp)
  mov    %r9d,%r15d
  jmp    L81
L82:
  mov    -0x18(%rsp),%r12
  mov    %r14d,%r11d
  movl   $0x0,-0x20(%rsp)
  mov    $0x0,%r10d
L69:
  mov    -0x20(%rsp),%r9d
  pxor   %xmm0,%xmm0
  cvtsi2ssl (%rdi,%r12,1),%xmm0
  mulss  0x0(%rbp),%xmm0
  cmpl   $0x1,-0x1c(%rsp)
  jle    L0
  movaps %xmm0,%xmm2
  mov    $0x1,%ecx
L26:
  movss  0x0(%rbp,%rcx,4),%xmm3
  lea    (%r8,%rcx,1),%eax
  cltd
  idiv   %esi
  add    %r9d,%edx
  add    %r13d,%edx
  movslq %edx,%rdx
  pxor   %xmm1,%xmm1
  cvtsi2ssl (%rdi,%rdx,4),%xmm1
  mulss  %xmm3,%xmm1
  addss  %xmm1,%xmm2
  mov    %ebx,%eax
  sub    %ecx,%eax
  cltd
  idiv   %esi
  add    %r9d,%edx
  add    %r13d,%edx
  movslq %edx,%rdx
  pxor   %xmm1,%xmm1
  cvtsi2ssl (%rdi,%rdx,4),%xmm1
  mulss  %xmm3,%xmm1
  addss  %xmm1,%xmm2
  lea    (%r10,%rcx,1),%eax
  cltd
  idiv   %r14d
  imul   %esi,%edx
  add    %r8d,%edx
  add    %r13d,%edx
  movslq %edx,%rdx
  pxor   %xmm1,%xmm1
  cvtsi2ssl (%rdi,%rdx,4),%xmm1
  mulss  %xmm3,%xmm1
  addss  %xmm1,%xmm0
  mov    %r11d,%eax
  sub    %ecx,%eax
  cltd
  idiv   %r14d
  imul   %esi,%edx
  add    %r8d,%edx
  add    %r13d,%edx
  movslq %edx,%rdx
  pxor   %xmm1,%xmm1
  cvtsi2ssl (%rdi,%rdx,4),%xmm1
  mulss  %xmm3,%xmm1
  addss  %xmm1,%xmm0
  add    $0x1,%rcx
  cmp    %r15,%rcx
  jne    L26
L109:
  mov    0x38(%rsp),%rax
  movss  %xmm0,(%rax,%r12,1)
  mov    0x40(%rsp),%rax
  movss  %xmm2,(%rax,%r12,1)
  add    $0x1,%r10d
  add    %esi,-0x20(%rsp)
  add    $0x1,%r11d
  add    -0x28(%rsp),%r12
  cmp    %r10d,%r14d
  jne    L69
L14:
  add    $0x1,%r8d
  addq   $0x4,-0x18(%rsp)
  add    $0x1,%ebx
  cmp    %r8d,%esi
  je     L115
L86:
  test   %r14d,%r14d
  jg     L82
  jmp    L14
L0:
  movaps %xmm0,%xmm2
  jmp    L109
L115:
  addl   $0x1,-0x10(%rsp)
  mov    -0x10(%rsp),%eax
  mov    -0x8(%rsp),%ecx
  add    %ecx,-0xc(%rsp)
  cmp    %eax,-0x4(%rsp)
  je     L21
L81:
  test   %esi,%esi
  jle    L115
  mov    -0xc(%rsp),%eax
  mov    %eax,%r13d
  cltq
  shl    $0x2,%rax
  mov    %rax,-0x18(%rsp)
  mov    %esi,%ebx
  mov    $0x0,%r8d
  jmp    L86
L21:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void separable_convolution(int *image, int width, int height, int depth, float *gau, int gau_width, float *smx, float *smy) {
  int i, j, k, d, I1, I2;
  float x, y;
  for (d = 0; d < depth; d++) {
    for (i = 0; i < height; i++) {
      for (j = 0; j < width; j++) {
        x = gau[0] * image[((height) * (j) + (i) + (width) * (height) * (d))];
        y = gau[0] * image[((height) * (j) + (i) + (width) * (height) * (d))];
        for (k = 1; k < gau_width; k++) {
          I1 = (i + k) % height;
          I2 = (i - k + height) % height;
          y = y + gau[k] * image[((height) * (j) + (I1) + (width) * (height) * (d))] + gau[k] * image[((height) * (j) + (I2) + (width) * (height) * (d))];
          I1 = (j + k) % width;
          I2 = (j - k + width) % width;
          x = x + gau[k] * image[((height) * (I1) + (i) + (width) * (height) * (d))] + gau[k] * image[((height) * (I2) + (i) + (width) * (height) * (d))];
        }
        smx[((height) * (j) + (i) + (width) * (height) * (d))] = x;
        smy[((height) * (j) + (i) + (width) * (height) * (d))] = y;
      }
    }
  }
}
```