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

```asm
<QLalgo2>:
  endbr64
  push   %r15
  movslq %edi,%rax
  mov    %rsi,%r15
  push   %r14
  mov    %rax,%r14
  push   %r13
  mov    %rcx,%r13
  push   %r12
  push   %rbp
  mov    %rdx,%rbp
  push   %rbx
  lea    -0x8(%rdx,%rax,8),%rbx
  sub    $0xa8,%rsp
  cmp    $0x1,%eax
  jle    L159
  lea    -0x1(%rax),%edx
  lea    0x8(%rbp),%rsi
  mov    %rbp,%rdi
  shl    $0x3,%rdx
  call   <memmove@plt>
  movq   $0x0,(%rbx)
L78:
  movq   $0x0,0x38(%rsp)
  mov    %r13,%rbx
  xor    %r12d,%r12d
  mov    %rbp,%r13
  movq   $0x0,0x98(%rsp)
L147:
  movsd  0x0(%r13,%r12,8),%xmm1
  movsd  (%r15,%r12,8),%xmm2
  mov    %r12d,%ecx
  mov    %r12d,%eax
  movsd  D16(%rip),%xmm7
  lea    0x0(%r13,%r12,8),%rdx
  movapd %xmm2,%xmm0
  movapd %xmm1,%xmm3
  andpd  D73(%rip),%xmm0
  andpd  D73(%rip),%xmm3
  addsd  %xmm3,%xmm0
  maxsd  0x98(%rsp),%xmm0
  movsd  %xmm0,0x98(%rsp)
  mulsd  %xmm0,%xmm7
  movapd %xmm1,%xmm0
  movsd  %xmm7,0x58(%rsp)
  movapd %xmm7,%xmm3
  jmp    L61
  nopl   0x0(%rax)
L130:
  movsd  (%rdx),%xmm0
L61:
  andpd  D73(%rip),%xmm0
  comisd %xmm0,%xmm3
  jae    L84
  add    $0x1,%eax
  add    $0x8,%rdx
  cmp    %eax,%r14d
  jg     L130
L84:
  mov    %r12d,0x64(%rsp)
  cmp    %r12d,%eax
  jg     L211
L210:
  addsd  0x38(%rsp),%xmm2
  movsd  %xmm2,(%r15,%r12,8)
  movq   $0x0,0x0(%r13,%r12,8)
  add    $0x1,%r12
  cmp    %r12d,%r14d
  jg     L147
  mov    %rbx,%r13
  cmp    $0x1,%r14d
  jle    L38
  lea    -0x1(%r14),%eax
  lea    -0x2(%r14),%ebp
  mov    $0x1,%r11d
  xor    %edi,%edi
  lea    0x8(%rbx,%rax,8),%rbx
L59:
  movsd  (%r15,%rdi,8),%xmm2
  mov    %edi,%edx
  lea    0x0(,%rdi,8),%rsi
  mov    %r11,%rax
  movapd %xmm2,%xmm0
  cmp    %r11d,%r14d
  jle    L10
  nopl   0x0(%rax)
L64:
  movsd  (%r15,%rax,8),%xmm1
  comisd %xmm1,%xmm0
  minsd  %xmm0,%xmm1
  cmova  %eax,%edx
  add    $0x1,%rax
  movapd %xmm1,%xmm0
  cmp    %eax,%r14d
  jg     L64
  cmp    %edi,%edx
  jne    L217
L10:
  lea    0x1(%rdi),%rax
  add    $0x1,%r11
  cmp    %rdi,%rbp
  je     L38
L53:
  mov    %rax,%rdi
  jmp    L59
L217:
  movslq %edx,%rdx
  movsd  %xmm2,(%r15,%rdx,8)
  lea    0x0(,%rdx,8),%r8
  mov    %r13,%rdx
  movsd  %xmm0,(%r15,%rdi,8)
  xchg   %ax,%ax
L16:
  mov    (%rdx),%rax
  add    $0x8,%rdx
  lea    (%rax,%rsi,1),%rcx
  add    %r8,%rax
  movsd  (%rcx),%xmm0
  movsd  (%rax),%xmm1
  movsd  %xmm1,(%rcx)
  movsd  %xmm0,(%rax)
  cmp    %rbx,%rdx
  jne    L16
  lea    0x1(%rdi),%rax
  add    $0x1,%r11
  cmp    %rdi,%rbp
  jne    L53
L38:
  add    $0xa8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L211:
  movslq %eax,%rsi
  lea    -0x1(%rax),%r9d
  sub    %ecx,%eax
  lea    (%r15,%rsi,8),%rdi
  sub    $0x1,%eax
  mov    %r9d,0x7c(%rsp)
  mov    %rdi,0x68(%rsp)
  lea    0x2(%rcx),%edi
  sub    %rax,%rsi
  movslq %edi,%rdx
  mov    %edi,0x78(%rsp)
  lea    -0x10(,%rsi,8),%rax
  lea    (%r15,%rdx,8),%rdi
  lea    -0x3(%r14),%edx
  mov    %rax,0x88(%rsp)
  sub    %ecx,%edx
  mov    %rdi,0x70(%rsp)
  add    %r12,%rdx
  lea    0x18(%r15,%rdx,8),%rdi
  movslq %r9d,%rdx
  mov    %rdi,0x80(%rsp)
  lea    0x0(,%rdx,8),%rdi
  mov    %rdi,0x90(%rsp)
  nopl   0x0(%rax,%rax,1)
L82:
  movsd  0x8(%r15,%r12,8),%xmm3
  addsd  %xmm1,%xmm1
  mov    D33(%rip),%rax
  movsd  %xmm2,0x10(%rsp)
  subsd  %xmm2,%xmm3
  divsd  %xmm1,%xmm3
  movq   %rax,%xmm1
  movapd %xmm3,%xmm0
  movsd  %xmm3,0x8(%rsp)
  call   <myhypot@plt>
  movsd  0x8(%rsp),%xmm3
  pxor   %xmm7,%xmm7
  movsd  0x10(%rsp),%xmm2
  comisd %xmm3,%xmm7
  jbe    L52
  xorpd  D94(%rip),%xmm0
L52:
  addsd  %xmm3,%xmm0
  movsd  0x0(%r13,%r12,8),%xmm3
  mov    0x70(%rsp),%rax
  mov    0x80(%rsp),%rdx
  divsd  %xmm0,%xmm3
  movsd  %xmm3,(%r15,%r12,8)
  mulsd  0x0(%r13,%r12,8),%xmm0
  subsd  %xmm3,%xmm2
  movapd %xmm2,%xmm1
  movsd  %xmm0,0x40(%rsp)
  movsd  %xmm0,0x8(%r15,%r12,8)
  cmp    %r14d,0x78(%rsp)
  jge    L123
  nopw   0x0(%rax,%rax,1)
L180:
  movsd  (%rax),%xmm0
  add    $0x8,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rax)
  cmp    %rax,%rdx
  jne    L180
L123:
  addsd  0x38(%rsp),%xmm1
  mov    0x68(%rsp),%rax
  movsd  0x8(%r13,%r12,8),%xmm7
  mov    0x7c(%rsp),%edi
  movsd  (%rax),%xmm5
  movsd  %xmm1,0x38(%rsp)
  movsd  %xmm7,0x48(%rsp)
  cmp    %edi,0x64(%rsp)
  jg     L185
  mov    D33(%rip),%rax
  mov    0x90(%rsp),%rbp
  mov    %r12,0x50(%rsp)
  mov    %r15,%r12
  mov    %r13,%r15
  mov    %rbx,%r13
  pxor   %xmm2,%xmm2
  mov    %rax,0x30(%rsp)
  mov    %rbp,%rbx
  mov    0x88(%rsp),%rbp
  movq   %rax,%xmm6
  nopl   0x0(%rax)
L37:
  movsd  (%r15,%rbx,1),%xmm1
  movapd %xmm5,%xmm7
  movapd %xmm5,%xmm0
  movsd  %xmm2,0x28(%rsp)
  mulsd  %xmm6,%xmm7
  movsd  %xmm6,0x18(%rsp)
  movapd %xmm1,%xmm4
  movsd  %xmm5,0x10(%rsp)
  mulsd  %xmm6,%xmm4
  movsd  %xmm7,0x8(%rsp)
  movsd  %xmm4,0x20(%rsp)
  call   <myhypot@plt>
  movsd  0x28(%rsp),%xmm2
  lea    0x8(%rbx),%rsi
  xor    %edx,%edx
  movapd %xmm0,%xmm1
  movsd  0x10(%rsp),%xmm5
  movsd  0x20(%rsp),%xmm4
  mulsd  %xmm2,%xmm1
  movapd %xmm2,%xmm7
  movsd  (%r15,%rbx,1),%xmm2
  movsd  0x18(%rsp),%xmm6
  divsd  %xmm0,%xmm5
  movsd  %xmm1,0x8(%r15,%rbx,1)
  movapd %xmm4,%xmm1
  divsd  %xmm0,%xmm2
  movsd  (%r12,%rbx,1),%xmm0
  movapd %xmm5,%xmm3
  movapd %xmm0,%xmm5
  mulsd  %xmm3,%xmm4
  mulsd  %xmm3,%xmm5
  mulsd  %xmm2,%xmm0
  mulsd  %xmm2,%xmm1
  addsd  %xmm0,%xmm4
  subsd  %xmm1,%xmm5
  mulsd  %xmm2,%xmm4
  addsd  0x8(%rsp),%xmm4
  movsd  %xmm4,0x8(%r12,%rbx,1)
  nopw   0x0(%rax,%rax,1)
L228:
  mov    0x0(%r13,%rdx,8),%rax
  movapd %xmm3,%xmm4
  add    $0x1,%rdx
  lea    (%rax,%rsi,1),%rcx
  add    %rbx,%rax
  movsd  (%rcx),%xmm0
  movsd  (%rax),%xmm1
  mulsd  %xmm0,%xmm4
  mulsd  %xmm2,%xmm1
  mulsd  %xmm2,%xmm0
  addsd  %xmm4,%xmm1
  movsd  %xmm1,(%rcx)
  movsd  (%rax),%xmm1
  mulsd  %xmm3,%xmm1
  subsd  %xmm0,%xmm1
  movsd  %xmm1,(%rax)
  cmp    %edx,%r14d
  jg     L228
  sub    $0x8,%rbx
  cmp    %rbx,%rbp
  je     L94
  movsd  %xmm6,0x30(%rsp)
  movapd %xmm3,%xmm6
  jmp    L37
  nopl   (%rax)
L94:
  movapd %xmm2,%xmm0
  xorpd  D94(%rip),%xmm0
  mov    %r13,%rbx
  mov    %r15,%r13
  mov    %r12,%r15
  mov    0x50(%rsp),%r12
  mulsd  %xmm7,%xmm0
  mulsd  0x30(%rsp),%xmm0
L99:
  movsd  0x48(%rsp),%xmm1
  mulsd  %xmm0,%xmm1
  mulsd  0x0(%r13,%r12,8),%xmm1
  divsd  0x40(%rsp),%xmm1
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0x0(%r13,%r12,8)
  movapd %xmm1,%xmm2
  mulsd  %xmm3,%xmm2
  movsd  %xmm2,(%r15,%r12,8)
  movsd  0x0(%r13,%r12,8),%xmm1
  movapd %xmm1,%xmm0
  andpd  D73(%rip),%xmm0
  comisd 0x58(%rsp),%xmm0
  ja     L82
  jmp    L210
L185:
  mov    D94(%rip),%rax
  pxor   %xmm2,%xmm2
  movq   %rax,%xmm0
  mov    D33(%rip),%rax
  movq   %rax,%xmm3
  jmp    L99
L159:
  movq   $0x0,(%rbx)
  jne    L38
  jmp    L78
```

