// File: /workspace/github_c_makefile_star2_O0/josch[P]stride/hbenergy.c
// Function: GRID_Energy
void GRID_Energy(float *CA2, float *C, float *O, float *H, float *N, COMMAND *Cmd, HBOND *HBond)
{

  float ProjH[3];

 /***** Distance dependence ( 8-6 potential ) ****/

  if( Cmd->Truncate && HBond->AccDonDist < RmGRID ) 
    HBond->AccDonDist = RmGRID;
  HBond->Er = CGRID/pow(HBond->AccDonDist,8.0) - DGRID/pow(HBond->AccDonDist,6.0);

 /************** Angular dependance ****************/

 /* Find projection of the hydrogen on the O-C-CA plane */
  Project4_123(O,C,CA2,H,ProjH); 


 /* Three angles determining the direction of the hydrogen bond */
  HBond->ti = fabs(180.0 - Ang(ProjH,O,C)); 
  HBond->to = Ang(H,O,ProjH);             
  HBond->p  = Ang(N,H,O);

 /* Calculate both angle-dependent HB energy components Et and Ep */ 
  if( HBond->ti >= 0.0 && HBond->ti < 90.0 )   
    HBond->Et = cos(RAD(HBond->to))*(0.9+0.1*sin(RAD(2*HBond->ti)));
  else
  if( HBond->ti >= 90.0 && HBond->ti < 110.0 ) 
    HBond->Et = K1GRID*cos(RAD(HBond->to))*
      (pow((K2GRID-pow(cos(RAD(HBond->ti)),2.0)),3.0));
  else
    HBond->Et = 0.0;

  if( HBond->p > 90.0 && HBond->p < 270.0 )
    HBond->Ep = pow(cos(RAD(HBond->p)),2.0);
  else
    HBond->Ep = 0.0;

    /******** Full hydrogen bond energy *********************/
  HBond->Energy = 1000.0*HBond->Er*HBond->Et*HBond->Ep;
}

