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

```asm
<convolutionGold>:
  endbr64
  push   %r15
  push   %r14
  mov    %esi,%r14d
  push   %r13
  movslq %r8d,%r13
  push   %r12
  mov    %rcx,%r12
  push   %rbp
  mov    %rdi,%rbp
  mov    %esi,%edi
  push   %rbx
  imul   %edx,%edi
  mov    %edx,%ebx
  sub    $0x48,%rsp
  shl    $0x2,%edi
  call   <malloc@plt>
  mov    %rax,0x28(%rsp)
  test   %rax,%rax
  je     L10
L26:
  mov    %r13d,%eax
  shr    $0x1f,%eax
  add    %r13d,%eax
  mov    %eax,%edi
  sar    %edi
  mov    %edi,0x14(%rsp)
  mov    %edi,%ecx
  test   %r14d,%r14d
  jle    L33
  lea    0x1(%rdi),%esi
  neg    %edi
  movslq %ecx,%rdx
  movl   $0x0,0x30(%rsp)
  lea    (%rsi,%r14,1),%eax
  mov    %esi,0x34(%rsp)
  shl    $0x2,%r13
  mov    %esi,%r15d
  mov    %eax,0x38(%rsp)
  movslq %edi,%rax
  lea    -0x1(%rbx),%r9d
  add    %rdx,%rax
  mov    %edi,0x10(%rsp)
  lea    (%r12,%rax,4),%rax
  mov    %rax,0x18(%rsp)
  lea    (%rbx,%rsi,1),%eax
  mov    %eax,0x20(%rsp)
  lea    -0x1(%rdi,%rdi,1),%eax
  mov    %eax,0x3c(%rsp)
  lea    -0x1(%r14),%eax
  mov    %eax,0x4(%rsp)
L29:
  test   %ebx,%ebx
  jle    L40
  movslq 0x30(%rsp),%rax
  mov    0x28(%rsp),%rdi
  mov    0x10(%rsp),%r12d
  mov    0x34(%rsp),%r8d
  lea    (%rdi,%rax,4),%rax
  mov    %rax,0x8(%rsp)
  mov    0x3c(%rsp),%eax
  add    %r15d,%eax
  mov    %eax,0x24(%rsp)
  nopl   0x0(%rax,%rax,1)
L36:
  mov    0x10(%rsp),%edi
  cmp    %edi,0x14(%rsp)
  jl     L34
  mov    0x18(%rsp),%r11
  mov    0x24(%rsp),%r10d
  xor    %esi,%esi
  nopw   0x0(%rax,%rax,1)
L50:
  xor    %edi,%edi
  test   %r10d,%r10d
  mov    %r11,%rcx
  mov    %r12d,%edx
  cmovns %r10d,%edi
  cmp    %edi,%r14d
  cmovle 0x4(%rsp),%edi
  imul   %ebx,%edi
  nopw   0x0(%rax,%rax,1)
L55:
  xor    %eax,%eax
  pxor   %xmm0,%xmm0
  pxor   %xmm1,%xmm1
  test   %edx,%edx
  cmovns %edx,%eax
  cvtsi2ss %esi,%xmm1
  cmp    %eax,%ebx
  cmovle %r9d,%eax
  add    $0x1,%edx
  add    $0x4,%rcx
  add    %edi,%eax
  cltq
  cvtsi2ssl 0x0(%rbp,%rax,4),%xmm0
  mulss  -0x4(%rcx),%xmm0
  addss  %xmm1,%xmm0
  cvttss2si %xmm0,%esi
  cmp    %r8d,%edx
  jne    L55
  add    $0x1,%r10d
  add    %r13,%r11
  cmp    %r15d,%r10d
  jne    L50
L69:
  mov    0x8(%rsp),%rax
  add    $0x1,%r8d
  add    $0x1,%r12d
  mov    %esi,(%rax)
  add    $0x4,%rax
  mov    %rax,0x8(%rsp)
  cmp    0x20(%rsp),%r8d
  jne    L36
L40:
  add    %ebx,0x30(%rsp)
  add    $0x1,%r15d
  cmp    0x38(%rsp),%r15d
  jne    L29
L33:
  mov    0x28(%rsp),%rax
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L34:
  xor    %esi,%esi
  jmp    L69
L10:
  mov    $0xffffffff,%edi
  call   <exit@plt>
  jmp    L26
```

assistant: ```c
int *convolutionGold(int *image, int rows, int cols, float *filter, int filterWidth) {
  int *outputImage;
  outputImage = (int *)malloc(rows * cols * sizeof(int));
  if (!outputImage) {
    exit(-1);
  }
  int halfFilterWidth = (int)filterWidth / 2;
  int i;
  int j;
  int k;
  int l;
  for (i = 0; i < rows; i++) {
    for (j = 0; j < cols; j++) {
      int sum = 0;
      for (k = -halfFilterWidth; k <= halfFilterWidth; k++) {
        for (l = -halfFilterWidth; l <= halfFilterWidth; l++) {
          int r = i + k;
          int c = j + l;
          r = (r < 0) ? 0 : r;
          c = (c < 0) ? 0 : c;
          r = (r >= rows) ? rows - 1 : r;
          c = (c >= cols) ? cols - 1 : c;
          sum += image[r * cols + c] * filter[(k + halfFilterWidth) * filterWidth + (l + halfFilterWidth)];
        }
      }
      outputImage[i * cols + j] = sum;
    }
  }
  return outputImage;
}
```