tool_calls: name: retrieve_data
arguments: label: D16
type: byte[8]
name: retrieve_data
arguments: label: D73
type: qword[2]
name: retrieve_data
arguments: label: D33
type: f64
name: retrieve_data
arguments: label: D94
type: f32[4]

tool: {"D16": [0x5e,0xad,0xe8,0xb7,0x5a,0xfe,0xaf,0x3c]}

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

tool: D33: 1

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

assistant: ```c
void QLalgo2(int n, double *d, double *e, double **V) {
  int i, k, l, m;
  double f = 0.0;
  double tst1 = 0.0;
  double eps = 2.22e-16;
  for (i = 1; i < n; i++) {
    e[i - 1] = e[i];
  }
  e[n - 1] = 0.0;
  for (l = 0; l < n; l++) {
    if (tst1 < fabs(d[l]) + fabs(e[l]))
      tst1 = fabs(d[l]) + fabs(e[l]);
    m = l;
    while (m < n) {
      if (fabs(e[m]) <= eps * tst1) {
        break;
      }
      m++;
    }
    if (m > l) {
      int iter = 0;
      do {
        double dl1, h;
        double g = d[l];
        double p = (d[l + 1] - g) / (2.0 * e[l]);
        double r = myhypot(p, 1.);
        iter = iter + 1;
        if (p < 0) {
          r = -r;
        }
        d[l] = e[l] / (p + r);
        d[l + 1] = e[l] * (p + r);
        dl1 = d[l + 1];
        h = g - d[l];
        for (i = l + 2; i < n; i++) {
          d[i] -= h;
        }
        f = f + h;
        p = d[m];
        {
          double c = 1.0;
          double c2 = c;
          double c3 = c;
          double el1 = e[l + 1];
          double s = 0.0;
          double s2 = 0.0;
          for (i = m - 1; i >= l; i--) {
            c3 = c2;
            c2 = c;
            s2 = s;
            g = c * e[i];
            h = c * p;
            r = myhypot(p, e[i]);
            e[i + 1] = s * r;
            s = e[i] / r;
            c = p / r;
            p = c * d[i] - s * g;
            d[i + 1] = h + s * (c * g + s * d[i]);
            for (k = 0; k < n; k++) {
              h = V[k][i + 1];
              V[k][i + 1] = s * V[k][i] + c * h;
              V[k][i] = c * V[k][i] - s * h;
            }
          }
          p = -s * s2 * c3 * el1 * e[l] / dl1;
          e[l] = s * p;
          d[l] = c * p;
        }
      } while (fabs(e[l]) > eps * tst1);
    }
    d[l] = d[l] + f;
    e[l] = 0.0;
  }
  {
    int j;
    double p;
    for (i = 0; i < n - 1; i++) {
      k = i;
      p = d[i];
      for (j = i + 1; j < n; j++) {
        if (d[j] < p) {
          k = j;
          p = d[j];
        }
      }
      if (k != i) {
        d[k] = d[i];
        d[i] = p;
        for (j = 0; j < n; j++) {
          p = V[j][i];
          V[j][i] = V[j][k];
          V[j][k] = p;
        }
      }
    }
  }
}
```