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

```asm
<echange>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %rbx
  sub    $0xf8,%rsp
  mov    %edi,-0xc4(%rbp)
  mov    %esi,-0xc8(%rbp)
  mov    %rdx,-0xd0(%rbp)
  mov    %rcx,-0xd8(%rbp)
  mov    %r8,-0xe0(%rbp)
  mov    %r9,-0xe8(%rbp)
  mov    0x10(%rbp),%rax
  mov    %rax,-0xf0(%rbp)
  mov    0x18(%rbp),%rax
  mov    %rax,-0xf8(%rbp)
  mov    0x20(%rbp),%rax
  mov    %rax,-0x100(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x18(%rbp)
  xor    %eax,%eax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  movsd  D73(%rip),%xmm0
  movsd  %xmm0,-0x48(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  mov    -0xc4(%rbp),%eax
  add    $0x1,%eax
  mov    %eax,-0x6c(%rbp)
  mov    -0x6c(%rbp),%eax
  add    $0x1,%eax
  mov    %eax,-0x68(%rbp)
  movl   $0x0,-0xc0(%rbp)
  movl   $0x0,-0xbc(%rbp)
  movl   $0x1,-0xb8(%rbp)
  mov    -0xc4(%rbp),%eax
  sub    $0x1,%eax
  movslq %eax,%rdx
  imul   $0xffffffff88888889,%rdx,%rdx
  shr    $0x20,%rdx
  add    %eax,%edx
  sar    $0x3,%edx
  sar    $0x1f,%eax
  mov    %eax,%ecx
  mov    %edx,%eax
  sub    %ecx,%eax
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  movsd  D20(%rip),%xmm0
  addsd  %xmm1,%xmm0
  cvttsd2si %xmm0,%eax
  mov    %eax,-0x64(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x58(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movl   $0x0,-0xb4(%rbp)
  movl   $0x0,-0xb0(%rbp)
  movl   $0x0,-0xac(%rbp)
  movl   $0x0,-0x60(%rbp)
  movl   $0x0,-0x5c(%rbp)
  movl   $0x0,-0xa8(%rbp)
  movl   $0x0,-0xa4(%rbp)
  movl   $0x1,-0xa0(%rbp)
  movl   $0x0,-0x9c(%rbp)
  movl   $0x0,-0x98(%rbp)
  jmp    L45
L827:
  addl   $0x1,-0xc0(%rbp)
  cmpl   $0xfa,-0xc0(%rbp)
  jg     L234
  movl   $0x0,-0x94(%rbp)
  jmp    L317
L325:
  mov    -0x94(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  movsd  D53(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  mov    -0x94(%rbp),%edx
  movslq %edx,%rdx
  lea    0x0(,%rdx,8),%rcx
  mov    -0xf0(%rbp),%rdx
  lea    (%rcx,%rdx,1),%rbx
  movq   %rax,%xmm0
  call   <cos@plt>
  movq   %xmm0,%rax
  mov    %rax,(%rbx)
  addl   $0x1,-0x94(%rbp)
L317:
  mov    -0x94(%rbp),%eax
  cmp    -0x6c(%rbp),%eax
  jl     L325
  movl   $0x0,-0x90(%rbp)
  jmp    L929
L429:
  mov    -0xf0(%rbp),%rcx
  mov    -0x64(%rbp),%edx
  mov    -0x6c(%rbp),%esi
  mov    -0x90(%rbp),%eax
  mov    %eax,%edi
  mov    $0x0,%eax
  call   <dd@plt>
  mov    %eax,%edx
  mov    -0x90(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0x100(%rbp),%rax
  add    %rcx,%rax
  pxor   %xmm0,%xmm0
  cvtsi2sd %edx,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x90(%rbp)
L929:
  mov    -0x90(%rbp),%eax
  cmp    -0x6c(%rbp),%eax
  jl     L429
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movl   $0x1,-0x8c(%rbp)
  movl   $0x0,-0x88(%rbp)
  jmp    L349
L150:
  mov    -0x88(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x100(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x88(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x30(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x8c(%rbp),%xmm1
  mov    -0x88(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x100(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  mov    -0x88(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xe8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  divsd  %xmm1,%xmm0
  movsd  -0x28(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  negl   -0x8c(%rbp)
  addl   $0x1,-0x88(%rbp)
L349:
  mov    -0x88(%rbp),%eax
  cmp    -0x6c(%rbp),%eax
  jl     L150
  movsd  -0x30(%rbp),%xmm0
  divsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x20(%rbp),%xmm0
  pxor   %xmm1,%xmm1
  comisd %xmm1,%xmm0
  jbe    L125
  mov    $0xffffffff,%eax
  jmp    L722
L125:
  mov    $0x1,%eax
L722:
  mov    %eax,-0x5c(%rbp)
  mov    -0x5c(%rbp),%eax
  neg    %eax
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  -0x20(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movl   $0x1,-0x84(%rbp)
  movl   $0x0,-0x80(%rbp)
  jmp    L617
L616:
  mov    -0x80(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  pxor   %xmm0,%xmm0
  cvtsi2sdl -0x5c(%rbp),%xmm0
  movapd %xmm0,%xmm2
  mulsd  -0x20(%rbp),%xmm2
  pxor   %xmm0,%xmm0
  cvtsi2sdl -0x84(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  mov    -0x80(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xe8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  divsd  %xmm2,%xmm0
  mov    -0x80(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xf8(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  negl   -0x84(%rbp)
  addl   $0x1,-0x80(%rbp)
L617:
  mov    -0x80(%rbp),%eax
  cmp    -0x6c(%rbp),%eax
  jl     L616
  movsd  -0x48(%rbp),%xmm0
  comisd -0x20(%rbp),%xmm0
  jb     L686
  lea    D50(%rip),%rax
  mov    %rax,%rdi
  call   <puts@plt>
  lea    D22(%rip),%rax
  mov    %rax,%rdi
  call   <puts@plt>
  mov    $0xffffffff,%eax
  jmp    L321
L686:
  movsd  -0x20(%rbp),%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movl   $0x0,-0xb8(%rbp)
  mov    -0xd0(%rbp),%rax
  mov    (%rax),%eax
  mov    %eax,-0xb4(%rbp)
  mov    -0x6c(%rbp),%eax
  cltq
  shl    $0x2,%rax
  lea    -0x4(%rax),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0xb0(%rbp)
  movl   $0xffffffff,-0xac(%rbp)
  mov    -0x5c(%rbp),%eax
  neg    %eax
  mov    %eax,-0xa8(%rbp)
  movl   $0x1,-0xa0(%rbp)
  jmp    L772
L124:
  mov    -0xa0(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0x60(%rbp)
  mov    -0xa0(%rbp),%eax
  cltq
  shl    $0x2,%rax
  lea    -0x4(%rax),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  add    $0x1,%eax
  mov    %eax,-0x98(%rbp)
  negl   -0xa8(%rbp)
  cmpl   $0x2,-0xa0(%rbp)
  jne    L441
  movsd  -0x50(%rbp),%xmm0
  movsd  %xmm0,-0x40(%rbp)
L441:
  movsd  -0x20(%rbp),%xmm0
  movsd  %xmm0,-0x50(%rbp)
  movl   $0x1,-0x9c(%rbp)
  mov    -0x98(%rbp),%eax
  cmp    -0x60(%rbp),%eax
  jge    L118
  mov    -0x50(%rbp),%rdi
  mov    -0xe8(%rbp),%r9
  mov    -0xe0(%rbp),%r10
  mov    -0xf8(%rbp),%rcx
  mov    -0xf0(%rbp),%rdx
  mov    -0x100(%rbp),%rsi
  mov    -0xd8(%rbp),%rax
  lea    -0x58(%rbp),%r8
  push   %r8
  mov    -0x6c(%rbp),%r8d
  push   %r8
  mov    -0x98(%rbp),%r8d
  push   %r8
  mov    -0xa8(%rbp),%r8d
  push   %r8
  movq   %rdi,%xmm0
  mov    %r10,%r8
  mov    %rax,%rdi
  mov    $0x1,%eax
  call   <gee@plt>
  add    $0x20,%rsp
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x38(%rbp),%xmm0
  pxor   %xmm1,%xmm1
  comisd %xmm1,%xmm0
  jbe    L118
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0xa8(%rbp),%xmm1
  movsd  -0x58(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  addl   $0x1,-0x98(%rbp)
  jmp    L90
L55:
  mov    -0x50(%rbp),%rdi
  mov    -0xe8(%rbp),%r9
  mov    -0xe0(%rbp),%r10
  mov    -0xf8(%rbp),%rcx
  mov    -0xf0(%rbp),%rdx
  mov    -0x100(%rbp),%rsi
  mov    -0xd8(%rbp),%rax
  lea    -0x58(%rbp),%r8
  push   %r8
  mov    -0x6c(%rbp),%r8d
  push   %r8
  mov    -0x98(%rbp),%r8d
  push   %r8
  mov    -0xa8(%rbp),%r8d
  push   %r8
  movq   %rdi,%xmm0
  mov    %r10,%r8
  mov    %rax,%rdi
  mov    $0x1,%eax
  call   <gee@plt>
  add    $0x20,%rsp
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  pxor   %xmm0,%xmm0
  comisd -0x38(%rbp),%xmm0
  jae    L224
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0xa8(%rbp),%xmm1
  movsd  -0x58(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  addl   $0x1,-0x98(%rbp)
L90:
  mov    -0x98(%rbp),%eax
  cmp    -0x60(%rbp),%eax
  jl     L55
  jmp    L164
L224:
  nop
L164:
  mov    -0xa0(%rbp),%eax
  cltq
  shl    $0x2,%rax
  lea    -0x4(%rax),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    -0x98(%rbp),%edx
  sub    $0x1,%edx
  mov    %edx,(%rax)
  addl   $0x1,-0xa0(%rbp)
  mov    -0x98(%rbp),%eax
  sub    $0x1,%eax
  mov    %eax,-0xac(%rbp)
  addl   $0x1,-0xb8(%rbp)
  movl   $0x0,-0x9c(%rbp)
L118:
  cmpl   $0x0,-0x9c(%rbp)
  je     L772
  subl   $0x2,-0x98(%rbp)
  jmp    L231
L152:
  mov    -0x50(%rbp),%rdi
  mov    -0xe8(%rbp),%r9
  mov    -0xe0(%rbp),%r10
  mov    -0xf8(%rbp),%rcx
  mov    -0xf0(%rbp),%rdx
  mov    -0x100(%rbp),%rsi
  mov    -0xd8(%rbp),%rax
  lea    -0x58(%rbp),%r8
  push   %r8
  mov    -0x6c(%rbp),%r8d
  push   %r8
  mov    -0x98(%rbp),%r8d
  push   %r8
  mov    -0xa8(%rbp),%r8d
  push   %r8
  movq   %rdi,%xmm0
  mov    %r10,%r8
  mov    %rax,%rdi
  mov    $0x1,%eax
  call   <gee@plt>
  add    $0x20,%rsp
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x38(%rbp),%xmm0
  pxor   %xmm1,%xmm1
  comisd %xmm1,%xmm0
  ja     L371
  cmpl   $0x0,-0xb8(%rbp)
  jg     L371
  subl   $0x1,-0x98(%rbp)
L231:
  mov    -0x98(%rbp),%eax
  cmp    -0xac(%rbp),%eax
  jg     L152
L371:
  mov    -0x98(%rbp),%eax
  cmp    -0xac(%rbp),%eax
  jg     L696
  mov    -0xa0(%rbp),%eax
  cltq
  shl    $0x2,%rax
  lea    -0x4(%rax),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  add    $0x1,%eax
  mov    %eax,-0x98(%rbp)
  cmpl   $0x0,-0xb8(%rbp)
  jle    L439
  mov    -0xa0(%rbp),%eax
  cltq
  shl    $0x2,%rax
  lea    -0x4(%rax),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    -0x98(%rbp),%edx
  sub    $0x1,%edx
  mov    %edx,(%rax)
  addl   $0x1,-0xa0(%rbp)
  mov    -0x98(%rbp),%eax
  sub    $0x1,%eax
  mov    %eax,-0xac(%rbp)
  addl   $0x1,-0xb8(%rbp)
  jmp    L772
L439:
  addl   $0x1,-0x98(%rbp)
  jmp    L850
L558:
  mov    -0x50(%rbp),%rdi
  mov    -0xe8(%rbp),%r9
  mov    -0xe0(%rbp),%r10
  mov    -0xf8(%rbp),%rcx
  mov    -0xf0(%rbp),%rdx
  mov    -0x100(%rbp),%rsi
  mov    -0xd8(%rbp),%rax
  lea    -0x58(%rbp),%r8
  push   %r8
  mov    -0x6c(%rbp),%r8d
  push   %r8
  mov    -0x98(%rbp),%r8d
  push   %r8
  mov    -0xa8(%rbp),%r8d
  push   %r8
  movq   %rdi,%xmm0
  mov    %r10,%r8
  mov    %rax,%rdi
  mov    $0x1,%eax
  call   <gee@plt>
  add    $0x20,%rsp
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x38(%rbp),%xmm0
  pxor   %xmm1,%xmm1
  comisd %xmm1,%xmm0
  ja     L810
  addl   $0x1,-0x98(%rbp)
L850:
  mov    -0x98(%rbp),%eax
  cmp    -0x60(%rbp),%eax
  jl     L558
  jmp    L600
L810:
  nop
L600:
  mov    -0x98(%rbp),%eax
  cmp    -0x60(%rbp),%eax
  jge    L344
  movsd  -0x38(%rbp),%xmm0
  pxor   %xmm1,%xmm1
  comisd %xmm1,%xmm0
  jbe    L344
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0xa8(%rbp),%xmm1
  movsd  -0x58(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  addl   $0x1,-0x98(%rbp)
  jmp    L922
L418:
  mov    -0x50(%rbp),%rdi
  mov    -0xe8(%rbp),%r9
  mov    -0xe0(%rbp),%r10
  mov    -0xf8(%rbp),%rcx
  mov    -0xf0(%rbp),%rdx
  mov    -0x100(%rbp),%rsi
  mov    -0xd8(%rbp),%rax
  lea    -0x58(%rbp),%r8
  push   %r8
  mov    -0x6c(%rbp),%r8d
  push   %r8
  mov    -0x98(%rbp),%r8d
  push   %r8
  mov    -0xa8(%rbp),%r8d
  push   %r8
  movq   %rdi,%xmm0
  mov    %r10,%r8
  mov    %rax,%rdi
  mov    $0x1,%eax
  call   <gee@plt>
  add    $0x20,%rsp
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  pxor   %xmm0,%xmm0
  comisd -0x38(%rbp),%xmm0
  jae    L774
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0xa8(%rbp),%xmm1
  movsd  -0x58(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  addl   $0x1,-0x98(%rbp)
L922:
  mov    -0x98(%rbp),%eax
  cmp    -0x60(%rbp),%eax
  jl     L418
  jmp    L695
L774:
  nop
L695:
  mov    -0xa0(%rbp),%eax
  cltq
  shl    $0x2,%rax
  lea    -0x4(%rax),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    -0x98(%rbp),%edx
  sub    $0x1,%edx
  mov    %edx,(%rax)
  addl   $0x1,-0xa0(%rbp)
  mov    -0x98(%rbp),%eax
  sub    $0x1,%eax
  mov    %eax,-0xac(%rbp)
  addl   $0x1,-0xb8(%rbp)
  jmp    L772
L344:
  mov    -0xa0(%rbp),%eax
  cltq
  shl    $0x2,%rax
  lea    -0x4(%rax),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0xac(%rbp)
  addl   $0x1,-0xa0(%rbp)
  jmp    L772
L696:
  movsd  -0x38(%rbp),%xmm0
  pxor   %xmm1,%xmm1
  comisd %xmm1,%xmm0
  jbe    L127
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0xa8(%rbp),%xmm1
  movsd  -0x58(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  subl   $0x1,-0x98(%rbp)
  jmp    L147
L739:
  mov    -0x50(%rbp),%rdi
  mov    -0xe8(%rbp),%r9
  mov    -0xe0(%rbp),%r10
  mov    -0xf8(%rbp),%rcx
  mov    -0xf0(%rbp),%rdx
  mov    -0x100(%rbp),%rsi
  mov    -0xd8(%rbp),%rax
  lea    -0x58(%rbp),%r8
  push   %r8
  mov    -0x6c(%rbp),%r8d
  push   %r8
  mov    -0x98(%rbp),%r8d
  push   %r8
  mov    -0xa8(%rbp),%r8d
  push   %r8
  movq   %rdi,%xmm0
  mov    %r10,%r8
  mov    %rax,%rdi
  mov    $0x1,%eax
  call   <gee@plt>
  add    $0x20,%rsp
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  pxor   %xmm0,%xmm0
  comisd -0x38(%rbp),%xmm0
  jae    L450
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0xa8(%rbp),%xmm1
  movsd  -0x58(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  subl   $0x1,-0x98(%rbp)
L147:
  mov    -0x98(%rbp),%eax
  cmp    -0xac(%rbp),%eax
  jg     L739
  jmp    L780
L450:
  nop
L780:
  mov    -0xa0(%rbp),%eax
  cltq
  shl    $0x2,%rax
  lea    -0x4(%rax),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0xac(%rbp)
  mov    -0xa0(%rbp),%eax
  cltq
  shl    $0x2,%rax
  lea    -0x4(%rax),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    -0x98(%rbp),%edx
  add    $0x1,%edx
  mov    %edx,(%rax)
  addl   $0x1,-0xa0(%rbp)
  addl   $0x1,-0xb8(%rbp)
  jmp    L772
L127:
  mov    -0xa0(%rbp),%eax
  cltq
  shl    $0x2,%rax
  lea    -0x4(%rax),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0xac(%rbp)
  addl   $0x1,-0xa0(%rbp)
L772:
  mov    -0xa0(%rbp),%eax
  cmp    -0x68(%rbp),%eax
  jl     L124
  jmp    L511
L27:
  mov    -0xd0(%rbp),%rax
  mov    (%rax),%eax
  mov    -0xb4(%rbp),%edx
  cmp    %eax,%edx
  cmovle %edx,%eax
  mov    %eax,-0xb4(%rbp)
  mov    -0x6c(%rbp),%eax
  cltq
  shl    $0x2,%rax
  lea    -0x4(%rax),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    -0xb0(%rbp),%edx
  cmp    %eax,%edx
  cmovge %edx,%eax
  mov    %eax,-0xb0(%rbp)
  mov    -0xa8(%rbp),%eax
  mov    %eax,-0xa4(%rbp)
  mov    -0x5c(%rbp),%eax
  neg    %eax
  mov    %eax,-0xa8(%rbp)
  movl   $0x0,-0x98(%rbp)
  mov    -0xb4(%rbp),%eax
  mov    %eax,-0x60(%rbp)
  movsd  -0x50(%rbp),%xmm1
  movsd  D79(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  movl   $0x1,-0xbc(%rbp)
  movl   $0x1,-0x9c(%rbp)
  jmp    L134
L870:
  mov    -0x50(%rbp),%rdi
  mov    -0xe8(%rbp),%r9
  mov    -0xe0(%rbp),%r10
  mov    -0xf8(%rbp),%rcx
  mov    -0xf0(%rbp),%rdx
  mov    -0x100(%rbp),%rsi
  mov    -0xd8(%rbp),%rax
  lea    -0x58(%rbp),%r8
  push   %r8
  mov    -0x6c(%rbp),%r8d
  push   %r8
  mov    -0x98(%rbp),%r8d
  push   %r8
  mov    -0xa8(%rbp),%r8d
  push   %r8
  movq   %rdi,%xmm0
  mov    %r10,%r8
  mov    %rax,%rdi
  mov    $0x1,%eax
  call   <gee@plt>
  add    $0x20,%rsp
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x38(%rbp),%xmm0
  pxor   %xmm1,%xmm1
  comisd %xmm1,%xmm0
  jbe    L308
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0xa8(%rbp),%xmm1
  movsd  -0x58(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  mov    -0x68(%rbp),%eax
  mov    %eax,-0xa0(%rbp)
  addl   $0x1,-0x98(%rbp)
  jmp    L461
L85:
  mov    -0x50(%rbp),%rdi
  mov    -0xe8(%rbp),%r9
  mov    -0xe0(%rbp),%r10
  mov    -0xf8(%rbp),%rcx
  mov    -0xf0(%rbp),%rdx
  mov    -0x100(%rbp),%rsi
  mov    -0xd8(%rbp),%rax
  lea    -0x58(%rbp),%r8
  push   %r8
  mov    -0x6c(%rbp),%r8d
  push   %r8
  mov    -0x98(%rbp),%r8d
  push   %r8
  mov    -0xa8(%rbp),%r8d
  push   %r8
  movq   %rdi,%xmm0
  mov    %r10,%r8
  mov    %rax,%rdi
  mov    $0x1,%eax
  call   <gee@plt>
  add    $0x20,%rsp
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  pxor   %xmm0,%xmm0
  comisd -0x38(%rbp),%xmm0
  jae    L172
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0xa8(%rbp),%xmm1
  movsd  -0x58(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  addl   $0x1,-0x98(%rbp)
L461:
  mov    -0x98(%rbp),%eax
  cmp    -0x60(%rbp),%eax
  jl     L85
  jmp    L497
L172:
  nop
L497:
  mov    -0xa0(%rbp),%eax
  cltq
  shl    $0x2,%rax
  lea    -0x4(%rax),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    -0x98(%rbp),%edx
  sub    $0x1,%edx
  mov    %edx,(%rax)
  addl   $0x1,-0xa0(%rbp)
  mov    -0x98(%rbp),%eax
  sub    $0x1,%eax
  mov    %eax,-0xac(%rbp)
  addl   $0x1,-0xb8(%rbp)
  movl   $0x0,-0x9c(%rbp)
  jmp    L841
L308:
  addl   $0x1,-0x98(%rbp)
L134:
  mov    -0x98(%rbp),%eax
  cmp    -0x60(%rbp),%eax
  jl     L870
L841:
  cmpl   $0x0,-0x9c(%rbp)
  je     L511
  movl   $0x6,-0xbc(%rbp)
  mov    -0xc8(%rbp),%eax
  mov    %eax,-0x98(%rbp)
  mov    -0xb0(%rbp),%eax
  mov    %eax,-0xac(%rbp)
  mov    -0xa4(%rbp),%eax
  neg    %eax
  mov    %eax,-0xa8(%rbp)
  movsd  -0x40(%rbp),%xmm1
  movsd  D79(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  subl   $0x1,-0x98(%rbp)
  jmp    L19
L391:
  mov    -0x50(%rbp),%rdi
  mov    -0xe8(%rbp),%r9
  mov    -0xe0(%rbp),%r10
  mov    -0xf8(%rbp),%rcx
  mov    -0xf0(%rbp),%rdx
  mov    -0x100(%rbp),%rsi
  mov    -0xd8(%rbp),%rax
  lea    -0x58(%rbp),%r8
  push   %r8
  mov    -0x6c(%rbp),%r8d
  push   %r8
  mov    -0x98(%rbp),%r8d
  push   %r8
  mov    -0xa8(%rbp),%r8d
  push   %r8
  movq   %rdi,%xmm0
  mov    %r10,%r8
  mov    %rax,%rdi
  mov    $0x1,%eax
  call   <gee@plt>
  add    $0x20,%rsp
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x38(%rbp),%xmm0
  pxor   %xmm1,%xmm1
  comisd %xmm1,%xmm0
  jbe    L336
  mov    -0x68(%rbp),%eax
  mov    %eax,-0xa0(%rbp)
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0xa8(%rbp),%xmm1
  movsd  -0x58(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  addl   $0xa,-0xbc(%rbp)
  subl   $0x1,-0x98(%rbp)
  jmp    L661
L437:
  mov    -0x50(%rbp),%rdi
  mov    -0xe8(%rbp),%r9
  mov    -0xe0(%rbp),%r10
  mov    -0xf8(%rbp),%rcx
  mov    -0xf0(%rbp),%rdx
  mov    -0x100(%rbp),%rsi
  mov    -0xd8(%rbp),%rax
  lea    -0x58(%rbp),%r8
  push   %r8
  mov    -0x6c(%rbp),%r8d
  push   %r8
  mov    -0x98(%rbp),%r8d
  push   %r8
  mov    -0xa8(%rbp),%r8d
  push   %r8
  movq   %rdi,%xmm0
  mov    %r10,%r8
  mov    %rax,%rdi
  mov    $0x1,%eax
  call   <gee@plt>
  add    $0x20,%rsp
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  pxor   %xmm0,%xmm0
  comisd -0x38(%rbp),%xmm0
  jae    L881
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0xa8(%rbp),%xmm1
  movsd  -0x58(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  subl   $0x1,-0x98(%rbp)
L661:
  mov    -0x98(%rbp),%eax
  cmp    -0xac(%rbp),%eax
  jg     L437
  jmp    L647
L881:
  nop
L647:
  mov    -0xa0(%rbp),%eax
  cltq
  shl    $0x2,%rax
  lea    -0x4(%rax),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0xac(%rbp)
  mov    -0xa0(%rbp),%eax
  cltq
  shl    $0x2,%rax
  lea    -0x4(%rax),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    -0x98(%rbp),%edx
  add    $0x1,%edx
  mov    %edx,(%rax)
  addl   $0x1,-0xa0(%rbp)
  addl   $0x1,-0xb8(%rbp)
  movl   $0x0,-0x9c(%rbp)
  jmp    L959
L336:
  subl   $0x1,-0x98(%rbp)
L19:
  mov    -0x98(%rbp),%eax
  cmp    -0xac(%rbp),%eax
  jg     L391
L959:
  cmpl   $0x0,-0x9c(%rbp)
  je     L511
  cmpl   $0x6,-0xbc(%rbp)
  je     L355
  mov    -0xc4(%rbp),%eax
  sub    $0x1,%eax
  mov    %eax,-0x7c(%rbp)
  jmp    L476
L583:
  mov    -0x7c(%rbp),%eax
  cltq
  shl    $0x2,%rax
  lea    -0x4(%rax),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    -0x7c(%rbp),%edx
  movslq %edx,%rdx
  lea    0x0(,%rdx,4),%rcx
  mov    -0xd0(%rbp),%rdx
  add    %rcx,%rdx
  mov    (%rax),%eax
  mov    %eax,(%rdx)
  subl   $0x1,-0x7c(%rbp)
L476:
  cmpl   $0x0,-0x7c(%rbp)
  jg     L583
  mov    -0xd0(%rbp),%rax
  mov    -0xb4(%rbp),%edx
  mov    %edx,(%rax)
  addl   $0x1,-0xb8(%rbp)
  jmp    L355
L511:
  mov    -0xa0(%rbp),%eax
  cmp    -0x68(%rbp),%eax
  je     L27
  jmp    L895
L355:
  nop
L895:
  mov    -0xa0(%rbp),%eax
  cmp    -0x68(%rbp),%eax
  jle    L45
  cmpl   $0x9,-0xbc(%rbp)
  jle    L728
  movl   $0x0,-0x78(%rbp)
  jmp    L909
L116:
  mov    -0x78(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    -0x78(%rbp),%edx
  movslq %edx,%rdx
  lea    0x0(,%rdx,4),%rcx
  mov    -0xd0(%rbp),%rdx
  add    %rcx,%rdx
  mov    (%rax),%eax
  mov    %eax,(%rdx)
  addl   $0x1,-0x78(%rbp)
L909:
  mov    -0x78(%rbp),%eax
  cmp    -0xc4(%rbp),%eax
  jl     L116
  mov    -0x6c(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    -0xc4(%rbp),%edx
  movslq %edx,%rdx
  lea    0x0(,%rdx,4),%rcx
  mov    -0xd0(%rbp),%rdx
  add    %rcx,%rdx
  mov    (%rax),%eax
  mov    %eax,(%rdx)
  addl   $0x1,-0xb8(%rbp)
  jmp    L45
L728:
  movsd  -0x40(%rbp),%xmm0
  comisd -0x50(%rbp),%xmm0
  jbe    L323
  movsd  -0x40(%rbp),%xmm0
  jmp    L785
L323:
  movsd  -0x50(%rbp),%xmm0
L785:
  movsd  %xmm0,-0x40(%rbp)
  mov    -0x68(%rbp),%eax
  cltq
  shl    $0x2,%rax
  lea    -0x4(%rax),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0xb4(%rbp)
  mov    -0xc8(%rbp),%eax
  mov    %eax,-0x98(%rbp)
  mov    -0xb0(%rbp),%eax
  mov    %eax,-0xac(%rbp)
  mov    -0xa4(%rbp),%eax
  neg    %eax
  mov    %eax,-0xa8(%rbp)
  movsd  -0x40(%rbp),%xmm1
  movsd  D79(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  subl   $0x1,-0x98(%rbp)
  jmp    L666
L587:
  mov    -0x50(%rbp),%rdi
  mov    -0xe8(%rbp),%r9
  mov    -0xe0(%rbp),%r10
  mov    -0xf8(%rbp),%rcx
  mov    -0xf0(%rbp),%rdx
  mov    -0x100(%rbp),%rsi
  mov    -0xd8(%rbp),%rax
  lea    -0x58(%rbp),%r8
  push   %r8
  mov    -0x6c(%rbp),%r8d
  push   %r8
  mov    -0x98(%rbp),%r8d
  push   %r8
  mov    -0xa8(%rbp),%r8d
  push   %r8
  movq   %rdi,%xmm0
  mov    %r10,%r8
  mov    %rax,%rdi
  mov    $0x1,%eax
  call   <gee@plt>
  add    $0x20,%rsp
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x38(%rbp),%xmm0
  pxor   %xmm1,%xmm1
  comisd %xmm1,%xmm0
  jbe    L465
  mov    -0x68(%rbp),%eax
  mov    %eax,-0xa0(%rbp)
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0xa8(%rbp),%xmm1
  movsd  -0x58(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  addl   $0xa,-0xbc(%rbp)
  subl   $0x1,-0x98(%rbp)
  jmp    L939
L94:
  mov    -0x50(%rbp),%rdi
  mov    -0xe8(%rbp),%r9
  mov    -0xe0(%rbp),%r10
  mov    -0xf8(%rbp),%rcx
  mov    -0xf0(%rbp),%rdx
  mov    -0x100(%rbp),%rsi
  mov    -0xd8(%rbp),%rax
  lea    -0x58(%rbp),%r8
  push   %r8
  mov    -0x6c(%rbp),%r8d
  push   %r8
  mov    -0x98(%rbp),%r8d
  push   %r8
  mov    -0xa8(%rbp),%r8d
  push   %r8
  movq   %rdi,%xmm0
  mov    %r10,%r8
  mov    %rax,%rdi
  mov    $0x1,%eax
  call   <gee@plt>
  add    $0x20,%rsp
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  pxor   %xmm0,%xmm0
  comisd -0x38(%rbp),%xmm0
  jae    L688
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0xa8(%rbp),%xmm1
  movsd  -0x58(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  subl   $0x1,-0x98(%rbp)
L939:
  mov    -0x98(%rbp),%eax
  cmp    -0xac(%rbp),%eax
  jg     L94
  jmp    L326
L688:
  nop
L326:
  mov    -0xa0(%rbp),%eax
  cltq
  shl    $0x2,%rax
  lea    -0x4(%rax),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0xac(%rbp)
  mov    -0xa0(%rbp),%eax
  cltq
  shl    $0x2,%rax
  lea    -0x4(%rax),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    -0x98(%rbp),%edx
  add    $0x1,%edx
  mov    %edx,(%rax)
  addl   $0x1,-0xa0(%rbp)
  addl   $0x1,-0xb8(%rbp)
  movl   $0x0,-0x74(%rbp)
  jmp    L212
L789:
  mov    -0x74(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    -0x74(%rbp),%edx
  movslq %edx,%rdx
  lea    0x0(,%rdx,4),%rcx
  mov    -0xd0(%rbp),%rdx
  add    %rcx,%rdx
  mov    (%rax),%eax
  mov    %eax,(%rdx)
  addl   $0x1,-0x74(%rbp)
L212:
  mov    -0x74(%rbp),%eax
  cmp    -0xc4(%rbp),%eax
  jl     L789
  mov    -0x6c(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    -0xc4(%rbp),%edx
  movslq %edx,%rdx
  lea    0x0(,%rdx,4),%rcx
  mov    -0xd0(%rbp),%rdx
  add    %rcx,%rdx
  mov    (%rax),%eax
  mov    %eax,(%rdx)
  jmp    L403
L465:
  subl   $0x1,-0x98(%rbp)
L666:
  mov    -0x98(%rbp),%eax
  cmp    -0xac(%rbp),%eax
  jg     L587
L403:
  cmpl   $0x6,-0xbc(%rbp)
  je     L45
  mov    -0xc4(%rbp),%eax
  mov    %eax,-0x70(%rbp)
  jmp    L318
L956:
  mov    -0x70(%rbp),%eax
  cltq
  shl    $0x2,%rax
  lea    -0x4(%rax),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    -0x70(%rbp),%edx
  movslq %edx,%rdx
  lea    0x0(,%rdx,4),%rcx
  mov    -0xd0(%rbp),%rdx
  add    %rcx,%rdx
  mov    (%rax),%eax
  mov    %eax,(%rdx)
  subl   $0x1,-0x70(%rbp)
L318:
  cmpl   $0x0,-0x70(%rbp)
  jg     L956
  mov    -0xd0(%rbp),%rax
  mov    -0xb4(%rbp),%edx
  mov    %edx,(%rax)
  addl   $0x1,-0xb8(%rbp)
L45:
  cmpl   $0x0,-0xb8(%rbp)
  jg     L827
  jmp    L522
L234:
  nop
L522:
  mov    $0x0,%eax
L321:
  mov    -0x18(%rbp),%rdx
  sub    %fs:0x28,%rdx
  je     L123
  call   <__stack_chk_fail@plt>
L123:
  mov    -0x8(%rbp),%rbx
  leave
  ret
```

tool_calls: name: read_data
arguments: data_tag: D73
datatype: f32[2]
name: read_data
arguments: data_tag: D20
datatype: f64
name: read_data
arguments: data_tag: D53
datatype: f64
name: read_data
arguments: data_tag: D50
datatype: string
name: read_data
arguments: data_tag: D22
datatype: string
name: read_data
arguments: data_tag: D79
datatype: f64

tool: {"D73": [0.0,-1.875]}

tool: D20: 1

tool: D53: 6.283185307179586

tool: D50: ERROR: remez exchange algorithm do not converge

tool: D22: ERROR: Filter designed will probably be incorrect

tool: D79: 1.00001

assistant: ```c
int echange(int nfcns, int taillegrilledense, int *iExt, double *FreqGrille, double *des, double *wt, double *x, double *y, double *ad) {
  double comp = 0.0;
  double deltal = -1;
  double delta = 0.0;
  int nz = nfcns + 1;
  int nzz = nz + 1;
  int iter = 0;
  int luck = 0;
  int jchnge = 1;
  int jet = floor((nfcns - 1) / 15) + 1;
  double y1 = 0.0;
  double E = 0.0;
  double dtemp = 0.0;
  int k1 = 0;
  int knz = 0;
  int klow = 0;
  int kup = 0;
  int nu = 0;
  int nut = 0;
  int nut1 = 0;
  int j = 1;
  int flag = 0;
  int l = 0;
  while (jchnge > 0) {
    iter++;
    if (iter > 250)
      break;
    for (int i = 0; i < nz; i++)
      x[i] = cos(6.2831853071795865 * FreqGrille[iExt[i]]);
    for (int nn = 0; nn < nz; nn++)
      ad[nn] = dd(nn, nz, jet, x);
    {
      double dnum = 0.0;
      double dden = 0.0;
      int add = 1;
      for (int i = 0; i < nz; i++) {
        dnum += ad[i] * des[iExt[i]];
        dden += add * ad[i] / wt[iExt[i]];
        add = -add;
      }
      delta = dnum / dden;
    }
    nu = (delta > 0 ? -1 : 1);
    delta *= -nu;
    {
      int add = 1;
      for (int i = 0; i < nz; i++) {
        y[i] = des[iExt[i]] + nu * delta * add / wt[iExt[i]];
        add = -add;
      }
    }
    if (delta <= deltal) {
      printf("ERROR: remez exchange algorithm do not converge\n");
      printf("ERROR: Filter designed will probably be incorrect\n");
      return -1;
    }
    deltal = delta;
    jchnge = 0;
    k1 = iExt[0];
    knz = iExt[nz - 1];
    klow = -1;
    nut = -nu;
    j = 1;
    while (j < nzz) {
      kup = iExt[j];
      l = iExt[j - 1] + 1;
      nut = -nut;
      if (j == 2)
        y1 = comp;
      comp = delta;
      flag = 1;
      if (l < kup) {
        dtemp = gee(FreqGrille, ad, x, y, des, wt, comp, nut, l, nz, &E);
        if (dtemp > 0) {
          comp = nut * E;
          l++;
          while (l < kup) {
            dtemp = gee(FreqGrille, ad, x, y, des, wt, comp, nut, l, nz, &E);
            if (dtemp <= 0)
              break;
            comp = nut * E;
            l++;
          }
          iExt[j - 1] = l - 1;
          j++;
          klow = l - 1;
          jchnge++;
          flag = 0;
        }
      }
      if (flag) {
        l -= 2;
        while (l > klow) {
          dtemp = gee(FreqGrille, ad, x, y, des, wt, comp, nut, l, nz, &E);
          if (dtemp > 0 || jchnge > 0)
            break;
          l--;
        }
        if (l <= klow) {
          l = iExt[j - 1] + 1;
          if (jchnge > 0) {
            iExt[j - 1] = l - 1;
            j++;
            klow = l - 1;
            jchnge++;
          } else {
            l++;
            while (l < kup) {
              dtemp = gee(FreqGrille, ad, x, y, des, wt, comp, nut, l, nz, &E);
              if (dtemp > 0)
                break;
              l++;
            }
            if (l < kup && dtemp > 0) {
              comp = nut * E;
              l++;
              while (l < kup) {
                dtemp = gee(FreqGrille, ad, x, y, des, wt, comp, nut, l, nz, &E);
                if (dtemp <= 0)
                  break;
                comp = nut * E;
                l++;
              }
              iExt[j - 1] = l - 1;
              j++;
              klow = l - 1;
              jchnge++;
            } else {
              klow = iExt[j - 1];
              j++;
            }
          }
        } else if (dtemp > 0) {
          comp = nut * E;
          l--;
          while (l > klow) {
            dtemp = gee(FreqGrille, ad, x, y, des, wt, comp, nut, l, nz, &E);
            if (dtemp <= 0)
              break;
            comp = nut * E;
            l--;
          }
          klow = iExt[j - 1];
          iExt[j - 1] = l + 1;
          j++;
          jchnge++;
        } else {
          klow = iExt[j - 1];
          j++;
        }
      }
    }
    while (j == nzz) {
      k1 = ((k1) < (iExt[0]) ? (k1) : (iExt[0]));
      knz = ((knz) > (iExt[nz - 1]) ? (knz) : (iExt[nz - 1]));
      nut1 = nut;
      nut = -nu;
      l = 0;
      kup = k1;
      comp = comp * 1.00001;
      luck = 1;
      flag = 1;
      while (l < kup) {
        dtemp = gee(FreqGrille, ad, x, y, des, wt, comp, nut, l, nz, &E);
        if (dtemp > 0) {
          comp = nut * E;
          j = nzz;
          l++;
          while (l < kup) {
            dtemp = gee(FreqGrille, ad, x, y, des, wt, comp, nut, l, nz, &E);
            if (dtemp <= 0)
              break;
            comp = nut * E;
            l++;
          }
          iExt[j - 1] = l - 1;
          j++;
          klow = l - 1;
          jchnge++;
          flag = 0;
          break;
        }
        l++;
      }
      if (flag) {
        luck = 6;
        l = taillegrilledense;
        klow = knz;
        nut = -nut1;
        comp = y1 * 1.00001;
        l--;
        while (l > klow) {
          dtemp = gee(FreqGrille, ad, x, y, des, wt, comp, nut, l, nz, &E);
          if (dtemp > 0) {
            j = nzz;
            comp = nut * E;
            luck += 10;
            l--;
            while (l > klow) {
              dtemp = gee(FreqGrille, ad, x, y, des, wt, comp, nut, l, nz, &E);
              if (dtemp <= 0)
                break;
              comp = nut * E;
              l--;
            }
            klow = iExt[j - 1];
            iExt[j - 1] = l + 1;
            j++;
            jchnge++;
            flag = 0;
            break;
          }
          l--;
        }
        if (flag) {
          if (luck != 6) {
            for (int i = nfcns - 1; i > 0; i--)
              iExt[i] = iExt[i - 1];
            iExt[0] = k1;
            jchnge++;
          }
          break;
        }
      }
    }
    if (j > nzz) {
      if (luck > 9) {
        for (int i = 0; i < nfcns; i++)
          iExt[i] = iExt[i + 1];
        iExt[nfcns] = iExt[nz];
        jchnge++;
      } else {
        y1 = ((y1) > (comp) ? (y1) : (comp));
        k1 = iExt[nzz - 1];
        l = taillegrilledense;
        klow = knz;
        nut = -nut1;
        comp = y1 * 1.00001;
        l--;
        while (l > klow) {
          dtemp = gee(FreqGrille, ad, x, y, des, wt, comp, nut, l, nz, &E);
          if (dtemp > 0) {
            j = nzz;
            comp = nut * E;
            luck += 10;
            l--;
            while (l > klow) {
              dtemp = gee(FreqGrille, ad, x, y, des, wt, comp, nut, l, nz, &E);
              if (dtemp <= 0)
                break;
              comp = nut * E;
              l--;
            }
            klow = iExt[j - 1];
            iExt[j - 1] = l + 1;
            j++;
            jchnge++;
            for (int i = 0; i < nfcns; i++)
              iExt[i] = iExt[i + 1];
            iExt[nfcns] = iExt[nz];
            break;
          }
          l--;
        }
        if (luck != 6) {
          for (int i = nfcns; i > 0; i--)
            iExt[i] = iExt[i - 1];
          iExt[0] = k1;
          jchnge++;
        }
      }
    }
  }
  return 0;
}
```