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

```asm
<Double_QR_Step>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x68,%rsp
  mov    %rdi,0x48(%rsp)
  mov    %esi,0x5c(%rsp)
  mov    %edx,0x3c(%rsp)
  mov    %ecx,0x38(%rsp)
  mov    %r8,0x50(%rsp)
  lea    -0x1(%rdx),%r12d
  mov    %ecx,%eax
  imul   %r9d,%eax
  cltq
  lea    (%rdi,%rax,8),%r14
  movslq %ecx,%rbx
  shl    $0x3,%rbx
  lea    (%r14,%rbx,1),%rdx
  movsd  (%rdx),%xmm3
  mov    %ecx,%edi
  lea    (%rcx,%r9,1),%eax
  movslq %eax,%rcx
  movapd %xmm3,%xmm2
  subsd  %xmm0,%xmm2
  mulsd  %xmm3,%xmm2
  addsd  %xmm2,%xmm1
  divsd  (%r14,%rcx,8),%xmm1
  addsd  0x8(%r14,%rbx,1),%xmm1
  addsd  0x8(%r14,%rcx,8),%xmm3
  subsd  %xmm0,%xmm3
  lea    (%r9,%r9,1),%esi
  mov    %esi,0x58(%rsp)
  add    %r9d,%eax
  cltq
  movsd  0x8(%r14,%rax,8),%xmm4
  movq   D3(%rip),%xmm2
  movapd %xmm1,%xmm0
  andpd  %xmm2,%xmm0
  movapd %xmm3,%xmm5
  andpd  %xmm2,%xmm5
  addsd  %xmm5,%xmm0
  andpd  %xmm4,%xmm2
  addsd  %xmm2,%xmm0
  divsd  %xmm0,%xmm1
  movsd  %xmm1,0x8(%rsp)
  divsd  %xmm0,%xmm3
  movsd  %xmm3,0x18(%rsp)
  divsd  %xmm0,%xmm4
  movsd  %xmm4,0x20(%rsp)
  cmp    %edi,%r12d
  jl     L215
  mov    %r9d,%r15d
  movslq %r9d,%r13
  lea    0x0(,%r13,8),%rax
  mov    %rax,0x28(%rsp)
  add    $0x8,%rax
  mov    %rax,0x40(%rsp)
  mov    %rdx,0x30(%rsp)
  add    $0x8,%rbx
  mov    %edi,%ebp
  jmp    L190
L85:
  pxor   %xmm3,%xmm3
  movsd  %xmm3,0x20(%rsp)
  jmp    L271
L24:
  call   <sqrt@plt>
  jmp    L25
L39:
  xorpd  D11(%rip),%xmm0
  jmp    L152
L149:
  mov    0x5c(%rsp),%esi
  cmp    %esi,0x38(%rsp)
  je     L287
  mov    0x30(%rsp),%rax
  movsd  -0x8(%rax),%xmm1
  xorpd  D11(%rip),%xmm1
  movsd  %xmm1,-0x8(%rax)
  jmp    L287
L83:
  movsd  (%rdx,%r13,8),%xmm0
  movsd  0x8(%rsp),%xmm6
  movapd %xmm6,%xmm1
  mulsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rdx,%r13,8)
  movsd  (%rdx),%xmm0
  mulsd  0x10(%rsp),%xmm6
  subsd  %xmm6,%xmm0
  movsd  %xmm0,(%rdx)
  add    $0x8,%rax
  add    $0x1,%ecx
  cmp    %esi,%ecx
  je     L236
L91:
  mov    %rax,%rdx
  movsd  0x18(%rsp),%xmm0
  mulsd  (%rax,%r13,8),%xmm0
  movapd %xmm0,%xmm4
  addsd  (%rax),%xmm4
  movsd  %xmm4,0x8(%rsp)
  cmp    %ebp,%r12d
  je     L83
  movslq %ecx,%rdi
  lea    (%r14,%rdi,8),%rdi
  movsd  (%rdi),%xmm0
  movsd  0x20(%rsp),%xmm1
  mulsd  %xmm0,%xmm1
  addsd  %xmm4,%xmm1
  movsd  %xmm1,0x8(%rsp)
  mulsd  %xmm3,%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rdi)
  jmp    L83
L236:
  lea    0x3(%rbp),%esi
  mov    0x3c(%rsp),%eax
  cmp    %eax,%esi
  cmovg  %eax,%esi
  test   %esi,%esi
  js     L233
  lea    -0x8(%rbx),%rcx
  lea    0x8(%rbx),%rdi
  add    $0x1,%esi
  mov    0x48(%rsp),%rax
  mov    $0x0,%edx
  jmp    L268
L171:
  movsd  (%rax,%rbx,1),%xmm0
  movsd  0x8(%rsp),%xmm7
  movapd %xmm7,%xmm1
  mulsd  0x18(%rsp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax,%rbx,1)
  movsd  (%rax,%rcx,1),%xmm0
  subsd  %xmm7,%xmm0
  movsd  %xmm0,(%rax,%rcx,1)
  add    $0x1,%edx
  add    0x28(%rsp),%rax
  cmp    %edx,%esi
  je     L233
L268:
  movsd  0x10(%rsp),%xmm1
  mulsd  (%rax,%rcx,1),%xmm1
  movapd %xmm2,%xmm0
  mulsd  (%rax,%rbx,1),%xmm0
  movapd %xmm1,%xmm7
  addsd  %xmm0,%xmm7
  movsd  %xmm7,0x8(%rsp)
  cmp    %ebp,%r12d
  je     L171
  movsd  (%rax,%rdi,1),%xmm0
  movapd %xmm0,%xmm1
  mulsd  %xmm3,%xmm1
  addsd  %xmm7,%xmm1
  movapd %xmm1,%xmm6
  movsd  %xmm1,0x8(%rsp)
  movsd  0x20(%rsp),%xmm1
  mulsd  %xmm6,%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax,%rdi,1)
  jmp    L171
L233:
  test   %r15d,%r15d
  jle    L170
  lea    -0x8(%rbx),%rcx
  lea    0x8(%rbx),%rsi
  mov    0x50(%rsp),%rax
  mov    $0x0,%edx
  jmp    L162
L134:
  movsd  (%rax,%rbx,1),%xmm0
  movsd  0x8(%rsp),%xmm5
  movapd %xmm5,%xmm1
  mulsd  0x18(%rsp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax,%rbx,1)
  movsd  (%rax,%rcx,1),%xmm0
  subsd  %xmm5,%xmm0
  movsd  %xmm0,(%rax,%rcx,1)
  add    0x28(%rsp),%rax
  add    $0x1,%edx
  cmp    %edx,%r15d
  je     L170
L162:
  movsd  0x10(%rsp),%xmm1
  mulsd  (%rax,%rcx,1),%xmm1
  movapd %xmm2,%xmm0
  mulsd  (%rax,%rbx,1),%xmm0
  movapd %xmm1,%xmm5
  addsd  %xmm0,%xmm5
  movsd  %xmm5,0x8(%rsp)
  cmp    %ebp,%r12d
  je     L134
  movsd  (%rax,%rsi,1),%xmm0
  movapd %xmm0,%xmm1
  mulsd  %xmm3,%xmm1
  addsd  %xmm5,%xmm1
  movapd %xmm1,%xmm4
  movsd  %xmm1,0x8(%rsp)
  movsd  0x20(%rsp),%xmm1
  mulsd  %xmm4,%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax,%rsi,1)
  jmp    L134
L170:
  add    $0x1,%ebp
  add    0x28(%rsp),%r14
  mov    0x40(%rsp),%rsi
  add    %rsi,0x30(%rsp)
  add    $0x8,%rbx
  cmp    %ebp,0x3c(%rsp)
  je     L215
L190:
  cmp    %ebp,0x38(%rsp)
  jge    L177
  cmp    %ebp,%r12d
  je     L85
  mov    0x58(%rsp),%eax
  add    %ebp,%eax
  cltq
  movsd  -0x8(%r14,%rax,8),%xmm2
  movsd  %xmm2,0x20(%rsp)
L271:
  mov    0x30(%rsp),%rax
  movsd  -0x8(%rax),%xmm2
  lea    (%r15,%rbp,1),%eax
  cltq
  movsd  -0x8(%r14,%rax,8),%xmm1
  movapd %xmm2,%xmm0
  andpd  D3(%rip),%xmm0
  movapd %xmm1,%xmm3
  andpd  D3(%rip),%xmm3
  addsd  %xmm3,%xmm0
  movsd  0x20(%rsp),%xmm3
  andpd  D3(%rip),%xmm3
  addsd  %xmm3,%xmm0
  movsd  %xmm0,0x10(%rsp)
  ucomisd D58(%rip),%xmm0
  jp     L186
  je     L170
L186:
  movsd  0x10(%rsp),%xmm3
  divsd  %xmm3,%xmm2
  movsd  %xmm2,0x8(%rsp)
  divsd  %xmm3,%xmm1
  movsd  %xmm1,0x18(%rsp)
  movsd  0x20(%rsp),%xmm2
  divsd  %xmm3,%xmm2
  movsd  %xmm2,0x20(%rsp)
L177:
  movsd  0x18(%rsp),%xmm2
  movapd %xmm2,%xmm1
  mulsd  %xmm2,%xmm1
  movsd  0x8(%rsp),%xmm2
  movapd %xmm2,%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  0x20(%rsp),%xmm2
  movapd %xmm2,%xmm1
  mulsd  %xmm2,%xmm1
  addsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm0,%xmm1
  ja     L24
  sqrtsd %xmm0,%xmm0
L25:
  pxor   %xmm1,%xmm1
  comisd 0x8(%rsp),%xmm1
  ja     L39
L152:
  cmp    %ebp,0x38(%rsp)
  jge    L149
  movapd %xmm0,%xmm1
  xorpd  D11(%rip),%xmm1
  mulsd  0x10(%rsp),%xmm1
  mov    0x30(%rsp),%rax
  movsd  %xmm1,-0x8(%rax)
L287:
  movsd  0x8(%rsp),%xmm7
  addsd  %xmm0,%xmm7
  movsd  %xmm7,0x8(%rsp)
  movapd %xmm7,%xmm3
  divsd  %xmm0,%xmm3
  movsd  %xmm3,0x10(%rsp)
  movsd  0x18(%rsp),%xmm5
  movapd %xmm5,%xmm2
  divsd  %xmm0,%xmm2
  movsd  0x20(%rsp),%xmm4
  movapd %xmm4,%xmm3
  divsd  %xmm0,%xmm3
  divsd  %xmm7,%xmm5
  movsd  %xmm5,0x18(%rsp)
  divsd  %xmm7,%xmm4
  movsd  %xmm4,0x20(%rsp)
  cmp    %ebp,%r15d
  jle    L236
  mov    0x58(%rsp),%eax
  lea    0x0(%rbp,%rax,1),%ecx
  lea    (%r15,%rax,1),%esi
  mov    0x30(%rsp),%rax
  jmp    L91
L215:
  add    $0x68,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: read_data
arguments: data_label: D3
type: qword[2]
name: read_data
arguments: data_label: D11
type: f32[4]
name: read_data
arguments: data_label: D58
type: f32[2]

tool: {"D3": [0x7fffffffffffffff,0x0000000000000000]}

tool: {"D11": [0.0,-0.0,0.0,0.0]}

tool: {"D58": [0.0,0.0]}

assistant: ```c
void Double_QR_Step(double *H, int min_row, int max_row, int min_col, double trace, double det, double *S, int n) {
  double s, x, y, z;
  double a, b, c;
  double *pH;
  double *tH;
  double *pS;
  int i, j, k;
  int last_test_row_col = max_row - 1;
  k = min_col;
  pH = H + min_col * n;
  a = (pH[k] * (pH[k] - trace) + det) / pH[n + k] + pH[k + 1];
  b = pH[k] + pH[n + k + 1] - trace;
  c = pH[n + n + k + 1];
  s = fabs(a) + fabs(b) + fabs(c);
  a /= s;
  b /= s;
  c /= s;
  for (; k <= last_test_row_col; k++, pH += n) {
    if (k > min_col) {
      c = (k == last_test_row_col) ? 0.0 : pH[n + n + k - 1];
      x = fabs(pH[k - 1]) + fabs(pH[n + k - 1]) + fabs(c);
      if (x == 0.0)
        continue;
      a = pH[k - 1] / x;
      b = pH[n + k - 1] / x;
      c /= x;
    }
    s = sqrt(a * a + b * b + c * c);
    if (a < 0.0)
      s = -s;
    if (k > min_col)
      pH[k - 1] = -s * x;
    else if (min_row != min_col)
      pH[k - 1] = -pH[k - 1];
    a += s;
    x = a / s;
    y = b / s;
    z = c / s;
    b /= a;
    c /= a;
    for (j = k; j < n; j++) {
      a = pH[j] + b * pH[n + j];
      if (k != last_test_row_col) {
        a += c * pH[n + n + j];
        pH[n + n + j] -= a * z;
      }
      pH[n + j] -= a * y;
      pH[j] -= a * x;
    }
    j = k + 3;
    if (j > max_row)
      j = max_row;
    for (i = 0, tH = H; i <= j; i++, tH += n) {
      a = x * tH[k] + y * tH[k + 1];
      if (k != last_test_row_col) {
        a += z * tH[k + 2];
        tH[k + 2] -= a * c;
      }
      tH[k + 1] -= a * b;
      tH[k] -= a;
    }
    for (i = 0, pS = S; i < n; pS += n, i++) {
      a = x * pS[k] + y * pS[k + 1];
      if (k != last_test_row_col) {
        a += z * pS[k + 2];
        pS[k + 2] -= a * c;
      }
      pS[k + 1] -= a * b;
      pS[k] -= a;
    }
  };
}
```