user: Please show me the C implementation that matches the assembly code listed below:

```asm
<get_evec>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x28,%rsp
  mov    %rsi,(%rsp)
  mov    %rdx,%rbx
  mov    %rcx,%rbp
  mov    %r8d,%r13d
  mov    %r9d,0x14(%rsp)
  movsd  %xmm0,0x8(%rsp)
  movapd %xmm1,%xmm3
  andpd  D9(%rip),%xmm3
  maxsd  %xmm3,%xmm0
  movapd %xmm0,%xmm3
  mov    0x68(%rsp),%rax
  movl   $0x0,(%rax)
  mov    %r8d,%edx
  imul   %r8d,%edx
  test   %edx,%edx
  jle    L124
  mov    %edx,%edx
  mov    $0x0,%eax
  mov    %rsi,%rcx
L35:
  movsd  (%rdi,%rax,8),%xmm0
  movsd  %xmm0,(%rcx,%rax,8)
  add    $0x1,%rax
  cmp    %rdx,%rax
  jne    L35
L124:
  test   %r13d,%r13d
  jle    L97
  movsd  0x8(%rsp),%xmm5
  movapd %xmm5,%xmm2
  mulsd  %xmm5,%xmm2
  mulsd  %xmm3,%xmm2
  movslq %r13d,%rax
  lea    0x8(,%rax,8),%rdx
  mov    (%rsp),%rax
  mov    $0x0,%esi
L205:
  movsd  (%rax),%xmm0
  subsd  %xmm1,%xmm0
  subsd  %xmm2,%xmm0
  movsd  %xmm0,(%rax)
  mov    %esi,%r14d
  add    $0x1,%esi
  add    %rdx,%rax
  cmp    %esi,%r13d
  jne    L205
  mov    0x68(%rsp),%r8
  mov    0x60(%rsp),%rcx
  mov    %esi,%edx
  mov    (%rsp),%rdi
  mov    $0x0,%eax
  call   <sgefa@plt>
  mov    0x68(%rsp),%rax
  mov    %r14d,%r14d
  mov    $0x0,%r12d
  cmpl   $0xffffffff,(%rax)
  jne    L108
L61:
  mov    $0x0,%eax
  call   <ndrand48@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  mulsd  D97(%rip),%xmm0
  addsd  D50(%rip),%xmm0
  movsd  %xmm0,(%rbx,%r12,8)
  movsd  %xmm0,0x0(%rbp,%r12,8)
  mov    %r12,%rax
  add    $0x1,%r12
  cmp    %rax,%r14
  jne    L61
L29:
  mov    0x68(%rsp),%rax
  movl   $0x0,(%rax)
  lea    -0x1(%r13),%eax
  mov    %rax,0x18(%rsp)
  mov    %rbx,%r14
  lea    0x8(%rbx,%rax,8),%r12
  mov    $0x0,%r15d
  jmp    L92
L93:
  movsd  (%rbx,%rdx,8),%xmm0
  andpd  D9(%rip),%xmm0
  comisd %xmm1,%xmm0
  maxsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  cmova  %ecx,%eax
  add    $0x1,%rdx
L171:
  mov    %edx,%ecx
  cmp    %rsi,%rdx
  jne    L93
  cltq
  movsd  (%rbx,%rax,8),%xmm1
  mov    %r14,%rax
L64:
  movsd  (%rax),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  add    $0x8,%rax
  cmp    %rax,%r12
  jne    L64
  mov    $0x0,%eax
  pxor   %xmm2,%xmm2
L31:
  movsd  (%rbx,%rax,8),%xmm1
  movapd %xmm1,%xmm0
  subsd  0x0(%rbp,%rax,8),%xmm0
  andpd  D9(%rip),%xmm0
  addsd  %xmm0,%xmm2
  movsd  %xmm1,0x0(%rbp,%rax,8)
  add    $0x1,%rax
  cmp    %rdx,%rax
  jne    L31
L189:
  movsd  0x8(%rsp),%xmm4
  comisd %xmm2,%xmm4
  ja     L212
  add    $0x1,%r15d
  cmp    0x14(%rsp),%r15d
  jg     L202
L92:
  mov    $0x0,%r9d
  mov    %rbx,%r8
  mov    0x60(%rsp),%rcx
  mov    %r13d,%edx
  mov    %r13d,%esi
  mov    (%rsp),%rdi
  mov    $0x0,%eax
  call   <sgesl@plt>
  movsd  (%rbx),%xmm1
  andpd  D9(%rip),%xmm1
  test   %r13d,%r13d
  jle    L153
  mov    0x18(%rsp),%rax
  lea    0x1(%rax),%rsi
  mov    $0x1,%edx
  mov    $0x0,%eax
  jmp    L171
L153:
  pxor   %xmm2,%xmm2
  jmp    L189
L202:
  lea    D86(%rip),%rdi
  mov    $0x0,%eax
  call   <plintf@plt>
  mov    0x68(%rsp),%rax
  movl   $0x1,(%rax)
  jmp    L148
L212:
  mov    0x68(%rsp),%rax
  mov    (%rax),%ecx
  test   %ecx,%ecx
  jne    L148
  movsd  (%rbx),%xmm1
  andpd  D9(%rip),%xmm1
  test   %r13d,%r13d
  jle    L148
  mov    0x18(%rsp),%rdx
  add    $0x1,%rdx
  mov    $0x1,%eax
  movq   D9(%rip),%xmm2
L121:
  mov    %eax,%esi
  cmp    %rax,%rdx
  je     L103
  movsd  (%rbx,%rax,8),%xmm0
  andpd  %xmm2,%xmm0
  comisd %xmm1,%xmm0
  maxsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  cmova  %esi,%ecx
  add    $0x1,%rax
  jmp    L121
L103:
  movslq %ecx,%rcx
  movsd  (%rbx,%rcx,8),%xmm1
L173:
  movsd  (%r14),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%r14)
  add    $0x8,%r14
  cmp    %r14,%r12
  jne    L173
  jmp    L148
L97:
  mov    0x68(%rsp),%r8
  mov    0x60(%rsp),%rcx
  mov    %r13d,%edx
  mov    %r13d,%esi
  mov    (%rsp),%rdi
  mov    $0x0,%eax
  call   <sgefa@plt>
  mov    0x68(%rsp),%rax
  cmpl   $0xffffffff,(%rax)
  je     L29
L108:
  lea    D94(%rip),%rdi
  mov    $0x0,%eax
  call   <plintf@plt>
L148:
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: extract_symbolic_data
arguments: data_tag: D9
data_type: qword[2]
name: extract_symbolic_data
arguments: data_tag: D97
data_type: byte[8]
name: extract_symbolic_data
arguments: data_tag: D50
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D86
data_type: string
name: extract_symbolic_data
arguments: data_tag: D94
data_type: string

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

tool: {"D97": [0x9a,0x99,0x99,0x99,0x99,0x99,0xb9,0x3f]}

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

tool: D86:  max iterates exceeded


tool: D94:  Pivot failed


assistant: ```c
void get_evec(a, anew, b, bp, n, maxit, err, ipivot, eval, ierr) double *a, *anew, *b, *bp, err, eval;
int n, maxit, *ipivot, *ierr;
{
  int j, iter, jmax;
  double temp;
  double zz = fabs(eval);
  if (zz < err)
    zz = err;
  *ierr = 0;
  for (j = 0; j < n * n; j++) {
    anew[j] = a[j];
  }
  for (j = 0; j < n; j++)
    anew[j * (1 + n)] = anew[j * (1 + n)] - eval - err * err * zz;
  sgefa(anew, n, n, ipivot, ierr);
  if (*ierr != -1) {
    plintf(" Pivot failed\n");
    return;
  }
  for (j = 0; j < n; j++) {
    b[j] = 1 + .1 * ndrand48();
    bp[j] = b[j];
  }
  iter = 0;
  *ierr = 0;
  while (1) {
    sgesl(anew, n, n, ipivot, b, 0);
    temp = fabs(b[0]);
    jmax = 0;
    for (j = 0; j < n; j++) {
      if (fabs(b[j]) > temp) {
        temp = fabs(b[j]);
        jmax = j;
      }
    }
    temp = b[jmax];
    for (j = 0; j < n; j++)
      b[j] = b[j] / temp;
    temp = 0.0;
    for (j = 0; j < n; j++) {
      temp = temp + fabs(b[j] - bp[j]);
      bp[j] = b[j];
    }
    if (temp < err)
      break;
    iter++;
    if (iter > maxit) {
      plintf(" max iterates exceeded\n");
      *ierr = 1;
      break;
    }
  }
  if (*ierr == 0) {
    temp = fabs(b[0]);
    jmax = 0;
    for (j = 0; j < n; j++) {
      if (fabs(b[j]) > temp) {
        temp = fabs(b[j]);
        jmax = j;
      }
    }
    temp = b[jmax];
    for (j = 0; j < n; j++)
      b[j] = b[j] / temp;
  }
  return;
}
```