user: 请将这段汇编代码片段转换成等效的C代码:

```asm
<cvJacobiEigens_64d>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x120,%rsp
  mov    %rdi,-0xf8(%rbp)
  mov    %rsi,-0x100(%rbp)
  mov    %rdx,-0x108(%rbp)
  mov    %ecx,-0x10c(%rbp)
  movsd  %xmm0,-0x118(%rbp)
  mov    -0xf8(%rbp),%rax
  mov    %rax,-0xd0(%rbp)
  mov    -0x100(%rbp),%rax
  mov    %rax,-0xc8(%rbp)
  mov    -0xf8(%rbp),%rax
  mov    %rax,-0xc0(%rbp)
  mov    -0x100(%rbp),%rax
  mov    %rax,-0xb8(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0xb0(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0xa8(%rbp)
  movsd  D76(%rip),%xmm0
  comisd -0x118(%rbp),%xmm0
  jbe    L217
  movsd  D76(%rip),%xmm0
  movsd  %xmm0,-0x118(%rbp)
L217:
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x10c(%rbp),%xmm1
  movsd  -0x118(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x90(%rbp)
  movl   $0x0,-0xe8(%rbp)
  jmp    L4
L19:
  movl   $0x0,-0xe4(%rbp)
  jmp    L158
L173:
  mov    -0xe4(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x8(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  -0xa8(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0xa8(%rbp)
  addl   $0x1,-0xe4(%rbp)
L158:
  mov    -0xe4(%rbp),%eax
  cmp    -0xe8(%rbp),%eax
  jl     L173
  movl   $0x0,-0xe4(%rbp)
  jmp    L151
L213:
  mov    -0xe4(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xc8(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0xe4(%rbp)
L151:
  mov    -0xe4(%rbp),%eax
  cmp    -0x10c(%rbp),%eax
  jl     L213
  mov    -0xe8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xc8(%rbp),%rax
  add    %rdx,%rax
  movsd  D73(%rip),%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0xe8(%rbp)
  mov    -0x10c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0xc8(%rbp)
  mov    -0x10c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0xd0(%rbp)
L4:
  mov    -0xe8(%rbp),%eax
  cmp    -0x10c(%rbp),%eax
  jl     L19
  movsd  -0xa8(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rax
  mov    %rax,-0xa8(%rbp)
  movsd  -0xa8(%rbp),%xmm0
  mulsd  -0x118(%rbp),%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x10c(%rbp),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x88(%rbp)
  movsd  -0xa8(%rbp),%xmm0
  movsd  %xmm0,-0xb0(%rbp)
  jmp    L15
L105:
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x10c(%rbp),%xmm1
  movsd  -0xb0(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0xb0(%rbp)
L32:
  movl   $0x0,-0xd4(%rbp)
  mov    -0xf8(%rbp),%rax
  mov    %rax,-0xd0(%rbp)
  mov    -0x100(%rbp),%rax
  mov    %rax,-0xc8(%rbp)
  movl   $0x0,-0xdc(%rbp)
  jmp    L222
L109:
  mov    -0xdc(%rbp),%eax
  add    $0x1,%eax
  imul   -0x10c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xf8(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0xc0(%rbp)
  mov    -0xdc(%rbp),%eax
  add    $0x1,%eax
  imul   -0x10c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x100(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0xb8(%rbp)
  mov    -0xdc(%rbp),%eax
  add    $0x1,%eax
  mov    %eax,-0xd8(%rbp)
  jmp    L52
L188:
  mov    -0xd8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movq   D62(%rip),%xmm1
  andpd  %xmm0,%xmm1
  movsd  -0xb0(%rbp),%xmm0
  comisd %xmm1,%xmm0
  ja     L6
  mov    -0xd8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x80(%rbp)
  movl   $0x1,-0xd4(%rbp)
  mov    -0xdc(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x78(%rbp)
  mov    -0xd8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xc0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x70(%rbp)
  movsd  -0x78(%rbp),%xmm0
  movapd %xmm0,%xmm1
  subsd  -0x70(%rbp),%xmm1
  movsd  D71(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x68(%rbp)
  movsd  -0x80(%rbp),%xmm0
  movq   D83(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,-0x120(%rbp)
  movsd  -0x80(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  %xmm0,%xmm1
  movsd  -0x68(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movapd %xmm0,%xmm1
  movsd  -0x120(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0xa0(%rbp)
  pxor   %xmm0,%xmm0
  comisd -0x68(%rbp),%xmm0
  jbe    L224
  movsd  -0xa0(%rbp),%xmm0
  movq   D83(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,-0xa0(%rbp)
L224:
  movsd  -0xa0(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  %xmm0,%xmm1
  movsd  D73(%rip),%xmm0
  subsd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movsd  D73(%rip),%xmm1
  addsd  %xmm1,%xmm0
  addsd  %xmm0,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movapd %xmm0,%xmm1
  movsd  -0xa0(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x60(%rbp)
  movsd  -0x60(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x58(%rbp)
  movsd  D73(%rip),%xmm0
  subsd  -0x58(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x50(%rbp)
  movsd  -0x50(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movsd  -0x80(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  mulsd  -0x50(%rbp),%xmm0
  movsd  -0x60(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  mov    -0xf8(%rbp),%rax
  mov    %rax,-0x98(%rbp)
  movl   $0x0,-0xe8(%rbp)
  jmp    L31
L74:
  mov    -0xdc(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x98(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  mov    -0xd8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x98(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  mov    -0xe8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xc8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0xe8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xb8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x28(%rbp),%xmm0
  mulsd  -0x50(%rbp),%xmm0
  movsd  -0x20(%rbp),%xmm1
  mulsd  -0x60(%rbp),%xmm1
  mov    -0xdc(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x98(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x20(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x50(%rbp),%xmm1
  movsd  -0x28(%rbp),%xmm0
  mulsd  -0x60(%rbp),%xmm0
  mov    -0xd8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x98(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x18(%rbp),%xmm0
  mulsd  -0x50(%rbp),%xmm0
  movsd  -0x10(%rbp),%xmm1
  mulsd  -0x60(%rbp),%xmm1
  mov    -0xe8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xc8(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x10(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x50(%rbp),%xmm1
  movsd  -0x18(%rbp),%xmm0
  mulsd  -0x60(%rbp),%xmm0
  mov    -0xe8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xb8(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0xe8(%rbp)
  mov    -0x10c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0x98(%rbp)
L31:
  mov    -0xe8(%rbp),%eax
  cmp    -0xdc(%rbp),%eax
  jl     L74
  jmp    L246
L17:
  mov    -0xe8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  mov    -0xd8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x98(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  mov    -0xe8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xc8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0xe8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xb8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x28(%rbp),%xmm0
  mulsd  -0x50(%rbp),%xmm0
  movsd  -0x20(%rbp),%xmm1
  mulsd  -0x60(%rbp),%xmm1
  mov    -0xe8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x20(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x50(%rbp),%xmm1
  movsd  -0x28(%rbp),%xmm0
  mulsd  -0x60(%rbp),%xmm0
  mov    -0xd8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x98(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x18(%rbp),%xmm0
  mulsd  -0x50(%rbp),%xmm0
  movsd  -0x10(%rbp),%xmm1
  mulsd  -0x60(%rbp),%xmm1
  mov    -0xe8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xc8(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x10(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x50(%rbp),%xmm1
  movsd  -0x18(%rbp),%xmm0
  mulsd  -0x60(%rbp),%xmm0
  mov    -0xe8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xb8(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0xe8(%rbp)
  mov    -0x10c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0x98(%rbp)
L246:
  mov    -0xe8(%rbp),%eax
  cmp    -0xd8(%rbp),%eax
  jl     L17
  jmp    L153
L7:
  mov    -0xe8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  mov    -0xe8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xc0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  mov    -0xe8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xc8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0xe8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xb8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x28(%rbp),%xmm0
  mulsd  -0x50(%rbp),%xmm0
  movsd  -0x20(%rbp),%xmm1
  mulsd  -0x60(%rbp),%xmm1
  mov    -0xe8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x20(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x50(%rbp),%xmm1
  movsd  -0x28(%rbp),%xmm0
  mulsd  -0x60(%rbp),%xmm0
  mov    -0xe8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xc0(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x18(%rbp),%xmm0
  mulsd  -0x50(%rbp),%xmm0
  movsd  -0x10(%rbp),%xmm1
  mulsd  -0x60(%rbp),%xmm1
  mov    -0xe8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xc8(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x10(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x50(%rbp),%xmm1
  movsd  -0x18(%rbp),%xmm0
  mulsd  -0x60(%rbp),%xmm0
  mov    -0xe8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xb8(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0xe8(%rbp)
L153:
  mov    -0xe8(%rbp),%eax
  cmp    -0x10c(%rbp),%eax
  jl     L7
  movsd  -0x78(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x48(%rbp),%xmm1
  movsd  -0x70(%rbp),%xmm0
  mulsd  -0x58(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  subsd  -0x40(%rbp),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x78(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x58(%rbp),%xmm1
  movsd  -0x70(%rbp),%xmm0
  mulsd  -0x48(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  -0x40(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  mov    -0xdc(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x38(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xd8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xc0(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x30(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xdc(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xc0(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xd8(%rbp),%edx
  movslq %edx,%rdx
  lea    0x0(,%rdx,8),%rcx
  mov    -0xd0(%rbp),%rdx
  add    %rcx,%rdx
  movsd  (%rax),%xmm0
  movsd  %xmm0,(%rdx)
  jmp    L101
L6:
  nop
L101:
  addl   $0x1,-0xd8(%rbp)
  mov    -0x10c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0xc0(%rbp)
  mov    -0x10c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0xb8(%rbp)
L52:
  mov    -0xd8(%rbp),%eax
  cmp    -0x10c(%rbp),%eax
  jl     L188
  addl   $0x1,-0xdc(%rbp)
  mov    -0x10c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0xd0(%rbp)
  mov    -0x10c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0xc8(%rbp)
L222:
  mov    -0x10c(%rbp),%eax
  sub    $0x1,%eax
  cmp    %eax,-0xdc(%rbp)
  jl     L109
  cmpl   $0x0,-0xd4(%rbp)
  jne    L32
L15:
  movsd  -0xb0(%rbp),%xmm0
  comisd -0x88(%rbp),%xmm0
  ja     L105
  movl   $0x0,-0xe8(%rbp)
  movl   $0x0,-0xe0(%rbp)
  jmp    L79
L218:
  mov    -0xe0(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xf8(%rbp),%rax
  add    %rax,%rdx
  mov    -0xe8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0x108(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0xe8(%rbp)
  mov    -0x10c(%rbp),%eax
  add    $0x1,%eax
  add    %eax,-0xe0(%rbp)
L79:
  mov    -0xe8(%rbp),%eax
  cmp    -0x10c(%rbp),%eax
  jl     L218
  mov    $0x0,%eax
  leave
  ret
```

tool_calls: name: extract_symbolic_data
arguments: tag: D76
data_type: byte[8]
name: extract_symbolic_data
arguments: tag: D73
data_type: f64
name: extract_symbolic_data
arguments: tag: D62
data_type: qword
name: extract_symbolic_data
arguments: tag: D71
data_type: f32[2]
name: extract_symbolic_data
arguments: tag: D83
data_type: f32[2]

tool: {"D76": [0x16,0x56,0xe7,0x9e,0xaf,0x03,0xd2,0x3c]}

tool: D73: 1

tool: D62: 9223372036854775807

tool: {"D71": [0.0,1.75]}

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

assistant: ```c
int cvJacobiEigens_64d(double *A, double *V, double *E, int n, double eps) {
  int i, j, k, p, q, ind;
  double *A1 = A, *V1 = V, *A2 = A, *V2 = V;
  double Amax = 0.0, anorm = 0.0, ax, deps;
  if (eps < 1.0e-15)
    eps = 1.0e-15;
  deps = eps / (double)n;
  for (i = 0; i < n; i++, V1 += n, A1 += n) {
    for (j = 0; j < i; j++) {
      double Am = A1[j];
      anorm += Am * Am;
    }
    for (j = 0; j < n; j++)
      V1[j] = 0.0;
    V1[i] = 1.0;
  }
  anorm = sqrt(anorm + anorm);
  ax = anorm * eps / n;
  Amax = anorm;
  while (Amax > ax) {
    Amax /= n;
    do {
      ind = 0;
      A1 = A;
      V1 = V;
      for (p = 0; p < n - 1; p++, A1 += n, V1 += n) {
        A2 = A + n * (p + 1);
        V2 = V + n * (p + 1);
        for (q = p + 1; q < n; q++, A2 += n, V2 += n) {
          double x, y, c, s, c2, s2, a;
          double *A3, Apq, App, Aqq, App2, Aqq2, Aip, Aiq, Vpi, Vqi;
          if (fabs(A1[q]) < Amax)
            continue;
          Apq = A1[q];
          ind = 1;
          App = A1[p];
          Aqq = A2[q];
          y = 5.0e-1 * (App - Aqq);
          x = -Apq / sqrt(Apq * Apq + y * y);
          if (y < 0.0)
            x = -x;
          s = x / sqrt(2.0 * (1.0 + sqrt(1.0 - x * x)));
          s2 = s * s;
          c = sqrt(1.0 - s2);
          c2 = c * c;
          a = 2.0 * Apq * c * s;
          A3 = A;
          for (i = 0; i < p; i++, A3 += n) {
            Aip = A3[p];
            Aiq = A3[q];
            Vpi = V1[i];
            Vqi = V2[i];
            A3[p] = Aip * c - Aiq * s;
            A3[q] = Aiq * c + Aip * s;
            V1[i] = Vpi * c - Vqi * s;
            V2[i] = Vqi * c + Vpi * s;
          }
          for (; i < q; i++, A3 += n) {
            Aip = A1[i];
            Aiq = A3[q];
            Vpi = V1[i];
            Vqi = V2[i];
            A1[i] = Aip * c - Aiq * s;
            A3[q] = Aiq * c + Aip * s;
            V1[i] = Vpi * c - Vqi * s;
            V2[i] = Vqi * c + Vpi * s;
          }
          for (; i < n; i++) {
            Aip = A1[i];
            Aiq = A2[i];
            Vpi = V1[i];
            Vqi = V2[i];
            A1[i] = Aip * c - Aiq * s;
            A2[i] = Aiq * c + Aip * s;
            V1[i] = Vpi * c - Vqi * s;
            V2[i] = Vqi * c + Vpi * s;
          }
          App2 = App * c2 + Aqq * s2 - a;
          Aqq2 = App * s2 + Aqq * c2 + a;
          A1[p] = App2;
          A2[q] = Aqq2;
          A1[q] = A2[p] = 0.0;
        }
      }
    } while (ind);
  }
  for (i = 0, k = 0; i < n; i++, k += n + 1)
    E[i] = A[k];
  return 0;
}
```