user: What is the c source code of the assembly code below?

```asm
<kernel_dsymv_4_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    L88
  mov    %edx,%edi
  cmpl   $0x1,0x58(%rsp)
  je     L30
  movsd  (%rcx),%xmm4
  movq   D8(%rip),%xmm0
  xorpd  %xmm0,%xmm4
  movsd  0x8(%rcx),%xmm5
  xorpd  %xmm0,%xmm5
  movsd  0x10(%rcx),%xmm6
  xorpd  %xmm0,%xmm6
  movsd  0x18(%rcx),%xmm7
  xorpd  %xmm0,%xmm7
L70:
  cmpl   $0x1,0x50(%rsp)
  je     L36
  pxor   %xmm1,%xmm1
  movapd %xmm1,%xmm3
  movapd %xmm1,%xmm2
  movapd %xmm1,%xmm10
  mov    $0x0,%ebp
L103:
  lea    -0x3(%r10),%eax
  cmp    %eax,%ebp
  jge    L49
  shl    $0x2,%edi
  movslq %edi,%rdi
  shl    $0x3,%rdi
  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    %rsi,%rax
L107:
  movsd  (%rcx),%xmm8
  movsd  (%rax),%xmm9
  movsd  0x20(%rax),%xmm12
  movsd  0x40(%rax),%xmm11
  movsd  0x60(%rax),%xmm0
  movapd %xmm4,%xmm14
  mulsd  %xmm9,%xmm14
  addsd  (%r8),%xmm14
  mulsd  %xmm8,%xmm9
  addsd  %xmm10,%xmm9
  movapd %xmm5,%xmm13
  mulsd  %xmm12,%xmm13
  addsd  %xmm14,%xmm13
  mulsd  %xmm8,%xmm12
  addsd  %xmm12,%xmm2
  movapd %xmm6,%xmm10
  mulsd  %xmm11,%xmm10
  addsd  %xmm13,%xmm10
  mulsd  %xmm8,%xmm11
  addsd  %xmm11,%xmm3
  mulsd  %xmm0,%xmm8
  addsd  %xmm8,%xmm1
  mulsd  %xmm7,%xmm0
  addsd  %xmm10,%xmm0
  movsd  %xmm0,(%rdx)
  movsd  0x8(%rcx),%xmm8
  movsd  0x8(%rax),%xmm10
  movsd  0x28(%rax),%xmm12
  movsd  0x48(%rax),%xmm11
  movsd  0x68(%rax),%xmm0
  movapd %xmm4,%xmm14
  mulsd  %xmm10,%xmm14
  addsd  0x8(%r8),%xmm14
  mulsd  %xmm8,%xmm10
  addsd  %xmm9,%xmm10
  movapd %xmm5,%xmm13
  mulsd  %xmm12,%xmm13
  addsd  %xmm14,%xmm13
  mulsd  %xmm8,%xmm12
  addsd  %xmm12,%xmm2
  movapd %xmm6,%xmm9
  mulsd  %xmm11,%xmm9
  addsd  %xmm13,%xmm9
  mulsd  %xmm8,%xmm11
  addsd  %xmm11,%xmm3
  mulsd  %xmm0,%xmm8
  addsd  %xmm1,%xmm8
  mulsd  %xmm7,%xmm0
  addsd  %xmm9,%xmm0
  movsd  %xmm0,0x8(%rdx)
  movsd  0x10(%rcx),%xmm1
  movsd  0x10(%rax),%xmm12
  movsd  0x30(%rax),%xmm11
  movsd  0x50(%rax),%xmm9
  movsd  0x70(%rax),%xmm0
  movapd %xmm4,%xmm13
  mulsd  %xmm12,%xmm13
  addsd  0x10(%r8),%xmm13
  mulsd  %xmm1,%xmm12
  addsd  %xmm12,%xmm10
  movapd %xmm5,%xmm12
  mulsd  %xmm11,%xmm12
  addsd  %xmm13,%xmm12
  mulsd  %xmm1,%xmm11
  addsd  %xmm2,%xmm11
  movapd %xmm6,%xmm2
  mulsd  %xmm9,%xmm2
  addsd  %xmm12,%xmm2
  mulsd  %xmm1,%xmm9
  addsd  %xmm3,%xmm9
  mulsd  %xmm0,%xmm1
  addsd  %xmm1,%xmm8
  mulsd  %xmm7,%xmm0
  addsd  %xmm2,%xmm0
  movsd  %xmm0,0x10(%rdx)
  movsd  0x18(%rcx),%xmm1
  movsd  0x18(%rax),%xmm12
  movsd  0x38(%rax),%xmm2
  movsd  0x58(%rax),%xmm3
  movsd  0x78(%rax),%xmm0
  movapd %xmm4,%xmm13
  mulsd  %xmm12,%xmm13
  addsd  0x18(%r8),%xmm13
  mulsd  %xmm1,%xmm12
  addsd  %xmm12,%xmm10
  movapd %xmm5,%xmm12
  mulsd  %xmm2,%xmm12
  addsd  %xmm13,%xmm12
  mulsd  %xmm1,%xmm2
  addsd  %xmm11,%xmm2
  movapd %xmm6,%xmm11
  mulsd  %xmm3,%xmm11
  addsd  %xmm12,%xmm11
  mulsd  %xmm1,%xmm3
  addsd  %xmm9,%xmm3
  mulsd  %xmm0,%xmm1
  addsd  %xmm8,%xmm1
  mulsd  %xmm7,%xmm0
  addsd  %xmm11,%xmm0
  movsd  %xmm0,0x18(%rdx)
  add    %rdi,%rax
  add    $0x20,%r8
  add    $0x20,%rdx
  add    $0x20,%rcx
  cmp    %r11,%r8
  jne    L107
  imul   %r15,%rdi
  add    %rdi,%rsi
  mov    -0x8(%rsp),%rax
  add    %rax,%r9
  add    %rax,%rbx
  lea    0x4(%rbp,%r12,4),%ebp
L101:
  cmp    %ebp,%r10d
  jle    L44
  sub    $0x1,%r10d
  sub    %ebp,%r10d
  mov    $0x0,%eax
L41:
  movsd  (%rbx,%rax,8),%xmm8
  movsd  (%rsi,%rax,8),%xmm12
  movsd  0x20(%rsi,%rax,8),%xmm11
  movsd  0x40(%rsi,%rax,8),%xmm9
  movsd  0x60(%rsi,%rax,8),%xmm0
  movapd %xmm4,%xmm13
  mulsd  %xmm12,%xmm13
  addsd  (%r11,%rax,8),%xmm13
  mulsd  %xmm8,%xmm12
  addsd  %xmm12,%xmm10
  movapd %xmm5,%xmm12
  mulsd  %xmm11,%xmm12
  addsd  %xmm13,%xmm12
  mulsd  %xmm8,%xmm11
  addsd  %xmm11,%xmm2
  movapd %xmm6,%xmm11
  mulsd  %xmm9,%xmm11
  addsd  %xmm12,%xmm11
  mulsd  %xmm8,%xmm9
  addsd  %xmm9,%xmm3
  mulsd  %xmm0,%xmm8
  addsd  %xmm8,%xmm1
  mulsd  %xmm7,%xmm0
  addsd  %xmm11,%xmm0
  movsd  %xmm0,(%r9,%rax,8)
  mov    %rax,%rdx
  add    $0x1,%rax
  cmp    %r10,%rdx
  jne    L41
L44:
  cmpl   $0x1,0x60(%rsp)
  je     L15
  movsd  (%r14),%xmm0
  subsd  %xmm10,%xmm0
  movsd  %xmm0,0x0(%r13)
  movsd  0x8(%r14),%xmm0
  subsd  %xmm2,%xmm0
  movsd  %xmm0,0x8(%r13)
  movsd  0x10(%r14),%xmm0
  subsd  %xmm3,%xmm0
  movsd  %xmm0,0x10(%r13)
  movsd  0x18(%r14),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,0x18(%r13)
L88:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L30:
  movsd  (%rcx),%xmm4
  movsd  0x8(%rcx),%xmm5
  movsd  0x10(%rcx),%xmm6
  movsd  0x18(%rcx),%xmm7
  jmp    L70
L36:
  movsd  (%r8),%xmm0
  movsd  (%rbx),%xmm1
  mulsd  (%rsi),%xmm1
  pxor   %xmm8,%xmm8
  addsd  %xmm8,%xmm1
  movsd  %xmm0,(%r9)
  movsd  0x8(%rbx),%xmm2
  movsd  0x8(%rsi),%xmm0
  movapd %xmm2,%xmm9
  mulsd  %xmm0,%xmm9
  addsd  %xmm1,%xmm9
  mulsd  0x28(%rsi),%xmm2
  addsd  %xmm8,%xmm2
  mulsd  %xmm4,%xmm0
  addsd  0x8(%r8),%xmm0
  movsd  %xmm0,0x8(%r9)
  movsd  0x10(%rbx),%xmm1
  movsd  0x10(%rsi),%xmm3
  movsd  0x30(%rsi),%xmm0
  movapd %xmm4,%xmm10
  mulsd  %xmm3,%xmm10
  addsd  0x10(%r8),%xmm10
  mulsd  %xmm1,%xmm3
  addsd  %xmm9,%xmm3
  movapd %xmm1,%xmm11
  mulsd  %xmm0,%xmm11
  addsd  %xmm2,%xmm11
  mulsd  0x50(%rsi),%xmm1
  addsd  %xmm8,%xmm1
  movapd %xmm1,%xmm9
  mulsd  %xmm5,%xmm0
  addsd  %xmm10,%xmm0
  movsd  %xmm0,0x10(%r9)
  movsd  0x18(%rbx),%xmm1
  movsd  0x18(%rsi),%xmm10
  movsd  0x38(%rsi),%xmm2
  movsd  0x58(%rsi),%xmm0
  movapd %xmm4,%xmm13
  mulsd  %xmm10,%xmm13
  addsd  0x18(%r8),%xmm13
  mulsd  %xmm1,%xmm10
  addsd  %xmm3,%xmm10
  movapd %xmm5,%xmm12
  mulsd  %xmm2,%xmm12
  addsd  %xmm13,%xmm12
  mulsd  %xmm1,%xmm2
  addsd  %xmm11,%xmm2
  movapd %xmm1,%xmm3
  mulsd  %xmm0,%xmm3
  addsd  %xmm9,%xmm3
  mulsd  0x78(%rsi),%xmm1
  addsd  %xmm8,%xmm1
  mulsd  %xmm6,%xmm0
  addsd  %xmm12,%xmm0
  movsd  %xmm0,0x18(%r9)
  add    $0x20,%r8
  add    $0x20,%r9
  add    $0x20,%rbx
  lea    -0x4(,%rdi,4),%eax
  cltq
  lea    0x20(%rsi,%rax,8),%rsi
  mov    $0x4,%ebp
  jmp    L103
L49:
  mov    %r8,%r11
  jmp    L101
L15:
  addsd  (%r14),%xmm10
  movsd  %xmm10,0x0(%r13)
  addsd  0x8(%r14),%xmm2
  movsd  %xmm2,0x8(%r13)
  addsd  0x10(%r14),%xmm3
  movsd  %xmm3,0x10(%r13)
  addsd  0x18(%r14),%xmm1
  movsd  %xmm1,0x18(%r13)
  jmp    L88
```

