user: 请将以下的汇编代码转换成C语言代码:

```asm
<qrevec>:
  endbr64
  push   %r15
  pxor   %xmm6,%xmm6
  mov    %rdx,%r15
  xor    %r8d,%r8d
  push   %r14
  movapd %xmm6,%xmm2
  push   %r13
  mov    %rdi,%r13
  push   %r12
  push   %rbp
  push   %rbx
  movslq %ecx,%rbx
  lea    -0x1(%rbx),%ecx
  mov    %rbx,%rdx
  lea    0x0(,%rbx,8),%r9d
  shl    $0x3,%rbx
  mov    %rcx,%r12
  add    $0x1,%rcx
  lea    0x0(,%rcx,8),%rdi
  sub    $0x88,%rsp
  movslq %r12d,%rbp
  movq   D107(%rip),%xmm4
  lea    (%rsi,%rdi,1),%rax
  movsd  D129(%rip),%xmm5
  movsd  D139(%rip),%xmm7
  mov    %rax,0x8(%rsp)
L158:
  shl    $0x3,%rbp
  jmp    L13
  nopl   0x0(%rax)
L68:
  movsd  0x0(%r13,%rbp,1),%xmm1
  lea    -0x1(%r12),%r10d
  lea    -0x8(%rbp),%rax
  mov    %r13,%r14
  movsd  -0x8(%r15,%rbp,1),%xmm3
  lea    0x0(%rbp,%r13,1),%r11
  movapd %xmm1,%xmm0
  andpd  %xmm4,%xmm0
  movapd %xmm3,%xmm8
  mulsd  %xmm5,%xmm0
  andpd  %xmm4,%xmm8
  comisd %xmm8,%xmm0
  jb     L160
  mov    %rax,%rbp
  mov    %r10d,%r12d
L13:
  test   %r12d,%r12d
  jg     L68
  xor    %eax,%eax
L244:
  add    $0x88,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax,%rax,1)
L160:
  add    %r13,%rax
  movapd %xmm3,%xmm0
  movsd  (%rax),%xmm8
  mulsd  %xmm3,%xmm0
  movapd %xmm8,%xmm9
  subsd  %xmm1,%xmm9
  mulsd  %xmm7,%xmm9
  movapd %xmm9,%xmm10
  mulsd  %xmm9,%xmm10
  addsd  %xmm0,%xmm10
  ucomisd %xmm10,%xmm6
  ja     L178
  sqrtsd %xmm10,%xmm10
L31:
  cmp    $0x1,%r12d
  je     L7
  andpd  %xmm4,%xmm8
  movsd  -0x10(%r15,%rbp,1),%xmm0
  mulsd  %xmm5,%xmm8
  andpd  %xmm4,%xmm0
  comisd %xmm8,%xmm0
  jbe    L7
  cmp    %r9d,%r8d
  jg     L208
  comisd %xmm6,%xmm9
  addsd  %xmm9,%xmm1
  jbe    L41
  subsd  %xmm10,%xmm1
  movapd %xmm1,%xmm9
L19:
  mov    %rcx,%r11
  mov    0x8(%rsp),%rbp
  xor    %r10d,%r10d
  movapd %xmm6,%xmm10
  movsd  0x0(%r13),%xmm8
  neg    %r11
  movsd  D13(%rip),%xmm3
  shl    $0x3,%r11
  subsd  %xmm9,%xmm8
  movsd  %xmm8,0x0(%r13)
  nopl   0x0(%rax,%rax,1)
L189:
  mulsd  %xmm3,%xmm8
  movsd  (%r15,%r10,8),%xmm0
  subsd  %xmm10,%xmm8
  movapd %xmm0,%xmm10
  mulsd  %xmm0,%xmm0
  mulsd  %xmm3,%xmm10
  movapd %xmm8,%xmm1
  mulsd  %xmm8,%xmm1
  addsd  %xmm0,%xmm1
  ucomisd %xmm1,%xmm6
  ja     L119
  sqrtsd %xmm1,%xmm1
L200:
  test   %r10,%r10
  je     L27
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,-0x8(%r15,%r10,8)
L27:
  mulsd  %xmm1,%xmm3
  lea    (%r11,%rbp,1),%rax
  movsd  %xmm3,(%r14)
  movsd  (%r15,%r10,8),%xmm2
  movapd %xmm3,%xmm0
  movapd %xmm8,%xmm3
  divsd  %xmm1,%xmm3
  divsd  %xmm1,%xmm2
  movsd  0x8(%r14),%xmm1
  subsd  %xmm9,%xmm1
  movapd %xmm1,%xmm8
  movsd  %xmm1,0x8(%r14)
  mulsd  %xmm2,%xmm10
  mulsd  %xmm2,%xmm8
  addsd  %xmm10,%xmm0
  mulsd  %xmm2,%xmm8
  mulsd  %xmm3,%xmm0
  addsd  %xmm8,%xmm0
  addsd  %xmm9,%xmm0
  movsd  %xmm0,(%r14)
  test   %edx,%edx
  jle    L233
  nopw   0x0(%rax,%rax,1)
L238:
  movsd  (%rax),%xmm0
  movsd  (%rax,%rbx,1),%xmm1
  movapd %xmm3,%xmm8
  add    $0x8,%rax
  mulsd  %xmm0,%xmm8
  mulsd  %xmm2,%xmm1
  mulsd  %xmm2,%xmm0
  addsd  %xmm8,%xmm1
  movsd  %xmm1,-0x8(%rax)
  movsd  -0x8(%rax,%rbx,1),%xmm1
  mulsd  %xmm3,%xmm1
  subsd  %xmm0,%xmm1
  movsd  %xmm1,-0x8(%rax,%rbx,1)
  cmp    %rbp,%rax
  jne    L238
  movsd  0x8(%r14),%xmm1
L233:
  add    $0x1,%r10
  movapd %xmm1,%xmm8
  add    $0x8,%r14
  add    %rbx,%rbp
  cmp    %r10d,%r12d
  jg     L189
  mulsd  %xmm3,%xmm1
  movslq %r12d,%rbp
  test   %r12d,%r12d
  mov    $0x8,%r10d
  lea    0x0(,%rbp,8),%rax
  cmovle %r10,%rax
  add    $0x1,%r8d
  lea    0x0(%r13,%rax,1),%r10
  subsd  %xmm10,%xmm1
  movsd  %xmm1,(%r10)
  mulsd  %xmm2,%xmm1
  movsd  %xmm1,-0x8(%r15,%rax,1)
  mulsd  (%r10),%xmm3
  addsd  %xmm9,%xmm3
  movsd  %xmm3,(%r10)
  jmp    L158
  nop
L7:
  movapd %xmm9,%xmm8
  movsd  D13(%rip),%xmm2
  divsd  %xmm10,%xmm8
  addsd  %xmm2,%xmm8
  mulsd  %xmm7,%xmm8
  ucomisd %xmm8,%xmm6
  ja     L100
  sqrtsd %xmm8,%xmm8
L64:
  ucomisd %xmm6,%xmm8
  jp     L59
  je     L104
L59:
  movapd %xmm8,%xmm0
  movapd %xmm3,%xmm2
  addsd  %xmm8,%xmm0
  mulsd  %xmm10,%xmm0
  divsd  %xmm0,%xmm2
L104:
  imul   %edx,%r10d
  addsd  %xmm1,%xmm9
  sub    $0x2,%r12d
  movapd %xmm9,%xmm0
  addsd  %xmm10,%xmm9
  movslq %r10d,%r10
  subsd  %xmm10,%xmm0
  movsd  %xmm0,(%r11)
  movsd  %xmm9,(%rax)
  lea    (%rsi,%r10,8),%rax
  lea    (%rax,%rdi,1),%r10
  test   %edx,%edx
  jle    L133
  nopl   0x0(%rax)
L217:
  movsd  (%rax),%xmm0
  movsd  (%rax,%rbx,1),%xmm1
  add    $0x8,%rax
  mulsd  %xmm2,%xmm1
  movapd %xmm0,%xmm3
  mulsd  %xmm8,%xmm3
  mulsd  %xmm2,%xmm0
  addsd  %xmm3,%xmm1
  movsd  %xmm1,-0x8(%rax)
  movsd  -0x8(%rax,%rbx,1),%xmm1
  mulsd  %xmm8,%xmm1
  subsd  %xmm0,%xmm1
  movsd  %xmm1,-0x8(%rax,%rbx,1)
  cmp    %r10,%rax
  jne    L217
L133:
  movslq %r12d,%rbp
  jmp    L158
L41:
  movapd %xmm10,%xmm9
  addsd  %xmm1,%xmm9
  jmp    L19
L208:
  mov    $0xffffffff,%eax
  jmp    L244
L119:
  movapd %xmm1,%xmm0
  mov    %rcx,0x70(%rsp)
  mov    %r11,0x68(%rsp)
  mov    %edx,0x60(%rsp)
  mov    %rsi,0x58(%rsp)
  mov    %rdi,0x50(%rsp)
  mov    %r9d,0x48(%rsp)
  mov    %r10,0x40(%rsp)
  mov    %r8d,0x38(%rsp)
  movsd  %xmm3,0x30(%rsp)
  movsd  %xmm10,0x28(%rsp)
  movsd  %xmm8,0x20(%rsp)
  movsd  %xmm9,0x18(%rsp)
  movsd  %xmm2,0x10(%rsp)
  call   <sqrt@plt>
  mov    D139(%rip),%rax
  mov    0x70(%rsp),%rcx
  pxor   %xmm6,%xmm6
  mov    0x68(%rsp),%r11
  mov    0x60(%rsp),%edx
  movapd %xmm0,%xmm1
  movq   %rax,%xmm7
  mov    0x58(%rsp),%rsi
  mov    0x50(%rsp),%rdi
  mov    D129(%rip),%rax
  mov    0x48(%rsp),%r9d
  mov    0x40(%rsp),%r10
  mov    0x38(%rsp),%r8d
  movq   D107(%rip),%xmm4
  movsd  0x30(%rsp),%xmm3
  movq   %rax,%xmm5
  movsd  0x28(%rsp),%xmm10
  movsd  0x20(%rsp),%xmm8
  movsd  0x18(%rsp),%xmm9
  movsd  0x10(%rsp),%xmm2
  jmp    L200
L100:
  movapd %xmm8,%xmm0
  mov    %rcx,0x70(%rsp)
  mov    %edx,0x68(%rsp)
  mov    %rsi,0x60(%rsp)
  mov    %rdi,0x58(%rsp)
  mov    %r9d,0x50(%rsp)
  mov    %r8d,0x48(%rsp)
  mov    %r11,0x40(%rsp)
  mov    %r10d,0x30(%rsp)
  mov    %rax,0x20(%rsp)
  movsd  %xmm9,0x38(%rsp)
  movsd  %xmm10,0x28(%rsp)
  movsd  %xmm1,0x18(%rsp)
  movsd  %xmm3,0x10(%rsp)
  call   <sqrt@plt>
  mov    D139(%rip),%rax
  mov    0x70(%rsp),%rcx
  pxor   %xmm6,%xmm6
  mov    0x68(%rsp),%edx
  mov    0x60(%rsp),%rsi
  movapd %xmm0,%xmm8
  movq   %rax,%xmm7
  mov    0x58(%rsp),%rdi
  mov    0x50(%rsp),%r9d
  mov    D129(%rip),%rax
  mov    0x48(%rsp),%r8d
  movq   D107(%rip),%xmm4
  mov    0x40(%rsp),%r11
  movq   %rax,%xmm5
  mov    D13(%rip),%rax
  movsd  0x38(%rsp),%xmm9
  mov    0x30(%rsp),%r10d
  movsd  0x28(%rsp),%xmm10
  movq   %rax,%xmm2
  movsd  0x18(%rsp),%xmm1
  mov    0x20(%rsp),%rax
  movsd  0x10(%rsp),%xmm3
  jmp    L64
L178:
  movapd %xmm10,%xmm0
  mov    %rcx,0x78(%rsp)
  mov    %edx,0x70(%rsp)
  mov    %rsi,0x68(%rsp)
  mov    %rdi,0x60(%rsp)
  mov    %r9d,0x58(%rsp)
  mov    %r8d,0x50(%rsp)
  mov    %r11,0x48(%rsp)
  mov    %r10d,0x38(%rsp)
  mov    %rax,0x20(%rsp)
  movsd  %xmm9,0x40(%rsp)
  movsd  %xmm2,0x30(%rsp)
  movsd  %xmm8,0x28(%rsp)
  movsd  %xmm1,0x18(%rsp)
  movsd  %xmm3,0x10(%rsp)
  call   <sqrt@plt>
  mov    D139(%rip),%rax
  mov    0x78(%rsp),%rcx
  pxor   %xmm6,%xmm6
  mov    0x70(%rsp),%edx
  mov    0x68(%rsp),%rsi
  movapd %xmm0,%xmm10
  movq   %rax,%xmm7
  mov    0x60(%rsp),%rdi
  mov    0x58(%rsp),%r9d
  mov    D129(%rip),%rax
  mov    0x50(%rsp),%r8d
  movq   D107(%rip),%xmm4
  mov    0x48(%rsp),%r11
  movq   %rax,%xmm5
  mov    0x38(%rsp),%r10d
  movsd  0x40(%rsp),%xmm9
  movsd  0x30(%rsp),%xmm2
  movsd  0x28(%rsp),%xmm8
  mov    0x20(%rsp),%rax
  movsd  0x18(%rsp),%xmm1
  movsd  0x10(%rsp),%xmm3
  jmp    L31
```

