user: 请将以下的汇编代码转换成C语言代码:

```asm
<kernel_dgemm_nt_4x2_vs_lib4>:
  endbr64
  push   %r15
  xor    %eax,%eax
  mov    %edx,%r15d
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %rcx,%rbx
  sub    $0x68,%rsp
  mov    0xa0(%rsp),%rcx
  mov    %edi,0x44(%rsp)
  mov    %esi,0x48(%rsp)
  mov    %rcx,0x28(%rsp)
  mov    0xa8(%rsp),%ecx
  mov    %edx,0x40(%rsp)
  mov    %ecx,0x34(%rsp)
  mov    0xb0(%rsp),%ecx
  mov    %r8,0x20(%rsp)
  mov    %ecx,0x5c(%rsp)
  mov    0xb8(%rsp),%ecx
  mov    %r9,0x38(%rsp)
  mov    %ecx,0x4c(%rsp)
  call   <_mm_setzero_pd@plt>
  mov    %eax,%r14d
  xor    %eax,%eax
  call   <_mm_setzero_pd@plt>
  mov    %eax,%r13d
  xor    %eax,%eax
  call   <_mm_setzero_pd@plt>
  mov    %eax,%r12d
  xor    %eax,%eax
  call   <_mm_setzero_pd@plt>
  mov    %eax,%ebp
  cmp    $0x3,%r15d
  jle    L215
  lea    -0x4(%r15),%eax
  mov    %r12d,0xc(%rsp)
  mov    %ebp,%r9d
  shr    $0x2,%eax
  mov    %r13d,0x30(%rsp)
  add    $0x1,%eax
  mov    %eax,0x58(%rsp)
  shl    $0x7,%rax
  lea    (%rbx,%rax,1),%rcx
  mov    %rax,0x50(%rsp)
  mov    0x20(%rsp),%rax
  mov    %rcx,0x18(%rsp)
  mov    %rax,(%rsp)
  nopl   0x0(%rax,%rax,1)
L156:
  mov    %rbx,%rdi
  mov    %r9d,0x10(%rsp)
  call   <_mm_load_pd@plt>
  lea    0x10(%rbx),%rdi
  mov    %eax,%r12d
  call   <_mm_load_pd@plt>
  mov    (%rsp),%rdi
  mov    %eax,%r15d
  call   <_mm_load_pd@plt>
  mov    $0x1,%edx
  mov    %eax,%r13d
  mov    %eax,%esi
  mov    %eax,%edi
  call   <_mm_shuffle_pd@plt>
  mov    %r13d,%esi
  mov    %r12d,%edi
  mov    %eax,%ebp
  call   <_mm_mul_pd@plt>
  mov    %r14d,%edi
  mov    %eax,%esi
  call   <_mm_add_pd@plt>
  mov    %r13d,%esi
  mov    %r15d,%edi
  mov    %eax,%r14d
  call   <_mm_mul_pd@plt>
  mov    0xc(%rsp),%edi
  mov    %eax,%esi
  call   <_mm_add_pd@plt>
  mov    %ebp,%esi
  mov    %r12d,%edi
  mov    %eax,0xc(%rsp)
  call   <_mm_mul_pd@plt>
  mov    0x30(%rsp),%edi
  mov    %eax,%esi
  call   <_mm_add_pd@plt>
  mov    %ebp,%esi
  mov    %r15d,%edi
  mov    %eax,0x30(%rsp)
  call   <_mm_mul_pd@plt>
  mov    0x10(%rsp),%edi
  mov    %eax,%esi
  call   <_mm_add_pd@plt>
  lea    0x20(%rbx),%rdi
  mov    %eax,0x10(%rsp)
  call   <_mm_load_pd@plt>
  lea    0x30(%rbx),%rdi
  mov    %eax,%r13d
  call   <_mm_load_pd@plt>
  mov    %eax,%ebp
  mov    (%rsp),%rax
  lea    0x20(%rax),%rdi
  call   <_mm_load_pd@plt>
  mov    $0x1,%edx
  mov    %eax,%r15d
  mov    %eax,%esi
  mov    %eax,%edi
  call   <_mm_shuffle_pd@plt>
  mov    %r15d,%esi
  mov    %r13d,%edi
  mov    %eax,%r12d
  call   <_mm_mul_pd@plt>
  mov    %r14d,%edi
  mov    %eax,%esi
  call   <_mm_add_pd@plt>
  mov    %r15d,%esi
  mov    %ebp,%edi
  mov    %eax,%r14d
  call   <_mm_mul_pd@plt>
  mov    0xc(%rsp),%edi
  mov    %eax,%esi
  call   <_mm_add_pd@plt>
  mov    %r12d,%esi
  mov    %r13d,%edi
  mov    %eax,0xc(%rsp)
  call   <_mm_mul_pd@plt>
  mov    0x30(%rsp),%edi
  mov    %eax,%esi
  call   <_mm_add_pd@plt>
  mov    %r12d,%esi
  mov    %ebp,%edi
  mov    %eax,0x30(%rsp)
  call   <_mm_mul_pd@plt>
  mov    0x10(%rsp),%edi
  mov    %eax,%esi
  call   <_mm_add_pd@plt>
  lea    0x40(%rbx),%rdi
  mov    %eax,0x10(%rsp)
  call   <_mm_load_pd@plt>
  lea    0x50(%rbx),%rdi
  mov    %eax,%r13d
  call   <_mm_load_pd@plt>
  mov    %eax,%ebp
  mov    (%rsp),%rax
  lea    0x40(%rax),%rdi
  call   <_mm_load_pd@plt>
  mov    $0x1,%edx
  mov    %eax,%r15d
  mov    %eax,%esi
  mov    %eax,%edi
  call   <_mm_shuffle_pd@plt>
  mov    %r15d,%esi
  mov    %r13d,%edi
  mov    %eax,%r12d
  call   <_mm_mul_pd@plt>
  mov    %r14d,%edi
  mov    %eax,%esi
  call   <_mm_add_pd@plt>
  mov    %r15d,%esi
  mov    %ebp,%edi
  mov    %eax,%r14d
  call   <_mm_mul_pd@plt>
  mov    0xc(%rsp),%edi
  mov    %eax,%esi
  call   <_mm_add_pd@plt>
  mov    %r12d,%esi
  mov    %r13d,%edi
  mov    %eax,0xc(%rsp)
  call   <_mm_mul_pd@plt>
  mov    0x30(%rsp),%edi
  mov    %eax,%esi
  call   <_mm_add_pd@plt>
  mov    %r12d,%esi
  mov    %ebp,%edi
  mov    %eax,0x30(%rsp)
  call   <_mm_mul_pd@plt>
  mov    0x10(%rsp),%edi
  mov    %eax,%esi
  call   <_mm_add_pd@plt>
  lea    0x60(%rbx),%rdi
  mov    %eax,0x10(%rsp)
  call   <_mm_load_pd@plt>
  lea    0x70(%rbx),%rdi
  sub    $0xffffffffffffff80,%rbx
  mov    %eax,%r13d
  call   <_mm_load_pd@plt>
  mov    %eax,%ebp
  mov    (%rsp),%rax
  lea    0x60(%rax),%rdi
  call   <_mm_load_pd@plt>
  mov    $0x1,%edx
  mov    %eax,%r12d
  mov    %eax,%esi
  mov    %eax,%edi
  call   <_mm_shuffle_pd@plt>
  mov    %r12d,%esi
  mov    %r13d,%edi
  mov    %eax,%r15d
  call   <_mm_mul_pd@plt>
  mov    %r14d,%edi
  mov    %eax,%esi
  call   <_mm_add_pd@plt>
  mov    %r12d,%esi
  mov    %ebp,%edi
  mov    %eax,%r14d
  call   <_mm_mul_pd@plt>
  mov    0xc(%rsp),%edi
  mov    %eax,%esi
  call   <_mm_add_pd@plt>
  mov    %r15d,%esi
  mov    %r13d,%edi
  mov    %eax,0xc(%rsp)
  call   <_mm_mul_pd@plt>
  mov    0x30(%rsp),%edi
  mov    %eax,%esi
  call   <_mm_add_pd@plt>
  mov    %r15d,%esi
  mov    %ebp,%edi
  mov    %eax,0x30(%rsp)
  call   <_mm_mul_pd@plt>
  mov    0x10(%rsp),%edi
  mov    %eax,%esi
  call   <_mm_add_pd@plt>
  subq   $0xffffffffffffff80,(%rsp)
  mov    %eax,%r9d
  cmp    0x18(%rsp),%rbx
  jne    L156
  mov    0x50(%rsp),%rcx
  add    %rcx,0x20(%rsp)
  mov    %eax,%ebp
  mov    0x58(%rsp),%eax
  mov    0xc(%rsp),%r12d
  mov    0x30(%rsp),%r13d
  shl    $0x2,%eax
L5:
  mov    0x40(%rsp),%ecx
  cmp    %eax,%ecx
  jle    L167
  mov    %ecx,%ebx
  mov    %r12d,(%rsp)
  mov    %ebp,%r8d
  sub    $0x1,%ebx
  mov    %r13d,0xc(%rsp)
  sub    %eax,%ebx
  add    $0x1,%rbx
  shl    $0x5,%rbx
  mov    %rbx,%rax
  mov    0x18(%rsp),%rbx
  add    %rbx,%rax
  mov    %rax,0x10(%rsp)
L186:
  mov    %rbx,%rdi
  mov    %r8d,0x30(%rsp)
  call   <_mm_load_pd@plt>
  lea    0x10(%rbx),%rdi
  add    $0x20,%rbx
  mov    %eax,%r13d
  call   <_mm_load_pd@plt>
  mov    0x20(%rsp),%rdi
  mov    %eax,%ebp
  call   <_mm_load_pd@plt>
  mov    $0x1,%edx
  mov    %eax,%r15d
  mov    %eax,%esi
  mov    %eax,%edi
  call   <_mm_shuffle_pd@plt>
  mov    %r15d,%esi
  mov    %r13d,%edi
  mov    %eax,%r12d
  call   <_mm_mul_pd@plt>
  mov    %r14d,%edi
  mov    %eax,%esi
  call   <_mm_add_pd@plt>
  mov    %r15d,%esi
  mov    %ebp,%edi
  mov    %eax,%r14d
  call   <_mm_mul_pd@plt>
  mov    (%rsp),%edi
  mov    %eax,%esi
  call   <_mm_add_pd@plt>
  mov    %r12d,%esi
  mov    %r13d,%edi
  mov    %eax,(%rsp)
  call   <_mm_mul_pd@plt>
  mov    0xc(%rsp),%edi
  mov    %eax,%esi
  call   <_mm_add_pd@plt>
  mov    %r12d,%esi
  mov    %ebp,%edi
  mov    %eax,0xc(%rsp)
  call   <_mm_mul_pd@plt>
  mov    0x30(%rsp),%edi
  mov    %eax,%esi
  call   <_mm_add_pd@plt>
  addq   $0x20,0x20(%rsp)
  mov    %eax,%r8d
  cmp    0x10(%rsp),%rbx
  jne    L186
  mov    (%rsp),%r12d
  mov    0xc(%rsp),%r13d
  mov    %eax,%ebp
L167:
  mov    0x34(%rsp),%edi
  test   %edi,%edi
  jne    L208
  mov    0x4c(%rsp),%esi
  test   %esi,%esi
  jne    L43
  mov    %r14d,%edi
  mov    $0x2,%edx
  mov    %r13d,%esi
  call   <_mm_shuffle_pd@plt>
  mov    %r14d,%esi
  mov    $0x2,%edx
  mov    %r13d,%edi
  mov    %eax,(%rsp)
  call   <_mm_shuffle_pd@plt>
  mov    %ebp,%esi
  mov    $0x2,%edx
  mov    %r12d,%edi
  mov    %eax,%r14d
  call   <_mm_shuffle_pd@plt>
  mov    %ebp,%edi
  mov    $0x2,%edx
  mov    %r12d,%esi
  mov    %eax,%r15d
  call   <_mm_shuffle_pd@plt>
  mov    (%rsp),%r8d
  mov    %eax,%ebp
L199:
  mov    0x28(%rsp),%rbx
  cmpl   $0x3,0x44(%rsp)
  mov    %r8d,%esi
  lea    0x10(%rbx),%r12
  jle    L88
  mov    %rbx,%rdi
  call   <_mm_store_pd@plt>
  mov    %r15d,%esi
  mov    %r12,%rdi
  call   <_mm_store_pd@plt>
  cmpl   $0x1,0x48(%rsp)
  jg     L57
L233:
  add    $0x68,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopw   0x0(%rax,%rax,1)
L208:
  mov    0x38(%rsp),%rbx
  mov    0x5c(%rsp),%ecx
  lea    0x20(%rbx),%r15
  test   %ecx,%ecx
  jne    L231
  mov    %r13d,%esi
  mov    %r14d,%edi
  mov    $0x2,%edx
  call   <_mm_shuffle_pd@plt>
  mov    %r14d,%esi
  mov    %r13d,%edi
  mov    $0x2,%edx
  mov    %eax,0xc(%rsp)
  call   <_mm_shuffle_pd@plt>
  mov    %ebp,%esi
  mov    %r12d,%edi
  mov    $0x2,%edx
  mov    %eax,%r14d
  call   <_mm_shuffle_pd@plt>
  mov    %r12d,%esi
  mov    $0x2,%edx
  mov    %ebp,%edi
  mov    %eax,%r13d
  call   <_mm_shuffle_pd@plt>
  mov    %rbx,%rdi
  mov    %eax,%ebp
  call   <_mm_load_pd@plt>
  lea    0x10(%rbx),%rdi
  mov    %eax,(%rsp)
  call   <_mm_load_pd@plt>
  mov    %r15,%rdi
  mov    %eax,%r12d
  call   <_mm_load_pd@plt>
  lea    0x30(%rbx),%rdi
  mov    %eax,%r15d
  call   <_mm_load_pd@plt>
  mov    (%rsp),%r9d
  cmpl   $0x1,0x34(%rsp)
  mov    0xc(%rsp),%esi
  mov    %eax,%ebx
  mov    %r9d,%edi
  je     L22
  call   <_mm_sub_pd@plt>
  mov    %r14d,%esi
  mov    %r15d,%edi
  mov    %eax,(%rsp)
  call   <_mm_sub_pd@plt>
  mov    %r13d,%esi
  mov    %r12d,%edi
  mov    %eax,%r14d
  call   <_mm_sub_pd@plt>
  mov    %ebp,%esi
  mov    %ebx,%edi
  mov    %eax,%r15d
  call   <_mm_sub_pd@plt>
  mov    (%rsp),%r8d
  mov    %eax,%ebp
L34:
  mov    0x4c(%rsp),%edx
  test   %edx,%edx
  je     L199
  mov    %r8d,%edi
  mov    %r14d,%esi
  mov    %r8d,(%rsp)
  call   <_mm_unpacklo_pd@plt>
  mov    (%rsp),%edi
  mov    %r14d,%esi
  mov    %eax,%ebx
  call   <_mm_unpackhi_pd@plt>
  mov    %ebp,%esi
  mov    %r15d,%edi
  mov    %eax,%r14d
  call   <_mm_unpacklo_pd@plt>
  mov    %ebp,%esi
  mov    %r15d,%edi
  mov    %eax,%r13d
  call   <_mm_unpackhi_pd@plt>
  mov    %eax,%ebp
  jmp    L2
  cs nopw 0x0(%rax,%rax,1)
L231:
  mov    %r13d,%esi
  mov    %r14d,%edi
  call   <_mm_unpacklo_pd@plt>
  mov    %r14d,%esi
  mov    %r13d,%edi
  mov    %eax,%ebx
  call   <_mm_unpackhi_pd@plt>
  mov    %ebp,%esi
  mov    %r12d,%edi
  mov    %eax,%r14d
  call   <_mm_unpacklo_pd@plt>
  mov    %r12d,%esi
  mov    %ebp,%edi
  mov    %eax,%r13d
  call   <_mm_unpackhi_pd@plt>
  mov    0x38(%rsp),%r12
  mov    %eax,%ebp
  mov    %r12,%rdi
  call   <_mm_load_pd@plt>
  mov    %r15,%rdi
  mov    %eax,(%rsp)
  call   <_mm_load_pd@plt>
  lea    0x40(%r12),%rdi
  mov    %r12,0x38(%rsp)
  mov    %eax,%r15d
  call   <_mm_load_pd@plt>
  mov    %eax,%r12d
  mov    0x38(%rsp),%rax
  lea    0x60(%rax),%rdi
  call   <_mm_load_pd@plt>
  mov    (%rsp),%r9d
  cmpl   $0x1,0x34(%rsp)
  mov    %ebx,%esi
  mov    %eax,(%rsp)
  mov    %r9d,%edi
  je     L87
  call   <_mm_sub_pd@plt>
  mov    %r14d,%esi
  mov    %r15d,%edi
  mov    %eax,%ebx
  call   <_mm_sub_pd@plt>
  mov    %r13d,%esi
  mov    %r12d,%edi
  mov    %eax,%r14d
  call   <_mm_sub_pd@plt>
  mov    (%rsp),%edi
  mov    %ebp,%esi
  mov    %eax,%r13d
  call   <_mm_sub_pd@plt>
  mov    %eax,%ebp
L206:
  mov    0x4c(%rsp),%eax
  test   %eax,%eax
  jne    L2
  mov    %r14d,%esi
  mov    %ebx,%edi
  call   <_mm_unpacklo_pd@plt>
  mov    %r14d,%esi
  mov    %ebx,%edi
  mov    %eax,(%rsp)
  call   <_mm_unpackhi_pd@plt>
  mov    %ebp,%esi
  mov    %r13d,%edi
  mov    %eax,%r14d
  call   <_mm_unpacklo_pd@plt>
  mov    %ebp,%esi
  mov    %r13d,%edi
  mov    %eax,%r15d
  call   <_mm_unpackhi_pd@plt>
  mov    (%rsp),%r8d
  mov    %eax,%ebp
  jmp    L199
  nopl   0x0(%rax)
L43:
  mov    %r13d,%esi
  mov    %r14d,%edi
  call   <_mm_unpacklo_pd@plt>
  mov    %r14d,%esi
  mov    %r13d,%edi
  mov    %eax,%ebx
  call   <_mm_unpackhi_pd@plt>
  mov    %ebp,%esi
  mov    %r12d,%edi
  mov    %eax,%r14d
  call   <_mm_unpacklo_pd@plt>
  mov    %ebp,%edi
  mov    %r12d,%esi
  mov    %eax,%r13d
  call   <_mm_unpackhi_pd@plt>
  mov    %eax,%ebp
L2:
  mov    0x28(%rsp),%rdi
  cmpl   $0x1,0x48(%rsp)
  mov    %ebx,%esi
  lea    0x20(%rdi),%r15
  lea    0x40(%rdi),%r12
  jle    L216
  call   <_mm_store_pd@plt>
  mov    %r14d,%esi
  mov    %r15,%rdi
  call   <_mm_store_pd@plt>
  mov    %r13d,%esi
  mov    %r12,%rdi
  call   <_mm_store_pd@plt>
  cmpl   $0x3,0x44(%rsp)
  jle    L233
  mov    0x28(%rsp),%rdi
  add    $0x60,%rdi
  jmp    L157
  nopw   0x0(%rax,%rax,1)
L88:
  mov    0x28(%rsp),%rbx
  mov    %rbx,%rdi
  call   <_mm_store_pd@plt>
  mov    %r15d,%esi
  mov    %r12,%rdi
  call   <_mm_store_sd@plt>
  cmpl   $0x1,0x48(%rsp)
  jle    L233
  lea    0x20(%rbx),%rdi
  mov    %r14d,%esi
  call   <_mm_store_pd@plt>
  lea    0x30(%rbx),%rdi
L178:
  add    $0x68,%rsp
  mov    %ebp,%esi
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    1150 <_mm_store_sd@plt>
L216:
  mov    0x28(%rsp),%rdi
  call   <_mm_store_sd@plt>
  mov    %r14d,%esi
  mov    %r15,%rdi
  call   <_mm_store_sd@plt>
  mov    %r13d,%esi
  mov    %r12,%rdi
  call   <_mm_store_sd@plt>
  cmpl   $0x3,0x44(%rsp)
  jle    L233
  mov    0x28(%rsp),%rdi
  add    $0x60,%rdi
  jmp    L178
  cs nopw 0x0(%rax,%rax,1)
L57:
  lea    0x20(%rbx),%rdi
  mov    %r14d,%esi
  call   <_mm_store_pd@plt>
  lea    0x30(%rbx),%rdi
L157:
  add    $0x68,%rsp
  mov    %ebp,%esi
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    1140 <_mm_store_pd@plt>
L87:
  call   <_mm_add_pd@plt>
  mov    %r14d,%esi
  mov    %r15d,%edi
  mov    %eax,%ebx
  call   <_mm_add_pd@plt>
  mov    %r13d,%esi
  mov    %r12d,%edi
  mov    %eax,%r14d
  call   <_mm_add_pd@plt>
  mov    (%rsp),%edi
  mov    %ebp,%esi
  mov    %eax,%r13d
  call   <_mm_add_pd@plt>
  mov    %eax,%ebp
  jmp    L206
L22:
  call   <_mm_add_pd@plt>
  mov    %r14d,%esi
  mov    %r15d,%edi
  mov    %eax,(%rsp)
  call   <_mm_add_pd@plt>
  mov    %r13d,%esi
  mov    %r12d,%edi
  mov    %eax,%r14d
  call   <_mm_add_pd@plt>
  mov    %ebp,%esi
  mov    %ebx,%edi
  mov    %eax,%r15d
  call   <_mm_add_pd@plt>
  mov    (%rsp),%r8d
  mov    %eax,%ebp
  jmp    L34
  nopl   0x0(%rax)
L215:
  mov    %rbx,0x18(%rsp)
  xor    %eax,%eax
  jmp    L5
```

