user: Please show me the C implementation that matches the assembly code listed below:

```asm
<convolutionGoldFloat>:
  endbr64
  push   %r15
  mov    %esi,%r15d
  push   %r14
  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     L49
L28:
  mov    %r13d,%eax
  shr    $0x1f,%eax
  add    %r13d,%eax
  sar    %eax
  mov    %eax,0x18(%rsp)
  mov    %eax,%edi
  test   %r15d,%r15d
  jle    L59
  mov    %eax,%esi
  neg    %esi
  mov    %esi,0x14(%rsp)
  mov    %esi,%eax
  test   %ebx,%ebx
  jle    L59
  movslq %edi,%rdx
  cltq
  shl    $0x2,%r13
  pxor   %xmm2,%xmm2
  lea    0x1(%rdx),%ecx
  movl   $0x0,0x34(%rsp)
  lea    -0x1(%rbx),%r14d
  lea    (%rcx,%r15,1),%edi
  mov    %ecx,0x38(%rsp)
  mov    %edi,0x3c(%rsp)
  mov    %rax,%rdi
  add    %rdx,%rax
  lea    (%r12,%rax,4),%rax
  mov    %ecx,%r12d
  mov    %rax,0x20(%rsp)
  lea    (%rbx,%rcx,1),%eax
  mov    %eax,0x1c(%rsp)
  lea    -0x1(%r15),%eax
  mov    %eax,0x10(%rsp)
  mov    %edi,%eax
  sub    %edx,%eax
  sub    $0x1,%eax
  mov    %eax,0x30(%rsp)
L76:
  mov    0x28(%rsp),%rdi
  movslq 0x34(%rsp),%rax
  mov    0x14(%rsp),%r11d
  lea    (%rdi,%rax,4),%rax
  mov    0x38(%rsp),%edi
  mov    %rax,0x8(%rsp)
  nopw   0x0(%rax,%rax,1)
L64:
  mov    0x14(%rsp),%esi
  cmp    %esi,0x18(%rsp)
  jl     L40
  mov    0x30(%rsp),%eax
  mov    0x20(%rsp),%r10
  movaps %xmm2,%xmm0
  lea    (%rax,%r12,1),%r9d
  xchg   %ax,%ax
L34:
  xor    %esi,%esi
  test   %r9d,%r9d
  mov    %r10,%rcx
  mov    %r11d,%edx
  cmovns %r9d,%esi
  cmp    %esi,%r15d
  cmovle 0x10(%rsp),%esi
  imul   %ebx,%esi
  lea    (%r14,%rsi,1),%r8d
  movslq %r8d,%r8
  jmp    L63
  nopl   0x0(%rax,%rax,1)
L110:
  movss  0x0(%rbp,%r8,4),%xmm1
L41:
  mulss  (%rcx),%xmm1
  add    $0x1,%edx
  add    $0x4,%rcx
  addss  %xmm1,%xmm0
  cmp    %edi,%edx
  je     L100
L63:
  xor    %eax,%eax
  test   %edx,%edx
  cmovns %edx,%eax
  cmp    %eax,%ebx
  jle    L110
  add    %esi,%eax
  cltq
  movss  0x0(%rbp,%rax,4),%xmm1
  jmp    L41
  nopl   0x0(%rax)
L100:
  add    $0x1,%r9d
  add    %r13,%r10
  cmp    %r12d,%r9d
  jne    L34
L102:
  mov    0x8(%rsp),%rax
  add    $0x1,%edi
  add    $0x1,%r11d
  movss  %xmm0,(%rax)
  add    $0x4,%rax
  mov    %rax,0x8(%rsp)
  cmp    0x1c(%rsp),%edi
  jne    L64
  add    %ebx,0x34(%rsp)
  add    $0x1,%r12d
  cmp    0x3c(%rsp),%r12d
  jne    L76
L59:
  mov    0x28(%rsp),%rax
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L40:
  movaps %xmm2,%xmm0
  jmp    L102
L49:
  mov    $0xffffffff,%edi
  call   <exit@plt>
  jmp    L28
```

assistant: ```c
float *convolutionGoldFloat(float *image, int rows, int cols, float *filter, int filterWidth) {
  float *outputImage;
  outputImage = (float *)malloc(rows * cols * sizeof(float));
  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++) {
      float 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;
}
```