user: Converte this assembly snippet into its equivalent C code:

```asm
<separable_convolution>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %ecx,-0x4(%rsp)
  mov    %r9d,-0x10(%rsp)
  test   %ecx,%ecx
  jle    L55
  mov    %rdi,%r10
  mov    %edx,%eax
  mov    %esi,%edi
  mov    %edx,%esi
  imul   %edi,%eax
  mov    %eax,-0x8(%rsp)
  test   %edx,%edx
  jle    L55
  test   %edi,%edi
  jle    L55
  movslq %edx,%rax
  movl   $0x0,-0xc(%rsp)
  mov    %r8,%r12
  movslq %r9d,%r15
  shl    $0x2,%rax
  xor    %r8d,%r8d
  mov    %rax,-0x20(%rsp)
L4:
  movslq %r8d,%rax
  mov    %esi,%r13d
  xor    %r9d,%r9d
  shl    $0x2,%rax
  mov    %rax,-0x18(%rsp)
  xchg   %ax,%ax
L72:
  mov    -0x18(%rsp),%r14
  mov    %edi,%ebp
  xor    %r11d,%r11d
  xor    %ebx,%ebx
  nopl   0x0(%rax)
L35:
  pxor   %xmm0,%xmm0
  cmpl   $0x1,-0x10(%rsp)
  cvtsi2ssl (%r10,%r14,1),%xmm0
  mulss  (%r12),%xmm0
  movaps %xmm0,%xmm2
  jle    L13
  mov    $0x1,%ecx
  nopl   0x0(%rax,%rax,1)
L3:
  lea    (%r9,%rcx,1),%eax
  pxor   %xmm1,%xmm1
  movss  (%r12,%rcx,4),%xmm3
  cltd
  idiv   %esi
  mov    %r13d,%eax
  sub    %ecx,%eax
  add    %r11d,%edx
  add    %r8d,%edx
  movslq %edx,%rdx
  cvtsi2ssl (%r10,%rdx,4),%xmm1
  cltd
  idiv   %esi
  mulss  %xmm3,%xmm1
  lea    (%rbx,%rcx,1),%eax
  addss  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  add    %r11d,%edx
  add    %r8d,%edx
  movslq %edx,%rdx
  cvtsi2ssl (%r10,%rdx,4),%xmm1
  cltd
  idiv   %edi
  mulss  %xmm3,%xmm1
  mov    %ebp,%eax
  sub    %ecx,%eax
  add    $0x1,%rcx
  addss  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  imul   %esi,%edx
  add    %r9d,%edx
  add    %r8d,%edx
  movslq %edx,%rdx
  cvtsi2ssl (%r10,%rdx,4),%xmm1
  cltd
  idiv   %edi
  mulss  %xmm3,%xmm1
  addss  %xmm1,%xmm2
  pxor   %xmm1,%xmm1
  imul   %esi,%edx
  add    %r9d,%edx
  add    %r8d,%edx
  movslq %edx,%rdx
  cvtsi2ssl (%r10,%rdx,4),%xmm1
  mulss  %xmm3,%xmm1
  addss  %xmm1,%xmm2
  cmp    %rcx,%r15
  jne    L3
L13:
  mov    0x38(%rsp),%rax
  add    $0x1,%ebx
  add    %esi,%r11d
  add    $0x1,%ebp
  movss  %xmm2,(%rax,%r14,1)
  mov    0x40(%rsp),%rax
  movss  %xmm0,(%rax,%r14,1)
  add    -0x20(%rsp),%r14
  cmp    %ebx,%edi
  jne    L35
  add    $0x1,%r9d
  addq   $0x4,-0x18(%rsp)
  add    $0x1,%r13d
  cmp    %r9d,%esi
  jne    L72
  addl   $0x1,-0xc(%rsp)
  add    -0x8(%rsp),%r8d
  mov    -0xc(%rsp),%eax
  cmp    %eax,-0x4(%rsp)
  jne    L4
L55:
  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;
      }
    }
  }
}
```