user: 请从给定的汇编指令推导出C源代码:

```asm
<merge_matrix_3fts>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x18,%rsp
  mov    %rsi,-0x20(%rsp)
  mov    %rdx,-0x18(%rsp)
  mov    %rcx,-0x10(%rsp)
  mov    %r8,-0x8(%rsp)
  mov    %r9,(%rsp)
  lea    (%rdi,%rdi,4),%esi
  lea    (%rdi,%rsi,2),%edx
  lea    (%rdi,%rdi,2),%ecx
  add    %ecx,%ecx
  test   %edi,%edi
  jle    L124
  mov    %edi,%eax
  mov    %esi,0xc(%rsp)
  lea    (%rsi,%rsi,4),%edi
  movslq %edi,%rbx
  mov    %rbx,-0x40(%rsp)
  imul   $0x79,%eax,%edi
  movslq %edi,%rbx
  mov    %rbx,-0x38(%rsp)
  imul   $0x42,%eax,%edi
  movslq %edi,%rbx
  mov    %rbx,-0x60(%rsp)
  imul   $0x1e,%eax,%edi
  movslq %edi,%rbx
  mov    %rbx,-0x30(%rsp)
  lea    (%rax,%rax,8),%edi
  shl    $0x2,%edi
  movslq %edi,%rbx
  mov    %rbx,-0x28(%rsp)
  mov    %ecx,0x10(%rsp)
  movslq %esi,%rsi
  lea    0x0(,%rsi,8),%r9
  movslq %edx,%rdx
  lea    0x0(,%rdx,8),%r10
  movslq %ecx,%rcx
  lea    0x0(,%rcx,8),%r15
  movq   $0x0,-0x48(%rsp)
  movq   $0x0,-0x50(%rsp)
  movl   $0x0,0x8(%rsp)
  movq   $0x0,-0x58(%rsp)
  movq   $0x0,-0x68(%rsp)
  mov    $0x0,%ebx
  lea    0x0(,%rdx,4),%rcx
  sub    %rcx,%rdx
  shl    $0x4,%rdx
  mov    %rdx,-0x70(%rsp)
  lea    (%rax,%rax,4),%eax
  shl    $0x2,%eax
  mov    %eax,0x14(%rsp)
  jmp    L93
L99:
  add    %r9,%rsi
  add    %r10,%rdx
  add    %r15,%rcx
  sub    $0x1,%r8d
  je     L29
L45:
  mov    $0x0,%eax
  mov    %rdx,%rdi
L141:
  movsd  (%rcx,%rax,1),%xmm0
  movsd  %xmm0,(%rdx,%rax,1)
  add    $0x8,%rax
  cmp    $0x30,%rax
  jne    L141
  mov    $0x0,%eax
L37:
  movsd  (%rsi,%rax,1),%xmm0
  movsd  %xmm0,0x30(%rdi,%rax,1)
  add    $0x8,%rax
  cmp    $0x28,%rax
  jne    L37
  jmp    L99
L66:
  add    %r9,%rsi
  add    %r10,%rdx
  add    %r15,%rcx
  sub    $0x1,%r8d
  je     L110
L92:
  mov    $0x0,%eax
  mov    %rdx,%rdi
L98:
  movsd  (%rcx,%rax,1),%xmm0
  movsd  %xmm0,(%rdx,%rax,1)
  add    $0x8,%rax
  cmp    $0x30,%rax
  jne    L98
  mov    $0x0,%eax
L22:
  movsd  (%rsi,%rax,1),%xmm0
  movsd  %xmm0,0x30(%rdi,%rax,1)
  add    $0x8,%rax
  cmp    $0x28,%rax
  jne    L22
  jmp    L66
L29:
  mov    %r13,%rcx
  mov    %r11,%rdx
  mov    %r14,%rsi
  mov    $0x5,%r8d
  jmp    L92
L110:
  add    $0x5,%ebx
  add    $0x28,%r14
  add    $0x58,%r11
  add    $0x30,%r13
  add    $0x28,%r12
  add    $0x30,%rbp
  cmp    -0x74(%rsp),%ebx
  je     L3
L67:
  mov    -0x70(%rsp),%rax
  lea    (%r11,%rax,1),%rdx
  mov    %rbp,%rcx
  mov    %r12,%rsi
  mov    $0x6,%r8d
  jmp    L45
L3:
  mov    -0x74(%rsp),%ebx
  add    0x14(%rsp),%ebx
  mov    -0x40(%rsp),%rdi
  add    %rdi,-0x68(%rsp)
  mov    -0x38(%rsp),%rdi
  add    %rdi,-0x60(%rsp)
  mov    -0x30(%rsp),%rdi
  add    %rdi,-0x58(%rsp)
  addl   $0x6,0x8(%rsp)
  mov    0x8(%rsp),%eax
  add    %rdi,-0x50(%rsp)
  mov    -0x28(%rsp),%rsi
  add    %rsi,-0x48(%rsp)
  cmp    0x10(%rsp),%eax
  je     L124
L93:
  mov    -0x8(%rsp),%rax
  mov    -0x68(%rsp),%rdi
  lea    (%rax,%rdi,8),%r14
  mov    (%rsp),%rax
  mov    -0x60(%rsp),%rdi
  lea    (%rax,%rdi,8),%r11
  mov    -0x10(%rsp),%rax
  mov    -0x58(%rsp),%rdi
  lea    (%rax,%rdi,8),%r13
  mov    -0x18(%rsp),%rax
  mov    -0x50(%rsp),%rdi
  lea    (%rax,%rdi,8),%r12
  mov    -0x20(%rsp),%rax
  mov    -0x48(%rsp),%rdi
  lea    (%rax,%rdi,8),%rbp
  mov    0xc(%rsp),%eax
  add    %ebx,%eax
  mov    %eax,-0x74(%rsp)
  jmp    L67
L124:
  add    $0x18,%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];
        }
      }
    }
  }
}
```