user: 下面的汇编指令对应的C源代码是什么:

```asm
<kernel_dgemm_diag_left_4_lib4>:
  endbr64
  mov    %rsi,%r10
  mov    %rcx,%rax
  mov    %rdx,%rsi
  test   %edi,%edi
  jle    L84
  movsd  (%r10),%xmm0
  movsd  0x8(%r10),%xmm1
  lea    -0x3(%rdi),%edx
  movsd  0x10(%r10),%xmm2
  movsd  0x18(%r10),%xmm3
  cmp    $0xffffffff,%r9d
  je     L45
  test   %r9d,%r9d
  jne    L7
  test   %edx,%edx
  jle    L88
  lea    -0x4(%rdi),%r9d
  mov    %rsi,%rax
  shr    $0x2,%r9d
  lea    0x1(%r9),%ecx
  mov    %rcx,%r9
  shl    $0x7,%rcx
  lea    (%r8,%rcx,1),%rdx
  nopl   0x0(%rax)
L50:
  movsd  0x8(%rax),%xmm6
  movsd  (%rax),%xmm7
  sub    $0xffffffffffffff80,%r8
  sub    $0xffffffffffffff80,%rax
  movsd  -0x70(%rax),%xmm5
  movsd  -0x68(%rax),%xmm4
  mulsd  %xmm1,%xmm6
  mulsd  %xmm2,%xmm5
  mulsd  %xmm3,%xmm4
  mulsd  %xmm0,%xmm7
  movsd  %xmm6,-0x78(%r8)
  movsd  %xmm5,-0x70(%r8)
  movsd  %xmm4,-0x68(%r8)
  movsd  %xmm7,-0x80(%r8)
  movsd  -0x58(%rax),%xmm6
  movsd  -0x50(%rax),%xmm5
  movsd  -0x48(%rax),%xmm4
  movsd  -0x60(%rax),%xmm7
  mulsd  %xmm1,%xmm6
  mulsd  %xmm2,%xmm5
  mulsd  %xmm3,%xmm4
  mulsd  %xmm0,%xmm7
  movsd  %xmm6,-0x58(%r8)
  movsd  %xmm5,-0x50(%r8)
  movsd  %xmm4,-0x48(%r8)
  movsd  %xmm7,-0x60(%r8)
  movsd  -0x38(%rax),%xmm6
  movsd  -0x30(%rax),%xmm5
  movsd  -0x28(%rax),%xmm4
  movsd  -0x40(%rax),%xmm7
  mulsd  %xmm1,%xmm6
  mulsd  %xmm2,%xmm5
  mulsd  %xmm3,%xmm4
  mulsd  %xmm0,%xmm7
  movsd  %xmm6,-0x38(%r8)
  movsd  %xmm5,-0x30(%r8)
  movsd  %xmm4,-0x28(%r8)
  movsd  %xmm7,-0x40(%r8)
  movsd  -0x18(%rax),%xmm6
  movsd  -0x10(%rax),%xmm5
  movsd  -0x8(%rax),%xmm4
  movsd  -0x20(%rax),%xmm7
  mulsd  %xmm1,%xmm6
  mulsd  %xmm2,%xmm5
  mulsd  %xmm3,%xmm4
  mulsd  %xmm0,%xmm7
  movsd  %xmm6,-0x18(%r8)
  movsd  %xmm5,-0x10(%r8)
  movsd  %xmm4,-0x8(%r8)
  movsd  %xmm7,-0x20(%r8)
  cmp    %r8,%rdx
  jne    L50
  shl    $0x2,%r9d
  add    %rcx,%rsi
  cmp    %r9d,%edi
  jle    L12
L89:
  movsd  (%rsi),%xmm7
  add    $0x1,%r9d
  add    $0x20,%rsi
  add    $0x20,%rdx
  movsd  -0x18(%rsi),%xmm6
  movsd  -0x10(%rsi),%xmm5
  movsd  -0x8(%rsi),%xmm4
  mulsd  %xmm0,%xmm7
  mulsd  %xmm1,%xmm6
  mulsd  %xmm2,%xmm5
  mulsd  %xmm3,%xmm4
  movsd  %xmm7,-0x20(%rdx)
  movsd  %xmm6,-0x18(%rdx)
  movsd  %xmm5,-0x10(%rdx)
  movsd  %xmm4,-0x8(%rdx)
  cmp    %r9d,%edi
  jg     L89
  ret
  nopl   0x0(%rax)
L45:
  movq   D2(%rip),%xmm4
  xorpd  %xmm4,%xmm0
  xorpd  %xmm4,%xmm1
  xorpd  %xmm4,%xmm2
  xorpd  %xmm4,%xmm3
L7:
  mov    %rax,%r9
  xor    %r10d,%r10d
  test   %edx,%edx
  jle    L55
  lea    -0x4(%rdi),%r10d
  mov    %r8,%rcx
  mov    %rsi,%rdx
  shr    $0x2,%r10d
  lea    0x1(%r10),%r11d
  mov    %r11,%r10
  shl    $0x7,%r11
  lea    (%rax,%r11,1),%r9
  nopl   0x0(%rax,%rax,1)
L0:
  movsd  (%rdx),%xmm7
  sub    $0xffffffffffffff80,%rax
  sub    $0xffffffffffffff80,%rdx
  sub    $0xffffffffffffff80,%rcx
  movsd  -0x78(%rdx),%xmm6
  movsd  -0x70(%rdx),%xmm5
  movsd  -0x68(%rdx),%xmm4
  mulsd  %xmm0,%xmm7
  mulsd  %xmm1,%xmm6
  addsd  -0x80(%rax),%xmm7
  mulsd  %xmm2,%xmm5
  addsd  -0x78(%rax),%xmm6
  mulsd  %xmm3,%xmm4
  addsd  -0x70(%rax),%xmm5
  addsd  -0x68(%rax),%xmm4
  movsd  %xmm7,-0x80(%rcx)
  movsd  %xmm6,-0x78(%rcx)
  movsd  %xmm5,-0x70(%rcx)
  movsd  %xmm4,-0x68(%rcx)
  movsd  -0x58(%rdx),%xmm6
  movsd  -0x50(%rdx),%xmm5
  movsd  -0x48(%rdx),%xmm4
  movsd  -0x60(%rdx),%xmm7
  mulsd  %xmm1,%xmm6
  mulsd  %xmm2,%xmm5
  addsd  -0x58(%rax),%xmm6
  mulsd  %xmm3,%xmm4
  addsd  -0x50(%rax),%xmm5
  mulsd  %xmm0,%xmm7
  addsd  -0x48(%rax),%xmm4
  addsd  -0x60(%rax),%xmm7
  movsd  %xmm6,-0x58(%rcx)
  movsd  %xmm5,-0x50(%rcx)
  movsd  %xmm4,-0x48(%rcx)
  movsd  %xmm7,-0x60(%rcx)
  movsd  -0x38(%rdx),%xmm6
  movsd  -0x30(%rdx),%xmm5
  movsd  -0x28(%rdx),%xmm4
  movsd  -0x40(%rdx),%xmm7
  mulsd  %xmm1,%xmm6
  mulsd  %xmm2,%xmm5
  addsd  -0x38(%rax),%xmm6
  mulsd  %xmm3,%xmm4
  addsd  -0x30(%rax),%xmm5
  mulsd  %xmm0,%xmm7
  addsd  -0x28(%rax),%xmm4
  addsd  -0x40(%rax),%xmm7
  movsd  %xmm7,-0x40(%rcx)
  movsd  %xmm6,-0x38(%rcx)
  movsd  %xmm5,-0x30(%rcx)
  movsd  %xmm4,-0x28(%rcx)
  movsd  -0x18(%rdx),%xmm6
  movsd  -0x10(%rdx),%xmm5
  movsd  -0x8(%rdx),%xmm4
  movsd  -0x20(%rdx),%xmm7
  mulsd  %xmm1,%xmm6
  mulsd  %xmm2,%xmm5
  addsd  -0x18(%rax),%xmm6
  mulsd  %xmm3,%xmm4
  addsd  -0x10(%rax),%xmm5
  mulsd  %xmm0,%xmm7
  addsd  -0x8(%rax),%xmm4
  addsd  -0x20(%rax),%xmm7
  movsd  %xmm6,-0x18(%rcx)
  movsd  %xmm5,-0x10(%rcx)
  movsd  %xmm4,-0x8(%rcx)
  movsd  %xmm7,-0x20(%rcx)
  cmp    %rax,%r9
  jne    L0
  shl    $0x2,%r10d
  add    %r11,%rsi
  add    %r11,%r8
  cmp    %r10d,%edi
  jle    L84
L55:
  movsd  (%rsi),%xmm7
  add    $0x1,%r10d
  add    $0x20,%rsi
  add    $0x20,%r9
  movsd  -0x18(%rsi),%xmm6
  movsd  -0x10(%rsi),%xmm5
  add    $0x20,%r8
  movsd  -0x8(%rsi),%xmm4
  mulsd  %xmm0,%xmm7
  mulsd  %xmm1,%xmm6
  addsd  -0x20(%r9),%xmm7
  mulsd  %xmm2,%xmm5
  addsd  -0x18(%r9),%xmm6
  mulsd  %xmm3,%xmm4
  addsd  -0x10(%r9),%xmm5
  addsd  -0x8(%r9),%xmm4
  movsd  %xmm7,-0x20(%r8)
  movsd  %xmm6,-0x18(%r8)
  movsd  %xmm5,-0x10(%r8)
  movsd  %xmm4,-0x8(%r8)
  cmp    %r10d,%edi
  jg     L55
L84:
  ret
  nopl   0x0(%rax)
L12:
  ret
  nopl   0x0(%rax)
L88:
  mov    %r8,%rdx
  jmp    L89
```