tool_calls: name: retrieve_data
arguments: label: D8
datatype: byte[8]

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

assistant: ```c
void kernel_dsymv_4_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_n, int alg_t) {
  if (kmax <= 0)
    return;
  const int bs = 4;
  int k;
  double a_00, a_01, a_02, a_03, x_n_0, x_n_1, x_n_2, x_n_3, y_n_0, x_t_0, y_t_0, y_t_1, y_t_2, y_t_3;
  if (alg_n == 1) {
    x_n_0 = x_n[0];
    x_n_1 = x_n[1];
    x_n_2 = x_n[2];
    x_n_3 = x_n[3];
  } else {
    x_n_0 = -x_n[0];
    x_n_1 = -x_n[1];
    x_n_2 = -x_n[2];
    x_n_3 = -x_n[3];
  }
  y_t_0 = 0;
  y_t_1 = 0;
  y_t_2 = 0;
  y_t_3 = 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;
    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];
    a_03 = A[3 + bs * 3];
    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;
    y_t_3 += a_03 * x_t_0;
    z_n[3] = y_n_0;
    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];
    a_03 = A[0 + bs * 3];
    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;
    y_n_0 += a_03 * x_n_3;
    y_t_3 += a_03 * 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];
    a_03 = A[1 + bs * 3];
    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;
    y_n_0 += a_03 * x_n_3;
    y_t_3 += a_03 * 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];
    a_03 = A[2 + bs * 3];
    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;
    y_n_0 += a_03 * x_n_3;
    y_t_3 += a_03 * 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];
    a_03 = A[3 + bs * 3];
    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;
    y_n_0 += a_03 * x_n_3;
    y_t_3 += a_03 * 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];
    a_03 = A[0 + bs * 3];
    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;
    y_n_0 += a_03 * x_n_3;
    y_t_3 += a_03 * x_t_0;
    z_n[0] = y_n_0;
    A += 1;
    y_n += 1;
    z_n += 1;
    x_t += 1;
  }
  if (alg_t == 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;
    z_t[3] = y_t[3] + y_t_3;
  } 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;
    z_t[3] = y_t[3] - y_t_3;
  }
}
```