// File: /workspace/github_c_makefile_star2_O3/ugobas[P]tnm/allostery.c
// Function: Ave_site
float Ave_site(double *sigprof, double *ave_pair,
	       double *sigpair, int *n_notfound,
	       int *site, int *atomres, int nr, int *atom_site,
	       float *profile, float **coupling, int Na)
{
  double sum_prof=0; 
  *sigprof=0; *ave_pair=0; *sigpair=0; *n_notfound=0;
  int npair=0;
  for(int i1=0; i1<nr; i1++){
    int j1=atomres[site[i1]];
    if((j1<0)||(j1>=Na)){
      if(atom_site)atom_site[j1]=-site[i1];
      (*n_notfound)++;
      continue;
    }
    if(atom_site)atom_site[j1]=1;
    sum_prof+=profile[j1];
    *sigprof+=profile[j1]*profile[j1];
    for(int i2=0; i2<i1; i2++){
      int j2=atomres[site[i2]];
      if((j2<0)||(j2>=Na)||(j2==j1))continue;
      float x=coupling[j1][j2];
      *ave_pair+=x; *sigpair+=x*x;
      npair++;
    }
  } // end of the site
  Ave_St_dev(&sum_prof, sigprof, nr-*n_notfound);
  Ave_St_dev(ave_pair, sigpair, npair);
  return(sum_prof);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x18, %rsp
movq %rdx, %rbx
movq %rsi, %r14
movl 0x50(%rsp), %r11d
xorl %eax, %eax
movq %rdi, 0x10(%rsp)
movq %rax, (%rdi)
movq %rax, (%rsi)
movq %rax, (%rdx)
movl $0x0, (%rcx)
testl %r11d, %r11d
jle 0x477f9
movl 0x70(%rsp), %edx
movq 0x68(%rsp), %rsi
movq 0x58(%rsp), %r10
movl %r11d, %r11d
xorpd %xmm3, %xmm3
xorl %r12d, %r12d
xorl %ebp, %ebp
movslq (%r8,%r12,4), %rax
movslq (%r9,%rax,4), %r13
testl %r13d, %r13d
setns %dil
cmpl %edx, %r13d
setl %r15b
testb %r15b, %dil
je 0x477d7
testq %r10, %r10
je 0x47759
movl $0x1, (%r10,%r13,4)
movq 0x60(%rsp), %rax
movss (%rax,%r13,4), %xmm0
xorps %xmm1, %xmm1
cvtss2sd %xmm0, %xmm1
addsd %xmm1, %xmm3
mulss %xmm0, %xmm0
cvtss2sd %xmm0, %xmm0
movq 0x10(%rsp), %rax
addsd (%rax), %xmm0
movsd %xmm0, (%rax)
testq %r12, %r12
je 0x477e4
xorl %eax, %eax
movslq (%r8,%rax,4), %rdi
movl (%r9,%rdi,4), %edi
cmpl %r13d, %edi
je 0x477cd
testl %edi, %edi
js 0x477cd
cmpl %edx, %edi
jge 0x477cd
movq (%rsi,%r13,8), %r15
movss (%r15,%rdi,4), %xmm0
xorps %xmm1, %xmm1
cvtss2sd %xmm0, %xmm1
addsd (%r14), %xmm1
movsd %xmm1, (%r14)
mulss %xmm0, %xmm0
cvtss2sd %xmm0, %xmm0
addsd (%rbx), %xmm0
movsd %xmm0, (%rbx)
incl %ebp
incq %rax
cmpq %rax, %r12
jne 0x4778b
jmp 0x477e4
testq %r10, %r10
je 0x477e2
negl %eax
movl %eax, (%r10,%r13,4)
incl (%rcx)
incq %r12
cmpq %r11, %r12
jne 0x4772d
movl (%rcx), %eax
movl 0x50(%rsp), %r11d
jmp 0x477ff
xorpd %xmm3, %xmm3
xorl %ebp, %ebp
subl %eax, %r11d
jle 0x47868
xorpd %xmm0, %xmm0
cmpl $0x1, %r11d
je 0x4785f
xorps %xmm1, %xmm1
cvtsi2sd %r11d, %xmm1
divsd %xmm1, %xmm3
movq 0x10(%rsp), %rax
movsd (%rax), %xmm0
mulsd %xmm3, %xmm1
mulsd %xmm3, %xmm1
subsd %xmm1, %xmm0
decl %r11d
xorps %xmm1, %xmm1
cvtsi2sd %r11d, %xmm1
divsd %xmm1, %xmm0
xorpd %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jb 0x4784e
sqrtsd %xmm0, %xmm0
jmp 0x4785f
movsd %xmm3, 0x8(%rsp)
callq 0x22a0
movsd 0x8(%rsp), %xmm3
movq 0x10(%rsp), %rax
movsd %xmm0, (%rax)
testl %ebp, %ebp
jle 0x478cc
xorpd %xmm0, %xmm0
cmpl $0x1, %ebp
je 0x478c8
xorps %xmm1, %xmm1
cvtsi2sd %ebp, %xmm1
movsd (%r14), %xmm2
divsd %xmm1, %xmm2
movsd %xmm2, (%r14)
movsd (%rbx), %xmm0
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
subsd %xmm1, %xmm0
decl %ebp
xorps %xmm1, %xmm1
cvtsi2sd %ebp, %xmm1
divsd %xmm1, %xmm0
xorpd %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jb 0x478b7
sqrtsd %xmm0, %xmm0
jmp 0x478c8
movsd %xmm3, 0x8(%rsp)
callq 0x22a0
movsd 0x8(%rsp), %xmm3
movsd %xmm0, (%rbx)
xorps %xmm0, %xmm0
cvtsd2ss %xmm3, %xmm0
addq $0x18, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq