user: Translate the below assembly code into C programming language:

```asm
<qrevec>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x88,%rsp
  mov    %rdi,%r12
  mov    %rsi,0x48(%rsp)
  mov    %rdx,0x20(%rsp)
  mov    %ecx,%eax
  mov    %ecx,0x40(%rsp)
  lea    0x0(,%rcx,8),%edi
  mov    %edi,0x70(%rsp)
  lea    -0x1(%rcx),%r15d
  cltq
  lea    0x0(,%rax,8),%r13
  mov    %r15d,%eax
  mov    %rax,0x60(%rsp)
  lea    0x8(,%rax,8),%rax
  mov    %rax,0x50(%rsp)
  add    %rax,%rsi
  mov    %rsi,0x58(%rsp)
  movl   $0x0,0x44(%rsp)
  pxor   %xmm7,%xmm7
  movsd  %xmm7,(%rsp)
  mov    %r15d,%r14d
L101:
  movslq %r14d,%rbx
  shl    $0x3,%rbx
  mov    0x20(%rsp),%rdx
  jmp    L220
L44:
  mov    %rax,%rbx
  mov    %r15d,%r14d
L220:
  test   %r14d,%r14d
  jle    L175
  lea    -0x1(%r14),%r15d
  lea    -0x8(%rbx),%rax
  movsd  -0x8(%rdx,%rbx,1),%xmm2
  lea    (%rbx,%r12,1),%rbp
  movsd  (%r12,%rbx,1),%xmm3
  movapd %xmm2,%xmm1
  andpd  D65(%rip),%xmm1
  movapd %xmm3,%xmm0
  andpd  D65(%rip),%xmm0
  mulsd  D28(%rip),%xmm0
  comisd %xmm1,%xmm0
  jae    L44
  movsd  %xmm2,0x28(%rsp)
  movsd  %xmm3,0x8(%rsp)
  mov    %r12,0x68(%rsp)
  add    %r12,%rax
  mov    %rax,0x38(%rsp)
  movsd  (%rax),%xmm6
  movsd  %xmm6,0x30(%rsp)
  movapd %xmm6,%xmm0
  subsd  %xmm3,%xmm0
  mulsd  D18(%rip),%xmm0
  movsd  %xmm0,0x18(%rsp)
  mulsd  %xmm0,%xmm0
  movapd %xmm2,%xmm1
  mulsd  %xmm2,%xmm1
  addsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm0,%xmm1
  ja     L24
  sqrtsd %xmm0,%xmm0
  movsd  %xmm0,0x10(%rsp)
L42:
  cmp    $0x1,%r14d
  jle    L0
  mov    0x20(%rsp),%rax
  movsd  -0x10(%rax,%rbx,1),%xmm1
  andpd  D65(%rip),%xmm1
  movsd  0x30(%rsp),%xmm0
  andpd  D65(%rip),%xmm0
  mulsd  D28(%rip),%xmm0
  comisd %xmm0,%xmm1
  ja     L35
L0:
  movsd  0x18(%rsp),%xmm2
  divsd  0x10(%rsp),%xmm2
  addsd  D125(%rip),%xmm2
  mulsd  D18(%rip),%xmm2
  pxor   %xmm0,%xmm0
  ucomisd %xmm2,%xmm0
  ja     L3
  sqrtsd %xmm2,%xmm2
L188:
  ucomisd D21(%rip),%xmm2
  jp     L269
  movsd  D125(%rip),%xmm6
  movsd  %xmm6,(%rsp)
  jne    L269
L199:
  movsd  0x8(%rsp),%xmm0
  addsd  0x18(%rsp),%xmm0
  movapd %xmm0,%xmm1
  movsd  0x10(%rsp),%xmm3
  subsd  %xmm3,%xmm1
  movsd  %xmm1,0x0(%rbp)
  sub    $0x2,%r14d
  addsd  %xmm3,%xmm0
  mov    0x38(%rsp),%rax
  movsd  %xmm0,(%rax)
  mov    0x40(%rsp),%edi
  imul   %edi,%r15d
  movslq %r15d,%r15
  mov    0x48(%rsp),%rax
  lea    (%rax,%r15,8),%rax
  test   %edi,%edi
  jle    L101
  mov    0x50(%rsp),%rdi
  lea    (%rax,%rdi,1),%rdx
L186:
  movsd  (%rax),%xmm0
  movsd  (%rsp),%xmm6
  movapd %xmm6,%xmm1
  mulsd  (%rax,%r13,1),%xmm1
  movapd %xmm0,%xmm3
  mulsd  %xmm2,%xmm3
  addsd  %xmm3,%xmm1
  movsd  %xmm1,(%rax)
  movapd %xmm2,%xmm1
  mulsd  (%rax,%r13,1),%xmm1
  mulsd  %xmm6,%xmm0
  subsd  %xmm0,%xmm1
  movsd  %xmm1,(%rax,%r13,1)
  add    $0x8,%rax
  cmp    %rdx,%rax
  jne    L186
  jmp    L101
L24:
  call   <sqrt@plt>
  movsd  %xmm0,0x10(%rsp)
  jmp    L42
L3:
  movapd %xmm2,%xmm0
  call   <sqrt@plt>
  movapd %xmm0,%xmm2
  jmp    L188
L269:
  movapd %xmm2,%xmm0
  addsd  %xmm2,%xmm0
  mulsd  0x10(%rsp),%xmm0
  movsd  0x28(%rsp),%xmm6
  divsd  %xmm0,%xmm6
  movsd  %xmm6,(%rsp)
  jmp    L199
L35:
  mov    0x70(%rsp),%edi
  cmp    %edi,0x44(%rsp)
  jg     L88
  movsd  0x18(%rsp),%xmm7
  comisd D21(%rip),%xmm7
  jbe    L59
  movsd  0x8(%rsp),%xmm0
  addsd  %xmm7,%xmm0
  subsd  0x10(%rsp),%xmm0
  movsd  %xmm0,0x28(%rsp)
L107:
  movsd  (%r12),%xmm0
  subsd  0x28(%rsp),%xmm0
  movsd  %xmm0,(%r12)
  mov    %r14d,%eax
  mov    %rax,0x30(%rsp)
  mov    0x58(%rsp),%rbx
  mov    $0x0,%ebp
  pxor   %xmm3,%xmm3
  movsd  D125(%rip),%xmm6
  movsd  %xmm6,0x8(%rsp)
  mov    0x60(%rsp),%rax
  not    %rax
  shl    $0x3,%rax
  mov    %rax,0x38(%rsp)
  mov    %r14d,0x74(%rsp)
  mov    %r12,0x78(%rsp)
  mov    0x68(%rsp),%r12
  mov    0x20(%rsp),%r15
  jmp    L25
L59:
  movsd  0x8(%rsp),%xmm0
  addsd  0x18(%rsp),%xmm0
  addsd  0x10(%rsp),%xmm0
  movsd  %xmm0,0x28(%rsp)
  jmp    L107
L10:
  movapd %xmm1,%xmm0
  call   <sqrt@plt>
  movapd %xmm0,%xmm1
  jmp    L33
L17:
  movsd  0x8(%rsp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,(%r14)
  movsd  0x10(%rsp),%xmm6
  divsd  %xmm1,%xmm6
  movsd  %xmm6,0x8(%rsp)
  movsd  (%r15,%rbp,8),%xmm2
  movapd %xmm2,%xmm7
  divsd  %xmm1,%xmm7
  movsd  %xmm7,(%rsp)
  movsd  0x8(%r14),%xmm1
  movsd  0x28(%rsp),%xmm5
  subsd  %xmm5,%xmm1
  movsd  %xmm1,0x8(%r14)
  movsd  0x18(%rsp),%xmm3
  mulsd  %xmm7,%xmm3
  addsd  %xmm3,%xmm0
  mulsd  %xmm6,%xmm0
  mulsd  %xmm7,%xmm1
  mulsd  %xmm7,%xmm1
  addsd  %xmm1,%xmm0
  addsd  %xmm5,%xmm0
  movsd  %xmm0,(%r14)
  mov    0x38(%rsp),%rax
  add    %rbx,%rax
  cmpl   $0x0,0x40(%rsp)
  jle    L135
L85:
  movsd  (%rax),%xmm0
  movsd  (%rsp),%xmm5
  movapd %xmm5,%xmm1
  mulsd  (%rax,%r13,1),%xmm1
  movsd  0x8(%rsp),%xmm4
  movapd %xmm4,%xmm2
  mulsd  %xmm0,%xmm2
  addsd  %xmm2,%xmm1
  movsd  %xmm1,(%rax)
  mulsd  (%rax,%r13,1),%xmm4
  movapd %xmm4,%xmm1
  mulsd  %xmm5,%xmm0
  subsd  %xmm0,%xmm1
  movsd  %xmm1,(%rax,%r13,1)
  add    $0x8,%rax
  cmp    %rbx,%rax
  jne    L85
L135:
  add    $0x1,%rbp
  add    $0x8,%r12
  add    %r13,%rbx
  cmp    0x30(%rsp),%rbp
  je     L187
L25:
  mov    %r12,%r14
  movsd  0x8(%rsp),%xmm7
  movapd %xmm7,%xmm0
  mulsd  (%r12),%xmm0
  movapd %xmm0,%xmm6
  subsd  %xmm3,%xmm6
  movsd  %xmm6,0x10(%rsp)
  movsd  (%r15,%rbp,8),%xmm0
  mulsd  %xmm0,%xmm7
  movsd  %xmm7,0x18(%rsp)
  movapd %xmm6,%xmm1
  mulsd  %xmm6,%xmm1
  mulsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm1
  pxor   %xmm3,%xmm3
  ucomisd %xmm1,%xmm3
  ja     L10
  sqrtsd %xmm1,%xmm1
L33:
  test   %ebp,%ebp
  jle    L17
  movsd  (%rsp),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,-0x8(%r15,%rbp,8)
  jmp    L17
L187:
  mov    0x74(%rsp),%r14d
  mov    0x78(%rsp),%r12
  movslq %r14d,%rdx
  shl    $0x3,%rdx
  lea    (%r12,%rdx,1),%rax
  movsd  0x8(%rsp),%xmm7
  movapd %xmm7,%xmm0
  mulsd  (%rax),%xmm0
  subsd  %xmm3,%xmm0
  movsd  %xmm0,(%rax)
  mulsd  (%rsp),%xmm0
  mov    0x20(%rsp),%rdi
  movsd  %xmm0,-0x8(%rdi,%rdx,1)
  mulsd  (%rax),%xmm7
  movapd %xmm7,%xmm0
  addsd  0x28(%rsp),%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,0x44(%rsp)
  jmp    L101
L175:
  mov    $0x0,%eax
L251:
  add    $0x88,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L88:
  mov    $0xffffffff,%eax
  jmp    L251
```

tool_calls: name: retrieve_data
arguments: data_tag: D65
type_infered: qword[2]
name: retrieve_data
arguments: data_tag: D28
type_infered: byte[8]
name: retrieve_data
arguments: data_tag: D18
type_infered: f32[2]
name: retrieve_data
arguments: data_tag: D125
type_infered: f64
name: retrieve_data
arguments: data_tag: D21
type_infered: f32

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

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

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

tool: D125: 1

tool: D21: 0

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;
}
```