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

```asm
<EigenTridagQLImplicit>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x68,%rsp
  mov    %rcx,0x10(%rsp)
  movslq %edx,%rcx
  lea    0x0(,%rcx,8),%rbp
  mov    %rdi,0x8(%rsp)
  lea    -0x8(%rsi,%rbp,1),%rbx
  cmp    $0x1,%edx
  jle    L153
  mov    %rsi,%r14
  mov    %edx,%r15d
  lea    -0x1(%rdx),%edx
  xor    %r13d,%r13d
  mov    %rdx,%r12
  lea    0x8(%rsi),%rsi
  shl    $0x3,%rdx
  mov    %r14,%rdi
  call   <memmove@plt>
  movl   $0x0,0x5c(%rsp)
  mov    %r12d,%esi
  movq   D20(%rip),%xmm6
  mov    %r13,%r12
  mov    0x8(%rsp),%r13
  movq   $0x0,(%rbx)
  mov    %r15d,%ebx
  mov    %rbp,%r15
  mov    %r14,%rbp
L128:
  movl   $0x1f,0x58(%rsp)
  mov    %ebx,%eax
  mov    %rbp,%r14
  mov    %r13,%rbx
  mov    %r12,%rbp
  mov    %eax,%r13d
  nopl   (%rax)
L54:
  movsd  (%rbx,%rbp,8),%xmm3
  mov    %rbp,%rax
  movapd %xmm3,%xmm2
  jmp    L42
  cs nopw 0x0(%rax,%rax,1)
L130:
  add    $0x1,%rax
  lea    0x1(%rdx),%r12d
  cmp    %eax,%esi
  jle    L85
L42:
  andpd  %xmm6,%xmm2
  mov    %eax,%edx
  mov    %eax,%r12d
  movapd %xmm2,%xmm0
  movsd  0x8(%rbx,%rax,8),%xmm2
  movapd %xmm2,%xmm1
  andpd  %xmm6,%xmm1
  addsd  %xmm0,%xmm1
  movsd  (%r14,%rax,8),%xmm0
  andpd  %xmm6,%xmm0
  addsd  %xmm1,%xmm0
  ucomisd %xmm1,%xmm0
  jp     L130
  jne    L130
L85:
  mov    %ebp,0x8(%rsp)
  cmp    %ebp,%r12d
  jne    L241
  mov    %r13d,%eax
  mov    %rbp,%r12
  mov    %rbx,%r13
  mov    %r14,%rbp
  mov    %eax,%ebx
L91:
  add    $0x1,%r12
  cmp    %r12d,%esi
  jg     L128
L218:
  mov    0x5c(%rsp),%eax
  add    $0x68,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax,%rax,1)
L241:
  subl   $0x1,0x58(%rsp)
  je     L131
  movsd  0x8(%rbx,%rbp,8),%xmm1
  movsd  (%r14,%rbp,8),%xmm0
  subsd  %xmm3,%xmm1
  addsd  %xmm0,%xmm0
  divsd  %xmm0,%xmm1
  movapd %xmm1,%xmm2
  movapd %xmm1,%xmm0
  mulsd  %xmm1,%xmm0
  andpd  %xmm6,%xmm2
  comisd D15(%rip),%xmm2
  movsd  %xmm2,0x18(%rsp)
  jbe    L27
  movsd  D15(%rip),%xmm3
  mov    %esi,0x28(%rsp)
  movsd  %xmm1,0x20(%rsp)
  divsd  %xmm0,%xmm3
  addsd  D15(%rip),%xmm3
  cvttsd2si %xmm3,%edi
  call   <sqrt@plt>
  movsd  0x18(%rsp),%xmm2
  movsd  0x20(%rsp),%xmm1
  pxor   %xmm7,%xmm7
  movq   D20(%rip),%xmm6
  mov    0x28(%rsp),%esi
  mulsd  %xmm2,%xmm0
L185:
  comisd %xmm7,%xmm1
  movslq %r12d,%rax
  movsd  (%r14,%rbp,8),%xmm2
  andpd  %xmm6,%xmm0
  movsd  (%rbx,%rax,8),%xmm9
  mov    %rax,0x40(%rsp)
  subsd  (%rbx,%rbp,8),%xmm9
  jae    L78
  xorpd  D40(%rip),%xmm0
L78:
  addsd  %xmm1,%xmm0
  movapd %xmm2,%xmm1
  sub    $0x1,%r12d
  divsd  %xmm0,%xmm1
  addsd  %xmm1,%xmm9
  cmp    0x8(%rsp),%r12d
  jl     L11
  mov    D15(%rip),%rax
  mov    %esi,0x54(%rsp)
  movslq %r12d,%r12
  movapd %xmm7,%xmm10
  mov    %rbp,0x48(%rsp)
  mov    %r13d,%ebp
  movq   %rax,%xmm4
  movq   %rax,%xmm3
  mov    %r14,%rax
  mov    %r12,%r14
  mov    %rbx,%r12
  mov    %rax,%rbx
  nopl   0x0(%rax,%rax,1)
L199:
  movsd  (%rbx,%r14,8),%xmm1
  mov    %r14d,%r8d
  mov    %r14d,%r13d
  mulsd  %xmm1,%xmm3
  mulsd  %xmm4,%xmm1
  movapd %xmm9,%xmm4
  andpd  %xmm6,%xmm4
  movapd %xmm3,%xmm2
  andpd  %xmm6,%xmm2
  comisd %xmm4,%xmm2
  ja     L279
  ucomisd %xmm7,%xmm9
  jp     L30
  jne    L30
  mov    %rbx,%rax
  mov    %r12,%rbx
  mov    %r14,%r12
  mov    %ebp,%r13d
  add    $0x1,%r12
  mov    0x48(%rsp),%rbp
  mov    0x54(%rsp),%esi
  mov    %r14d,%r10d
  movq   $0x0,(%rax,%r12,8)
  mov    %rax,%r14
  lea    0x0(,%r12,8),%rdi
L154:
  add    %rbx,%rdi
  mov    0x40(%rsp),%rax
  movsd  (%rdi),%xmm0
  subsd  %xmm10,%xmm0
  movsd  %xmm0,(%rdi)
  movq   $0x0,(%r14,%rax,8)
  mov    $0x1,%eax
L100:
  cmp    0x8(%rsp),%r10d
  jl     L2
  test   %al,%al
  jne    L54
L2:
  movsd  (%rbx,%rbp,8),%xmm0
  mov    0x40(%rsp),%rax
  subsd  %xmm10,%xmm0
  movsd  %xmm0,(%rbx,%rbp,8)
  movsd  %xmm9,(%r14,%rbp,8)
  movq   $0x0,(%r14,%rax,8)
  jmp    L54
  nopw   0x0(%rax,%rax,1)
L279:
  movapd %xmm9,%xmm0
  movsd  %xmm2,0x38(%rsp)
  divsd  %xmm3,%xmm0
  movsd  %xmm1,0x30(%rsp)
  movsd  %xmm10,0x28(%rsp)
  movsd  %xmm9,0x20(%rsp)
  movsd  %xmm3,0x18(%rsp)
  mov    %r14d,0x50(%rsp)
  mulsd  %xmm0,%xmm0
  addsd  D15(%rip),%xmm0
  cvttsd2si %xmm0,%edi
  call   <sqrt@plt>
  movsd  0x38(%rsp),%xmm2
  pxor   %xmm7,%xmm7
  movsd  0x18(%rsp),%xmm3
  movsd  0x20(%rsp),%xmm9
  movsd  0x28(%rsp),%xmm10
  mulsd  %xmm2,%xmm0
  movsd  0x30(%rsp),%xmm1
  mov    0x50(%rsp),%r8d
  movq   D20(%rip),%xmm6
L0:
  ucomisd %xmm7,%xmm0
  lea    0x0(,%r14,8),%rax
  movsd  %xmm0,0x8(%rbx,%r14,8)
  lea    0x8(%rax),%rdi
  jp     L49
  je     L86
L49:
  divsd  %xmm0,%xmm3
  movapd %xmm9,%xmm4
  add    0x10(%rsp),%rax
  xor    %edi,%edi
  movsd  (%r12,%r14,8),%xmm8
  divsd  %xmm0,%xmm4
  movsd  0x8(%r12,%r14,8),%xmm0
  subsd  %xmm10,%xmm0
  subsd  %xmm0,%xmm8
  mulsd  %xmm3,%xmm8
  movapd %xmm3,%xmm10
  movapd %xmm4,%xmm2
  movapd %xmm4,%xmm9
  addsd  %xmm4,%xmm2
  mulsd  %xmm1,%xmm2
  addsd  %xmm2,%xmm8
  mulsd  %xmm8,%xmm10
  mulsd  %xmm8,%xmm9
  addsd  %xmm10,%xmm0
  subsd  %xmm1,%xmm9
  movsd  %xmm0,0x8(%r12,%r14,8)
  nopl   (%rax)
L170:
  movsd  (%rax),%xmm0
  movsd  0x8(%rax),%xmm1
  movapd %xmm4,%xmm5
  add    $0x1,%edi
  mulsd  %xmm1,%xmm5
  movapd %xmm0,%xmm2
  mulsd  %xmm3,%xmm2
  mulsd  %xmm4,%xmm0
  mulsd  %xmm3,%xmm1
  addsd  %xmm5,%xmm2
  subsd  %xmm1,%xmm0
  movsd  %xmm2,0x8(%rax)
  movsd  %xmm0,(%rax)
  add    %r15,%rax
  cmp    %edi,%ebp
  jne    L170
  lea    -0x1(%r8),%r13d
  sub    $0x1,%r14
  cmp    %r14d,0x8(%rsp)
  jle    L199
  ucomisd %xmm7,%xmm8
  mov    $0x0,%edi
  mov    %r13d,%r10d
  mov    %rbx,%r14
  mov    %ebp,%r13d
  mov    0x54(%rsp),%esi
  mov    0x48(%rsp),%rbp
  mov    %r12,%rbx
  setnp  %al
  cmovne %edi,%eax
  jmp    L100
  nopl   0x0(%rax)
L30:
  movapd %xmm3,%xmm0
  mov    %r8d,0x50(%rsp)
  divsd  %xmm9,%xmm0
  movsd  %xmm4,0x38(%rsp)
  movsd  %xmm1,0x30(%rsp)
  movsd  %xmm10,0x28(%rsp)
  movsd  %xmm3,0x20(%rsp)
  movsd  %xmm9,0x18(%rsp)
  mulsd  %xmm0,%xmm0
  addsd  D15(%rip),%xmm0
  cvttsd2si %xmm0,%edi
  call   <sqrt@plt>
  movsd  0x38(%rsp),%xmm4
  mov    0x50(%rsp),%r8d
  pxor   %xmm7,%xmm7
  movq   D20(%rip),%xmm6
  movsd  0x30(%rsp),%xmm1
  mulsd  %xmm4,%xmm0
  movsd  0x28(%rsp),%xmm10
  movsd  0x20(%rsp),%xmm3
  movsd  0x18(%rsp),%xmm9
  jmp    L0
L27:
  addsd  D15(%rip),%xmm0
  mov    %esi,0x20(%rsp)
  movsd  %xmm1,0x18(%rsp)
  cvttsd2si %xmm0,%edi
  call   <sqrt@plt>
  mov    0x20(%rsp),%esi
  movsd  0x18(%rsp),%xmm1
  pxor   %xmm7,%xmm7
  movq   D20(%rip),%xmm6
  jmp    L185
L131:
  mov    %r13d,%eax
  mov    %rbp,%r12
  mov    %rbx,%r13
  mov    %r14,%rbp
  movl   $0xffffffff,0x5c(%rsp)
  mov    %eax,%ebx
  jmp    L91
L11:
  movapd %xmm7,%xmm10
  jmp    L2
L153:
  movq   $0x0,(%rbx)
  movl   $0x0,0x5c(%rsp)
  jmp    L218
L86:
  mov    %r13d,%r10d
  mov    %rbx,%r14
  mov    %ebp,%r13d
  mov    0x54(%rsp),%esi
  mov    0x48(%rsp),%rbp
  mov    %r12,%rbx
  jmp    L154
```

