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

```asm
<qrsolve_hh>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x78,%rsp
  mov    %rdi,0x58(%rsp)
  mov    %esi,0x68(%rsp)
  mov    %rdx,0x60(%rsp)
  mov    %rcx,0x50(%rsp)
  mov    %r9d,0x6c(%rsp)
  mov    0xb0(%rsp),%r14d
  test   %r9d,%r9d
  jle    L65
  mov    %r8d,%eax
  movslq %esi,%rdx
  lea    0x8(,%rdx,8),%rsi
  mov    %rsi,0x38(%rsp)
  mov    %rdi,0x28(%rsp)
  lea    0x0(,%rdx,8),%rsi
  mov    %rsi,0x40(%rsp)
  mov    %r9d,%ecx
  mov    %rcx,0x48(%rsp)
  cltq
  shl    $0x3,%rax
  mov    %rax,0x20(%rsp)
  mov    %rdi,0x30(%rsp)
  mov    $0x0,%r15d
L99:
  mov    %r15d,0x1c(%rsp)
  mov    0x30(%rsp),%rbp
  lea    0x0(,%r15,8),%rax
  mov    %rax,0x10(%rsp)
  mov    0x28(%rsp),%rax
  mov    %rax,(%rsp)
  movsd  (%rax),%xmm3
  movsd  %xmm3,0x8(%rsp)
  cmpl   $0x0,0xb8(%rsp)
  jle    L125
  mov    0x50(%rsp),%r12
  mov    $0x0,%r13d
L103:
  mov    %r12,%rbx
  mov    %r14d,%edx
  sub    0x1c(%rsp),%edx
  mov    0x10(%rsp),%rax
  lea    (%r12,%rax,1),%rsi
  mov    (%rsp),%rdi
  call   <vdot@plt>
  xorpd  D5(%rip),%xmm0
  divsd  0x8(%rsp),%xmm0
  cmp    %r15d,%r14d
  jle    L59
  mov    %r15,%rax
L67:
  movapd %xmm0,%xmm1
  mulsd  0x0(%rbp,%rax,8),%xmm1
  addsd  (%rbx,%rax,8),%xmm1
  movsd  %xmm1,(%rbx,%rax,8)
  add    $0x1,%rax
  cmp    %eax,%r14d
  jg     L67
L59:
  add    $0x1,%r13d
  add    0x20(%rsp),%r12
  cmp    %r13d,0xb8(%rsp)
  jne    L103
L125:
  add    $0x1,%r15
  mov    0x38(%rsp),%rcx
  add    %rcx,0x28(%rsp)
  mov    0x40(%rsp),%rcx
  add    %rcx,0x30(%rsp)
  cmp    0x48(%rsp),%r15
  jne    L99
  mov    0x6c(%rsp),%r10d
  sub    $0x1,%r10d
  mov    0x68(%rsp),%eax
  mov    %eax,%r11d
  neg    %r11d
  movslq %r11d,%r11
  shl    $0x3,%r11
  imul   %r10d,%eax
  cltq
  mov    0x58(%rsp),%rcx
  lea    (%rcx,%rax,8),%rsi
  movslq %r10d,%rbp
  shl    $0x3,%rbp
  mov    0x20(%rsp),%r12
  mov    0x60(%rsp),%r14
  mov    0x50(%rsp),%r15
  mov    0xb8(%rsp),%r13d
  jmp    L110
L91:
  cmp    %r13d,%ebx
  jge    L65
L0:
  add    %r12,%r8
  add    $0x1,%ebx
  add    %r12,%r9
L92:
  mov    %r9,%rdx
  mov    %r8,%rcx
  movsd  (%r8),%xmm0
  divsd  (%r14,%rbp,1),%xmm0
  movsd  %xmm0,(%r8)
  test   %r10d,%r10d
  jle    L91
  mov    $0x0,%eax
L27:
  movsd  (%rcx),%xmm1
  mulsd  (%rsi,%rax,8),%xmm1
  movsd  (%rdx,%rax,8),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rdx,%rax,8)
  add    $0x1,%rax
  cmp    %rdi,%rax
  jne    L27
  cmp    %ebx,%r13d
  jg     L0
L129:
  sub    $0x1,%r10d
  add    %r11,%rsi
  sub    $0x8,%rbp
  cmp    $0xffffffff,%r10d
  je     L65
L110:
  test   %r13d,%r13d
  jle    L129
  mov    %r15,%r9
  lea    0x0(%rbp,%r15,1),%r8
  mov    $0x1,%ebx
  mov    %r10d,%edi
  jmp    L92
L65:
  add    $0x78,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: retrieve_data
arguments: data_label: D5
datatype: byte[16]

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

assistant: ```c
void qrsolve_hh(double *qrt, int qrt_stride, double *d, double *bbt, int bbt_stride, int n, int m, int l) {
  int i;
  int j;
  int k;
  for (k = 0; k < n; k++) {
    double *qrtk;
    double e;
    qrtk = qrt + k * qrt_stride;
    e = qrtk[k];
    for (j = 0; j < l; j++) {
      double *bbtj;
      double s;
      bbtj = bbt + j * bbt_stride;
      s = -vdot(qrtk + k, bbtj + k, m - k) / e;
      for (i = k; i < m; i++)
        bbtj[i] += s * qrtk[i];
    }
  }
  for (k = n; k-- > 0;) {
    double *qrtk;
    qrtk = qrt + k * qrt_stride;
    for (j = 0; j < l; j++) {
      double *bbtj;
      bbtj = bbt + j * bbt_stride;
      bbtj[k] /= d[k];
      for (i = 0; i < k; i++)
        bbtj[i] -= bbtj[k] * qrtk[i];
    }
  }
}
```