user: 请从给定的汇编指令推导出C源代码:

```asm
<aan_mascara_canal>:
  endbr64
  push   %r15
  mov    %rdi,%r15
  mov    $0x18,%edi
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  mov    %ecx,%ebp
  push   %rbx
  mov    %edx,%ebx
  sub    $0xd8,%rsp
  mov    %rsi,0xa0(%rsp)
  mov    %edx,0x60(%rsp)
  mov    %ecx,0x9c(%rsp)
  mov    %r8,0x18(%rsp)
  call   <malloc@plt>
  mov    $0x24,%edi
  mov    %rax,%r14
  call   <malloc@plt>
  mov    %rax,%r13
  add    $0xc,%rax
  movq   %rax,%xmm1
  movq   %r13,%xmm0
  lea    0x18(%r13),%rax
  punpcklqdq %xmm1,%xmm0
  mov    %rax,0x10(%r14)
  movups %xmm0,(%r14)
  test   %ebp,%ebp
  jle    L3
  mov    %ebx,%eax
  movl   $0x0,0x6c(%rsp)
  imul   %ebp,%eax
  movl   $0x0,0x64(%rsp)
  movl   $0x0,0x4(%rsp)
  cltq
  shl    $0x2,%rax
  lea    -0x4(%r15,%rax,1),%rdx
  lea    -0x8(%r15,%rax,1),%rax
  mov    %rax,0xc8(%rsp)
  movslq %ebx,%rax
  shl    $0x2,%rax
  mov    %rdx,0xc0(%rsp)
  lea    -0x4(%r15,%rax,1),%rdx
  mov    %rdx,0xa8(%rsp)
  lea    -0x8(%r15,%rax,1),%rdx
  add    %r15,%rax
  mov    %rax,0x30(%rsp)
  lea    0x1(%rbx),%eax
  cltq
  mov    %rdx,0xb0(%rsp)
  lea    (%rbx,%rbx,1),%edx
  lea    (%r15,%rax,4),%rax
  mov    %edx,0x20(%rsp)
  movslq %edx,%rdx
  mov    %rax,0x38(%rsp)
  mov    %ebx,%eax
  lea    -0x4(%r15,%rdx,4),%rdx
  neg    %eax
  mov    %rdx,0xb8(%rsp)
  mov    %eax,0x68(%rsp)
  lea    -0x1(%rbx),%eax
  mov    %eax,0xc(%rsp)
  mov    %eax,0x98(%rsp)
  lea    -0x1(%rbp),%eax
  mov    %eax,0x24(%rsp)
  nopl   0x0(%rax)
L216:
  mov    0x60(%rsp),%eax
  test   %eax,%eax
  jle    L111
  mov    0x68(%rsp),%ebp
  movslq 0x6c(%rsp),%rdx
  mov    0x98(%rsp),%r8d
  movslq %ebp,%rax
  mov    %rdx,%rsi
  shl    $0x2,%rdx
  mov    %rax,%rbx
  shl    $0x2,%rax
  lea    (%r15,%rax,1),%rcx
  lea    0x4(%r15,%rax,1),%rax
  mov    %rcx,0x40(%rsp)
  lea    (%r15,%rdx,1),%rcx
  lea    -0x4(%r15,%rdx,1),%rdx
  mov    %rcx,0x48(%rsp)
  lea    0x1(%r8),%ecx
  movslq %ecx,%rcx
  mov    %rdx,0x90(%rsp)
  mov    %ebp,%edx
  lea    (%r15,%rcx,4),%rcx
  mov    %rax,0x88(%rsp)
  movslq 0x64(%rsp),%rax
  mov    %rcx,0x70(%rsp)
  lea    0x1(%rbx),%ecx
  movslq %ecx,%rcx
  lea    (%r15,%rcx,4),%rcx
  mov    %rcx,0x50(%rsp)
  lea    0x1(%rsi),%ecx
  movslq %ecx,%rcx
  lea    (%r15,%rcx,4),%rcx
  mov    %rcx,0x58(%rsp)
  lea    0x2(%r8),%ecx
  movslq %ecx,%rcx
  lea    (%r15,%rcx,4),%rcx
  mov    %rcx,0x78(%rsp)
  mov    0x60(%rsp),%ecx
  sub    %ecx,%edx
  movslq %edx,%rdx
  lea    (%r15,%rdx,4),%rbx
  mov    %rbx,0x80(%rsp)
  mov    0xa0(%rsp),%rbx
  lea    (%rbx,%rax,4),%rax
  xor    %ebx,%ebx
  mov    %rax,0x10(%rsp)
  jmp    L210
  nopl   0x0(%rax)
L218:
  mov    %ebx,%edi
  cmp    %ebx,0xc(%rsp)
  jne    L122
  mov    0x4(%rsp),%r9d
  test   %r9d,%r9d
  je     L156
L122:
  test   %rbx,%rbx
  je     L10
  cmp    %edi,0xc(%rsp)
  je     L154
  mov    0x4(%rsp),%edi
  test   %edi,%edi
  je     L45
  lea    -0x1(%rbp),%edi
  mov    0x4(%rsp),%r10d
  movslq %edi,%rdi
  lea    (%r15,%rdi,4),%r12
  movslq %ebp,%rdi
  add    $0x1,%ebp
  lea    (%r15,%rdi,4),%r11
  movslq %ebp,%rdi
  lea    (%r15,%rdi,4),%rdi
  mov    %rdi,0x28(%rsp)
  lea    -0x1(%rsi),%edi
  movslq %edi,%rdi
  lea    (%r15,%rdi,4),%r9
  movslq %esi,%rdi
  add    $0x1,%esi
  mov    %esi,0x8(%rsp)
  movslq %esi,%rsi
  lea    (%r15,%rdi,4),%rdi
  lea    (%r15,%rsi,4),%rsi
  cmp    %r10d,0x24(%rsp)
  je     L86
  movss  (%r12),%xmm0
  lea    0x1(%r8),%r12d
  movss  %xmm0,0x0(%r13)
  movss  (%r11),%xmm0
  mov    0x28(%rsp),%r11
  movss  %xmm0,0x4(%r13)
  movss  (%r11),%xmm0
  movss  %xmm0,0x8(%r13)
  movss  (%r9),%xmm0
  lea    0x1(%rcx),%r13d
  movss  %xmm0,(%rax)
  movss  (%rdi),%xmm0
  movss  %xmm0,0x4(%rax)
  movss  (%rsi),%xmm0
  movss  %xmm0,0x8(%rax)
  movslq %r8d,%rax
  movss  (%r15,%rax,4),%xmm0
  movslq %r12d,%rax
  movss  %xmm0,(%rdx)
  movss  (%r15,%rax,4),%xmm0
  lea    0x2(%r8),%eax
  cltq
  movss  %xmm0,0x4(%rdx)
  movss  (%r15,%rax,4),%xmm0
  movss  %xmm0,0x8(%rdx)
L23:
  mov    0x18(%rsp),%rdi
  mov    %r14,%rsi
  xor    %eax,%eax
  call   <pixel_resultante@plt>
  pxor   %xmm0,%xmm0
  mov    %r13d,%ecx
  mov    0x8(%rsp),%esi
  cvtsi2ss %eax,%xmm0
  mov    0x10(%rsp),%rax
  mov    (%r14),%r13
  mov    %r12d,%r8d
  movss  %xmm0,(%rax,%rbx,4)
  cmp    0x20(%rsp),%ecx
  je     L111
  add    $0x1,%rbx
L210:
  mov    0x4(%rsp),%edi
  mov    0x8(%r14),%rax
  mov    0x10(%r14),%rdx
  or     %ebx,%edi
  jne    L218
  movq   (%r15),%xmm0
  mov    0x30(%rsp),%rdi
  shufps $0xe0,%xmm0,%xmm0
  movlps %xmm0,0x0(%r13)
  movss  0x4(%r15),%xmm0
  movss  %xmm0,0x8(%r13)
  movss  (%r15),%xmm0
  movss  %xmm0,(%rax)
  movss  (%rdi),%xmm0
  movss  %xmm0,(%rdx)
  movss  (%r15),%xmm0
  movss  %xmm0,0x4(%rax)
  movss  0x4(%r15),%xmm0
  movss  %xmm0,0x8(%rax)
  mov    0x38(%rsp),%rax
  movss  (%rdi),%xmm0
  movss  %xmm0,0x4(%rdx)
  movss  (%rax),%xmm0
  movss  %xmm0,0x8(%rdx)
L47:
  lea    0x1(%rsi),%eax
  lea    0x1(%r8),%r12d
  add    $0x1,%ebp
  mov    %eax,0x8(%rsp)
  lea    0x1(%rcx),%r13d
  jmp    L23
  nopl   (%rax)
L10:
  mov    0x40(%rsp),%rdi
  mov    0x4(%rsp),%r11d
  movss  (%rdi),%xmm0
  cmp    %r11d,0x24(%rsp)
  je     L134
  mov    0x48(%rsp),%r9
  movss  %xmm0,0x0(%r13)
  mov    0x70(%rsp),%r11
  movss  (%r9),%xmm0
  movss  %xmm0,(%rax)
  movss  (%r11),%xmm0
  movss  %xmm0,(%rdx)
  movss  (%rdi),%xmm0
  mov    0x50(%rsp),%rdi
  movss  %xmm0,0x4(%r13)
  movss  (%rdi),%xmm0
  mov    0x58(%rsp),%rdi
  movss  %xmm0,0x8(%r13)
  movss  (%r9),%xmm0
  movss  %xmm0,0x4(%rax)
  movss  (%rdi),%xmm0
  movss  %xmm0,0x8(%rax)
  mov    0x78(%rsp),%rax
  movss  (%r11),%xmm0
  movss  %xmm0,0x4(%rdx)
  movss  (%rax),%xmm0
  movss  %xmm0,0x8(%rdx)
  jmp    L47
  nopl   0x0(%rax)
L111:
  mov    0x60(%rsp),%eax
  addl   $0x1,0x4(%rsp)
  add    %eax,0x64(%rsp)
  mov    0x4(%rsp),%ecx
  add    %eax,0x68(%rsp)
  add    %eax,0x6c(%rsp)
  add    %eax,0x98(%rsp)
  cmp    %ecx,0x9c(%rsp)
  jne    L216
L3:
  mov    %r13,%rdi
  call   <free@plt>
  add    $0xd8,%rsp
  mov    %r14,%rdi
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    1070 <free@plt>
  xchg   %ax,%ax
L86:
  movss  (%r9),%xmm0
  movss  %xmm0,(%rdx)
  movss  (%rdi),%xmm0
  movss  %xmm0,0x4(%rdx)
  movss  (%rsi),%xmm0
  movss  %xmm0,0x8(%rdx)
  movss  (%r12),%xmm0
  mov    0x28(%rsp),%rdx
  lea    0x1(%r8),%r12d
  movss  %xmm0,0x0(%r13)
  movss  (%r11),%xmm0
  movss  %xmm0,0x4(%r13)
  movss  (%rdx),%xmm0
  movss  %xmm0,0x8(%r13)
  movss  (%r9),%xmm0
  lea    0x1(%rcx),%r13d
  movss  %xmm0,(%rax)
  movss  (%rdi),%xmm0
  movss  %xmm0,0x4(%rax)
  movss  (%rsi),%xmm0
  movss  %xmm0,0x8(%rax)
  jmp    L23
L154:
  mov    0x4(%rsp),%r10d
  cmp    %r10d,0x24(%rsp)
  je     L147
  movslq %ebp,%rdi
  lea    0x1(%r8),%r12d
  lea    (%r15,%rdi,4),%r11
  movslq %esi,%rdi
  movss  (%r11),%xmm0
  lea    (%r15,%rdi,4),%r10
  movslq %r12d,%rdi
  lea    (%r15,%rdi,4),%r9
  lea    -0x1(%rbp),%edi
  add    $0x1,%ebp
  movss  %xmm0,0x8(%r13)
  movss  (%r10),%xmm0
  movslq %edi,%rdi
  movss  %xmm0,0x8(%rax)
  movss  (%r9),%xmm0
  movss  %xmm0,0x8(%rdx)
  movss  (%r15,%rdi,4),%xmm0
  lea    -0x1(%rsi),%edi
  movslq %edi,%rdi
  movss  %xmm0,0x0(%r13)
  movss  (%r11),%xmm0
  movss  %xmm0,0x4(%r13)
  movss  (%r15,%rdi,4),%xmm0
  lea    0x1(%rcx),%r13d
  movss  %xmm0,(%rax)
  movss  (%r10),%xmm0
  movss  %xmm0,0x4(%rax)
  movslq %r8d,%rax
  movss  (%r15,%rax,4),%xmm0
  lea    0x1(%rsi),%eax
  mov    %eax,0x8(%rsp)
  movss  %xmm0,(%rdx)
  movss  (%r9),%xmm0
  movss  %xmm0,0x4(%rdx)
  jmp    L23
  nopw   0x0(%rax,%rax,1)
L45:
  mov    (%r14),%rdi
  lea    0x1(%rcx),%r13d
  lea    0x1(%r8),%r12d
  add    $0x1,%ebp
  movss  -0x4(%r15,%rbx,4),%xmm0
  movss  %xmm0,(%rdi)
  movss  (%r15,%rbx,4),%xmm0
  movss  %xmm0,0x4(%rdi)
  movss  0x4(%r15,%rbx,4),%xmm0
  movss  %xmm0,0x8(%rdi)
  movss  -0x4(%r15,%rbx,4),%xmm0
  movss  %xmm0,(%rax)
  movss  (%r15,%rbx,4),%xmm0
  movss  %xmm0,0x4(%rax)
  movss  0x4(%r15,%rbx,4),%xmm0
  movss  %xmm0,0x8(%rax)
  lea    -0x1(%rcx),%eax
  cltq
  movss  (%r15,%rax,4),%xmm0
  movslq %ecx,%rax
  movss  %xmm0,(%rdx)
  movss  (%r15,%rax,4),%xmm0
  movslq %r13d,%rax
  movss  %xmm0,0x4(%rdx)
  movss  (%r15,%rax,4),%xmm0
  lea    0x1(%rsi),%eax
  mov    %eax,0x8(%rsp)
  movss  %xmm0,0x8(%rdx)
  jmp    L23
  nopl   0x0(%rax)
L156:
  mov    0xa8(%rsp),%rdi
  mov    0xb0(%rsp),%r9
  movss  (%rdi),%xmm0
  movss  %xmm0,0x8(%r13)
  movss  %xmm0,0x4(%r13)
  movss  (%r9),%xmm0
  movss  %xmm0,0x0(%r13)
  movss  (%rdi),%xmm0
  mov    0xb8(%rsp),%rdi
  movss  %xmm0,0x8(%rax)
  movss  (%rdi),%xmm0
  movss  %xmm0,0x8(%rdx)
  movss  -0x4(%r15,%rbx,4),%xmm0
  movss  %xmm0,(%rax)
  movss  (%r15,%rbx,4),%xmm0
  movss  %xmm0,0x4(%rax)
  lea    -0x1(%rcx),%eax
  cltq
  movss  (%r15,%rax,4),%xmm0
  movslq %ecx,%rax
  movss  %xmm0,(%rdx)
  movss  (%r15,%rax,4),%xmm0
  movss  %xmm0,0x4(%rdx)
  jmp    L47
  nopl   (%rax)
L134:
  movss  %xmm0,0x8(%r13)
  mov    0x80(%rsp),%r10
  mov    0x88(%rsp),%r9
  movss  %xmm0,0x4(%r13)
  movss  (%rdi),%xmm0
  movss  %xmm0,0x8(%rax)
  movss  (%r10),%xmm0
  movss  %xmm0,0x0(%r13)
  movss  (%r9),%xmm0
  movss  %xmm0,0x8(%rdx)
  mov    0x50(%rsp),%rdx
  movss  (%rdi),%xmm0
  movss  %xmm0,0x4(%r13)
  movss  (%rdx),%xmm0
  mov    0x48(%rsp),%rdx
  movss  %xmm0,0x8(%r13)
  movss  (%rdx),%xmm0
  mov    0x58(%rsp),%rdx
  movss  %xmm0,0x4(%rax)
  movss  (%rdx),%xmm0
  movss  %xmm0,0x8(%rax)
  jmp    L47
  nopl   0x0(%rax)
L147:
  mov    0xc0(%rsp),%rdi
  mov    0x90(%rsp),%r10
  movss  (%rdi),%xmm0
  movss  %xmm0,0x8(%rdx)
  movss  %xmm0,0x8(%rax)
  movss  (%r10),%xmm0
  movss  %xmm0,0x8(%r13)
  movss  (%rdi),%xmm0
  mov    0xc8(%rsp),%rdi
  movss  %xmm0,0x4(%rdx)
  movss  (%rdi),%xmm0
  movss  %xmm0,(%rdx)
  lea    -0x1(%rbp),%edx
  movslq %edx,%rdx
  movss  (%r15,%rdx,4),%xmm0
  movslq %ebp,%rdx
  movss  %xmm0,0x0(%r13)
  movss  (%r15,%rdx,4),%xmm0
  lea    -0x1(%rsi),%edx
  movslq %edx,%rdx
  movss  %xmm0,0x4(%r13)
  movss  (%r15,%rdx,4),%xmm0
  movslq %esi,%rdx
  movss  %xmm0,(%rax)
  movss  (%r15,%rdx,4),%xmm0
  movss  %xmm0,0x4(%rax)
  jmp    L47
```

