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

```asm
<kernel_dsymv_3_lib4>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %edi,%r10d
  mov    0x38(%rsp),%rbx
  mov    0x40(%rsp),%r14
  mov    0x48(%rsp),%r13
  test   %edi,%edi
  jle    L54
  mov    %rsi,%rdi
  mov    %edx,%esi
  cmpl   $0x1,0x58(%rsp)
  je     L108
  movsd  (%rcx),%xmm2
  movq   D0(%rip),%xmm0
  xorpd  %xmm0,%xmm2
  movsd  0x8(%rcx),%xmm3
  xorpd  %xmm0,%xmm3
  movsd  0x10(%rcx),%xmm4
  xorpd  %xmm0,%xmm4
L133:
  cmpl   $0x1,0x50(%rsp)
  je     L111
  pxor   %xmm1,%xmm1
  movapd %xmm1,%xmm5
  movapd %xmm1,%xmm8
  mov    $0x0,%ebp
L0:
  lea    -0x3(%r10),%eax
  cmp    %eax,%ebp
  jge    L123
  shl    $0x2,%esi
  movslq %esi,%rsi
  shl    $0x3,%rsi
  lea    -0x4(%r10),%r12d
  sub    %ebp,%r12d
  shr    $0x2,%r12d
  lea    0x1(%r12),%r15d
  mov    %r15d,%r15d
  mov    %r15,%rax
  shl    $0x5,%rax
  mov    %rax,-0x8(%rsp)
  lea    (%r8,%rax,1),%r11
  mov    %rbx,%rcx
  mov    %r9,%rdx
  mov    %rdi,%rax
L64:
  movsd  (%rcx),%xmm7
  movsd  (%rax),%xmm6
  movsd  0x20(%rax),%xmm9
  movsd  0x40(%rax),%xmm0
  movapd %xmm2,%xmm10
  mulsd  %xmm6,%xmm10
  addsd  (%r8),%xmm10
  mulsd  %xmm7,%xmm6
  addsd  %xmm8,%xmm6
  movapd %xmm3,%xmm8
  mulsd  %xmm9,%xmm8
  addsd  %xmm10,%xmm8
  mulsd  %xmm7,%xmm9
  addsd  %xmm9,%xmm5
  mulsd  %xmm0,%xmm7
  addsd  %xmm7,%xmm1
  mulsd  %xmm4,%xmm0
  addsd  %xmm8,%xmm0
  movsd  %xmm0,(%rdx)
  movsd  0x8(%rcx),%xmm7
  movsd  0x8(%rax),%xmm8
  movsd  0x28(%rax),%xmm9
  movsd  0x48(%rax),%xmm0
  movapd %xmm2,%xmm10
  mulsd  %xmm8,%xmm10
  addsd  0x8(%r8),%xmm10
  mulsd  %xmm7,%xmm8
  addsd  %xmm6,%xmm8
  movapd %xmm3,%xmm6
  mulsd  %xmm9,%xmm6
  addsd  %xmm10,%xmm6
  mulsd  %xmm7,%xmm9
  addsd  %xmm9,%xmm5
  mulsd  %xmm0,%xmm7
  addsd  %xmm1,%xmm7
  mulsd  %xmm4,%xmm0
  addsd  %xmm6,%xmm0
  movsd  %xmm0,0x8(%rdx)
  movsd  0x10(%rcx),%xmm6
  movsd  0x10(%rax),%xmm10
  movsd  0x30(%rax),%xmm9
  movsd  0x50(%rax),%xmm0
  movapd %xmm2,%xmm11
  mulsd  %xmm10,%xmm11
  addsd  0x10(%r8),%xmm11
  mulsd  %xmm6,%xmm10
  addsd  %xmm8,%xmm10
  movapd %xmm3,%xmm1
  mulsd  %xmm9,%xmm1
  addsd  %xmm11,%xmm1
  mulsd  %xmm6,%xmm9
  addsd  %xmm5,%xmm9
  mulsd  %xmm0,%xmm6
  addsd  %xmm7,%xmm6
  mulsd  %xmm4,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,0x10(%rdx)
  movsd  0x18(%rcx),%xmm1
  movsd  0x18(%rax),%xmm8
  movsd  0x38(%rax),%xmm5
  movsd  0x58(%rax),%xmm0
  movapd %xmm2,%xmm11
  mulsd  %xmm8,%xmm11
  addsd  0x18(%r8),%xmm11
  mulsd  %xmm1,%xmm8
  addsd  %xmm10,%xmm8
  movapd %xmm3,%xmm7
  mulsd  %xmm5,%xmm7
  addsd  %xmm11,%xmm7
  mulsd  %xmm1,%xmm5
  addsd  %xmm9,%xmm5
  mulsd  %xmm0,%xmm1
  addsd  %xmm6,%xmm1
  mulsd  %xmm4,%xmm0
  addsd  %xmm7,%xmm0
  movsd  %xmm0,0x18(%rdx)
  add    %rsi,%rax
  add    $0x20,%r8
  add    $0x20,%rdx
  add    $0x20,%rcx
  cmp    %r11,%r8
  jne    L64
  imul   %r15,%rsi
  add    %rsi,%rdi
  mov    -0x8(%rsp),%rax
  add    %rax,%r9
  add    %rax,%rbx
  lea    0x4(%rbp,%r12,4),%ebp
L5:
  cmp    %ebp,%r10d
  jle    L9
  sub    $0x1,%r10d
  sub    %ebp,%r10d
  mov    $0x0,%eax
L39:
  movsd  (%rbx,%rax,8),%xmm6
  movsd  (%rdi,%rax,8),%xmm9
  movsd  0x20(%rdi,%rax,8),%xmm7
  movsd  0x40(%rdi,%rax,8),%xmm0
  movapd %xmm2,%xmm10
  mulsd  %xmm9,%xmm10
  addsd  (%r11,%rax,8),%xmm10
  mulsd  %xmm6,%xmm9
  addsd  %xmm9,%xmm8
  movapd %xmm3,%xmm9
  mulsd  %xmm7,%xmm9
  addsd  %xmm10,%xmm9
  mulsd  %xmm6,%xmm7
  addsd  %xmm7,%xmm5
  mulsd  %xmm0,%xmm6
  addsd  %xmm6,%xmm1
  mulsd  %xmm4,%xmm0
  addsd  %xmm9,%xmm0
  movsd  %xmm0,(%r9,%rax,8)
  mov    %rax,%rdx
  add    $0x1,%rax
  cmp    %r10,%rdx
  jne    L39
L9:
  cmpl   $0x1,0x58(%rsp)
  je     L3
  movsd  (%r14),%xmm0
  subsd  %xmm8,%xmm0
  movsd  %xmm0,0x0(%r13)
  movsd  0x8(%r14),%xmm0
  subsd  %xmm5,%xmm0
  movsd  %xmm0,0x8(%r13)
  movsd  0x10(%r14),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,0x10(%r13)
L54:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L108:
  movsd  (%rcx),%xmm2
  movsd  0x8(%rcx),%xmm3
  movsd  0x10(%rcx),%xmm4
  jmp    L133
L111:
  movsd  (%r8),%xmm0
  movsd  (%rbx),%xmm1
  mulsd  (%rdi),%xmm1
  pxor   %xmm6,%xmm6
  addsd  %xmm6,%xmm1
  movsd  %xmm0,(%r9)
  movsd  0x8(%rbx),%xmm5
  movsd  0x8(%rdi),%xmm0
  movapd %xmm5,%xmm9
  mulsd  %xmm0,%xmm9
  addsd  %xmm1,%xmm9
  mulsd  0x28(%rdi),%xmm5
  addsd  %xmm6,%xmm5
  mulsd  %xmm2,%xmm0
  addsd  0x8(%r8),%xmm0
  movsd  %xmm0,0x8(%r9)
  movsd  0x10(%rbx),%xmm1
  movsd  0x10(%rdi),%xmm8
  movsd  0x30(%rdi),%xmm0
  movapd %xmm2,%xmm7
  mulsd  %xmm8,%xmm7
  addsd  0x10(%r8),%xmm7
  mulsd  %xmm1,%xmm8
  addsd  %xmm9,%xmm8
  movapd %xmm1,%xmm9
  mulsd  %xmm0,%xmm9
  addsd  %xmm9,%xmm5
  mulsd  0x50(%rdi),%xmm1
  addsd  %xmm6,%xmm1
  mulsd  %xmm3,%xmm0
  addsd  %xmm7,%xmm0
  movsd  %xmm0,0x10(%r9)
  cmp    $0x3,%r10d
  je     L9
  movsd  0x18(%rbx),%xmm6
  movsd  0x18(%rdi),%xmm9
  movsd  0x38(%rdi),%xmm7
  movsd  0x58(%rdi),%xmm0
  movapd %xmm2,%xmm10
  mulsd  %xmm9,%xmm10
  addsd  0x18(%r8),%xmm10
  mulsd  %xmm6,%xmm9
  addsd  %xmm9,%xmm8
  movapd %xmm3,%xmm9
  mulsd  %xmm7,%xmm9
  addsd  %xmm10,%xmm9
  mulsd  %xmm6,%xmm7
  addsd  %xmm7,%xmm5
  mulsd  %xmm0,%xmm6
  addsd  %xmm6,%xmm1
  mulsd  %xmm4,%xmm0
  addsd  %xmm9,%xmm0
  movsd  %xmm0,0x18(%r9)
  cmp    $0x4,%r10d
  je     L9
  add    $0x20,%r8
  add    $0x20,%r9
  add    $0x20,%rbx
  lea    -0x4(,%rsi,4),%eax
  cltq
  lea    0x20(%rdi,%rax,8),%rdi
  mov    $0x4,%ebp
  jmp    L0
L123:
  mov    %r8,%r11
  jmp    L5
L3:
  addsd  (%r14),%xmm8
  movsd  %xmm8,0x0(%r13)
  addsd  0x8(%r14),%xmm5
  movsd  %xmm5,0x8(%r13)
  addsd  0x10(%r14),%xmm1
  movsd  %xmm1,0x10(%r13)
  jmp    L54
```