assistant: ```c
void kernel_dgemm_nt_4x2_vs_lib4(int km, int kn, int kmax, double *A, double *B, double *C, double *D, int alg, int tc, int td) {
  const int ldc = 4;
  int k;
  __m128d c_00_11, c_01_10, c_20_31, c_21_30, a_01, a_23, b_01, b_10, b_temp_0, b_temp_1;
  c_00_11 = _mm_setzero_pd();
  c_01_10 = _mm_setzero_pd();
  c_20_31 = _mm_setzero_pd();
  c_21_30 = _mm_setzero_pd();
  k = 0;
  for (; k < kmax - 3; k += 4) {
    a_01 = _mm_load_pd(&A[0]);
    a_23 = _mm_load_pd(&A[2]);
    b_01 = _mm_load_pd(&B[0]);
    b_10 = _mm_shuffle_pd(b_01, b_01, 1);
    b_temp_0 = b_01;
    b_01 = _mm_mul_pd(a_01, b_01);
    c_00_11 = _mm_add_pd(c_00_11, b_01);
    b_temp_0 = _mm_mul_pd(a_23, b_temp_0);
    c_20_31 = _mm_add_pd(c_20_31, b_temp_0);
    b_temp_1 = b_10;
    b_10 = _mm_mul_pd(a_01, b_10);
    c_01_10 = _mm_add_pd(c_01_10, b_10);
    b_temp_1 = _mm_mul_pd(a_23, b_temp_1);
    c_21_30 = _mm_add_pd(c_21_30, b_temp_1);
    a_01 = _mm_load_pd(&A[4]);
    a_23 = _mm_load_pd(&A[6]);
    b_01 = _mm_load_pd(&B[4]);
    b_10 = _mm_shuffle_pd(b_01, b_01, 1);
    b_temp_0 = b_01;
    b_01 = _mm_mul_pd(a_01, b_01);
    c_00_11 = _mm_add_pd(c_00_11, b_01);
    b_temp_0 = _mm_mul_pd(a_23, b_temp_0);
    c_20_31 = _mm_add_pd(c_20_31, b_temp_0);
    b_temp_1 = b_10;
    b_10 = _mm_mul_pd(a_01, b_10);
    c_01_10 = _mm_add_pd(c_01_10, b_10);
    b_temp_1 = _mm_mul_pd(a_23, b_temp_1);
    c_21_30 = _mm_add_pd(c_21_30, b_temp_1);
    a_01 = _mm_load_pd(&A[8]);
    a_23 = _mm_load_pd(&A[10]);
    b_01 = _mm_load_pd(&B[8]);
    b_10 = _mm_shuffle_pd(b_01, b_01, 1);
    b_temp_0 = b_01;
    b_01 = _mm_mul_pd(a_01, b_01);
    c_00_11 = _mm_add_pd(c_00_11, b_01);
    b_temp_0 = _mm_mul_pd(a_23, b_temp_0);
    c_20_31 = _mm_add_pd(c_20_31, b_temp_0);
    b_temp_1 = b_10;
    b_10 = _mm_mul_pd(a_01, b_10);
    c_01_10 = _mm_add_pd(c_01_10, b_10);
    b_temp_1 = _mm_mul_pd(a_23, b_temp_1);
    c_21_30 = _mm_add_pd(c_21_30, b_temp_1);
    a_01 = _mm_load_pd(&A[12]);
    a_23 = _mm_load_pd(&A[14]);
    b_01 = _mm_load_pd(&B[12]);
    b_10 = _mm_shuffle_pd(b_01, b_01, 1);
    b_temp_0 = b_01;
    b_01 = _mm_mul_pd(a_01, b_01);
    c_00_11 = _mm_add_pd(c_00_11, b_01);
    b_temp_0 = _mm_mul_pd(a_23, b_temp_0);
    c_20_31 = _mm_add_pd(c_20_31, b_temp_0);
    b_temp_1 = b_10;
    b_10 = _mm_mul_pd(a_01, b_10);
    c_01_10 = _mm_add_pd(c_01_10, b_10);
    b_temp_1 = _mm_mul_pd(a_23, b_temp_1);
    c_21_30 = _mm_add_pd(c_21_30, b_temp_1);
    A += 16;
    B += 16;
  }
  for (; k < kmax; k++) {
    a_01 = _mm_load_pd(&A[0]);
    a_23 = _mm_load_pd(&A[2]);
    b_01 = _mm_load_pd(&B[0]);
    b_10 = _mm_shuffle_pd(b_01, b_01, 1);
    b_temp_0 = b_01;
    b_01 = _mm_mul_pd(a_01, b_01);
    c_00_11 = _mm_add_pd(c_00_11, b_01);
    b_temp_0 = _mm_mul_pd(a_23, b_temp_0);
    c_20_31 = _mm_add_pd(c_20_31, b_temp_0);
    b_temp_1 = b_10;
    b_10 = _mm_mul_pd(a_01, b_10);
    c_01_10 = _mm_add_pd(c_01_10, b_10);
    b_temp_1 = _mm_mul_pd(a_23, b_temp_1);
    c_21_30 = _mm_add_pd(c_21_30, b_temp_1);
    A += 4;
    B += 4;
  }
  __m128d t_0, c_0, c_1, c_2, c_3, d_0, d_1, d_2, d_3;
  if (alg == 0) {
    if (td == 0) {
      d_0 = _mm_shuffle_pd(c_00_11, c_01_10, 0x2);
      d_1 = _mm_shuffle_pd(c_01_10, c_00_11, 0x2);
      d_2 = _mm_shuffle_pd(c_20_31, c_21_30, 0x2);
      d_3 = _mm_shuffle_pd(c_21_30, c_20_31, 0x2);
      goto store_n;
    } else {
      d_0 = _mm_unpacklo_pd(c_00_11, c_01_10);
      d_1 = _mm_unpackhi_pd(c_01_10, c_00_11);
      d_2 = _mm_unpacklo_pd(c_20_31, c_21_30);
      d_3 = _mm_unpackhi_pd(c_21_30, c_20_31);
      goto store_t;
    }
  } else {
    if (tc == 0) {
      c_0 = _mm_shuffle_pd(c_00_11, c_01_10, 0x2);
      c_1 = _mm_shuffle_pd(c_01_10, c_00_11, 0x2);
      c_2 = _mm_shuffle_pd(c_20_31, c_21_30, 0x2);
      c_3 = _mm_shuffle_pd(c_21_30, c_20_31, 0x2);
      d_0 = _mm_load_pd(&C[0 + ldc * 0]);
      d_2 = _mm_load_pd(&C[2 + ldc * 0]);
      d_1 = _mm_load_pd(&C[0 + ldc * 1]);
      d_3 = _mm_load_pd(&C[2 + ldc * 1]);
      if (alg == 1) {
        d_0 = _mm_add_pd(d_0, c_0);
        d_1 = _mm_add_pd(d_1, c_1);
        d_2 = _mm_add_pd(d_2, c_2);
        d_3 = _mm_add_pd(d_3, c_3);
      } else {
        d_0 = _mm_sub_pd(d_0, c_0);
        d_1 = _mm_sub_pd(d_1, c_1);
        d_2 = _mm_sub_pd(d_2, c_2);
        d_3 = _mm_sub_pd(d_3, c_3);
      }
      if (td == 0) {
        goto store_n;
      } else {
        t_0 = d_0;
        d_0 = _mm_unpacklo_pd(d_0, d_1);
        d_1 = _mm_unpackhi_pd(t_0, d_1);
        t_0 = d_2;
        d_2 = _mm_unpacklo_pd(d_2, d_3);
        d_3 = _mm_unpackhi_pd(t_0, d_3);
        goto store_t;
      }
    } else {
      c_0 = _mm_unpacklo_pd(c_00_11, c_01_10);
      c_1 = _mm_unpackhi_pd(c_01_10, c_00_11);
      c_2 = _mm_unpacklo_pd(c_20_31, c_21_30);
      c_3 = _mm_unpackhi_pd(c_21_30, c_20_31);
      d_0 = _mm_load_pd(&C[0 + ldc * 0]);
      d_1 = _mm_load_pd(&C[0 + ldc * 1]);
      d_2 = _mm_load_pd(&C[0 + ldc * 2]);
      d_3 = _mm_load_pd(&C[0 + ldc * 3]);
      if (alg == 1) {
        d_0 = _mm_add_pd(d_0, c_0);
        d_1 = _mm_add_pd(d_1, c_1);
        d_2 = _mm_add_pd(d_2, c_2);
        d_3 = _mm_add_pd(d_3, c_3);
      } else {
        d_0 = _mm_sub_pd(d_0, c_0);
        d_1 = _mm_sub_pd(d_1, c_1);
        d_2 = _mm_sub_pd(d_2, c_2);
        d_3 = _mm_sub_pd(d_3, c_3);
      }
      if (td == 0) {
        t_0 = d_0;
        d_0 = _mm_unpacklo_pd(d_0, d_1);
        d_1 = _mm_unpackhi_pd(t_0, d_1);
        t_0 = d_2;
        d_2 = _mm_unpacklo_pd(d_2, d_3);
        d_3 = _mm_unpackhi_pd(t_0, d_3);
        goto store_n;
      } else {
        goto store_t;
      }
    }
  }
store_n:
  if (km >= 4) {
    _mm_store_pd(&D[0 + ldc * 0], d_0);
    _mm_store_pd(&D[2 + ldc * 0], d_2);
    if (kn >= 2) {
      _mm_store_pd(&D[0 + ldc * 1], d_1);
      _mm_store_pd(&D[2 + ldc * 1], d_3);
    }
  } else {
    _mm_store_pd(&D[0 + ldc * 0], d_0);
    _mm_store_sd(&D[2 + ldc * 0], d_2);
    if (kn >= 2) {
      _mm_store_pd(&D[0 + ldc * 1], d_1);
      _mm_store_sd(&D[2 + ldc * 1], d_3);
    }
  }
  return;
store_t:
  if (kn >= 2) {
    _mm_store_pd(&D[0 + ldc * 0], d_0);
    _mm_store_pd(&D[0 + ldc * 1], d_1);
    _mm_store_pd(&D[0 + ldc * 2], d_2);
    if (km >= 4)
      _mm_store_pd(&D[0 + ldc * 3], d_3);
  } else {
    _mm_store_sd(&D[0 + ldc * 0], d_0);
    _mm_store_sd(&D[0 + ldc * 1], d_1);
    _mm_store_sd(&D[0 + ldc * 2], d_2);
    if (km >= 4)
      _mm_store_sd(&D[0 + ldc * 3], d_3);
  }
  return;
}
```