// File: /workspace/github_c_makefile_star2_O0/ugobas[P]tnm/mutation.c
// Function: Compute_mut_def
void Compute_mut_def(float *pred_mut_Cart, int N_Cart,
		     float *pred_mut_Tors, int N_axes,
		     float *DE, float *Force_coeff, float *Force,
		     struct Normal_Mode NM, float *sigma2, float mass)
{
  int a, i;
  for(i=0; i<N_Cart; i++)pred_mut_Cart[i]=0;
  for(i=0; i<N_axes; i++)pred_mut_Tors[i]=0;

  /*double F=0;
  for(i=0; i<N_Cart; i++){F+=Force[i]*Force[i];}
  printf("RMS of mutational force: %.3f\n", sqrt(3*F/N_Cart)); */
  double DE_tmp=0;

  // Compute predicted deformation of kinetic atoms (Cart)
  for(a=0; a<NM.N; a++){
    if((NM.Cart_coll[a]<COLL_THR_MUT)||(sigma2[a]<=0)){
      Force_coeff[a]=0; continue;
    }
    float *f=Force, *x=NM.Cart[a];
    double xf=0; // Projection of force on normal mode
    for(i=0; i<N_Cart; i++){xf+=(*x)*(*f); x++; f++;}
    Force_coeff[a]=xf;
    //(*DE)+=xf*xf/sigma2[a];
    DE_tmp+=xf*xf*sigma2[a]; // OK
    xf*=sigma2[a];
    float *m=pred_mut_Cart; x=NM.Cart[a];
    for(i=0; i<N_Cart; i++){
      (*m)+=(*x)*xf; x++; m++;
    }
    m=pred_mut_Tors; x=NM.Tors[a];
    for(i=0; i<N_axes; i++){
      (*m)+=(*x)*xf; x++; m++;
    }
  }
  *DE=DE_tmp;
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
movq 0xe8(%rbp), %rax
leaq 0x18(%rbp), %rax
movq %rax, -0x70(%rbp)
movq 0x10(%rbp), %rax
movq %rdi, -0x8(%rbp)
movl %esi, -0xc(%rbp)
movq %rdx, -0x18(%rbp)
movl %ecx, -0x1c(%rbp)
movq %r8, -0x28(%rbp)
movq %r9, -0x30(%rbp)
movss %xmm0, -0x34(%rbp)
movl $0x0, -0x3c(%rbp)
movl -0x3c(%rbp), %eax
cmpl -0xc(%rbp), %eax
jge 0x59f0c
movq -0x8(%rbp), %rax
movslq -0x3c(%rbp), %rcx
xorps %xmm0, %xmm0
movss %xmm0, (%rax,%rcx,4)
movl -0x3c(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x3c(%rbp)
jmp 0x59ee9
movl $0x0, -0x3c(%rbp)
movl -0x3c(%rbp), %eax
cmpl -0x1c(%rbp), %eax
jge 0x59f36
movq -0x18(%rbp), %rax
movslq -0x3c(%rbp), %rcx
xorps %xmm0, %xmm0
movss %xmm0, (%rax,%rcx,4)
movl -0x3c(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x3c(%rbp)
jmp 0x59f13
xorps %xmm0, %xmm0
movsd %xmm0, -0x48(%rbp)
movl $0x0, -0x38(%rbp)
movq -0x70(%rbp), %rcx
movl -0x38(%rbp), %eax
cmpl (%rcx), %eax
jge 0x5a179
movq -0x70(%rbp), %rax
movq 0x68(%rax), %rax
movslq -0x38(%rbp), %rcx
movss (%rax,%rcx,4), %xmm1
movss 0x31f9f(%rip), %xmm0    # 0x8bf0c
ucomiss %xmm1, %xmm0
ja 0x59f86
movq 0xe8(%rbp), %rax
movslq -0x38(%rbp), %rcx
xorps %xmm0, %xmm0
ucomiss (%rax,%rcx,4), %xmm0
jb 0x59f9b
movq -0x30(%rbp), %rax
movslq -0x38(%rbp), %rcx
xorps %xmm0, %xmm0
movss %xmm0, (%rax,%rcx,4)
jmp 0x5a16b
movq -0x70(%rbp), %rax
movq 0x10(%rbp), %rcx
movq %rcx, -0x50(%rbp)
movq 0x38(%rax), %rax
movslq -0x38(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movq %rax, -0x58(%rbp)
xorps %xmm0, %xmm0
movsd %xmm0, -0x60(%rbp)
movl $0x0, -0x3c(%rbp)
movl -0x3c(%rbp), %eax
cmpl -0xc(%rbp), %eax
jge 0x5a00f
movq -0x58(%rbp), %rax
movss (%rax), %xmm0
movq -0x50(%rbp), %rax
mulss (%rax), %xmm0
cvtss2sd %xmm0, %xmm0
addsd -0x60(%rbp), %xmm0
movsd %xmm0, -0x60(%rbp)
movq -0x58(%rbp), %rax
addq $0x4, %rax
movq %rax, -0x58(%rbp)
movq -0x50(%rbp), %rax
addq $0x4, %rax
movq %rax, -0x50(%rbp)
movl -0x3c(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x3c(%rbp)
jmp 0x59fc6
movq -0x70(%rbp), %rax
movsd -0x60(%rbp), %xmm0
cvtsd2ss %xmm0, %xmm0
movq -0x30(%rbp), %rcx
movslq -0x38(%rbp), %rdx
movss %xmm0, (%rcx,%rdx,4)
movsd -0x60(%rbp), %xmm0
mulsd -0x60(%rbp), %xmm0
movq 0xe8(%rbp), %rcx
movslq -0x38(%rbp), %rdx
movss (%rcx,%rdx,4), %xmm1
cvtss2sd %xmm1, %xmm2
movsd -0x48(%rbp), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x48(%rbp)
movq 0xe8(%rbp), %rcx
movslq -0x38(%rbp), %rdx
movss (%rcx,%rdx,4), %xmm0
cvtss2sd %xmm0, %xmm0
mulsd -0x60(%rbp), %xmm0
movsd %xmm0, -0x60(%rbp)
movq -0x8(%rbp), %rcx
movq %rcx, -0x68(%rbp)
movq 0x38(%rax), %rax
movslq -0x38(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movq %rax, -0x58(%rbp)
movl $0x0, -0x3c(%rbp)
movl -0x3c(%rbp), %eax
cmpl -0xc(%rbp), %eax
jge 0x5a0ee
movq -0x58(%rbp), %rax
movss (%rax), %xmm0
cvtss2sd %xmm0, %xmm0
movsd -0x60(%rbp), %xmm2
movq -0x68(%rbp), %rax
movss (%rax), %xmm1
cvtss2sd %xmm1, %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
cvtsd2ss %xmm0, %xmm0
movss %xmm0, (%rax)
movq -0x58(%rbp), %rax
addq $0x4, %rax
movq %rax, -0x58(%rbp)
movq -0x68(%rbp), %rax
addq $0x4, %rax
movq %rax, -0x68(%rbp)
movl -0x3c(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x3c(%rbp)
jmp 0x5a096
movq -0x70(%rbp), %rax
movq -0x18(%rbp), %rcx
movq %rcx, -0x68(%rbp)
movq 0x40(%rax), %rax
movslq -0x38(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movq %rax, -0x58(%rbp)
movl $0x0, -0x3c(%rbp)
movl -0x3c(%rbp), %eax
cmpl -0x1c(%rbp), %eax
jge 0x5a169
movq -0x58(%rbp), %rax
movss (%rax), %xmm0
cvtss2sd %xmm0, %xmm0
movsd -0x60(%rbp), %xmm2
movq -0x68(%rbp), %rax
movss (%rax), %xmm1
cvtss2sd %xmm1, %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
cvtsd2ss %xmm0, %xmm0
movss %xmm0, (%rax)
movq -0x58(%rbp), %rax
addq $0x4, %rax
movq %rax, -0x58(%rbp)
movq -0x68(%rbp), %rax
addq $0x4, %rax
movq %rax, -0x68(%rbp)
movl -0x3c(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x3c(%rbp)
jmp 0x5a111
jmp 0x5a16b
movl -0x38(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x38(%rbp)
jmp 0x59f45
movsd -0x48(%rbp), %xmm0
cvtsd2ss %xmm0, %xmm0
movq -0x28(%rbp), %rax
movss %xmm0, (%rax)
popq %rbp
retq
nopl (%rax)