user: 请显示与下面列出的汇编代码相匹配的C语言实现:

```asm
<vec_subtract_half>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x60,%rsp
  mov    %rdi,-0x48(%rbp)
  mov    %rsi,-0x50(%rbp)
  mov    %rdx,-0x58(%rbp)
  mov    -0x48(%rbp),%rax
  mov    $0x40,%esi
  mov    %rax,%rdi
  call   <__assume_aligned@plt>
  mov    -0x50(%rbp),%rax
  mov    $0x40,%esi
  mov    %rax,%rdi
  call   <__assume_aligned@plt>
  movq   $0x40,-0x10(%rbp)
  movq   $0x0,-0x38(%rbp)
  jmp    L78
L15:
  movq   $0x0,-0x30(%rbp)
  jmp    L45
L56:
  movq   $0x0,-0x28(%rbp)
  jmp    L16
L99:
  mov    -0x30(%rbp),%rax
  imul   -0x58(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x30(%rbp),%rax
  imul   -0x58(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  movsd  D1(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  mov    -0x30(%rbp),%rax
  imul   -0x58(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addq   $0x1,-0x28(%rbp)
L16:
  mov    -0x28(%rbp),%rax
  cmp    -0x10(%rbp),%rax
  jb     L99
  addq   $0x1,-0x30(%rbp)
L45:
  mov    -0x30(%rbp),%rax
  cmp    -0x58(%rbp),%rax
  jb     L56
  mov    -0x10(%rbp),%rax
  add    %rax,-0x38(%rbp)
L78:
  mov    -0x58(%rbp),%rax
  sub    -0x10(%rbp),%rax
  cmp    %rax,-0x38(%rbp)
  jb     L15
  mov    -0x58(%rbp),%rax
  mov    $0x0,%edx
  divq   -0x10(%rbp)
  mov    -0x58(%rbp),%rax
  sub    %rdx,%rax
  mov    %rax,-0x8(%rbp)
  movq   $0x0,-0x20(%rbp)
  jmp    L27
L77:
  mov    -0x8(%rbp),%rax
  mov    %rax,-0x18(%rbp)
  jmp    L82
L36:
  mov    -0x20(%rbp),%rax
  imul   -0x58(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x20(%rbp),%rax
  imul   -0x58(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  movsd  D1(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  mov    -0x20(%rbp),%rax
  imul   -0x58(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addq   $0x1,-0x18(%rbp)
L82:
  mov    -0x18(%rbp),%rax
  cmp    -0x58(%rbp),%rax
  jb     L36
  addq   $0x1,-0x20(%rbp)
L27:
  mov    -0x20(%rbp),%rax
  cmp    -0x58(%rbp),%rax
  jb     L77
  nop
  nop
  leave
  ret
```

tool_calls: name: extract_symbolic_data
arguments: tag: D1
type_guessed: f64

tool: D1: 0.5

assistant: ```c
void vec_subtract_half(double *restrict u, const double *restrict v, size_t n) {
  __assume_aligned(u, 64);
  __assume_aligned(v, 64);
  size_t block = 64;
  for (size_t i = 0; i < n - block; i += block)
    for (size_t j = 0; j < n; ++j) {
      for (size_t b = 0; b < block; ++b) {
        u[j * n + b] -= 0.5 * v[j * n + b];
      }
    }
  size_t start = n - n % block;
  for (size_t j = 0; j < n; ++j)
    for (size_t b = start; b < n; ++b) {
      u[j * n + b] -= 0.5 * v[j * n + b];
    }
}
```