// File: /workspace/github_c_makefile_star2_O3/ugobas[P]tnm/mutation.c
// Function: Normalize_parameters
void Normalize_parameters()
{
  // a.a. 20 = gap
  double E1=0, E2=0, D1=0, D2=0, na1=0, na2=0;
  int a, b, n=21;
  for(a=0; a<=20; a++){
    na1+=atom_numb[a];
    na2+=atom_numb[a]*atom_numb[a];
  }
  na1/=n; na2=sqrt(na2/n-na1*na1);
  for(a=0; a<=20; a++){
    natm_norm[a]=(atom_numb[a]-na1)/na2;
  }
  // E(gap)=0; D(gap)=min D_Res
  n=230;
  double D_min=1000;
  for(a=0; a<20; a++){
    for(b=0; b<=a; b++){
      E1+=Econt[a][b];
      E2+=Econt[a][b]*Econt[a][b];
      D1+=D_Res[a][b];
      D2+=D_Res[a][b]*D_Res[a][b];
      if(D_Res[a][b]<D_min)D_min=D_Res[a][b];
    }
  }
  E1/=n; E2=sqrt(E2/n-E1*E1);
  D1+=D_min*20; D2+=D_min*D_min*20;
  D1/=n; D2=sqrt(D2/n-D1*D1);
  float E_gap=-E1/E2;
  float D_gap=(D_min-D1)/D2;
  for(a=0; a<20; a++){
    for(b=0; b<20; b++){
      Econt_norm[a][b]=(Econt[a][b]-E1)/E2;
      D_Res_norm[a][b]=(D_Res[a][b]-D1)/D2;
    }
    Econt_norm[a][20]=E_gap;
    Econt_norm[20][a]=E_gap;
    D_Res_norm[a][20]=D_gap;
    D_Res_norm[20][a]=D_gap;
  }
  printf("Parameters of mutation model type, ave, s.d.:\n");
  printf("SIZE natom: %.3g %.3g\n\n", na1, na2);
  printf("DIST D_Res: %.3g %.3g\n", D1, D2);
  printf("STAB Econt: %.3g %.3g\n", E1, E2);
}

// Assembly:
pushq %rbx
subq $0x60, %rsp
xorpd %xmm0, %xmm0
xorl %eax, %eax
leaq 0x2985a(%rip), %rbx     # 0x68ef0
xorpd %xmm3, %xmm3
movl (%rax,%rbx), %ecx
xorps %xmm1, %xmm1
cvtsi2sd %ecx, %xmm1
imull %ecx, %ecx
xorps %xmm2, %xmm2
cvtsi2sd %ecx, %xmm2
addsd %xmm1, %xmm3
addsd %xmm2, %xmm0
addq $0x4, %rax
cmpq $0x54, %rax
jne 0x3f69a
movsd 0x1c8e0(%rip), %xmm1    # 0x5bfa8
divsd %xmm1, %xmm3
divsd %xmm1, %xmm0
movapd %xmm3, %xmm1
mulsd %xmm3, %xmm1
subsd %xmm1, %xmm0
xorpd %xmm1, %xmm1
ucomisd %xmm1, %xmm0
movsd %xmm3, 0x8(%rsp)
jb 0x3f6f2
sqrtsd %xmm0, %xmm0
jmp 0x3f6fd
callq 0x22a0
movsd 0x8(%rsp), %xmm3
xorl %eax, %eax
leaq 0x3a4da(%rip), %rcx     # 0x79be0
xorps %xmm1, %xmm1
cvtsi2sdl (%rax,%rbx), %xmm1
subsd %xmm3, %xmm1
divsd %xmm0, %xmm1
cvtsd2ss %xmm1, %xmm1
movss %xmm1, (%rax,%rcx)
addq $0x4, %rax
cmpq $0x54, %rax
jne 0x3f706
movsd %xmm0, 0x48(%rsp)
leaq 0x29a4a(%rip), %rbx     # 0x69180
movq (%rbx), %rax
movsd 0x1b61f(%rip), %xmm3    # 0x5ad60
xorpd %xmm4, %xmm4
movl $0x1, %ecx
leaq 0x2a07f(%rip), %rdx     # 0x697d0
xorl %esi, %esi
xorpd %xmm0, %xmm0
xorpd %xmm7, %xmm7
xorpd %xmm8, %xmm8
movq (%rax,%rsi,8), %rdi
xorl %r8d, %r8d
movapd %xmm3, %xmm1
movss (%rdi,%r8,4), %xmm2
xorps %xmm3, %xmm3
cvtss2sd %xmm2, %xmm3
addsd %xmm3, %xmm4
mulss %xmm2, %xmm2
cvtss2sd %xmm2, %xmm2
addsd %xmm2, %xmm0
movss (%rdx,%r8,4), %xmm2
xorps %xmm3, %xmm3
cvtss2sd %xmm2, %xmm3
addsd %xmm3, %xmm8
mulss %xmm2, %xmm2
cvtss2sd %xmm2, %xmm2
addsd %xmm2, %xmm7
minsd %xmm1, %xmm3
incq %r8
cmpq %r8, %rcx
jne 0x3f767
incq %rsi
incq %rcx
addq $0x50, %rdx
cmpq $0x14, %rsi
jne 0x3f760
movsd 0x1c7e6(%rip), %xmm2    # 0x5bfb0
divsd %xmm2, %xmm4
divsd %xmm2, %xmm0
movapd 0x118b6(%rip), %xmm1    # 0x51090
xorpd %xmm4, %xmm1
movapd %xmm1, 0x50(%rsp)
movapd %xmm4, %xmm1
mulsd %xmm4, %xmm1
subsd %xmm1, %xmm0
xorpd %xmm9, %xmm9
ucomisd %xmm9, %xmm0
movapd %xmm4, 0x30(%rsp)
movsd %xmm3, (%rsp)
jb 0x3f80d
sqrtsd %xmm0, %xmm5
jmp 0x3f842
movapd %xmm7, 0x20(%rsp)
movapd %xmm8, 0x10(%rsp)
callq 0x22a0
xorpd %xmm9, %xmm9
movsd 0x1c784(%rip), %xmm2    # 0x5bfb0
movapd 0x10(%rsp), %xmm8
movapd 0x20(%rsp), %xmm7
movsd (%rsp), %xmm3
movapd %xmm0, %xmm5
movsd 0x17b16(%rip), %xmm1    # 0x57360
movapd %xmm3, %xmm6
mulsd %xmm1, %xmm6
addsd %xmm8, %xmm6
movapd %xmm3, %xmm0
mulsd %xmm3, %xmm0
mulsd %xmm1, %xmm0
addsd %xmm7, %xmm0
divsd %xmm2, %xmm6
divsd %xmm2, %xmm0
movapd %xmm6, %xmm1
mulsd %xmm6, %xmm1
subsd %xmm1, %xmm0
ucomisd %xmm9, %xmm0
movsd %xmm5, 0x20(%rsp)
movsd %xmm6, 0x10(%rsp)
jb 0x3f894
sqrtsd %xmm0, %xmm0
jmp 0x3f8aa
callq 0x22a0
movsd 0x10(%rsp), %xmm6
movsd 0x20(%rsp), %xmm5
movsd (%rsp), %xmm3
movapd 0x30(%rsp), %xmm4
movapd 0x50(%rsp), %xmm1
divsd %xmm5, %xmm1
xorps %xmm7, %xmm7
cvtsd2ss %xmm1, %xmm7
subsd %xmm6, %xmm3
divsd %xmm0, %xmm3
xorps %xmm1, %xmm1
cvtsd2ss %xmm3, %xmm1
movq (%rbx), %rax
leaq 0x29ef6(%rip), %rcx     # 0x697d0
xorl %edx, %edx
leaq 0x3a35d(%rip), %rsi     # 0x79c40
leaq 0x3aa46(%rip), %rdi     # 0x7a330
xorl %r8d, %r8d
movq (%rax,%r8,8), %r9
xorl %r10d, %r10d
xorps %xmm2, %xmm2
cvtss2sd (%r9,%r10), %xmm2
subsd %xmm4, %xmm2
divsd %xmm5, %xmm2
cvtsd2ss %xmm2, %xmm2
leaq (%rdx,%r10), %r11
xorps %xmm3, %xmm3
cvtss2sd (%rcx,%r10), %xmm3
subsd %xmm6, %xmm3
divsd %xmm0, %xmm3
cvtsd2ss %xmm3, %xmm3
movss %xmm2, (%rsi,%r11)
movss %xmm3, (%rdi,%r11)
addq $0x4, %r10
cmpq $0x50, %r10
jne 0x3f8f4
imulq $0x54, %r8, %r9
movss %xmm7, 0x50(%r9,%rsi)
movss %xmm7, 0x690(%rsi,%r8,4)
movss %xmm1, 0x50(%r9,%rdi)
movss %xmm1, 0x690(%rdi,%r8,4)
incq %r8
addq $0x54, %rdx
addq $0x50, %rcx
cmpq $0x14, %r8
jne 0x3f8ed
leaq 0x1d24e(%rip), %rdi     # 0x5cbc8
movsd %xmm0, (%rsp)
callq 0x2090
leaq 0x1cf9b(%rip), %rdi     # 0x5c926
movsd 0x8(%rsp), %xmm0
movsd 0x48(%rsp), %xmm1
movb $0x2, %al
callq 0x20f0
leaq 0x1cf99(%rip), %rdi     # 0x5c93e
movsd 0x10(%rsp), %xmm0
movsd (%rsp), %xmm1
movb $0x2, %al
callq 0x20f0
leaq 0x1cf97(%rip), %rdi     # 0x5c955
movaps 0x30(%rsp), %xmm0
movsd 0x20(%rsp), %xmm1
movb $0x2, %al
addq $0x60, %rsp
popq %rbx
jmp 0x20f0