// Assembly:
subq $0x78, %rsp
movq 0x80(%rsp), %rax
movq %rdi, 0x70(%rsp)
movq %rsi, 0x68(%rsp)
movq %rdx, 0x60(%rsp)
movq %rcx, 0x58(%rsp)
movq %r8, 0x50(%rsp)
movq %r9, 0x48(%rsp)
movq 0x48(%rsp), %rax
cmpl $0x0, 0x1d64(%rax)
je 0x6dcc
movq 0x80(%rsp), %rax
movss 0x38(%rax), %xmm0
cvtss2sd %xmm0, %xmm1
movsd 0x24cd7(%rip), %xmm0    # 0x2ba88
ucomisd %xmm1, %xmm0
jbe 0x6dcc
movq 0x80(%rsp), %rax
movss 0x24d29(%rip), %xmm0    # 0x2baf0
movss %xmm0, 0x38(%rax)
movsd 0x24cb4(%rip), %xmm0    # 0x2ba88
movsd 0x24cc4(%rip), %xmm1    # 0x2baa0
callq 0x2090
movaps %xmm0, %xmm1
movsd 0x24cbc(%rip), %xmm0    # 0x2baa8
mulsd %xmm1, %xmm0
movsd %xmm0, 0x20(%rsp)
movq 0x80(%rsp), %rax
movss 0x38(%rax), %xmm0
cvtss2sd %xmm0, %xmm0
movsd 0x24c91(%rip), %xmm1    # 0x2baa0
callq 0x2090
movaps %xmm0, %xmm1
movsd 0x20(%rsp), %xmm0
divsd %xmm1, %xmm0
movsd %xmm0, 0x30(%rsp)
movsd 0x24c59(%rip), %xmm0    # 0x2ba88
movsd 0x24c59(%rip), %xmm1    # 0x2ba90
callq 0x2090
movaps %xmm0, %xmm1
movsd 0x24c51(%rip), %xmm0    # 0x2ba98
mulsd %xmm1, %xmm0
movsd %xmm0, 0x28(%rsp)
movq 0x80(%rsp), %rax
movss 0x38(%rax), %xmm0
cvtss2sd %xmm0, %xmm0
movsd 0x24c26(%rip), %xmm1    # 0x2ba90
callq 0x2090
movsd 0x28(%rsp), %xmm1
movaps %xmm0, %xmm2
movsd 0x30(%rsp), %xmm0
divsd %xmm2, %xmm1
subsd %xmm1, %xmm0
cvtsd2ss %xmm0, %xmm0
movq 0x80(%rsp), %rax
movss %xmm0, 0x20(%rax)
movq 0x60(%rsp), %rdi
movq 0x68(%rsp), %rsi
movq 0x70(%rsp), %rdx
movq 0x58(%rsp), %rcx
leaq 0x3c(%rsp), %r8
callq 0x55e0
leaq 0x3c(%rsp), %rdi
movq 0x60(%rsp), %rsi
movq 0x68(%rsp), %rdx
callq 0x4ad0
cvtss2sd %xmm0, %xmm1
movsd 0x2493b(%rip), %xmm0    # 0x2b810
subsd %xmm1, %xmm0
movaps 0x24970(%rip), %xmm1    # 0x2b850
pand %xmm1, %xmm0
cvtsd2ss %xmm0, %xmm0
movq 0x80(%rsp), %rax
movss %xmm0, 0x2c(%rax)
movq 0x58(%rsp), %rdi
movq 0x60(%rsp), %rsi
leaq 0x3c(%rsp), %rdx
callq 0x4ad0
movq 0x80(%rsp), %rax
movss %xmm0, 0x30(%rax)
movq 0x50(%rsp), %rdi
movq 0x58(%rsp), %rsi
movq 0x60(%rsp), %rdx
callq 0x4ad0
movq 0x80(%rsp), %rax
movss %xmm0, 0x34(%rax)
movq 0x80(%rsp), %rax
movss 0x2c(%rax), %xmm0
cvtss2sd %xmm0, %xmm0
xorps %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jb 0x701d
movq 0x80(%rsp), %rax
movss 0x2c(%rax), %xmm0
cvtss2sd %xmm0, %xmm1
movsd 0x248b2(%rip), %xmm0    # 0x2b820
ucomisd %xmm1, %xmm0
jbe 0x701d
movq 0x80(%rsp), %rax
movss 0x30(%rax), %xmm0
cvtss2sd %xmm0, %xmm0
movsd 0x24887(%rip), %xmm1    # 0x2b818
mulsd %xmm1, %xmm0
movsd 0x24873(%rip), %xmm1    # 0x2b810
divsd %xmm1, %xmm0
callq 0x20f0
movsd %xmm0, 0x18(%rsp)
movq 0x80(%rsp), %rax
movss 0x24b38(%rip), %xmm0    # 0x2baf4
mulss 0x2c(%rax), %xmm0
cvtss2sd %xmm0, %xmm0
movsd 0x2484b(%rip), %xmm1    # 0x2b818
mulsd %xmm1, %xmm0
movsd 0x24837(%rip), %xmm1    # 0x2b810
divsd %xmm1, %xmm0
callq 0x21b0
movaps %xmm0, %xmm1
movsd 0x18(%rsp), %xmm0
movsd 0x24add(%rip), %xmm2    # 0x2bad0
mulsd %xmm2, %xmm1
movsd 0x24ac9(%rip), %xmm2    # 0x2bac8
addsd %xmm2, %xmm1
mulsd %xmm1, %xmm0
cvtsd2ss %xmm0, %xmm0
movq 0x80(%rsp), %rax
movss %xmm0, 0x24(%rax)
jmp 0x7178
movq 0x80(%rsp), %rax
movss 0x2c(%rax), %xmm0
cvtss2sd %xmm0, %xmm0
movsd 0x247ea(%rip), %xmm1    # 0x2b820
ucomisd %xmm1, %xmm0
jb 0x7166
movq 0x80(%rsp), %rax
movss 0x2c(%rax), %xmm0
cvtss2sd %xmm0, %xmm1
movsd 0x24a57(%rip), %xmm0    # 0x2bab0
ucomisd %xmm1, %xmm0
jbe 0x7166
movsd 0x24a55(%rip), %xmm0    # 0x2bac0
callq 0x20f0
movsd 0x24a18(%rip), %xmm1    # 0x2ba90
callq 0x2090
movaps %xmm0, %xmm1
movsd 0x24a40(%rip), %xmm0    # 0x2bac8
divsd %xmm1, %xmm0
movsd %xmm0, (%rsp)
movq 0x80(%rsp), %rax
movss 0x30(%rax), %xmm0
cvtss2sd %xmm0, %xmm0
movsd 0x2476e(%rip), %xmm1    # 0x2b818
mulsd %xmm1, %xmm0
movsd 0x2475a(%rip), %xmm1    # 0x2b810
divsd %xmm1, %xmm0
callq 0x20f0
movaps %xmm0, %xmm1
movsd (%rsp), %xmm0
mulsd %xmm1, %xmm0
movsd %xmm0, 0x10(%rsp)
movsd 0x249e7(%rip), %xmm0    # 0x2bac0
callq 0x20f0
movsd 0x249d2(%rip), %xmm1    # 0x2bab8
callq 0x2090
movsd %xmm0, 0x8(%rsp)
movq 0x80(%rsp), %rax
movss 0x2c(%rax), %xmm0
cvtss2sd %xmm0, %xmm0
movsd 0x2470e(%rip), %xmm1    # 0x2b818
mulsd %xmm1, %xmm0
movsd 0x246fa(%rip), %xmm1    # 0x2b810
divsd %xmm1, %xmm0
callq 0x20f0
movsd 0x24991(%rip), %xmm1    # 0x2bab8
callq 0x2090
movaps %xmm0, %xmm1
movsd 0x8(%rsp), %xmm0
subsd %xmm1, %xmm0
movsd 0x24947(%rip), %xmm1    # 0x2ba88
callq 0x2090
movaps %xmm0, %xmm1
movsd 0x10(%rsp), %xmm0
mulsd %xmm1, %xmm0
cvtsd2ss %xmm0, %xmm0
movq 0x80(%rsp), %rax
movss %xmm0, 0x24(%rax)
jmp 0x7176
movq 0x80(%rsp), %rax
xorps %xmm0, %xmm0
movss %xmm0, 0x24(%rax)
jmp 0x7178
movq 0x80(%rsp), %rax
movss 0x34(%rax), %xmm0
cvtss2sd %xmm0, %xmm0
movsd 0x2468f(%rip), %xmm1    # 0x2b820
ucomisd %xmm1, %xmm0
jbe 0x7204
movq 0x80(%rsp), %rax
movss 0x34(%rax), %xmm0
cvtss2sd %xmm0, %xmm1
movsd 0x24928(%rip), %xmm0    # 0x2bad8
ucomisd %xmm1, %xmm0
jbe 0x7204
movq 0x80(%rsp), %rax
movss 0x34(%rax), %xmm0
cvtss2sd %xmm0, %xmm0
movsd 0x24649(%rip), %xmm1    # 0x2b818
mulsd %xmm1, %xmm0
movsd 0x24635(%rip), %xmm1    # 0x2b810
divsd %xmm1, %xmm0
callq 0x20f0
movsd 0x248cc(%rip), %xmm1    # 0x2bab8
callq 0x2090
cvtsd2ss %xmm0, %xmm0
movq 0x80(%rsp), %rax
movss %xmm0, 0x28(%rax)
jmp 0x7214
movq 0x80(%rsp), %rax
xorps %xmm0, %xmm0
movss %xmm0, 0x28(%rax)
movq 0x80(%rsp), %rax
movss 0x20(%rax), %xmm0
cvtss2sd %xmm0, %xmm1
movsd 0x248b3(%rip), %xmm0    # 0x2bae0
mulsd %xmm1, %xmm0
movq 0x80(%rsp), %rax
movss 0x24(%rax), %xmm1
cvtss2sd %xmm1, %xmm1
mulsd %xmm1, %xmm0
movq 0x80(%rsp), %rax
movss 0x28(%rax), %xmm1
cvtss2sd %xmm1, %xmm1
mulsd %xmm1, %xmm0
cvtsd2ss %xmm0, %xmm0
movq 0x80(%rsp), %rax
movss %xmm0, 0x1c(%rax)
addq $0x78, %rsp
retq
nopw %cs:(%rax,%rax)