assistant: ```c
void aan_mascara_canal(float *canal_input, float *canal_output, int width, int height, float **m) {
  int i, j, k, l;
  float **area;
  {
    int ml;
    area = (float **)malloc(sizeof(float *) * (3));
    area[0] = (float *)malloc(sizeof(float) * (3) * (3));
    for (ml = 0; ml < (3); ml++)
      area[ml] = &(area[0][ml * (3)]);
  };
  for (j = 0; j < height; j++) {
    for (i = 0; i < width; i++) {
      if (i == 0 && j == 0) {
        area[0][0] = canal_input[0];
        area[0][1] = canal_input[0];
        area[0][2] = canal_input[1];
        area[1][0] = canal_input[0];
        area[2][0] = canal_input[width];
        for (k = 1; k < 3; k++)
          for (l = 1; l < 3; l++)
            area[k][l] = canal_input[(width * (j - 1 + k)) + (i - 1 + l)];
      } else if (i == (width - 1) && j == 0) {
        area[0][2] = canal_input[width - 1];
        area[0][1] = canal_input[width - 1];
        area[0][0] = canal_input[width - 2];
        area[1][2] = canal_input[width - 1];
        area[2][2] = canal_input[width * 2 - 1];
        for (k = 1; k < 3; k++)
          for (l = 0; l < 2; l++)
            area[k][l] = canal_input[(width * (j - 1 + k)) + (i - 1 + l)];
      } else if (i == 0 && j == (height - 1)) {
        area[0][2] = canal_input[width * (j - 1)];
        area[0][1] = canal_input[width * (j - 1)];
        area[1][2] = canal_input[width * (j - 1)];
        area[0][0] = canal_input[width * (j - 2)];
        area[2][2] = canal_input[(width * (j - 1)) + 1];
        for (k = 0; k < 2; k++)
          for (l = 1; l < 3; l++)
            area[k][l] = canal_input[(width * (j - 1 + k)) + (i - 1 + l)];
      } else if (i == (width - 1) && j == (height - 1)) {
        area[2][2] = canal_input[width * height - 1];
        area[1][2] = canal_input[width * height - 1];
        area[0][2] = canal_input[width * (height - 1) - 1];
        area[2][1] = canal_input[width * height - 1];
        area[2][0] = canal_input[width * height - 2];
        for (k = 0; k < 2; k++)
          for (l = 0; l < 2; l++)
            area[k][l] = canal_input[(width * (j - 1 + k)) + (i - 1 + l)];
      } else if (i == 0) {
        area[0][0] = canal_input[(width * (j - 1 + 0)) + i];
        area[1][0] = canal_input[(width * (j - 1 + 1)) + i];
        area[2][0] = canal_input[(width * (j - 1 + 2)) + i];
        for (k = 0; k < 3; k++)
          for (l = 1; l < 3; l++)
            area[k][l] = canal_input[(width * (j - 1 + k)) + (i - 1 + l)];
      } else if (i == (width - 1)) {
        area[0][2] = canal_input[(width * (j - 1 + 0)) + i];
        area[1][2] = canal_input[(width * (j - 1 + 1)) + i];
        area[2][2] = canal_input[(width * (j - 1 + 2)) + i];
        for (k = 0; k < 3; k++)
          for (l = 0; l < 2; l++)
            area[k][l] = canal_input[(width * (j - 1 + k)) + (i - 1 + l)];
      } else if (j == 0) {
        area[0][0] = canal_input[(width * j) + (i - 1 + 0)];
        area[0][1] = canal_input[(width * j) + (i - 1 + 1)];
        area[0][2] = canal_input[(width * j) + (i - 1 + 2)];
        for (k = 1; k < 3; k++)
          for (l = 0; l < 3; l++)
            area[k][l] = canal_input[(width * (j - 1 + k)) + (i - 1 + l)];
      } else if (j == (height - 1)) {
        area[2][0] = canal_input[(width * j) + (i - 1 + 0)];
        area[2][1] = canal_input[(width * j) + (i - 1 + 1)];
        area[2][2] = canal_input[(width * j) + (i - 1 + 2)];
        for (k = 0; k < 2; k++)
          for (l = 0; l < 3; l++)
            area[k][l] = canal_input[(width * (j - 1 + k)) + (i - 1 + l)];
      } else {
        for (k = 0; k < 3; k++)
          for (l = 0; l < 3; l++)
            area[k][l] = canal_input[width * (j - 1 + k) + (i - 1 + l)];
      }
      canal_output[width * j + i] = pixel_resultante(m, area);
    }
  }
  {
    free(area[0]);
    free(area);
  };
}
```