tool_calls: name: parse_data
arguments: label: D2
type_infered: byte[8]

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

assistant: ```c
void kernel_dgemm_diag_left_4_lib4(int kmax, double *A, double *B, double *C, double *D, int alg) {
  if (kmax <= 0)
    return;
  const int bs = 4;
  int k;
  double a_0, a_1, a_2, a_3, b_0, b_1, b_2, b_3, c_0, c_1, c_2, c_3;
  if (alg == -1) {
    a_0 = -A[0];
    a_1 = -A[1];
    a_2 = -A[2];
    a_3 = -A[3];
  } else {
    a_0 = A[0];
    a_1 = A[1];
    a_2 = A[2];
    a_3 = A[3];
  }
  if (alg == 0) {
    for (k = 0; k < kmax - 3; k += 4) {
      b_0 = B[0 + bs * 0];
      b_1 = B[1 + bs * 0];
      b_2 = B[2 + bs * 0];
      b_3 = B[3 + bs * 0];
      c_0 = a_0 * b_0;
      c_1 = a_1 * b_1;
      c_2 = a_2 * b_2;
      c_3 = a_3 * b_3;
      D[0 + bs * 0] = c_0;
      D[1 + bs * 0] = c_1;
      D[2 + bs * 0] = c_2;
      D[3 + bs * 0] = c_3;
      b_0 = B[0 + bs * 1];
      b_1 = B[1 + bs * 1];
      b_2 = B[2 + bs * 1];
      b_3 = B[3 + bs * 1];
      c_0 = a_0 * b_0;
      c_1 = a_1 * b_1;
      c_2 = a_2 * b_2;
      c_3 = a_3 * b_3;
      D[0 + bs * 1] = c_0;
      D[1 + bs * 1] = c_1;
      D[2 + bs * 1] = c_2;
      D[3 + bs * 1] = c_3;
      b_0 = B[0 + bs * 2];
      b_1 = B[1 + bs * 2];
      b_2 = B[2 + bs * 2];
      b_3 = B[3 + bs * 2];
      c_0 = a_0 * b_0;
      c_1 = a_1 * b_1;
      c_2 = a_2 * b_2;
      c_3 = a_3 * b_3;
      D[0 + bs * 2] = c_0;
      D[1 + bs * 2] = c_1;
      D[2 + bs * 2] = c_2;
      D[3 + bs * 2] = c_3;
      b_0 = B[0 + bs * 3];
      b_1 = B[1 + bs * 3];
      b_2 = B[2 + bs * 3];
      b_3 = B[3 + bs * 3];
      c_0 = a_0 * b_0;
      c_1 = a_1 * b_1;
      c_2 = a_2 * b_2;
      c_3 = a_3 * b_3;
      D[0 + bs * 3] = c_0;
      D[1 + bs * 3] = c_1;
      D[2 + bs * 3] = c_2;
      D[3 + bs * 3] = c_3;
      B += 16;
      D += 16;
    }
    for (; k < kmax; k++) {
      b_0 = B[0 + bs * 0];
      b_1 = B[1 + bs * 0];
      b_2 = B[2 + bs * 0];
      b_3 = B[3 + bs * 0];
      c_0 = a_0 * b_0;
      c_1 = a_1 * b_1;
      c_2 = a_2 * b_2;
      c_3 = a_3 * b_3;
      D[0 + bs * 0] = c_0;
      D[1 + bs * 0] = c_1;
      D[2 + bs * 0] = c_2;
      D[3 + bs * 0] = c_3;
      B += 4;
      D += 4;
    }
  } else {
    for (k = 0; k < kmax - 3; k += 4) {
      b_0 = B[0 + bs * 0];
      b_1 = B[1 + bs * 0];
      b_2 = B[2 + bs * 0];
      b_3 = B[3 + bs * 0];
      c_0 = C[0 + bs * 0] + a_0 * b_0;
      c_1 = C[1 + bs * 0] + a_1 * b_1;
      c_2 = C[2 + bs * 0] + a_2 * b_2;
      c_3 = C[3 + bs * 0] + a_3 * b_3;
      D[0 + bs * 0] = c_0;
      D[1 + bs * 0] = c_1;
      D[2 + bs * 0] = c_2;
      D[3 + bs * 0] = c_3;
      b_0 = B[0 + bs * 1];
      b_1 = B[1 + bs * 1];
      b_2 = B[2 + bs * 1];
      b_3 = B[3 + bs * 1];
      c_0 = C[0 + bs * 1] + a_0 * b_0;
      c_1 = C[1 + bs * 1] + a_1 * b_1;
      c_2 = C[2 + bs * 1] + a_2 * b_2;
      c_3 = C[3 + bs * 1] + a_3 * b_3;
      D[0 + bs * 1] = c_0;
      D[1 + bs * 1] = c_1;
      D[2 + bs * 1] = c_2;
      D[3 + bs * 1] = c_3;
      b_0 = B[0 + bs * 2];
      b_1 = B[1 + bs * 2];
      b_2 = B[2 + bs * 2];
      b_3 = B[3 + bs * 2];
      c_0 = C[0 + bs * 2] + a_0 * b_0;
      c_1 = C[1 + bs * 2] + a_1 * b_1;
      c_2 = C[2 + bs * 2] + a_2 * b_2;
      c_3 = C[3 + bs * 2] + a_3 * b_3;
      D[0 + bs * 2] = c_0;
      D[1 + bs * 2] = c_1;
      D[2 + bs * 2] = c_2;
      D[3 + bs * 2] = c_3;
      b_0 = B[0 + bs * 3];
      b_1 = B[1 + bs * 3];
      b_2 = B[2 + bs * 3];
      b_3 = B[3 + bs * 3];
      c_0 = C[0 + bs * 3] + a_0 * b_0;
      c_1 = C[1 + bs * 3] + a_1 * b_1;
      c_2 = C[2 + bs * 3] + a_2 * b_2;
      c_3 = C[3 + bs * 3] + a_3 * b_3;
      D[0 + bs * 3] = c_0;
      D[1 + bs * 3] = c_1;
      D[2 + bs * 3] = c_2;
      D[3 + bs * 3] = c_3;
      B += 16;
      C += 16;
      D += 16;
    }
    for (; k < kmax; k++) {
      b_0 = B[0 + bs * 0];
      b_1 = B[1 + bs * 0];
      b_2 = B[2 + bs * 0];
      b_3 = B[3 + bs * 0];
      c_0 = C[0 + bs * 0] + a_0 * b_0;
      c_1 = C[1 + bs * 0] + a_1 * b_1;
      c_2 = C[2 + bs * 0] + a_2 * b_2;
      c_3 = C[3 + bs * 0] + a_3 * b_3;
      D[0 + bs * 0] = c_0;
      D[1 + bs * 0] = c_1;
      D[2 + bs * 0] = c_2;
      D[3 + bs * 0] = c_3;
      B += 4;
      C += 4;
      D += 4;
    }
  }
}
```