user: What is the c source code of the assembly code below:

```asm
<eigen>:
  endbr64
  push   %r15
  mov    %rsi,%r15
  push   %r14
  mov    %edi,%r14d
  push   %r13
  mov    %r8,%r13
  push   %r12
  mov    %rcx,%r12
  push   %rbp
  push   %rbx
  mov    %edx,%ebx
  sub    $0x48,%rsp
  movsd  D14(%rip),%xmm1
  movsd  D20(%rip),%xmm0
  mov    0x80(%rsp),%rax
  mov    %r9,0x8(%rsp)
  mov    0x88(%rsp),%rbp
  mov    %rax,0x10(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x38(%rsp)
  xor    %eax,%eax
  call   <pow@plt>
  mov    %eax,%edi
  call   <sqrt@plt>
  lea    0x34(%rsp),%rcx
  mov    %rbp,%r8
  mov    %ebx,%esi
  lea    0x30(%rsp),%rdx
  mov    %r15,%rdi
  movsd  %xmm0,0x28(%rsp)
  mov    %rbp,0x18(%rsp)
  call   <balance@plt>
  mov    0x18(%rsp),%rdi
  movslq %ebx,%r10
  mov    %ebx,%edx
  lea    0x0(,%r10,8),%rbp
  mov    %r10,0x20(%rsp)
  mov    %r15,%rsi
  lea    (%rdi,%rbp,1),%rax
  mov    %r14d,%edi
  push   %rax
  push   0x18(%rsp)
  mov    0x18(%rsp),%r9
  mov    0x44(%rsp),%r8d
  mov    0x40(%rsp),%ecx
  call   <elemhess@plt>
  pop    %rax
  push   0x18(%rsp)
  mov    %r12,%r9
  push   0x18(%rsp)
  mov    %ebx,%edx
  mov    %r15,%rsi
  mov    %r14d,%edi
  push   %r13
  mov    0x54(%rsp),%r8d
  mov    0x50(%rsp),%ecx
  call   <realeig@plt>
  add    $0x20,%rsp
  mov    %eax,%r9d
  cmp    $0xffffffff,%eax
  je     L61
  test   %r14d,%r14d
  mov    0x20(%rsp),%r10
  movsd  0x28(%rsp),%xmm0
  jne    L78
  test   %ebx,%ebx
  jle    L6
L109:
  mov    %r12,%r11
  mov    $0x1,%r14d
  xor    %r8d,%r8d
  xor    %r9d,%r9d
  movq   D25(%rip),%xmm4
  mov    $0x1,%r15d
  xchg   %ax,%ax
L38:
  movsd  (%r11),%xmm3
  movslq %r8d,%rdi
  lea    0x0(,%r8,8),%rax
  cmp    %r14d,%ebx
  jle    L72
  mov    %r14,%rax
  movapd %xmm3,%xmm1
L14:
  movsd  (%r12,%rax,8),%xmm2
  comisd %xmm1,%xmm2
  maxsd  %xmm1,%xmm2
  cmova  %eax,%edi
  add    $0x1,%rax
  movapd %xmm2,%xmm1
  cmp    %eax,%ebx
  jg     L14
  movslq %edi,%rdi
  lea    0x0(,%rdi,8),%rax
  lea    (%r12,%rax,1),%rdx
L52:
  movsd  %xmm3,(%rdx)
  add    %r13,%rax
  mov    0x8(%rsp),%rsi
  mov    %r8,%rcx
  movsd  %xmm1,(%r11)
  movsd  0x0(%r13,%r8,8),%xmm2
  movsd  (%rax),%xmm1
  sub    %rdi,%rcx
  movsd  %xmm2,(%rax)
  lea    0x0(,%rdi,8),%rax
  lea    (%rsi,%rax,1),%rdx
  movsd  %xmm1,0x0(%r13,%r8,8)
  add    0x10(%rsp),%rax
  xor    %esi,%esi
  nopl   0x0(%rax)
L59:
  movsd  (%rdx),%xmm1
  movsd  (%rdx,%rcx,8),%xmm2
  add    $0x1,%esi
  movsd  %xmm2,(%rdx)
  movsd  %xmm1,(%rdx,%rcx,8)
  movsd  (%rax),%xmm1
  movsd  (%rax,%rcx,8),%xmm2
  add    %rbp,%rdx
  movsd  %xmm2,(%rax)
  movsd  %xmm1,(%rax,%rcx,8)
  add    %rbp,%rax
  cmp    %esi,%ebx
  jne    L59
  movsd  0x0(%r13,%r8,8),%xmm1
  andpd  %xmm4,%xmm1
  ucomisd %xmm0,%xmm1
  cmova  %r15d,%r9d
  add    $0x1,%r8
  add    $0x1,%r14
  add    $0x8,%r11
  cmp    %r10,%r8
  jne    L38
L61:
  mov    0x38(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L12
  add    $0x48,%rsp
  mov    %r9d,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax,%rax,1)
L72:
  mov    %r11,%rdx
  movapd %xmm3,%xmm1
  jmp    L52
L78:
  mov    0x18(%rsp),%r9
  mov    0x30(%rsp),%ecx
  mov    %ebx,%edi
  mov    %r10,0x28(%rsp)
  mov    0x34(%rsp),%r8d
  mov    0x10(%rsp),%rdx
  movsd  %xmm0,0x20(%rsp)
  mov    0x8(%rsp),%rsi
  call   <unbalance@plt>
  mov    0x28(%rsp),%r10
  movsd  0x20(%rsp),%xmm0
  test   %ebx,%ebx
  jg     L109
L6:
  xor    %r9d,%r9d
  jmp    L61
L12:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: read_data
arguments: label: D14
datatype: byte[8]
name: read_data
arguments: label: D20
datatype: f64
name: read_data
arguments: label: D25
datatype: qword

tool: {"D14": [0x00,0x00,0x00,0x00,0x00,0x80,0x43,0xc0]}

tool: D20: 2

tool: D25: 9223372036854775807

assistant: ```c
int eigen(int job, double *A, int n, double *rr, double *ri, double *vr, double *vi, double *work) {
  int low, hi, i, j, k, it, istate = 0;
  double tiny = sqrt(pow((double)2, (double)(1 - 40))), t;
  balance(A, n, &low, &hi, work);
  elemhess(job, A, n, low, hi, vr, vi, (int *)(work + n));
  if (-1 == realeig(job, A, n, low, hi, rr, ri, vr, vi))
    return (-1);
  if (job)
    unbalance(n, vr, vi, low, hi, work);
  for (i = 0; i < n; i++) {
    for (j = i + 1, it = i, t = rr[i]; j < n; j++)
      if (t < rr[j]) {
        t = rr[j];
        it = j;
      }
    rr[it] = rr[i];
    rr[i] = t;
    t = ri[it];
    ri[it] = ri[i];
    ri[i] = t;
    for (k = 0; k < n; k++) {
      t = vr[k * n + it];
      vr[k * n + it] = vr[k * n + i];
      vr[k * n + i] = t;
      t = vi[k * n + it];
      vi[k * n + it] = vi[k * n + i];
      vi[k * n + i] = t;
    }
    if (fabs(ri[i]) > tiny)
      istate = 1;
  }
  return (istate);
}
```