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

```asm
<merge_matrix_3fts>:
  endbr64
  push   %r15
  lea    (%rdi,%rdi,4),%eax
  mov    %edi,%r11d
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x8,%rsp
  mov    %r8,-0x8(%rsp)
  lea    (%rdi,%rax,2),%r8d
  lea    (%rdi,%rdi,2),%edi
  mov    %rsi,-0x20(%rsp)
  add    %edi,%edi
  mov    %rdx,-0x18(%rsp)
  mov    %rcx,-0x10(%rsp)
  mov    %r9,(%rsp)
  test   %r11d,%r11d
  jle    L39
  imul   $0x1e,%r11d,%edx
  movslq %edi,%rdi
  movslq %r8d,%r8
  movq   $0x0,-0x70(%rsp)
  movq   $0x0,-0x68(%rsp)
  shl    $0x3,%rdi
  shl    $0x3,%r8
  movq   $0x0,-0x60(%rsp)
  movslq %edx,%rbx
  imul   $0x79,%r11d,%edx
  movq   $0x0,-0x58(%rsp)
  mov    %rbx,-0x28(%rsp)
  movq   $0x0,-0x48(%rsp)
  movslq %edx,%rbx
  imul   $0x42,%r11d,%edx
  movl   $0x0,-0x74(%rsp)
  mov    %rbx,-0x30(%rsp)
  movslq %edx,%rbx
  lea    (%rax,%rax,4),%edx
  cltq
  mov    %rbx,-0x50(%rsp)
  movslq %edx,%rbx
  lea    (%r11,%r11,8),%edx
  lea    0x0(,%rax,8),%rsi
  shl    $0x2,%edx
  mov    %rbx,-0x38(%rsp)
  movslq %edx,%rbx
  mov    %rbx,-0x40(%rsp)
L31:
  mov    -0x10(%rsp),%rax
  mov    -0x48(%rsp),%rbx
  xor    %r13d,%r13d
  mov    -0x58(%rsp),%rcx
  mov    -0x68(%rsp),%r14
  lea    (%rax,%rbx,8),%r12
  mov    (%rsp),%rbx
  mov    -0x50(%rsp),%rax
  mov    -0x70(%rsp),%rdx
  lea    (%rbx,%rax,8),%r9
  mov    -0x8(%rsp),%rbx
  sub    %rax,%r14
  shl    $0x3,%r14
  lea    (%rbx,%rcx,8),%rbp
  mov    -0x60(%rsp),%rcx
  mov    -0x20(%rsp),%rbx
  lea    (%rbx,%rcx,8),%rbx
  mov    -0x18(%rsp),%rcx
  lea    (%rcx,%rdx,8),%r10
  nopl   0x0(%rax)
L42:
  lea    (%r14,%r9,1),%rax
  mov    %r10,%rcx
  mov    %rbx,%rdx
  mov    $0x6,%r15d
  nopl   0x0(%rax,%rax,1)
L10:
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  movsd  0x8(%rdx),%xmm0
  movsd  %xmm0,0x8(%rax)
  movsd  0x10(%rdx),%xmm0
  movsd  %xmm0,0x10(%rax)
  movsd  0x18(%rdx),%xmm0
  movsd  %xmm0,0x18(%rax)
  movsd  0x20(%rdx),%xmm0
  movsd  %xmm0,0x20(%rax)
  movsd  0x28(%rdx),%xmm0
  add    %rdi,%rdx
  movsd  %xmm0,0x28(%rax)
  movsd  (%rcx),%xmm0
  movsd  %xmm0,0x30(%rax)
  movsd  0x8(%rcx),%xmm0
  movsd  %xmm0,0x38(%rax)
  movsd  0x10(%rcx),%xmm0
  movsd  %xmm0,0x40(%rax)
  movsd  0x18(%rcx),%xmm0
  movsd  %xmm0,0x48(%rax)
  movsd  0x20(%rcx),%xmm0
  add    %rsi,%rcx
  movsd  %xmm0,0x50(%rax)
  add    %r8,%rax
  sub    $0x1,%r15d
  jne    L10
  mov    %rbp,%rcx
  mov    %r9,%rax
  mov    %r12,%rdx
  mov    $0x5,%r15d
L57:
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  movsd  0x8(%rdx),%xmm0
  movsd  %xmm0,0x8(%rax)
  movsd  0x10(%rdx),%xmm0
  movsd  %xmm0,0x10(%rax)
  movsd  0x18(%rdx),%xmm0
  movsd  %xmm0,0x18(%rax)
  movsd  0x20(%rdx),%xmm0
  movsd  %xmm0,0x20(%rax)
  movsd  0x28(%rdx),%xmm0
  add    %rdi,%rdx
  movsd  %xmm0,0x28(%rax)
  movsd  (%rcx),%xmm0
  movsd  %xmm0,0x30(%rax)
  movsd  0x8(%rcx),%xmm0
  movsd  %xmm0,0x38(%rax)
  movsd  0x10(%rcx),%xmm0
  movsd  %xmm0,0x40(%rax)
  movsd  0x18(%rcx),%xmm0
  movsd  %xmm0,0x48(%rax)
  movsd  0x20(%rcx),%xmm0
  add    %rsi,%rcx
  movsd  %xmm0,0x50(%rax)
  add    %r8,%rax
  sub    $0x1,%r15d
  jne    L57
  lea    0x1(%r13),%eax
  add    $0x30,%r12
  add    $0x58,%r9
  add    $0x28,%rbp
  add    $0x30,%rbx
  add    $0x28,%r10
  cmp    %eax,%r11d
  je     L58
  mov    %eax,%r13d
  jmp    L42
L58:
  mov    -0x74(%rsp),%ebx
  mov    -0x28(%rsp),%rdx
  mov    -0x30(%rsp),%r14
  mov    -0x38(%rsp),%r10
  add    %rdx,-0x48(%rsp)
  lea    0x1(%rbx),%eax
  add    %r10,-0x58(%rsp)
  mov    -0x40(%rsp),%r10
  add    %r14,-0x50(%rsp)
  add    %r10,-0x60(%rsp)
  add    %r14,-0x68(%rsp)
  add    %rdx,-0x70(%rsp)
  cmp    %ebx,%r13d
  je     L39
  mov    %eax,-0x74(%rsp)
  jmp    L31
L39:
  add    $0x8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void merge_matrix_3fts(int np, const double *mat_ll, const double *mat_lh, const double *mat_hl, const double *mat_hh, double *mat) {
  int i, j;
  int ii, jj;
  int i11, i6, i5;
  int j11, j6, j5;
  int n11, n6, n5;
  n11 = np * 11;
  n6 = np * 6;
  n5 = np * 5;
  for (i = 0; i < np; ++i) {
    i11 = i * 11;
    i6 = i * 6;
    i5 = i * 5;
    for (j = 0; j < np; ++j) {
      j11 = j * 11;
      j6 = j * 6;
      j5 = j * 5;
      for (ii = 0; ii < 6; ++ii) {
        for (jj = 0; jj < 6; ++jj) {
          mat[(i11 + ii) * n11 + j11 + jj] = mat_ll[(i6 + ii) * n6 + j6 + jj];
        }
        for (jj = 0; jj < 5; ++jj) {
          mat[(i11 + ii) * n11 + j11 + 6 + jj] = mat_lh[(i6 + ii) * n5 + j5 + jj];
        }
      }
      for (ii = 0; ii < 5; ++ii) {
        for (jj = 0; jj < 6; ++jj) {
          mat[(i11 + 6 + ii) * n11 + j11 + jj] = mat_hl[(i5 + ii) * n6 + j6 + jj];
        }
        for (jj = 0; jj < 5; ++jj) {
          mat[(i11 + 6 + ii) * n11 + j11 + 6 + jj] = mat_hh[(i5 + ii) * n5 + j5 + jj];
        }
      }
    }
  }
}
```