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

```asm
<sqfun>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %rbx
  sub    $0x3a8,%rsp
  mov    %edi,-0x394(%rbp)
  mov    %esi,-0x398(%rbp)
  mov    %rdx,-0x3a0(%rbp)
  movsd  D356(%rip),%xmm0
  movsd  %xmm0,-0x318(%rbp)
  movl   $0x28,-0x37c(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x370(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x368(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x360(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x358(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x350(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x348(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x340(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x338(%rbp)
  mov    -0x3a0(%rbp),%rax
  movsd  0x80(%rax),%xmm0
  movsd  %xmm0,-0x310(%rbp)
  movsd  -0x310(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x308(%rbp)
  movsd  -0x308(%rbp),%xmm0
  mulsd  -0x310(%rbp),%xmm0
  movsd  %xmm0,-0x300(%rbp)
  movsd  -0x310(%rbp),%xmm1
  movsd  D170(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x2f8(%rbp)
  movsd  -0x2f8(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  movsd  %xmm0,-0x2f0(%rbp)
  mov    -0x3a0(%rbp),%rax
  add    $0x20,%rax
  movsd  (%rax),%xmm0
  mov    -0x3a0(%rbp),%rax
  sub    $0xffffffffffffff80,%rax
  movsd  (%rax),%xmm1
  divsd  %xmm1,%xmm0
  movq   %xmm0,%rax
  mov    -0x3a0(%rbp),%rdx
  lea    0x68(%rdx),%rbx
  movsd  D890(%rip),%xmm0
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  call   <pow@plt>
  movq   %xmm0,%rax
  mov    %rax,(%rbx)
  mov    -0x3a0(%rbp),%rax
  add    $0x30,%rax
  movsd  (%rax),%xmm0
  mov    -0x3a0(%rbp),%rax
  add    $0x68,%rax
  movsd  (%rax),%xmm1
  mov    -0x3a0(%rbp),%rax
  add    $0x60,%rax
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movl   $0x0,-0x380(%rbp)
  mov    -0x3a0(%rbp),%rax
  add    $0x60,%rax
  movsd  (%rax),%xmm0
  comisd D88(%rip),%xmm0
  jbe    L46
  cmpl   $0x1,-0x394(%rbp)
  jne    L46
  movl   $0x1,-0x380(%rbp)
L46:
  mov    -0x3a0(%rbp),%rax
  add    $0x28,%rax
  movsd  (%rax),%xmm1
  mov    -0x3a0(%rbp),%rax
  add    $0x68,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  %xmm1,-0x3a8(%rbp)
  mov    -0x3a0(%rbp),%rax
  add    $0x30,%rax
  movsd  (%rax),%xmm0
  mov    -0x3a0(%rbp),%rax
  add    $0x60,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <exp@plt>
  mov    -0x3a0(%rbp),%rax
  add    $0x58,%rax
  mulsd  -0x3a8(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x3a0(%rbp),%rax
  movsd  0x58(%rax),%xmm0
  movsd  %xmm0,-0x2e8(%rbp)
  mov    -0x3a0(%rbp),%rax
  movsd  0x60(%rax),%xmm0
  movsd  %xmm0,-0x2e0(%rbp)
  movsd  -0x2e0(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x2d8(%rbp)
  movsd  -0x2e0(%rbp),%xmm1
  movsd  D131(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x2d0(%rbp)
  movsd  D131(%rip),%xmm0
  divsd  -0x2d8(%rbp),%xmm0
  movsd  %xmm0,-0x2c8(%rbp)
  movsd  -0x2c8(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x2c0(%rbp)
  movsd  D131(%rip),%xmm0
  subsd  -0x310(%rbp),%xmm0
  movsd  %xmm0,-0x2b8(%rbp)
  movsd  -0x2b8(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x2b0(%rbp)
  movsd  -0x2b8(%rbp),%xmm0
  divsd  -0x2e0(%rbp),%xmm0
  movsd  %xmm0,-0x2a8(%rbp)
  movsd  D131(%rip),%xmm0
  divsd  -0x2b0(%rbp),%xmm0
  movsd  %xmm0,-0x2a0(%rbp)
  movsd  -0x2a0(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x298(%rbp)
  movsd  -0x298(%rbp),%xmm1
  movsd  D65(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x290(%rbp)
  movsd  -0x310(%rbp),%xmm1
  movsd  D636(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x288(%rbp)
  movsd  -0x288(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  movsd  %xmm0,-0x280(%rbp)
  movsd  -0x300(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  movsd  %xmm0,-0x278(%rbp)
  movsd  -0x310(%rbp),%xmm1
  movsd  D336(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x270(%rbp)
  movsd  -0x270(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x268(%rbp)
  movsd  -0x310(%rbp),%xmm0
  movapd %xmm0,%xmm1
  addsd  %xmm0,%xmm1
  movsd  D131(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x260(%rbp)
  movsd  -0x260(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x258(%rbp)
  movsd  -0x260(%rbp),%xmm0
  movq   D538(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  -0x2a8(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x250(%rbp)
  movsd  -0x308(%rbp),%xmm1
  movsd  D147(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x310(%rbp),%xmm2
  movsd  D209(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  movsd  D131(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movsd  -0x2c8(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x248(%rbp)
  movsd  -0x308(%rbp),%xmm1
  movsd  D648(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x2c0(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x240(%rbp)
  movsd  -0x310(%rbp),%xmm1
  movsd  D140(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  addsd  -0x308(%rbp),%xmm1
  movsd  D131(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movq   D538(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  -0x2a8(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x238(%rbp)
  movsd  -0x310(%rbp),%xmm1
  movsd  D108(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x310(%rbp),%xmm2
  movsd  D209(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  -0x308(%rbp),%xmm0
  movsd  D336(%rip),%xmm2
  subsd  %xmm2,%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x2c8(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x230(%rbp)
  movsd  -0x310(%rbp),%xmm1
  movsd  D170(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x308(%rbp),%xmm0
  movapd %xmm0,%xmm2
  addsd  %xmm0,%xmm2
  movsd  -0x310(%rbp),%xmm3
  movsd  D4(%rip),%xmm0
  mulsd  %xmm3,%xmm0
  addsd  %xmm2,%xmm0
  movsd  D131(%rip),%xmm2
  subsd  %xmm2,%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x2c0(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x228(%rbp)
  movsd  -0x308(%rbp),%xmm1
  movsd  D636(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  addsd  -0x300(%rbp),%xmm1
  movsd  -0x310(%rbp),%xmm2
  movsd  D667(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movsd  D573(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movq   D538(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  -0x2a8(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x220(%rbp)
  movsd  -0x308(%rbp),%xmm1
  movsd  D636(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  addsd  -0x278(%rbp),%xmm1
  movsd  -0x310(%rbp),%xmm2
  movsd  D825(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  D37(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movsd  -0x2c8(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x218(%rbp)
  movsd  -0x310(%rbp),%xmm1
  movsd  D480(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  -0x278(%rbp),%xmm0
  movsd  D573(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movsd  -0x2c0(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x210(%rbp)
  movsd  -0x2f0(%rbp),%xmm0
  mulsd  -0x2e0(%rbp),%xmm0
  mulsd  -0x210(%rbp),%xmm0
  movsd  -0x220(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x208(%rbp)
  movsd  -0x210(%rbp),%xmm1
  movsd  D209(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  -0x218(%rbp),%xmm0
  movsd  -0x280(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x200(%rbp)
  movsd  -0x280(%rbp),%xmm0
  divsd  -0x2e0(%rbp),%xmm0
  movsd  %xmm0,-0x1f8(%rbp)
  movsd  -0x2f8(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x2c8(%rbp),%xmm1
  movsd  -0x2b8(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x1f0(%rbp)
  movsd  -0x310(%rbp),%xmm1
  movsd  D573(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  -0x2e0(%rbp),%xmm2
  movsd  D573(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x1e8(%rbp)
  movsd  -0x270(%rbp),%xmm0
  mulsd  -0x2c8(%rbp),%xmm0
  movsd  %xmm0,-0x1e0(%rbp)
  movsd  -0x2f8(%rbp),%xmm0
  movq   D538(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x2e8(%rbp),%xmm1
  movsd  -0x1e8(%rbp),%xmm0
  addsd  -0x1e0(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x1d8(%rbp)
  movsd  -0x288(%rbp),%xmm0
  movapd %xmm0,%xmm2
  mulsd  -0x2d8(%rbp),%xmm2
  movsd  -0x1e8(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  %xmm0,%xmm1
  movsd  -0x1e0(%rbp),%xmm0
  movapd %xmm0,%xmm3
  mulsd  %xmm0,%xmm3
  subsd  %xmm3,%xmm1
  movapd %xmm1,%xmm0
  mulsd  %xmm2,%xmm0
  movsd  %xmm0,-0x1d0(%rbp)
  movsd  -0x310(%rbp),%xmm1
  movsd  D4(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x288(%rbp),%xmm1
  movsd  D545(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x258(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  movapd %xmm0,%xmm2
  mulsd  -0x2c8(%rbp),%xmm2
  subsd  %xmm2,%xmm1
  movapd %xmm1,%xmm0
  movsd  %xmm0,-0x1c8(%rbp)
  mov    -0x2e0(%rbp),%rax
  movq   %rax,%xmm0
  call   <exp@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x1c0(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x378(%rbp)
  mov    -0x3a0(%rbp),%rax
  add    $0x60,%rax
  movsd  (%rax),%xmm1
  movsd  D37(%rip),%xmm0
  comisd %xmm1,%xmm0
  jbe    L142
  movsd  -0x2e0(%rbp),%xmm0
  movq   D538(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <exp@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x378(%rbp)
L142:
  movsd  -0x1c0(%rbp),%xmm0
  movapd %xmm0,%xmm1
  subsd  -0x378(%rbp),%xmm1
  movsd  D231(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x1b8(%rbp)
  movsd  -0x1c0(%rbp),%xmm0
  movapd %xmm0,%xmm1
  addsd  -0x378(%rbp),%xmm1
  movsd  D231(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x1b0(%rbp)
  movsd  -0x1b0(%rbp),%xmm0
  movsd  D131(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movsd  -0x2e0(%rbp),%xmm1
  mulsd  -0x1b8(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x1a8(%rbp)
  movsd  -0x2e0(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x1b0(%rbp),%xmm1
  movsd  -0x1b8(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x1a0(%rbp)
  movsd  -0x2f0(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x2e8(%rbp),%xmm1
  movsd  -0x250(%rbp),%xmm0
  movapd %xmm0,%xmm2
  addsd  -0x248(%rbp),%xmm2
  movsd  -0x2d0(%rbp),%xmm0
  mulsd  -0x240(%rbp),%xmm0
  addsd  %xmm2,%xmm0
  mulsd  %xmm1,%xmm0
  subsd  -0x258(%rbp),%xmm0
  movsd  -0x298(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x198(%rbp)
  cmpl   $0x0,-0x380(%rbp)
  jne    L98
  movsd  -0x240(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x1a0(%rbp),%xmm1
  movsd  -0x248(%rbp),%xmm0
  mulsd  -0x1b8(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  -0x250(%rbp),%xmm1
  mulsd  -0x1b0(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  mulsd  -0x2f0(%rbp),%xmm0
  movsd  -0x298(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x370(%rbp)
  movsd  -0x258(%rbp),%xmm0
  mulsd  -0x2c8(%rbp),%xmm0
  movsd  -0x2b0(%rbp),%xmm2
  movsd  D231(%rip),%xmm1
  mulsd  %xmm1,%xmm2
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  movsd  -0x240(%rbp),%xmm0
  mulsd  -0x1a8(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  -0x250(%rbp),%xmm1
  movapd %xmm1,%xmm2
  mulsd  -0x1b8(%rbp),%xmm2
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  movsd  -0x248(%rbp),%xmm0
  mulsd  -0x1b0(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  mulsd  -0x2f0(%rbp),%xmm0
  movsd  -0x298(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x368(%rbp)
L98:
  movsd  -0x310(%rbp),%xmm1
  movsd  D417(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  mulsd  -0x268(%rbp),%xmm0
  movsd  -0x2f8(%rbp),%xmm1
  movapd %xmm1,%xmm2
  mulsd  -0x2e8(%rbp),%xmm2
  movsd  -0x238(%rbp),%xmm1
  movapd %xmm1,%xmm3
  addsd  -0x230(%rbp),%xmm3
  movsd  -0x2d0(%rbp),%xmm1
  mulsd  -0x228(%rbp),%xmm1
  addsd  %xmm3,%xmm1
  mulsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  movsd  -0x298(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x190(%rbp)
  cmpl   $0x0,-0x380(%rbp)
  jne    L137
  movsd  -0x228(%rbp),%xmm0
  movq   D538(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  mulsd  -0x1a0(%rbp),%xmm0
  movsd  -0x230(%rbp),%xmm1
  movapd %xmm1,%xmm2
  mulsd  -0x1b8(%rbp),%xmm2
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  movsd  -0x238(%rbp),%xmm0
  mulsd  -0x1b0(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  mulsd  -0x2f8(%rbp),%xmm0
  movsd  -0x298(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x360(%rbp)
  movsd  -0x2b0(%rbp),%xmm1
  movsd  D231(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  mulsd  -0x270(%rbp),%xmm0
  movsd  -0x288(%rbp),%xmm1
  mulsd  -0x268(%rbp),%xmm1
  mulsd  -0x2c8(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  -0x228(%rbp),%xmm1
  movapd %xmm1,%xmm2
  mulsd  -0x1a8(%rbp),%xmm2
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  movsd  -0x238(%rbp),%xmm0
  mulsd  -0x1b8(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  -0x230(%rbp),%xmm1
  mulsd  -0x1b0(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  mulsd  -0x2f8(%rbp),%xmm0
  movsd  -0x298(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x358(%rbp)
L137:
  movsd  -0x310(%rbp),%xmm0
  movapd %xmm0,%xmm2
  addsd  %xmm0,%xmm2
  movsd  -0x308(%rbp),%xmm0
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  movsd  D517(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x260(%rbp),%xmm1
  movsd  D15(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x208(%rbp),%xmm1
  addsd  -0x200(%rbp),%xmm1
  mulsd  -0x2e8(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  -0x290(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x188(%rbp)
  cmpl   $0x0,-0x380(%rbp)
  jne    L194
  movsd  -0x208(%rbp),%xmm0
  mulsd  -0x1b0(%rbp),%xmm0
  movsd  -0x200(%rbp),%xmm1
  mulsd  -0x1b8(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  -0x290(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x350(%rbp)
  movsd  -0x308(%rbp),%xmm1
  movsd  D897(%rip),%xmm0
  movapd %xmm1,%xmm2
  mulsd  %xmm0,%xmm2
  movsd  -0x300(%rbp),%xmm0
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  movsd  D573(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x2b8(%rbp),%xmm1
  movsd  -0x300(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  movsd  -0x308(%rbp),%xmm3
  movsd  D636(%rip),%xmm2
  mulsd  %xmm2,%xmm3
  movapd %xmm0,%xmm2
  subsd  %xmm3,%xmm2
  movsd  -0x310(%rbp),%xmm3
  movsd  D2(%rip),%xmm0
  mulsd  %xmm3,%xmm0
  addsd  %xmm0,%xmm2
  movsd  D517(%rip),%xmm0
  addsd  %xmm2,%xmm0
  mulsd  -0x280(%rbp),%xmm0
  mulsd  -0x2c8(%rbp),%xmm0
  subsd  %xmm0,%xmm1
  movsd  -0x2f0(%rbp),%xmm0
  mulsd  -0x210(%rbp),%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x208(%rbp),%xmm0
  mulsd  -0x1b8(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  -0x200(%rbp),%xmm1
  mulsd  -0x1b0(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  -0x290(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x348(%rbp)
L194:
  movsd  -0x1f8(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x180(%rbp)
  movsd  -0x1f0(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x178(%rbp)
  movsd  -0x180(%rbp),%xmm0
  addsd  -0x178(%rbp),%xmm0
  movsd  %xmm0,-0x170(%rbp)
  movsd  -0x1f8(%rbp),%xmm0
  mulsd  -0x1f0(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  movsd  %xmm0,-0x168(%rbp)
  movsd  -0x180(%rbp),%xmm0
  addsd  -0x178(%rbp),%xmm0
  subsd  -0x168(%rbp),%xmm0
  mulsd  -0x2e8(%rbp),%xmm0
  movsd  -0x270(%rbp),%xmm2
  movsd  D231(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  movsd  -0x2a0(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x160(%rbp)
  cmpl   $0x0,-0x380(%rbp)
  jne    L26
  movsd  -0x170(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x1b8(%rbp),%xmm1
  movsd  -0x168(%rbp),%xmm0
  mulsd  -0x1b0(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  -0x2a0(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x340(%rbp)
  movsd  -0x170(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x1b0(%rbp),%xmm1
  movsd  -0x168(%rbp),%xmm0
  mulsd  -0x1b8(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  addsd  -0x180(%rbp),%xmm0
  subsd  -0x178(%rbp),%xmm0
  movsd  -0x2a0(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x338(%rbp)
L26:
  movsd  -0x1d0(%rbp),%xmm0
  mulsd  -0x198(%rbp),%xmm0
  movapd %xmm0,%xmm1
  addsd  -0x1d8(%rbp),%xmm1
  movsd  -0x1c8(%rbp),%xmm0
  mulsd  -0x190(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x158(%rbp)
  cmpl   $0x0,-0x380(%rbp)
  je     L40
  movsd  -0x308(%rbp),%xmm1
  movsd  D897(%rip),%xmm0
  movapd %xmm1,%xmm2
  mulsd  %xmm0,%xmm2
  movsd  -0x300(%rbp),%xmm0
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  movsd  D573(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x2b8(%rbp),%xmm1
  movsd  -0x300(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  movsd  -0x308(%rbp),%xmm3
  movsd  D636(%rip),%xmm2
  mulsd  %xmm2,%xmm3
  movapd %xmm0,%xmm2
  subsd  %xmm3,%xmm2
  movsd  -0x310(%rbp),%xmm3
  movsd  D2(%rip),%xmm0
  mulsd  %xmm3,%xmm0
  addsd  %xmm0,%xmm2
  movsd  D517(%rip),%xmm0
  addsd  %xmm2,%xmm0
  mulsd  -0x280(%rbp),%xmm0
  mulsd  -0x2c8(%rbp),%xmm0
  subsd  %xmm0,%xmm1
  movsd  -0x2f0(%rbp),%xmm0
  mulsd  -0x210(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  -0x290(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x348(%rbp)
  movsd  -0x1d0(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x368(%rbp),%xmm1
  movsd  -0x1c8(%rbp),%xmm0
  mulsd  -0x358(%rbp),%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x2f8(%rbp),%xmm0
  mulsd  -0x1e0(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  -0x310(%rbp),%xmm2
  movsd  D626(%rip),%xmm1
  mulsd  %xmm1,%xmm2
  movsd  -0x310(%rbp),%xmm3
  movsd  D517(%rip),%xmm1
  addsd  %xmm3,%xmm1
  mulsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  movsd  D131(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x148(%rbp)
  movsd  -0x180(%rbp),%xmm0
  subsd  -0x178(%rbp),%xmm0
  movsd  -0x2a0(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x338(%rbp)
  movsd  -0x2b0(%rbp),%xmm1
  movsd  D231(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  mulsd  -0x270(%rbp),%xmm0
  movsd  -0x288(%rbp),%xmm1
  mulsd  -0x268(%rbp),%xmm1
  mulsd  -0x2c8(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  addsd  -0x228(%rbp),%xmm0
  mulsd  -0x2f8(%rbp),%xmm0
  movsd  -0x298(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x358(%rbp)
  movsd  -0x258(%rbp),%xmm0
  mulsd  -0x2c8(%rbp),%xmm0
  movsd  -0x2b0(%rbp),%xmm2
  movsd  D231(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  subsd  -0x240(%rbp),%xmm0
  mulsd  -0x2f0(%rbp),%xmm0
  movsd  -0x298(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x368(%rbp)
  movsd  -0x148(%rbp),%xmm0
  mulsd  -0x368(%rbp),%xmm0
  movsd  -0x2f8(%rbp),%xmm1
  mulsd  -0x348(%rbp),%xmm1
  mulsd  -0x348(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x118(%rbp)
  movsd  -0x158(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x368(%rbp),%xmm1
  movsd  -0x198(%rbp),%xmm0
  mulsd  -0x148(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  -0x2f0(%rbp),%xmm1
  mulsd  -0x188(%rbp),%xmm1
  mulsd  -0x348(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x120(%rbp)
  movsd  -0x158(%rbp),%xmm0
  mulsd  -0x198(%rbp),%xmm0
  movsd  -0x2f8(%rbp),%xmm1
  mulsd  -0x188(%rbp),%xmm1
  mulsd  -0x188(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x128(%rbp)
  movsd  -0x2f8(%rbp),%xmm0
  mulsd  -0x338(%rbp),%xmm0
  movsd  -0x338(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x100(%rbp)
  movsd  -0x2f0(%rbp),%xmm0
  mulsd  -0x160(%rbp),%xmm0
  mulsd  -0x338(%rbp),%xmm0
  subsd  -0x358(%rbp),%xmm0
  subsd  -0x358(%rbp),%xmm0
  subsd  -0x2d8(%rbp),%xmm0
  movsd  %xmm0,-0x108(%rbp)
  movsd  -0x2f8(%rbp),%xmm0
  mulsd  -0x160(%rbp),%xmm0
  mulsd  -0x160(%rbp),%xmm0
  subsd  -0x190(%rbp),%xmm0
  subsd  -0x190(%rbp),%xmm0
  movsd  %xmm0,-0x110(%rbp)
  movsd  -0x120(%rbp),%xmm0
  mulsd  -0x100(%rbp),%xmm0
  movsd  -0x108(%rbp),%xmm1
  mulsd  -0x118(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x90(%rbp)
  movsd  -0x128(%rbp),%xmm0
  mulsd  -0x100(%rbp),%xmm0
  movsd  -0x110(%rbp),%xmm1
  mulsd  -0x118(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x98(%rbp)
  movsd  -0x98(%rbp),%xmm0
  movq   D538(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  divsd  -0x90(%rbp),%xmm0
  movsd  %xmm0,-0x330(%rbp)
  movsd  -0x330(%rbp),%xmm0
  movq   D538(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,-0x320(%rbp)
  mov    -0x3a0(%rbp),%rax
  add    $0x18,%rax
  movsd  -0x330(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x3a0(%rbp),%rax
  add    $0x10,%rax
  movsd  -0x320(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x358(%rbp),%xmm0
  mulsd  -0x330(%rbp),%xmm0
  mov    -0x3a0(%rbp),%rax
  add    $0x8,%rax
  addsd  -0x190(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x368(%rbp),%xmm0
  mulsd  -0x330(%rbp),%xmm0
  addsd  -0x198(%rbp),%xmm0
  mov    -0x3a0(%rbp),%rax
  movsd  %xmm0,(%rax)
  movsd  -0x348(%rbp),%xmm0
  mulsd  -0x330(%rbp),%xmm0
  mov    -0x3a0(%rbp),%rax
  add    $0x40,%rax
  addsd  -0x188(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x338(%rbp),%xmm0
  mulsd  -0x330(%rbp),%xmm0
  addsd  -0x160(%rbp),%xmm0
  mov    -0x3a0(%rbp),%rax
  add    $0x70,%rax
  movq   D538(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x3a0(%rbp),%rax
  lea    0x78(%rax),%rdx
  mov    -0x3a0(%rbp),%rax
  movsd  0x70(%rax),%xmm0
  movsd  %xmm0,(%rdx)
  mov    -0x3a0(%rbp),%rax
  add    $0x78,%rax
  movsd  (%rax),%xmm0
  movq   D502(%rip),%xmm1
  andpd  %xmm0,%xmm1
  movsd  D516(%rip),%xmm0
  comisd %xmm1,%xmm0
  jbe    L16
  mov    -0x3a0(%rbp),%rax
  add    $0x78,%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
L16:
  mov    -0x3a0(%rbp),%rax
  lea    0x50(%rax),%rdx
  mov    -0x3a0(%rbp),%rax
  movsd  0x80(%rax),%xmm0
  movsd  %xmm0,(%rdx)
  jmp    L62
L40:
  movsd  -0x1d0(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x370(%rbp),%xmm1
  movsd  -0x1c8(%rbp),%xmm0
  mulsd  -0x360(%rbp),%xmm0
  movapd %xmm1,%xmm2
  addsd  %xmm0,%xmm2
  movsd  -0x1e8(%rbp),%xmm0
  mulsd  -0x1b0(%rbp),%xmm0
  movsd  -0x1e0(%rbp),%xmm1
  mulsd  -0x1b8(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  mulsd  -0x2f8(%rbp),%xmm0
  addsd  %xmm2,%xmm0
  movsd  %xmm0,-0x150(%rbp)
  movsd  -0x1d0(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x368(%rbp),%xmm1
  movsd  -0x1c8(%rbp),%xmm0
  mulsd  -0x358(%rbp),%xmm0
  movapd %xmm1,%xmm2
  addsd  %xmm0,%xmm2
  movsd  -0x1e8(%rbp),%xmm0
  mulsd  -0x1b8(%rbp),%xmm0
  movsd  -0x1b0(%rbp),%xmm1
  movsd  D131(%rip),%xmm3
  subsd  %xmm3,%xmm1
  mulsd  -0x1e0(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  mulsd  -0x2f8(%rbp),%xmm0
  addsd  %xmm2,%xmm0
  movsd  -0x310(%rbp),%xmm2
  movsd  D626(%rip),%xmm1
  mulsd  %xmm1,%xmm2
  movsd  -0x310(%rbp),%xmm3
  movsd  D517(%rip),%xmm1
  addsd  %xmm3,%xmm1
  mulsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  movsd  D131(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x148(%rbp)
  movsd  -0x150(%rbp),%xmm0
  mulsd  -0x370(%rbp),%xmm0
  movsd  -0x2f8(%rbp),%xmm1
  mulsd  -0x350(%rbp),%xmm1
  mulsd  -0x350(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x140(%rbp)
  movsd  -0x158(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x370(%rbp),%xmm1
  movsd  -0x150(%rbp),%xmm0
  mulsd  -0x198(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  -0x2f0(%rbp),%xmm1
  mulsd  -0x188(%rbp),%xmm1
  mulsd  -0x350(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x138(%rbp)
  movsd  -0x150(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x368(%rbp),%xmm1
  movsd  -0x148(%rbp),%xmm0
  mulsd  -0x370(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  -0x2f0(%rbp),%xmm1
  mulsd  -0x350(%rbp),%xmm1
  mulsd  -0x348(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x130(%rbp)
  movsd  -0x158(%rbp),%xmm0
  mulsd  -0x198(%rbp),%xmm0
  movsd  -0x2f8(%rbp),%xmm1
  mulsd  -0x188(%rbp),%xmm1
  mulsd  -0x188(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x128(%rbp)
  movsd  -0x158(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x368(%rbp),%xmm1
  movsd  -0x148(%rbp),%xmm0
  mulsd  -0x198(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  -0x2f0(%rbp),%xmm1
  mulsd  -0x188(%rbp),%xmm1
  mulsd  -0x348(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x120(%rbp)
  movsd  -0x148(%rbp),%xmm0
  mulsd  -0x368(%rbp),%xmm0
  movsd  -0x2f8(%rbp),%xmm1
  mulsd  -0x348(%rbp),%xmm1
  mulsd  -0x348(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x118(%rbp)
  cmpl   $0x1,-0x394(%rbp)
  je     L214
  cmpl   $0x3,-0x394(%rbp)
  jne    L229
L214:
  movsd  -0x2f8(%rbp),%xmm0
  mulsd  -0x340(%rbp),%xmm0
  movsd  -0x340(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x250(%rbp)
  movsd  -0x2f0(%rbp),%xmm0
  mulsd  -0x160(%rbp),%xmm0
  mulsd  -0x340(%rbp),%xmm0
  subsd  -0x360(%rbp),%xmm0
  subsd  -0x360(%rbp),%xmm0
  movsd  %xmm0,-0x248(%rbp)
  movsd  -0x2f0(%rbp),%xmm0
  mulsd  -0x340(%rbp),%xmm0
  movsd  -0x338(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x240(%rbp)
  movsd  -0x2f8(%rbp),%xmm0
  mulsd  -0x160(%rbp),%xmm0
  mulsd  -0x160(%rbp),%xmm0
  subsd  -0x190(%rbp),%xmm0
  subsd  -0x190(%rbp),%xmm0
  movsd  %xmm0,-0x110(%rbp)
  movsd  -0x2f0(%rbp),%xmm0
  mulsd  -0x160(%rbp),%xmm0
  mulsd  -0x338(%rbp),%xmm0
  subsd  -0x358(%rbp),%xmm0
  subsd  -0x358(%rbp),%xmm0
  subsd  -0x2d8(%rbp),%xmm0
  movsd  %xmm0,-0x108(%rbp)
  movsd  -0x2f8(%rbp),%xmm0
  mulsd  -0x338(%rbp),%xmm0
  movsd  -0x338(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x100(%rbp)
  movsd  -0x140(%rbp),%xmm0
  mulsd  -0x100(%rbp),%xmm0
  movsd  -0x250(%rbp),%xmm1
  mulsd  -0x118(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0xf8(%rbp)
  movsd  -0x140(%rbp),%xmm0
  mulsd  -0x108(%rbp),%xmm0
  movsd  -0x250(%rbp),%xmm1
  mulsd  -0x120(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0xf0(%rbp)
  movsd  -0x140(%rbp),%xmm0
  mulsd  -0x110(%rbp),%xmm0
  movsd  -0x250(%rbp),%xmm1
  mulsd  -0x128(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0xe8(%rbp)
  movsd  -0x140(%rbp),%xmm0
  mulsd  -0x240(%rbp),%xmm0
  movsd  -0x250(%rbp),%xmm1
  mulsd  -0x130(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0xe0(%rbp)
  movsd  -0x140(%rbp),%xmm0
  mulsd  -0x248(%rbp),%xmm0
  movsd  -0x250(%rbp),%xmm1
  mulsd  -0x138(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0xd8(%rbp)
  movsd  -0x138(%rbp),%xmm0
  mulsd  -0x100(%rbp),%xmm0
  movsd  -0x248(%rbp),%xmm1
  mulsd  -0x118(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0xd0(%rbp)
  movsd  -0x138(%rbp),%xmm0
  mulsd  -0x108(%rbp),%xmm0
  movsd  -0x248(%rbp),%xmm1
  mulsd  -0x120(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0xc8(%rbp)
  movsd  -0x138(%rbp),%xmm0
  mulsd  -0x110(%rbp),%xmm0
  movsd  -0x248(%rbp),%xmm1
  mulsd  -0x128(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0xc0(%rbp)
  movsd  -0x130(%rbp),%xmm0
  mulsd  -0x100(%rbp),%xmm0
  movsd  -0x240(%rbp),%xmm1
  mulsd  -0x118(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0xb8(%rbp)
  movsd  -0x130(%rbp),%xmm0
  mulsd  -0x108(%rbp),%xmm0
  movsd  -0x240(%rbp),%xmm1
  mulsd  -0x120(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0xb0(%rbp)
  movsd  -0x130(%rbp),%xmm0
  mulsd  -0x110(%rbp),%xmm0
  movsd  -0x240(%rbp),%xmm1
  mulsd  -0x128(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0xa8(%rbp)
  movsd  -0x130(%rbp),%xmm0
  mulsd  -0x248(%rbp),%xmm0
  movsd  -0x240(%rbp),%xmm1
  mulsd  -0x138(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0xa0(%rbp)
  movsd  -0x128(%rbp),%xmm0
  mulsd  -0x100(%rbp),%xmm0
  movsd  -0x110(%rbp),%xmm1
  mulsd  -0x118(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x98(%rbp)
  movsd  -0x120(%rbp),%xmm0
  mulsd  -0x100(%rbp),%xmm0
  movsd  -0x108(%rbp),%xmm1
  mulsd  -0x118(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x90(%rbp)
  movsd  -0xb0(%rbp),%xmm0
  addsd  -0xd0(%rbp),%xmm0
  movsd  %xmm0,-0x88(%rbp)
  movsd  -0xa8(%rbp),%xmm0
  addsd  -0xc8(%rbp),%xmm0
  movsd  %xmm0,-0x80(%rbp)
  movsd  -0xf8(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  -0xe0(%rbp),%xmm1
  mulsd  -0xb8(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x78(%rbp)
  movsd  -0xf8(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  mulsd  -0xf0(%rbp),%xmm0
  movsd  -0xe0(%rbp),%xmm1
  mulsd  -0x88(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  -0xd8(%rbp),%xmm1
  mulsd  -0xb8(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x70(%rbp)
  movsd  -0xf0(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  %xmm0,%xmm1
  movsd  -0xf8(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  mulsd  -0xe8(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  -0xe0(%rbp),%xmm1
  mulsd  -0x80(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  -0xd8(%rbp),%xmm1
  mulsd  -0x88(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x68(%rbp)
  movsd  -0xf0(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  mulsd  -0xe8(%rbp),%xmm0
  movsd  -0xe0(%rbp),%xmm1
  mulsd  -0xc0(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  -0xd8(%rbp),%xmm1
  mulsd  -0x80(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x60(%rbp)
  movsd  -0xe8(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  -0xd8(%rbp),%xmm1
  mulsd  -0xc0(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x58(%rbp)
  cmpl   $0x1,-0x394(%rbp)
  jne    L76
  movsd  -0xe8(%rbp),%xmm0
  subsd  -0xa8(%rbp),%xmm0
  subsd  -0x98(%rbp),%xmm0
  movsd  -0xd8(%rbp),%xmm1
  subsd  -0xf0(%rbp),%xmm1
  addsd  -0xb0(%rbp),%xmm1
  subsd  -0xd0(%rbp),%xmm1
  addsd  -0x90(%rbp),%xmm1
  subsd  -0xa0(%rbp),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x328(%rbp)
  jmp    L66
L76:
  movsd  -0x270(%rbp),%xmm1
  movsd  D231(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  mulsd  -0x2a0(%rbp),%xmm0
  movsd  %xmm0,-0x3a8(%rbp)
  movsd  -0x2e8(%rbp),%xmm0
  movq   D538(%rip),%xmm1
  movapd %xmm0,%xmm5
  xorpd  %xmm1,%xmm5
  movq   %xmm5,%rax
  movq   %rax,%xmm0
  call   <exp@plt>
  mov    -0x3a0(%rbp),%rax
  add    $0x70,%rax
  mulsd  -0x3a8(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x3a0(%rbp),%rax
  add    $0x58,%rax
  movsd  (%rax),%xmm1
  movsd  D336(%rip),%xmm0
  comisd %xmm1,%xmm0
  jb     L20
  mov    -0x3a0(%rbp),%rax
  add    $0x58,%rax
  movsd  (%rax),%xmm0
  pxor   %xmm1,%xmm1
  comisd %xmm1,%xmm0
  jb     L20
  mov    -0x3a0(%rbp),%rax
  add    $0x60,%rax
  movsd  (%rax),%xmm1
  movsd  D131(%rip),%xmm0
  comisd %xmm1,%xmm0
  jb     L20
  movsd  -0x2f0(%rbp),%xmm0
  mulsd  -0x2e8(%rbp),%xmm0
  movsd  %xmm0,-0x3a8(%rbp)
  mov    -0x2e0(%rbp),%rax
  movq   %rax,%xmm0
  call   <exp@plt>
  mulsd  -0x3a8(%rbp),%xmm0
  movsd  %xmm0,-0x50(%rbp)
  mov    -0x50(%rbp),%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x48(%rbp)
  movsd  -0x2b0(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  mulsd  -0x2b8(%rbp),%xmm0
  mulsd  -0x48(%rbp),%xmm0
  movapd %xmm0,%xmm1
  divsd  -0x258(%rbp),%xmm1
  movsd  D131(%rip),%xmm0
  addsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movapd %xmm0,%xmm1
  movsd  D131(%rip),%xmm0
  subsd  %xmm1,%xmm0
  mulsd  -0x260(%rbp),%xmm0
  divsd  -0x2b8(%rbp),%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movsd  -0x40(%rbp),%xmm0
  movq   D538(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  mulsd  -0x40(%rbp),%xmm0
  movapd %xmm0,%xmm1
  divsd  -0x2f0(%rbp),%xmm1
  movsd  -0x270(%rbp),%xmm2
  movsd  D231(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  mulsd  -0x2a0(%rbp),%xmm0
  mov    -0x3a0(%rbp),%rax
  add    $0x70,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
L20:
  mov    -0x3a0(%rbp),%rax
  add    $0x70,%rax
  movsd  (%rax),%xmm0
  movsd  -0x160(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x2d8(%rbp),%xmm0
  movapd %xmm0,%xmm2
  mulsd  -0x38(%rbp),%xmm2
  movsd  -0x358(%rbp),%xmm0
  mulsd  -0x38(%rbp),%xmm0
  movsd  -0x190(%rbp),%xmm1
  mulsd  -0x338(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  addsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm2
  movapd %xmm2,%xmm1
  mov    -0x3a0(%rbp),%rax
  add    $0x70,%rax
  movsd  (%rax),%xmm0
  movapd %xmm0,%xmm2
  mulsd  -0x2f8(%rbp),%xmm2
  mov    -0x3a0(%rbp),%rax
  add    $0x70,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm2,%xmm0
  mulsd  -0x338(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x320(%rbp)
  movsd  -0x320(%rbp),%xmm0
  movq   D538(%rip),%xmm1
  xorpd  %xmm0,%xmm1
  movsd  -0x2d8(%rbp),%xmm0
  movapd %xmm0,%xmm3
  mulsd  -0x340(%rbp),%xmm3
  movsd  -0x358(%rbp),%xmm0
  mulsd  -0x340(%rbp),%xmm0
  movsd  -0x360(%rbp),%xmm2
  mulsd  -0x338(%rbp),%xmm2
  subsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm3
  movapd %xmm3,%xmm2
  divsd  %xmm2,%xmm1
  movapd %xmm1,%xmm0
  movsd  %xmm0,-0x320(%rbp)
  movsd  -0x340(%rbp),%xmm0
  mulsd  -0x320(%rbp),%xmm0
  addsd  -0x38(%rbp),%xmm0
  movq   D538(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  divsd  -0x338(%rbp),%xmm0
  movsd  %xmm0,-0x328(%rbp)
L66:
  movl   $0x0,-0x384(%rbp)
L82:
  addl   $0x1,-0x384(%rbp)
  mov    -0x384(%rbp),%eax
  cmp    -0x37c(%rbp),%eax
  jle    L151
  movl   $0xfffffffe,-0x398(%rbp)
  mov    -0x398(%rbp),%eax
  jmp    L56
L151:
  movsd  -0x328(%rbp),%xmm0
  movsd  %xmm0,-0x330(%rbp)
  movsd  -0x78(%rbp),%xmm0
  mulsd  -0x330(%rbp),%xmm0
  addsd  -0x70(%rbp),%xmm0
  mulsd  -0x330(%rbp),%xmm0
  addsd  -0x68(%rbp),%xmm0
  mulsd  -0x330(%rbp),%xmm0
  addsd  -0x60(%rbp),%xmm0
  mulsd  -0x330(%rbp),%xmm0
  movsd  -0x58(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  -0x68(%rbp),%xmm0
  movapd %xmm0,%xmm1
  addsd  %xmm0,%xmm1
  movsd  -0x70(%rbp),%xmm2
  movsd  D636(%rip),%xmm0
  mulsd  %xmm0,%xmm2
  movsd  -0x78(%rbp),%xmm3
  movsd  D209(%rip),%xmm0
  mulsd  %xmm3,%xmm0
  mulsd  -0x330(%rbp),%xmm0
  addsd  %xmm2,%xmm0
  mulsd  -0x330(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  mulsd  -0x330(%rbp),%xmm0
  movsd  -0x60(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x30(%rbp),%xmm0
  movapd %xmm0,%xmm1
  divsd  -0x28(%rbp),%xmm1
  movsd  -0x330(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x328(%rbp)
  movsd  -0x328(%rbp),%xmm0
  subsd  -0x330(%rbp),%xmm0
  divsd  -0x330(%rbp),%xmm0
  movq   D502(%rip),%xmm1
  andpd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x20(%rbp),%xmm0
  comisd -0x318(%rbp),%xmm0
  ja     L82
  mov    -0x384(%rbp),%eax
  add    %eax,-0x398(%rbp)
  movsd  -0x328(%rbp),%xmm0
  movsd  %xmm0,-0x330(%rbp)
  movsd  -0xf8(%rbp),%xmm0
  mulsd  -0x330(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x330(%rbp),%xmm1
  movsd  -0xf0(%rbp),%xmm0
  mulsd  -0x330(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  addsd  -0xe8(%rbp),%xmm0
  movq   D538(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  -0xe0(%rbp),%xmm1
  mulsd  -0x330(%rbp),%xmm1
  addsd  -0xd8(%rbp),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x320(%rbp)
  movsd  -0x340(%rbp),%xmm0
  mulsd  -0x320(%rbp),%xmm0
  movapd %xmm0,%xmm1
  addsd  -0x160(%rbp),%xmm1
  movsd  -0x338(%rbp),%xmm0
  mulsd  -0x330(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  mov    -0x3a0(%rbp),%rax
  add    $0x70,%rax
  movq   D538(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x3a0(%rbp),%rax
  lea    0x78(%rax),%rdx
  mov    -0x3a0(%rbp),%rax
  movsd  0x70(%rax),%xmm0
  movsd  %xmm0,(%rdx)
  mov    -0x3a0(%rbp),%rax
  add    $0x78,%rax
  movsd  (%rax),%xmm0
  movq   D502(%rip),%xmm1
  andpd  %xmm0,%xmm1
  movsd  D516(%rip),%xmm0
  comisd %xmm1,%xmm0
  jbe    L136
  mov    -0x3a0(%rbp),%rax
  add    $0x78,%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
L136:
  mov    -0x3a0(%rbp),%rax
  lea    0x50(%rax),%rdx
  mov    -0x3a0(%rbp),%rax
  movsd  0x80(%rax),%xmm0
  movsd  %xmm0,(%rdx)
  jmp    L179
L229:
  movsd  -0x2d8(%rbp),%xmm0
  movapd %xmm0,%xmm2
  mulsd  -0x160(%rbp),%xmm2
  movsd  -0x358(%rbp),%xmm0
  mulsd  -0x160(%rbp),%xmm0
  movsd  -0x190(%rbp),%xmm1
  mulsd  -0x338(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  addsd  %xmm0,%xmm0
  addsd  %xmm2,%xmm0
  movsd  %xmm0,-0x320(%rbp)
  movsd  -0x320(%rbp),%xmm0
  movq   D538(%rip),%xmm1
  xorpd  %xmm0,%xmm1
  movsd  -0x2d8(%rbp),%xmm0
  movapd %xmm0,%xmm3
  mulsd  -0x340(%rbp),%xmm3
  movsd  -0x358(%rbp),%xmm0
  mulsd  -0x340(%rbp),%xmm0
  movsd  -0x360(%rbp),%xmm2
  mulsd  -0x338(%rbp),%xmm2
  subsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm3
  movapd %xmm3,%xmm2
  divsd  %xmm2,%xmm1
  movapd %xmm1,%xmm0
  movsd  %xmm0,-0x320(%rbp)
  movsd  -0x340(%rbp),%xmm0
  mulsd  -0x320(%rbp),%xmm0
  addsd  -0x160(%rbp),%xmm0
  movq   D538(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  divsd  -0x338(%rbp),%xmm0
  movsd  %xmm0,-0x330(%rbp)
  cmpl   $0x2,-0x394(%rbp)
  jne    L161
  movsd  -0x140(%rbp),%xmm0
  mulsd  -0x320(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x320(%rbp),%xmm1
  movsd  -0x130(%rbp),%xmm0
  mulsd  -0x330(%rbp),%xmm0
  addsd  -0x138(%rbp),%xmm0
  mulsd  -0x320(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  addsd  -0x128(%rbp),%xmm1
  movsd  -0x120(%rbp),%xmm0
  mulsd  -0x330(%rbp),%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x118(%rbp),%xmm0
  mulsd  -0x330(%rbp),%xmm0
  mulsd  -0x330(%rbp),%xmm0
  mov    -0x3a0(%rbp),%rax
  add    $0x78,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
L161:
  cmpl   $0x4,-0x394(%rbp)
  jne    L179
  movsd  -0x340(%rbp),%xmm0
  mulsd  -0x320(%rbp),%xmm0
  movapd %xmm0,%xmm1
  addsd  -0x160(%rbp),%xmm1
  movsd  -0x338(%rbp),%xmm0
  mulsd  -0x330(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  mov    -0x3a0(%rbp),%rax
  add    $0x78,%rax
  movq   D538(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
L179:
  mov    -0x3a0(%rbp),%rax
  add    $0x18,%rax
  movsd  -0x330(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x3a0(%rbp),%rax
  add    $0x10,%rax
  movsd  -0x320(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x360(%rbp),%xmm0
  mulsd  -0x320(%rbp),%xmm0
  movapd %xmm0,%xmm1
  addsd  -0x190(%rbp),%xmm1
  movsd  -0x358(%rbp),%xmm0
  mulsd  -0x330(%rbp),%xmm0
  mov    -0x3a0(%rbp),%rax
  add    $0x8,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x370(%rbp),%xmm0
  mulsd  -0x320(%rbp),%xmm0
  movapd %xmm0,%xmm1
  addsd  -0x198(%rbp),%xmm1
  movsd  -0x368(%rbp),%xmm0
  mulsd  -0x330(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  mov    -0x3a0(%rbp),%rax
  movsd  %xmm0,(%rax)
  movsd  -0x350(%rbp),%xmm0
  mulsd  -0x320(%rbp),%xmm0
  movapd %xmm0,%xmm1
  addsd  -0x188(%rbp),%xmm1
  movsd  -0x348(%rbp),%xmm0
  mulsd  -0x330(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  mov    -0x3a0(%rbp),%rax
  movsd  (%rax),%xmm1
  mov    -0x3a0(%rbp),%rax
  add    $0x40,%rax
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
L62:
  movsd  -0x2f0(%rbp),%xmm0
  mulsd  -0x2e8(%rbp),%xmm0
  movsd  -0x1c0(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x2e0(%rbp),%xmm0
  mulsd  -0x2d8(%rbp),%xmm0
  mulsd  -0x320(%rbp),%xmm0
  subsd  -0x18(%rbp),%xmm0
  movsd  -0x2d8(%rbp),%xmm1
  mulsd  -0x18(%rbp),%xmm1
  mov    -0x3a0(%rbp),%rax
  add    $0x38,%rax
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x398(%rbp),%eax
L56:
  mov    -0x8(%rbp),%rbx
  leave
  ret
```

