user: Demonstrate the C code that generates the assembly output below:

```asm
<EigenTridagQLImplicit>:
  endbr64
  push   %r15
  mov    %edx,%r15d
  push   %r14
  push   %r13
  push   %r12
  mov    %rsi,%r12
  push   %rbp
  mov    %rdi,%rbp
  push   %rbx
  sub    $0x68,%rsp
  mov    %rcx,0x8(%rsp)
  movslq %edx,%rcx
  lea    0x0(,%rcx,8),%r14
  lea    -0x8(%rsi,%r14,1),%rbx
  cmp    $0x1,%edx
  jle    L286
  lea    -0x1(%rdx),%edx
  lea    0x8(%rsi),%rsi
  mov    %r12,%rdi
  mov    %edx,0x10(%rsp)
  shl    $0x3,%rdx
  call   <memmove@plt>
  movq   $0x0,(%rbx)
  mov    0x10(%rsp),%r9d
L163:
  movl   $0x0,0x5c(%rsp)
  xor    %r13d,%r13d
  mov    %r12,%rax
  movq   D26(%rip),%xmm6
  mov    %r13,%r12
  mov    %r14,%r13
  mov    %rbp,%r14
  mov    %rax,%rbp
L260:
  movl   $0x1f,0x58(%rsp)
  cs nopw 0x0(%rax,%rax,1)
L147:
  mov    %r12d,%esi
  cmp    %r9d,%r12d
  jge    L126
  mov    %r12,%rax
  jmp    L77
  nopl   (%rax)
L41:
  add    $0x1,%rax
  lea    0x1(%rdx),%ebx
  cmp    %eax,%r9d
  jle    L256
L77:
  movsd  (%r14,%rax,8),%xmm0
  movsd  0x8(%r14,%rax,8),%xmm1
  mov    %eax,%edx
  mov    %eax,%ebx
  andpd  %xmm6,%xmm1
  andpd  %xmm6,%xmm0
  addsd  %xmm1,%xmm0
  movsd  0x0(%rbp,%rax,8),%xmm1
  andpd  %xmm6,%xmm1
  addsd  %xmm0,%xmm1
  ucomisd %xmm0,%xmm1
  jp     L41
  jne    L41
L256:
  cmp    %esi,%ebx
  jne    L235
L126:
  add    $0x1,%r12
  cmp    %r12d,%r15d
  jg     L260
L1:
  mov    0x5c(%rsp),%eax
  add    $0x68,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nop
L235:
  subl   $0x1,0x58(%rsp)
  je     L27
  movsd  0x0(%rbp,%r12,8),%xmm0
  movsd  0x8(%r14,%r12,8),%xmm1
  subsd  (%r14,%r12,8),%xmm1
  addsd  %xmm0,%xmm0
  divsd  %xmm0,%xmm1
  movapd %xmm1,%xmm2
  movapd %xmm1,%xmm0
  mulsd  %xmm1,%xmm0
  andpd  %xmm6,%xmm2
  comisd D14(%rip),%xmm2
  movsd  %xmm2,0x10(%rsp)
  jbe    L155
  movsd  D14(%rip),%xmm3
  mov    %r9d,0x28(%rsp)
  mov    %esi,0x20(%rsp)
  divsd  %xmm0,%xmm3
  movsd  %xmm1,0x18(%rsp)
  addsd  D14(%rip),%xmm3
  cvttsd2si %xmm3,%edi
  call   <sqrt@plt>
  movsd  0x10(%rsp),%xmm2
  movsd  0x18(%rsp),%xmm1
  pxor   %xmm7,%xmm7
  mov    0x20(%rsp),%esi
  movq   D26(%rip),%xmm6
  mulsd  %xmm2,%xmm0
  mov    0x28(%rsp),%r9d
L287:
  comisd %xmm7,%xmm1
  movslq %ebx,%rax
  movsd  0x0(%rbp,%r12,8),%xmm2
  andpd  %xmm6,%xmm0
  movsd  (%r14,%rax,8),%xmm9
  mov    %rax,0x40(%rsp)
  subsd  (%r14,%r12,8),%xmm9
  jae    L67
  xorpd  D98(%rip),%xmm0
L67:
  addsd  %xmm1,%xmm0
  sub    $0x1,%ebx
  divsd  %xmm0,%xmm2
  addsd  %xmm2,%xmm9
  cmp    %esi,%ebx
  jl     L232
  mov    D14(%rip),%rax
  mov    %r9d,0x54(%rsp)
  movslq %ebx,%rbx
  movapd %xmm7,%xmm10
  mov    %r12,0x48(%rsp)
  mov    %rbp,%r12
  mov    %r14,%rbp
  mov    %r13,%r14
  movq   %rax,%xmm4
  mov    %rbx,%r13
  movq   %rax,%xmm3
  mov    %esi,%ebx
  nopw   0x0(%rax,%rax,1)
L25:
  movsd  (%r12,%r13,8),%xmm1
  mov    %r13d,%r8d
  mov    %r13d,%r10d
  mulsd  %xmm1,%xmm3
  mulsd  %xmm4,%xmm1
  movapd %xmm9,%xmm4
  andpd  %xmm6,%xmm4
  movapd %xmm3,%xmm2
  andpd  %xmm6,%xmm2
  comisd %xmm4,%xmm2
  ja     L277
  ucomisd %xmm7,%xmm9
  jp     L226
  jne    L226
  mov    %ebx,%esi
  mov    %r13,%rbx
  mov    %r14,%r13
  mov    %rbp,%r14
  add    $0x1,%rbx
  mov    %r12,%rbp
  mov    0x54(%rsp),%r9d
  mov    0x48(%rsp),%r12
  movq   $0x0,0x0(%rbp,%rbx,8)
  lea    0x0(,%rbx,8),%rdi
L289:
  add    %r14,%rdi
  mov    0x40(%rsp),%rax
  movsd  (%rdi),%xmm0
  subsd  %xmm10,%xmm0
  movsd  %xmm0,(%rdi)
  movq   $0x0,0x0(%rbp,%rax,8)
  mov    $0x1,%eax
L210:
  cmp    %esi,%r10d
  jl     L106
  test   %al,%al
  jne    L147
L106:
  movsd  (%r14,%r12,8),%xmm0
  mov    0x40(%rsp),%rax
  subsd  %xmm10,%xmm0
  movsd  %xmm0,(%r14,%r12,8)
  movsd  %xmm9,0x0(%rbp,%r12,8)
  movq   $0x0,0x0(%rbp,%rax,8)
  jmp    L147
L277:
  movapd %xmm9,%xmm0
  movsd  %xmm2,0x30(%rsp)
  divsd  %xmm3,%xmm0
  movsd  %xmm1,0x28(%rsp)
  movsd  %xmm10,0x20(%rsp)
  movsd  %xmm9,0x18(%rsp)
  movsd  %xmm3,0x10(%rsp)
  mov    %r13d,0x50(%rsp)
  mov    %r13d,0x3c(%rsp)
  mulsd  %xmm0,%xmm0
  addsd  D14(%rip),%xmm0
  cvttsd2si %xmm0,%edi
  call   <sqrt@plt>
  movsd  0x30(%rsp),%xmm2
  pxor   %xmm7,%xmm7
  movsd  0x10(%rsp),%xmm3
  movsd  0x18(%rsp),%xmm9
  movsd  0x20(%rsp),%xmm10
  mulsd  %xmm2,%xmm0
  movsd  0x28(%rsp),%xmm1
  mov    0x3c(%rsp),%r8d
  movq   D26(%rip),%xmm6
  mov    0x50(%rsp),%r10d
L117:
  ucomisd %xmm7,%xmm0
  lea    0x0(,%r13,8),%rax
  movsd  %xmm0,0x8(%r12,%r13,8)
  lea    0x8(%rax),%rdi
  jp     L193
  je     L276
L193:
  divsd  %xmm0,%xmm3
  movapd %xmm9,%xmm4
  add    0x8(%rsp),%rax
  xor    %edi,%edi
  movsd  0x0(%rbp,%r13,8),%xmm8
  divsd  %xmm0,%xmm4
  movsd  0x8(%rbp,%r13,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(%rbp,%r13,8)
  nopl   0x0(%rax,%rax,1)
L139:
  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    %r14,%rax
  cmp    %edi,%r15d
  jne    L139
  sub    $0x1,%r13
  lea    -0x1(%r8),%r10d
  cmp    %r13d,%ebx
  jle    L25
  ucomisd %xmm7,%xmm8
  mov    $0x0,%ecx
  mov    %r14,%r13
  mov    %ebx,%esi
  mov    %rbp,%r14
  mov    0x54(%rsp),%r9d
  mov    %r12,%rbp
  mov    0x48(%rsp),%r12
  setnp  %al
  cmovne %ecx,%eax
  jmp    L210
  nopl   0x0(%rax,%rax,1)
L226:
  movapd %xmm3,%xmm0
  mov    %r10d,0x50(%rsp)
  divsd  %xmm9,%xmm0
  mov    %r8d,0x3c(%rsp)
  movsd  %xmm4,0x30(%rsp)
  movsd  %xmm1,0x28(%rsp)
  movsd  %xmm10,0x20(%rsp)
  movsd  %xmm3,0x18(%rsp)
  movsd  %xmm9,0x10(%rsp)
  mulsd  %xmm0,%xmm0
  addsd  D14(%rip),%xmm0
  cvttsd2si %xmm0,%edi
  call   <sqrt@plt>
  movsd  0x30(%rsp),%xmm4
  mov    0x50(%rsp),%r10d
  pxor   %xmm7,%xmm7
  movq   D26(%rip),%xmm6
  mov    0x3c(%rsp),%r8d
  mulsd  %xmm4,%xmm0
  movsd  0x28(%rsp),%xmm1
  movsd  0x20(%rsp),%xmm10
  movsd  0x18(%rsp),%xmm3
  movsd  0x10(%rsp),%xmm9
  jmp    L117
L155:
  addsd  D14(%rip),%xmm0
  mov    %r9d,0x20(%rsp)
  mov    %esi,0x18(%rsp)
  movsd  %xmm1,0x10(%rsp)
  cvttsd2si %xmm0,%edi
  call   <sqrt@plt>
  mov    0x20(%rsp),%r9d
  mov    0x18(%rsp),%esi
  pxor   %xmm7,%xmm7
  movq   D26(%rip),%xmm6
  movsd  0x10(%rsp),%xmm1
  jmp    L287
L27:
  movl   $0xffffffff,0x5c(%rsp)
  jmp    L126
L232:
  movapd %xmm7,%xmm10
  jmp    L106
L286:
  movq   $0x0,(%rbx)
  je     L231
  movl   $0x0,0x5c(%rsp)
  jmp    L1
L231:
  xor    %r9d,%r9d
  jmp    L163
L276:
  mov    %r14,%r13
  mov    0x54(%rsp),%r9d
  mov    %rbp,%r14
  mov    %ebx,%esi
  mov    %r12,%rbp
  mov    0x48(%rsp),%r12
  jmp    L289
```

tool_calls: name: retrieve_data
arguments: data_tag: D26
data_type: qword
name: retrieve_data
arguments: data_tag: D14
data_type: f32[4]
name: retrieve_data
arguments: data_tag: D98
data_type: f32[4]

tool: D26: 9223372036854775807

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

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