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

```asm
<laplacian>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %rsi,-0x10(%rsp)
  mov    %edx,%r11d
  mov    %ecx,%r13d
  lea    -0x1(%rcx),%eax
  mov    %eax,-0x8(%rsp)
  cmp    $0x1,%eax
  jle    L49
  mov    %rdi,%r12
  mov    %eax,%r14d
  movslq %edx,%r9
  mov    %r9,%r8
  neg    %r8
  shl    $0x2,%r8
  mov    $0x0,%ebx
  mov    $0x1,%ebp
  lea    -0x3(%rdx),%eax
  lea    0x4(%rdi),%r15
  movsd  D26(%rip),%xmm4
  movsd  D4(%rip),%xmm3
  movsd  D3(%rip),%xmm2
  mov    %ecx,-0x4(%rsp)
  mov    %rax,%r13
  jmp    L60
L75:
  lea    (%rbx,%r11,1),%r10d
  lea    (%r11,%r10,1),%ecx
  movslq %r10d,%r10
  lea    0x0(,%r10,4),%rdx
  lea    (%r12,%rdx,1),%rax
  add    -0x10(%rsp),%rdx
  lea    0x0(%r13,%r10,1),%rsi
  lea    (%r15,%rsi,4),%rdi
  movslq %ebx,%rsi
  sub    %r10,%rsi
  movslq %ecx,%rcx
  sub    %r10,%rcx
L36:
  movss  0x8(%rax),%xmm0
  addss  0x4(%rax,%r9,4),%xmm0
  addss  (%rax),%xmm0
  addss  0x4(%rax,%r8,1),%xmm0
  cvtss2sd %xmm0,%xmm0
  mulsd  %xmm4,%xmm0
  movss  (%rax,%rsi,4),%xmm1
  addss  0x8(%rax,%rsi,4),%xmm1
  addss  0x8(%rax,%rcx,4),%xmm1
  addss  (%rax,%rcx,4),%xmm1
  cvtss2sd %xmm1,%xmm1
  addsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd 0x4(%rax),%xmm1
  mulsd  %xmm3,%xmm1
  addsd  %xmm1,%xmm0
  divsd  %xmm2,%xmm0
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,0x4(%rdx)
  add    $0x4,%rax
  add    $0x4,%rdx
  cmp    %rdi,%rax
  jne    L36
L28:
  add    $0x1,%ebp
  add    %r11d,%ebx
  cmp    %r14d,%ebp
  je     L11
L60:
  cmp    $0x2,%r11d
  jg     L75
  jmp    L28
L11:
  mov    -0x4(%rsp),%r13d
L49:
  test   %r11d,%r11d
  jle    L53
  mov    -0x8(%rsp),%edx
  imul   %r11d,%edx
  movslq %edx,%rdx
  mov    -0x10(%rsp),%rbx
  lea    (%rbx,%rdx,4),%rax
  lea    -0x1(%r11),%ecx
  add    %rdx,%rcx
  lea    0x4(%rbx,%rcx,4),%rcx
  neg    %rdx
  shl    $0x2,%rdx
L76:
  movl   $0x0,(%rax)
  movl   $0x0,(%rax,%rdx,1)
  add    $0x4,%rax
  cmp    %rcx,%rax
  jne    L76
L53:
  test   %r13d,%r13d
  jle    L69
  movslq %r11d,%r11
  lea    0x0(,%r11,4),%rcx
  mov    -0x10(%rsp),%rax
  mov    $0x0,%edx
L32:
  movl   $0x0,-0x4(%rax,%r11,4)
  movl   $0x0,(%rax)
  add    $0x1,%edx
  add    %rcx,%rax
  cmp    %edx,%r13d
  jne    L32
L69:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: unpack_data
arguments: label: D26
type: f64
name: unpack_data
arguments: label: D4
type: f32[2]
name: unpack_data
arguments: label: D3
type: f64

tool: D26: 4

tool: {"D4": [0.0,-2.8125]}

tool: D3: 6

assistant: ```c
void laplacian(image, lap, nx, ny) float *lap, *image;
int nx, ny;
{
  int index, i, j;
  for (i = 1; i < ny - 1; i++) {
    for (j = 1; j < nx - 1; j++) {
      index = (i * nx) + j;
      lap[index] = (image[(i - 1) * nx + (j - 1)] + image[(i - 1) * nx + (j + 1)] + image[(i + 1) * nx + (j + 1)] + image[(i + 1) * nx + (j - 1)] + (4.0 * (image[index + 1] + image[index + nx] + image[index - 1] + image[index - nx])) + (-20.0 * image[index])) / 6.0;
    }
  }
  for (j = 0; j < nx; j++) {
    lap[(ny - 1) * nx + j] = 0.0;
    lap[j] = 0.0;
  }
  for (i = 0; i < ny; i++) {
    lap[(i * nx) + nx - 1] = 0.0;
    lap[i * nx] = 0.0;
  }
}
```