tool_calls: name: disassemble_data
arguments: data_label: D356
type: byte[8]
name: disassemble_data
arguments: data_label: D170
type: f64
name: disassemble_data
arguments: data_label: D890
type: byte[8]
name: disassemble_data
arguments: data_label: D88
type: f64[2]
name: disassemble_data
arguments: data_label: D131
type: f64
name: disassemble_data
arguments: data_label: D65
type: byte[8]
name: disassemble_data
arguments: data_label: D636
type: f64
name: disassemble_data
arguments: data_label: D336
type: f64
name: disassemble_data
arguments: data_label: D538
type: f32[2]
name: disassemble_data
arguments: data_label: D147
type: f64
name: disassemble_data
arguments: data_label: D209
type: f64
name: disassemble_data
arguments: data_label: D648
type: f64
name: disassemble_data
arguments: data_label: D140
type: f64
name: disassemble_data
arguments: data_label: D108
type: f64
name: disassemble_data
arguments: data_label: D4
type: f64
name: disassemble_data
arguments: data_label: D667
type: f64
name: disassemble_data
arguments: data_label: D573
type: f64
name: disassemble_data
arguments: data_label: D825
type: f64
name: disassemble_data
arguments: data_label: D37
type: f64
name: disassemble_data
arguments: data_label: D480
type: f64
name: disassemble_data
arguments: data_label: D545
type: byte[8]
name: disassemble_data
arguments: data_label: D231
type: f64
name: disassemble_data
arguments: data_label: D417
type: f64
name: disassemble_data
arguments: data_label: D517
type: f64
name: disassemble_data
arguments: data_label: D15
type: f32[2]
name: disassemble_data
arguments: data_label: D897
type: f64
name: disassemble_data
arguments: data_label: D2
type: f64
name: disassemble_data
arguments: data_label: D626
type: byte[8]
name: disassemble_data
arguments: data_label: D502
type: qword
name: disassemble_data
arguments: data_label: D516
type: byte[8]

