// File: /workspace/github_c_makefile_star2_O2/ugobas[P]tnm/anharmonic_tnm.c
// Function: Compute_V
float Compute_V(struct interaction *Int, float r)
{
  if(Int==NULL)return(0);
  double V=0, f;
  if(r<=Int->r1){
    // Knowledge based zone for SB or structure based if r0^SB < r0^KB
    if(POW){f=pow(r,-lambda);}else{f=exp(-lambda*r);}
    if(Int->A0)V=Int->A0;
    if(Int->A1)V+=Int->A1*f; 
    f*=f; if(Int->A2)V+=Int->A2*f;
    f*=f; if(Int->A4)V+=Int->A4*f;
  }else if(r<=Int->r2){
    double d=(r-Int->r2), d2=d*d, d3=d2*d;
    V=Int->a0+Int->a1*d+Int->a2*d2+Int->a3*d3+Int->a4*d3*d+Int->a5*d3*d2;
  }else if(r<=Int->r3){
    if(POW){f=pow(r,-lambda);}else{f=exp(-lambda*r);}
    if(Int->B0)V=Int->B0;
    if(Int->B1)V+=Int->B1*f; 
    f*=f; if(Int->B2)V+=Int->B2*f;
    f*=f; if(Int->B4)V+=Int->B4*f;
  }else if(r<C_THR){
    // Final interpolation region
    double d=(r-C_THR), d2=d*d, d3=d2*d;
    V=Int->b3*d3+Int->b4*d3*d+Int->b5*d3*d2;
  }
  return(V);
}

// Assembly:
testq %rdi, %rdi
je 0x1acae
pushq %rbx
movq %rdi, %rbx
movss 0x10(%rdi), %xmm1
ucomiss %xmm0, %xmm1
jae 0x1acb2
movss 0x14(%rbx), %xmm1
ucomiss %xmm0, %xmm1
jae 0x1ace6
movss 0x18(%rbx), %xmm1
ucomiss %xmm0, %xmm1
jae 0x1adf8
leaq 0x49416(%rip), %rax     # 0x64068
movss (%rax), %xmm2
xorps %xmm1, %xmm1
ucomiss %xmm0, %xmm2
jbe 0x1adf2
subss %xmm2, %xmm0
cvtss2sd %xmm0, %xmm0
movaps %xmm0, %xmm2
mulsd %xmm0, %xmm2
movapd %xmm2, %xmm1
mulsd %xmm0, %xmm1
movsd 0xa8(%rbx), %xmm3
mulsd %xmm1, %xmm3
mulsd %xmm0, %xmm3
movsd 0xa0(%rbx), %xmm0
mulsd %xmm1, %xmm0
addsd %xmm3, %xmm0
mulsd 0xb0(%rbx), %xmm1
mulsd %xmm2, %xmm1
addsd %xmm0, %xmm1
jmp 0x1adf2
xorps %xmm0, %xmm0
retq
leaq 0x48933(%rip), %rax     # 0x635ec
cmpl $0x0, (%rax)
je 0x1ad54
cvtss2sd %xmm0, %xmm0
leaq 0x49407(%rip), %rax     # 0x640d4
movss (%rax), %xmm1
xorps 0x35b88(%rip), %xmm1    # 0x50860
cvtss2sd %xmm1, %xmm1
callq 0x20a0
jmp 0x1ad73
subss %xmm1, %xmm0
cvtss2sd %xmm0, %xmm0
movaps %xmm0, %xmm2
mulsd %xmm0, %xmm2
movapd %xmm2, %xmm1
mulsd %xmm0, %xmm1
movsd 0x78(%rbx), %xmm3
mulsd %xmm0, %xmm3
addsd 0x70(%rbx), %xmm3
movsd 0x80(%rbx), %xmm4
mulsd %xmm2, %xmm4
addsd %xmm3, %xmm4
movsd 0x88(%rbx), %xmm3
mulsd %xmm1, %xmm3
addsd %xmm4, %xmm3
movsd 0x90(%rbx), %xmm4
mulsd %xmm1, %xmm4
mulsd %xmm0, %xmm4
addsd %xmm3, %xmm4
mulsd 0x98(%rbx), %xmm1
mulsd %xmm2, %xmm1
addsd %xmm4, %xmm1
jmp 0x1adf2
leaq 0x49379(%rip), %rax     # 0x640d4
movss (%rax), %xmm1
xorps 0x35afa(%rip), %xmm1    # 0x50860
mulss %xmm0, %xmm1
cvtss2sd %xmm1, %xmm0
callq 0x22a0
movsd 0x30(%rbx), %xmm1
movsd 0x38(%rbx), %xmm3
xorpd %xmm2, %xmm2
movapd %xmm1, %xmm4
cmpeqsd %xmm2, %xmm4
andnpd %xmm1, %xmm4
movapd %xmm3, %xmm1
mulsd %xmm0, %xmm1
addsd %xmm4, %xmm1
cmpneqsd %xmm2, %xmm3
movapd %xmm3, %xmm5
andnpd %xmm4, %xmm5
andpd %xmm1, %xmm3
orpd %xmm5, %xmm3
mulsd %xmm0, %xmm0
movsd 0x40(%rbx), %xmm1
movapd %xmm1, %xmm4
mulsd %xmm0, %xmm4
addsd %xmm3, %xmm4
cmpneqsd %xmm2, %xmm1
movapd %xmm1, %xmm5
andnpd %xmm3, %xmm5
andpd %xmm4, %xmm1
orpd %xmm5, %xmm1
movsd 0x48(%rbx), %xmm3
ucomisd %xmm2, %xmm3
jne 0x1ade6
jnp 0x1adf2
mulsd %xmm0, %xmm0
mulsd %xmm0, %xmm3
addsd %xmm3, %xmm1
cvtsd2ss %xmm1, %xmm0
popq %rbx
retq
leaq 0x487ed(%rip), %rax     # 0x635ec
cmpl $0x0, (%rax)
je 0x1ae25
cvtss2sd %xmm0, %xmm0
leaq 0x492c5(%rip), %rax     # 0x640d4
movss (%rax), %xmm1
xorps 0x35a46(%rip), %xmm1    # 0x50860
cvtss2sd %xmm1, %xmm1
callq 0x20a0
jmp 0x1ae44
leaq 0x492a8(%rip), %rax     # 0x640d4
movss (%rax), %xmm1
xorps 0x35a29(%rip), %xmm1    # 0x50860
mulss %xmm0, %xmm1
cvtss2sd %xmm1, %xmm0
callq 0x22a0
movsd 0x50(%rbx), %xmm1
movsd 0x58(%rbx), %xmm3
xorpd %xmm2, %xmm2
movapd %xmm1, %xmm4
cmpeqsd %xmm2, %xmm4
andnpd %xmm1, %xmm4
movapd %xmm3, %xmm1
mulsd %xmm0, %xmm1
addsd %xmm4, %xmm1
cmpneqsd %xmm2, %xmm3
movapd %xmm3, %xmm5
andnpd %xmm4, %xmm5
andpd %xmm1, %xmm3
orpd %xmm5, %xmm3
mulsd %xmm0, %xmm0
movsd 0x60(%rbx), %xmm1
movapd %xmm1, %xmm4
mulsd %xmm0, %xmm4
addsd %xmm3, %xmm4
cmpneqsd %xmm2, %xmm1
movapd %xmm1, %xmm5
andnpd %xmm3, %xmm5
andpd %xmm4, %xmm1
orpd %xmm5, %xmm1
movsd 0x68(%rbx), %xmm3
jmp 0x1adde