tool_calls: name: retrieve_data
arguments: tag: D20
type: qword
name: retrieve_data
arguments: tag: D15
type: f32[4]
name: retrieve_data
arguments: tag: D40
type: f32[4]

tool: D20: 9223372036854775807

tool: {"D15": [0.0,1.875,0.0,0.0]}

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

assistant: ```c
int EigenTridagQLImplicit(double d[], double e[], int n, double z[]) {
  int m, j, iter, niter = 30, status = 0, i, k;
  double s, r, p, g, f, dd, c, b, aa, bb;
  for (i = 1; i < n; i++)
    e[i - 1] = e[i];
  e[n - 1] = 0;
  for (j = 0; j < n; j++) {
    iter = 0;
    do {
      for (m = j; m < n - 1; m++) {
        dd = fabs(d[m]) + fabs(d[m + 1]);
        if (fabs(e[m]) + dd == dd)
          break;
      }
      if (m != j) {
        if (iter++ == niter) {
          status = -1;
          break;
        }
        g = (d[j + 1] - d[j]) / (2 * e[j]);
        if ((aa = fabs(g)) > 1)
          r = aa * sqrt(1 + 1 / (g * g));
        else
          r = sqrt(1 + g * g);
        g = d[m] - d[j] + e[j] / (g + ((g) >= 0.0 ? fabs(r) : -fabs(r)));
        s = c = 1;
        p = 0;
        for (i = m - 1; i >= j; i--) {
          f = s * e[i];
          b = c * e[i];
          aa = fabs(f);
          bb = fabs(g);
          if (aa > bb) {
            bb /= aa;
            r = aa * sqrt(1 + bb * bb);
          } else if (bb == 0)
            r = 0;
          else {
            aa /= bb;
            r = bb * sqrt(1 + aa * aa);
          }
          e[i + 1] = r;
          if (r == 0) {
            d[i + 1] -= p;
            e[m] = 0;
            break;
          }
          s = f / r;
          c = g / r;
          g = d[i + 1] - p;
          r = (d[i] - g) * s + 2 * c * b;
          d[i + 1] = g + (p = s * r);
          g = c * r - b;
          for (k = 0; k < n; k++) {
            f = z[k * n + i + 1];
            z[k * n + i + 1] = s * z[k * n + i] + c * f;
            z[k * n + i] = c * z[k * n + i] - s * f;
          }
        }
        if (r == 0 && i >= j)
          continue;
        d[j] -= p;
        e[j] = g;
        e[m] = 0;
      }
    } while (m != j);
  }
  return (status);
}
```