tool: {"D356": [0x8d,0xed,0xb5,0xa0,0xf7,0xc6,0xb0,0x3e]}

tool: D170: 12

tool: {"D890": [0x55,0x55,0x55,0x55,0x55,0x55,0xd5,0x3f]}

tool: {"D88": [15.0,1.0]}

tool: D131: 1

tool: {"D65": [0x9a,0x99,0x99,0x99,0x99,0x99,0xc9,0x3f]}

tool: D636: 3

tool: D336: 2

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

tool: D147: 14

tool: D209: 4

tool: D648: 36

tool: D140: 7

tool: D108: 9

tool: D4: 8

tool: D667: 45

tool: D573: 5

tool: D825: 42

tool: D37: 20

tool: D480: 30

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

tool: D231: 0.5

tool: D417: 1.5

tool: D517: 10

tool: {"D15": [0.0,1.625]}

tool: D897: 6

tool: D2: 18

tool: {"D626": [0x9a,0x99,0x99,0x99,0x99,0x99,0xd9,0x3f]}

tool: D502: 9223372036854775807

tool: {"D516": [0xfc,0xa9,0xf1,0xd2,0x4d,0x62,0x50,0x3f]}

assistant: ```c
int sqfun(int ix, int ir, double gMSAWave[]) {
  double acc = 1.0e-6;
  double reta, eta2, eta21, eta22, eta3, eta32, eta2d, eta2d2, eta3d, eta6d, e12, e24, rgek;
  double rak, ak1, ak2, dak, dak2, dak4, d, d2, dd2, dd4, dd45, ex1, ex2, sk, ck, ckma, skma;
  double al1, al2, al3, al4, al5, al6, be1, be2, be3, vu1, vu2, vu3, vu4, vu5, ph1, ph2, ta1, ta2, ta3, ta4, ta5;
  double a1, a2, a3, b1, b2, b3, v1, v2, v3, p1, p2, p3, pp, pp1, pp2, p1p2, t1, t2, t3, um1, um2, um3, um4, um5, um6;
  double w0, w1, w2, w3, w4, w12, w13, w14, w15, w16, w24, w25, w26, w32, w34, w3425, w35, w3526, w36, w46, w56;
  double fa, fap, ca, e24g, pwk, qpw, pg, del, fun, fund, g24;
  int ii, ibig, itm = 40;
  a2 = 0;
  a3 = 0;
  b2 = 0;
  b3 = 0;
  v2 = 0;
  v3 = 0;
  p2 = 0;
  p3 = 0;
  reta = gMSAWave[16];
  eta2 = reta * reta;
  eta3 = eta2 * reta;
  e12 = 12.0 * reta;
  e24 = e12 + e12;
  gMSAWave[13] = pow((gMSAWave[4] / gMSAWave[16]), (1.0 / 3.0));
  gMSAWave[12] = gMSAWave[6] / gMSAWave[13];
  ibig = 0;
  if ((gMSAWave[12] > 15.0) && (ix == 1)) {
    ibig = 1;
  }
  gMSAWave[11] = gMSAWave[5] * gMSAWave[13] * exp(gMSAWave[6] - gMSAWave[12]);
  rgek = gMSAWave[11];
  rak = gMSAWave[12];
  ak2 = rak * rak;
  ak1 = 1.0 + rak;
  dak2 = 1.0 / ak2;
  dak4 = dak2 * dak2;
  d = 1.0 - reta;
  d2 = d * d;
  dak = d / rak;
  dd2 = 1.0 / d2;
  dd4 = dd2 * dd2;
  dd45 = dd4 * 2.0e-1;
  eta3d = 3.0 * reta;
  eta6d = eta3d + eta3d;
  eta32 = eta3 + eta3;
  eta2d = reta + 2.0;
  eta2d2 = eta2d * eta2d;
  eta21 = 2.0 * reta + 1.0;
  eta22 = eta21 * eta21;
  al1 = -eta21 * dak;
  al2 = (14.0 * eta2 - 4.0 * reta - 1.0) * dak2;
  al3 = 36.0 * eta2 * dak4;
  be1 = -(eta2 + 7.0 * reta + 1.0) * dak;
  be2 = 9.0 * reta * (eta2 + 4.0 * reta - 2.0) * dak2;
  be3 = 12.0 * reta * (2.0 * eta2 + 8.0 * reta - 1.0) * dak4;
  vu1 = -(eta3 + 3.0 * eta2 + 45.0 * reta + 5.0) * dak;
  vu2 = (eta32 + 3.0 * eta2 + 42.0 * reta - 2.0e1) * dak2;
  vu3 = (eta32 + 3.0e1 * reta - 5.0) * dak4;
  vu4 = vu1 + e24 * rak * vu3;
  vu5 = eta6d * (vu2 + 4.0 * vu3);
  ph1 = eta6d / rak;
  ph2 = d - e12 * dak2;
  ta1 = (reta + 5.0) / (5.0 * rak);
  ta2 = eta2d * dak2;
  ta3 = -e12 * rgek * (ta1 + ta2);
  ta4 = eta3d * ak2 * (ta1 * ta1 - ta2 * ta2);
  ta5 = eta3d * (reta + 8.0) * 1.0e-1 - 2.0 * eta22 * dak2;
  ex1 = exp(rak);
  ex2 = 0.0;
  if (gMSAWave[12] < 20.0) {
    ex2 = exp(-rak);
  }
  sk = 0.5 * (ex1 - ex2);
  ck = 0.5 * (ex1 + ex2);
  ckma = ck - 1.0 - rak * sk;
  skma = sk - rak * ck;
  a1 = (e24 * rgek * (al1 + al2 + ak1 * al3) - eta22) * dd4;
  if (ibig == 0) {
    a2 = e24 * (al3 * skma + al2 * sk - al1 * ck) * dd4;
    a3 = e24 * (eta22 * dak2 - 0.5 * d2 + al3 * ckma - al1 * sk + al2 * ck) * dd4;
  }
  b1 = (1.5 * reta * eta2d2 - e12 * rgek * (be1 + be2 + ak1 * be3)) * dd4;
  if (ibig == 0) {
    b2 = e12 * (-be3 * skma - be2 * sk + be1 * ck) * dd4;
    b3 = e12 * (0.5 * d2 * eta2d - eta3d * eta2d2 * dak2 - be3 * ckma + be1 * sk - be2 * ck) * dd4;
  }
  v1 = (eta21 * (eta2 - 2.0 * reta + 1.0e1) * 2.5e-1 - rgek * (vu4 + vu5)) * dd45;
  if (ibig == 0) {
    v2 = (vu4 * ck - vu5 * sk) * dd45;
    v3 = ((eta3 - 6.0 * eta2 + 5.0) * d - eta6d * (2.0 * eta3 - 3.0 * eta2 + 18.0 * reta + 1.0e1) * dak2 + e24 * vu3 + vu4 * sk - vu5 * ck) * dd45;
  }
  pp1 = ph1 * ph1;
  pp2 = ph2 * ph2;
  pp = pp1 + pp2;
  p1p2 = ph1 * ph2 * 2.0;
  p1 = (rgek * (pp1 + pp2 - p1p2) - 0.5 * eta2d) * dd2;
  if (ibig == 0) {
    p2 = (pp * sk + p1p2 * ck) * dd2;
    p3 = (pp * ck + p1p2 * sk + pp1 - pp2) * dd2;
  }
  t1 = ta3 + ta4 * a1 + ta5 * b1;
  if (ibig != 0) {
    v3 = ((eta3 - 6.0 * eta2 + 5.0) * d - eta6d * (2.0 * eta3 - 3.0 * eta2 + 18.0 * reta + 1.0e1) * dak2 + e24 * vu3) * dd45;
    t3 = ta4 * a3 + ta5 * b3 + e12 * ta2 - 4.0e-1 * reta * (reta + 1.0e1) - 1.0;
    p3 = (pp1 - pp2) * dd2;
    b3 = e12 * (0.5 * d2 * eta2d - eta3d * eta2d2 * dak2 + be3) * dd4;
    a3 = e24 * (eta22 * dak2 - 0.5 * d2 - al3) * dd4;
    um6 = t3 * a3 - e12 * v3 * v3;
    um5 = t1 * a3 + a1 * t3 - e24 * v1 * v3;
    um4 = t1 * a1 - e12 * v1 * v1;
    al6 = e12 * p3 * p3;
    al5 = e24 * p1 * p3 - b3 - b3 - ak2;
    al4 = e12 * p1 * p1 - b1 - b1;
    w56 = um5 * al6 - al5 * um6;
    w46 = um4 * al6 - al4 * um6;
    fa = -w46 / w56;
    ca = -fa;
    gMSAWave[3] = fa;
    gMSAWave[2] = ca;
    gMSAWave[1] = b1 + b3 * fa;
    gMSAWave[0] = a1 + a3 * fa;
    gMSAWave[8] = v1 + v3 * fa;
    gMSAWave[14] = -(p1 + p3 * fa);
    gMSAWave[15] = gMSAWave[14];
    if (fabs(gMSAWave[15]) < 1.0e-3) {
      gMSAWave[15] = 0.0;
    }
    gMSAWave[10] = gMSAWave[16];
  } else {
    t2 = ta4 * a2 + ta5 * b2 + e12 * (ta1 * ck - ta2 * sk);
    t3 = ta4 * a3 + ta5 * b3 + e12 * (ta1 * sk - ta2 * (ck - 1.0)) - 4.0e-1 * reta * (reta + 1.0e1) - 1.0;
    um1 = t2 * a2 - e12 * v2 * v2;
    um2 = t1 * a2 + t2 * a1 - e24 * v1 * v2;
    um3 = t2 * a3 + t3 * a2 - e24 * v2 * v3;
    um4 = t1 * a1 - e12 * v1 * v1;
    um5 = t1 * a3 + t3 * a1 - e24 * v1 * v3;
    um6 = t3 * a3 - e12 * v3 * v3;
    if ((ix == 1) || (ix == 3)) {
      al1 = e12 * p2 * p2;
      al2 = e24 * p1 * p2 - b2 - b2;
      al3 = e24 * p2 * p3;
      al4 = e12 * p1 * p1 - b1 - b1;
      al5 = e24 * p1 * p3 - b3 - b3 - ak2;
      al6 = e12 * p3 * p3;
      w16 = um1 * al6 - al1 * um6;
      w15 = um1 * al5 - al1 * um5;
      w14 = um1 * al4 - al1 * um4;
      w13 = um1 * al3 - al1 * um3;
      w12 = um1 * al2 - al1 * um2;
      w26 = um2 * al6 - al2 * um6;
      w25 = um2 * al5 - al2 * um5;
      w24 = um2 * al4 - al2 * um4;
      w36 = um3 * al6 - al3 * um6;
      w35 = um3 * al5 - al3 * um5;
      w34 = um3 * al4 - al3 * um4;
      w32 = um3 * al2 - al3 * um2;
      w46 = um4 * al6 - al4 * um6;
      w56 = um5 * al6 - al5 * um6;
      w3526 = w35 + w26;
      w3425 = w34 + w25;
      w4 = w16 * w16 - w13 * w36;
      w3 = 2.0 * w16 * w15 - w13 * w3526 - w12 * w36;
      w2 = w15 * w15 + 2.0 * w16 * w14 - w13 * w3425 - w12 * w3526;
      w1 = 2.0 * w15 * w14 - w13 * w24 - w12 * w3425;
      w0 = w14 * w14 - w12 * w24;
      if (ix == 1) {
        fap = (w14 - w34 - w46) / (w12 - w15 + w35 - w26 + w56 - w32);
      } else {
        gMSAWave[14] = 0.5 * eta2d * dd2 * exp(-rgek);
        if ((gMSAWave[11] <= 2.0) && (gMSAWave[11] >= 0.0) && (gMSAWave[12] <= 1.0)) {
          e24g = e24 * rgek * exp(rak);
          pwk = sqrt(e24g);
          qpw = (1.0 - sqrt(1.0 + 2.0 * d2 * d * pwk / eta22)) * eta21 / d;
          gMSAWave[14] = -qpw * qpw / e24 + 0.5 * eta2d * dd2;
        }
        pg = p1 + gMSAWave[14];
        ca = ak2 * pg + 2.0 * (b3 * pg - b1 * p3) + e12 * gMSAWave[14] * gMSAWave[14] * p3;
        ca = -ca / (ak2 * p2 + 2.0 * (b3 * p2 - b2 * p3));
        fap = -(pg + p2 * ca) / p3;
      }
      ii = 0;
      do {
        ii = ii + 1;
        if (ii > itm) {
          ir = -2;
          return (ir);
        }
        fa = fap;
        fun = w0 + (w1 + (w2 + (w3 + w4 * fa) * fa) * fa) * fa;
        fund = w1 + (2.0 * w2 + (3.0 * w3 + 4.0 * w4 * fa) * fa) * fa;
        fap = fa - fun / fund;
        del = fabs((fap - fa) / fa);
      } while (del > acc);
      ir = ir + ii;
      fa = fap;
      ca = -(w16 * fa * fa + w15 * fa + w14) / (w13 * fa + w12);
      gMSAWave[14] = -(p1 + p2 * ca + p3 * fa);
      gMSAWave[15] = gMSAWave[14];
      if (fabs(gMSAWave[15]) < 1.0e-3) {
        gMSAWave[15] = 0.0;
      }
      gMSAWave[10] = gMSAWave[16];
    } else {
      ca = ak2 * p1 + 2.0 * (b3 * p1 - b1 * p3);
      ca = -ca / (ak2 * p2 + 2.0 * (b3 * p2 - b2 * p3));
      fa = -(p1 + p2 * ca) / p3;
      if (ix == 2) {
        gMSAWave[15] = um1 * ca * ca + (um2 + um3 * fa) * ca + um4 + um5 * fa + um6 * fa * fa;
      }
      if (ix == 4) {
        gMSAWave[15] = -(p1 + p2 * ca + p3 * fa);
      }
    }
    gMSAWave[3] = fa;
    gMSAWave[2] = ca;
    gMSAWave[1] = b1 + b2 * ca + b3 * fa;
    gMSAWave[0] = a1 + a2 * ca + a3 * fa;
    gMSAWave[8] = (v1 + v2 * ca + v3 * fa) / gMSAWave[0];
  }
  g24 = e24 * rgek * ex1;
  gMSAWave[7] = (rak * ak2 * ca - g24) / (ak2 * g24);
  return (ir);
}
```