tool_calls: name: disassemble_data
arguments: tag: D0
type: byte[8]

tool: {"D0": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80]}

assistant: ```c
void kernel_dsymv_3_lib4(int kmax, double *A, int sda, double *x_n, double *y_n, double *z_n, double *x_t, double *y_t, double *z_t, int tri, int alg) {
  if (kmax <= 0)
    return;
  const int bs = 4;
  int k;
  double a_00, a_01, a_02, x_n_0, x_n_1, x_n_2, y_n_0, x_t_0, y_t_0, y_t_1, y_t_2;
  if (alg == 1) {
    x_n_0 = x_n[0];
    x_n_1 = x_n[1];
    x_n_2 = x_n[2];
  } else {
    x_n_0 = -x_n[0];
    x_n_1 = -x_n[1];
    x_n_2 = -x_n[2];
  }
  y_t_0 = 0;
  y_t_1 = 0;
  y_t_2 = 0;
  k = 0;
  if (tri == 1) {
    y_n_0 = y_n[0];
    x_t_0 = x_t[0];
    a_00 = A[0 + bs * 0];
    y_t_0 += a_00 * x_t_0;
    z_n[0] = y_n_0;
    y_n_0 = y_n[1];
    x_t_0 = x_t[1];
    a_00 = A[1 + bs * 0];
    a_01 = A[1 + bs * 1];
    y_n_0 += a_00 * x_n_0;
    y_t_0 += a_00 * x_t_0;
    y_t_1 += a_01 * x_t_0;
    z_n[1] = y_n_0;
    y_n_0 = y_n[2];
    x_t_0 = x_t[2];
    a_00 = A[2 + bs * 0];
    a_01 = A[2 + bs * 1];
    a_02 = A[2 + bs * 2];
    y_n_0 += a_00 * x_n_0;
    y_t_0 += a_00 * x_t_0;
    y_n_0 += a_01 * x_n_1;
    y_t_1 += a_01 * x_t_0;
    y_t_2 += a_02 * x_t_0;
    z_n[2] = y_n_0;
    if (kmax == 3)
      goto STORE_3;
    y_n_0 = y_n[3];
    x_t_0 = x_t[3];
    a_00 = A[3 + bs * 0];
    a_01 = A[3 + bs * 1];
    a_02 = A[3 + bs * 2];
    y_n_0 += a_00 * x_n_0;
    y_t_0 += a_00 * x_t_0;
    y_n_0 += a_01 * x_n_1;
    y_t_1 += a_01 * x_t_0;
    y_n_0 += a_02 * x_n_2;
    y_t_2 += a_02 * x_t_0;
    z_n[3] = y_n_0;
    if (kmax == 4)
      goto STORE_3;
    A += 4;
    y_n += 4;
    z_n += 4;
    x_t += 4;
    k += 4;
    A += (sda - 1) * bs;
  }
  for (; k < kmax - 3; k += bs) {
    y_n_0 = y_n[0];
    x_t_0 = x_t[0];
    a_00 = A[0 + bs * 0];
    a_01 = A[0 + bs * 1];
    a_02 = A[0 + bs * 2];
    y_n_0 += a_00 * x_n_0;
    y_t_0 += a_00 * x_t_0;
    y_n_0 += a_01 * x_n_1;
    y_t_1 += a_01 * x_t_0;
    y_n_0 += a_02 * x_n_2;
    y_t_2 += a_02 * x_t_0;
    z_n[0] = y_n_0;
    y_n_0 = y_n[1];
    x_t_0 = x_t[1];
    a_00 = A[1 + bs * 0];
    a_01 = A[1 + bs * 1];
    a_02 = A[1 + bs * 2];
    y_n_0 += a_00 * x_n_0;
    y_t_0 += a_00 * x_t_0;
    y_n_0 += a_01 * x_n_1;
    y_t_1 += a_01 * x_t_0;
    y_n_0 += a_02 * x_n_2;
    y_t_2 += a_02 * x_t_0;
    z_n[1] = y_n_0;
    y_n_0 = y_n[2];
    x_t_0 = x_t[2];
    a_00 = A[2 + bs * 0];
    a_01 = A[2 + bs * 1];
    a_02 = A[2 + bs * 2];
    y_n_0 += a_00 * x_n_0;
    y_t_0 += a_00 * x_t_0;
    y_n_0 += a_01 * x_n_1;
    y_t_1 += a_01 * x_t_0;
    y_n_0 += a_02 * x_n_2;
    y_t_2 += a_02 * x_t_0;
    z_n[2] = y_n_0;
    y_n_0 = y_n[3];
    x_t_0 = x_t[3];
    a_00 = A[3 + bs * 0];
    a_01 = A[3 + bs * 1];
    a_02 = A[3 + bs * 2];
    y_n_0 += a_00 * x_n_0;
    y_t_0 += a_00 * x_t_0;
    y_n_0 += a_01 * x_n_1;
    y_t_1 += a_01 * x_t_0;
    y_n_0 += a_02 * x_n_2;
    y_t_2 += a_02 * x_t_0;
    z_n[3] = y_n_0;
    A += sda * bs;
    y_n += 4;
    z_n += 4;
    x_t += 4;
  }
  for (; k < kmax; k++) {
    y_n_0 = y_n[0];
    x_t_0 = x_t[0];
    a_00 = A[0 + bs * 0];
    a_01 = A[0 + bs * 1];
    a_02 = A[0 + bs * 2];
    y_n_0 += a_00 * x_n_0;
    y_t_0 += a_00 * x_t_0;
    y_n_0 += a_01 * x_n_1;
    y_t_1 += a_01 * x_t_0;
    y_n_0 += a_02 * x_n_2;
    y_t_2 += a_02 * x_t_0;
    z_n[0] = y_n_0;
    A += 1;
    y_n += 1;
    z_n += 1;
    x_t += 1;
  }
STORE_3:
  if (alg == 1) {
    z_t[0] = y_t[0] + y_t_0;
    z_t[1] = y_t[1] + y_t_1;
    z_t[2] = y_t[2] + y_t_2;
  } else {
    z_t[0] = y_t[0] - y_t_0;
    z_t[1] = y_t[1] - y_t_1;
    z_t[2] = y_t[2] - y_t_2;
  }
}
```