tool_calls: name: disassemble_data
arguments: tag: D107
type_guessed: qword
name: disassemble_data
arguments: tag: D129
type_guessed: byte[8]
name: disassemble_data
arguments: tag: D139
type_guessed: f32[2]
name: disassemble_data
arguments: tag: D13
type_guessed: f64

tool: D107: 9223372036854775807

tool: {"D129": [0x16,0x56,0xe7,0x9e,0xaf,0x03,0xd2,0x3c]}

tool: {"D139": [0.0,1.75]}

tool: D13: 1

assistant: ```c
int qrevec(double *ev, double *evec, double *dp, int n) {
  double cc, sc = 0.0, d, x, y, h, tzr = 1.e-15;
  int i, j, k, m, mqr = 8 * n;
  double *p;
  for (j = 0, m = n - 1;; ++j) {
    while (1) {
      if (m < 1)
        return 0;
      k = m - 1;
      if (fabs(dp[k]) <= fabs(ev[m]) * tzr)
        --m;
      else {
        x = (ev[k] - ev[m]) / 2.;
        h = sqrt(x * x + dp[k] * dp[k]);
        if (m > 1 && fabs(dp[m - 2]) > fabs(ev[k]) * tzr)
          break;
        if ((cc = sqrt((1. + x / h) / 2.)) != 0.)
          sc = dp[k] / (2. * cc * h);
        else
          sc = 1.;
        x += ev[m];
        ev[m--] = x - h;
        ev[m--] = x + h;
        for (i = 0, p = evec + n * (m + 1); i < n; ++i, ++p) {
          h = p[0];
          p[0] = cc * h + sc * p[n];
          p[n] = cc * p[n] - sc * h;
        }
      }
    }
    if (j > mqr)
      return -1;
    if (x > 0.)
      d = ev[m] + x - h;
    else
      d = ev[m] + x + h;
    cc = 1.;
    y = 0.;
    ev[0] -= d;
    for (k = 0; k < m; ++k) {
      x = ev[k] * cc - y;
      y = dp[k] * cc;
      h = sqrt(x * x + dp[k] * dp[k]);
      if (k > 0)
        dp[k - 1] = sc * h;
      ev[k] = cc * h;
      cc = x / h;
      sc = dp[k] / h;
      ev[k + 1] -= d;
      y *= sc;
      ev[k] = cc * (ev[k] + y) + ev[k + 1] * sc * sc + d;
      for (i = 0, p = evec + n * k; i < n; ++i, ++p) {
        h = p[0];
        p[0] = cc * h + sc * p[n];
        p[n] = cc * p[n] - sc * h;
      }
    }
    ev[k] = ev[k] * cc - y;
    dp[k - 1] = ev[k] * sc;
    ev[k] = ev[k] * cc + d;
  }
  return 0;
}
```