user: Converte this assembly snippet into its equivalent C code:

```asm
<Double_QR_Step>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x90,%rsp
  mov    %rdi,-0x68(%rbp)
  mov    %esi,-0x6c(%rbp)
  mov    %edx,-0x70(%rbp)
  mov    %ecx,-0x74(%rbp)
  movsd  %xmm0,-0x80(%rbp)
  movsd  %xmm1,-0x88(%rbp)
  mov    %r8,-0x90(%rbp)
  mov    %r9d,-0x78(%rbp)
  mov    -0x70(%rbp),%eax
  sub    $0x1,%eax
  mov    %eax,-0x54(%rbp)
  mov    -0x74(%rbp),%eax
  mov    %eax,-0x58(%rbp)
  mov    -0x74(%rbp),%eax
  imul   -0x78(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x28(%rbp)
  mov    -0x58(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x58(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  subsd  -0x80(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  -0x88(%rbp),%xmm0
  mov    -0x78(%rbp),%edx
  mov    -0x58(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  movapd %xmm0,%xmm1
  divsd  %xmm2,%xmm1
  mov    -0x58(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  mov    -0x58(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x78(%rbp),%edx
  mov    -0x58(%rbp),%eax
  add    %edx,%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  addsd  %xmm1,%xmm0
  subsd  -0x80(%rbp),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  mov    -0x78(%rbp),%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x58(%rbp),%eax
  add    %edx,%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  -0x40(%rbp),%xmm0
  movq   D3(%rip),%xmm1
  andpd  %xmm0,%xmm1
  movsd  -0x38(%rbp),%xmm0
  movq   D3(%rip),%xmm2
  andpd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x30(%rbp),%xmm0
  movq   D3(%rip),%xmm2
  andpd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  movsd  -0x40(%rbp),%xmm0
  divsd  -0x50(%rbp),%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movsd  -0x38(%rbp),%xmm0
  divsd  -0x50(%rbp),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x30(%rbp),%xmm0
  divsd  -0x50(%rbp),%xmm0
  movsd  %xmm0,-0x30(%rbp)
  jmp    L115
L33:
  mov    -0x58(%rbp),%eax
  cmp    -0x74(%rbp),%eax
  jle    L126
  mov    -0x58(%rbp),%eax
  cmp    -0x54(%rbp),%eax
  je     L19
  mov    -0x78(%rbp),%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x58(%rbp),%eax
  add    %edx,%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  jmp    L176
L19:
  pxor   %xmm0,%xmm0
L176:
  movsd  %xmm0,-0x30(%rbp)
  mov    -0x58(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movq   D3(%rip),%xmm1
  andpd  %xmm0,%xmm1
  mov    -0x78(%rbp),%edx
  mov    -0x58(%rbp),%eax
  add    %edx,%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movq   D3(%rip),%xmm2
  andpd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x30(%rbp),%xmm0
  movq   D3(%rip),%xmm2
  andpd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  pxor   %xmm0,%xmm0
  ucomisd -0x48(%rbp),%xmm0
  jp     L181
  pxor   %xmm0,%xmm0
  ucomisd -0x48(%rbp),%xmm0
  je     L83
L181:
  mov    -0x58(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  divsd  -0x48(%rbp),%xmm0
  movsd  %xmm0,-0x40(%rbp)
  mov    -0x78(%rbp),%edx
  mov    -0x58(%rbp),%eax
  add    %edx,%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  divsd  -0x48(%rbp),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x30(%rbp),%xmm0
  divsd  -0x48(%rbp),%xmm0
  movsd  %xmm0,-0x30(%rbp)
L126:
  movsd  -0x40(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  %xmm0,%xmm1
  movsd  -0x38(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x30(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x50(%rbp)
  pxor   %xmm0,%xmm0
  comisd -0x40(%rbp),%xmm0
  jbe    L11
  movsd  -0x50(%rbp),%xmm0
  movq   D74(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
L11:
  mov    -0x58(%rbp),%eax
  cmp    -0x74(%rbp),%eax
  jle    L128
  movsd  -0x50(%rbp),%xmm0
  movq   D74(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  mov    -0x58(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  mulsd  -0x48(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  jmp    L140
L128:
  mov    -0x6c(%rbp),%eax
  cmp    -0x74(%rbp),%eax
  je     L140
  mov    -0x58(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x58(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movq   D74(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
L140:
  movsd  -0x40(%rbp),%xmm0
  addsd  -0x50(%rbp),%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movsd  -0x40(%rbp),%xmm0
  divsd  -0x50(%rbp),%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movsd  -0x38(%rbp),%xmm0
  divsd  -0x50(%rbp),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x30(%rbp),%xmm0
  divsd  -0x50(%rbp),%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x38(%rbp),%xmm0
  divsd  -0x40(%rbp),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x30(%rbp),%xmm0
  divsd  -0x40(%rbp),%xmm0
  movsd  %xmm0,-0x30(%rbp)
  mov    -0x58(%rbp),%eax
  mov    %eax,-0x5c(%rbp)
  jmp    L192
L182:
  mov    -0x5c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x78(%rbp),%edx
  mov    -0x5c(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x38(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  mov    -0x58(%rbp),%eax
  cmp    -0x54(%rbp),%eax
  je     L96
  mov    -0x78(%rbp),%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x5c(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x30(%rbp),%xmm0
  movsd  -0x40(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  mov    -0x78(%rbp),%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x5c(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x40(%rbp),%xmm1
  mulsd  -0x8(%rbp),%xmm1
  mov    -0x78(%rbp),%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x5c(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
L96:
  mov    -0x78(%rbp),%edx
  mov    -0x5c(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x40(%rbp),%xmm1
  mulsd  -0x10(%rbp),%xmm1
  mov    -0x78(%rbp),%edx
  mov    -0x5c(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x5c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x40(%rbp),%xmm1
  mulsd  -0x48(%rbp),%xmm1
  mov    -0x5c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x5c(%rbp)
L192:
  mov    -0x5c(%rbp),%eax
  cmp    -0x78(%rbp),%eax
  jl     L182
  mov    -0x58(%rbp),%eax
  add    $0x3,%eax
  mov    %eax,-0x5c(%rbp)
  mov    -0x5c(%rbp),%eax
  cmp    -0x70(%rbp),%eax
  jle    L106
  mov    -0x70(%rbp),%eax
  mov    %eax,-0x5c(%rbp)
L106:
  movl   $0x0,-0x60(%rbp)
  mov    -0x68(%rbp),%rax
  mov    %rax,-0x20(%rbp)
  jmp    L25
L196:
  mov    -0x58(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x48(%rbp),%xmm1
  mov    -0x58(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x10(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  mov    -0x58(%rbp),%eax
  cmp    -0x54(%rbp),%eax
  je     L163
  mov    -0x58(%rbp),%eax
  cltq
  add    $0x2,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x8(%rbp),%xmm0
  movsd  -0x40(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  mov    -0x58(%rbp),%eax
  cltq
  add    $0x2,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x40(%rbp),%xmm1
  mulsd  -0x30(%rbp),%xmm1
  mov    -0x58(%rbp),%eax
  cltq
  add    $0x2,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
L163:
  mov    -0x58(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x40(%rbp),%xmm1
  mulsd  -0x38(%rbp),%xmm1
  mov    -0x58(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x58(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x58(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  subsd  -0x40(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x60(%rbp)
  mov    -0x78(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0x20(%rbp)
L25:
  mov    -0x60(%rbp),%eax
  cmp    -0x5c(%rbp),%eax
  jle    L196
  movl   $0x0,-0x60(%rbp)
  mov    -0x90(%rbp),%rax
  mov    %rax,-0x18(%rbp)
  jmp    L74
L216:
  mov    -0x58(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x48(%rbp),%xmm1
  mov    -0x58(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x10(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  mov    -0x58(%rbp),%eax
  cmp    -0x54(%rbp),%eax
  je     L218
  mov    -0x58(%rbp),%eax
  cltq
  add    $0x2,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x8(%rbp),%xmm0
  movsd  -0x40(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  mov    -0x58(%rbp),%eax
  cltq
  add    $0x2,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x40(%rbp),%xmm1
  mulsd  -0x30(%rbp),%xmm1
  mov    -0x58(%rbp),%eax
  cltq
  add    $0x2,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
L218:
  mov    -0x58(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x40(%rbp),%xmm1
  mulsd  -0x38(%rbp),%xmm1
  mov    -0x58(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x58(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x58(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  subsd  -0x40(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x78(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0x18(%rbp)
  addl   $0x1,-0x60(%rbp)
L74:
  mov    -0x60(%rbp),%eax
  cmp    -0x78(%rbp),%eax
  jl     L216
  jmp    L82
L83:
  nop
L82:
  addl   $0x1,-0x58(%rbp)
  mov    -0x78(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0x28(%rbp)
L115:
  mov    -0x58(%rbp),%eax
  cmp    -0x54(%rbp),%eax
  jle    L33
  nop
  nop
  leave
  ret
```

tool_calls: name: retrieve_data
arguments: data_label: D3
data_type: qword
name: retrieve_data
arguments: data_label: D74
data_type: f32[2]

tool: D3: 9223372036854775807

tool: {"D74": [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;
    }
  };
}
```