// File: /workspace/github_c_makefile_star2_O0/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 %rbp
movq %rsp, %rbp
subq $0x50, %rsp
xorps %xmm0, %xmm0
movsd %xmm0, -0x8(%rbp)
xorps %xmm0, %xmm0
movsd %xmm0, -0x10(%rbp)
xorps %xmm0, %xmm0
movsd %xmm0, -0x18(%rbp)
xorps %xmm0, %xmm0
movsd %xmm0, -0x20(%rbp)
xorps %xmm0, %xmm0
movsd %xmm0, -0x28(%rbp)
xorps %xmm0, %xmm0
movsd %xmm0, -0x30(%rbp)
movl $0x15, -0x3c(%rbp)
movl $0x0, -0x34(%rbp)
cmpl $0x14, -0x34(%rbp)
jg 0x5a50c
movslq -0x34(%rbp), %rcx
leaq 0x319e9(%rip), %rax     # 0x8beb0
cvtsi2sdl (%rax,%rcx,4), %xmm0
addsd -0x28(%rbp), %xmm0
movsd %xmm0, -0x28(%rbp)
movslq -0x34(%rbp), %rcx
leaq 0x319cf(%rip), %rax     # 0x8beb0
movl (%rax,%rcx,4), %eax
movslq -0x34(%rbp), %rdx
leaq 0x319c1(%rip), %rcx     # 0x8beb0
imull (%rcx,%rdx,4), %eax
cvtsi2sd %eax, %xmm0
addsd -0x30(%rbp), %xmm0
movsd %xmm0, -0x30(%rbp)
movl -0x34(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x34(%rbp)
jmp 0x5a4b6
cvtsi2sdl -0x3c(%rbp), %xmm1
movsd -0x28(%rbp), %xmm0
divsd %xmm1, %xmm0
movsd %xmm0, -0x28(%rbp)
movsd -0x30(%rbp), %xmm1
cvtsi2sdl -0x3c(%rbp), %xmm0
divsd %xmm0, %xmm1
movsd -0x28(%rbp), %xmm0
movsd -0x28(%rbp), %xmm2
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
callq 0x1220
movsd %xmm0, -0x30(%rbp)
movl $0x0, -0x34(%rbp)
cmpl $0x14, -0x34(%rbp)
jg 0x5a5a6
movslq -0x34(%rbp), %rcx
leaq 0x31938(%rip), %rax     # 0x8beb0
cvtsi2sdl (%rax,%rcx,4), %xmm0
subsd -0x28(%rbp), %xmm0
divsd -0x30(%rbp), %xmm0
cvtsd2ss %xmm0, %xmm0
movslq -0x34(%rbp), %rcx
leaq 0x4261a(%rip), %rax     # 0x9cbb0
movss %xmm0, (%rax,%rcx,4)
movl -0x34(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x34(%rbp)
jmp 0x5a567
movl $0xe6, -0x3c(%rbp)
movsd 0x26d63(%rip), %xmm0    # 0x81318
movsd %xmm0, -0x48(%rbp)
movl $0x0, -0x34(%rbp)
cmpl $0x14, -0x34(%rbp)
jge 0x5a726
movl $0x0, -0x38(%rbp)
movl -0x38(%rbp), %eax
cmpl -0x34(%rbp), %eax
jg 0x5a716
leaq 0x31b5b(%rip), %rax     # 0x8c140
movq (%rax), %rax
movslq -0x34(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movslq -0x38(%rbp), %rcx
movss (%rax,%rcx,4), %xmm0
cvtss2sd %xmm0, %xmm0
addsd -0x8(%rbp), %xmm0
movsd %xmm0, -0x8(%rbp)
leaq 0x31b32(%rip), %rax     # 0x8c140
movq (%rax), %rax
movslq -0x34(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movslq -0x38(%rbp), %rcx
movss (%rax,%rcx,4), %xmm0
leaq 0x31b17(%rip), %rax     # 0x8c140
movq (%rax), %rax
movslq -0x34(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movslq -0x38(%rbp), %rcx
mulss (%rax,%rcx,4), %xmm0
cvtss2sd %xmm0, %xmm0
addsd -0x10(%rbp), %xmm0
movsd %xmm0, -0x10(%rbp)
movslq -0x34(%rbp), %rcx
leaq 0x3213a(%rip), %rax     # 0x8c790
imulq $0x50, %rcx, %rcx
addq %rcx, %rax
movslq -0x38(%rbp), %rcx
movss (%rax,%rcx,4), %xmm0
cvtss2sd %xmm0, %xmm0
addsd -0x18(%rbp), %xmm0
movsd %xmm0, -0x18(%rbp)
movslq -0x34(%rbp), %rcx
leaq 0x32111(%rip), %rax     # 0x8c790
imulq $0x50, %rcx, %rcx
addq %rcx, %rax
movslq -0x38(%rbp), %rcx
movss (%rax,%rcx,4), %xmm0
movslq -0x34(%rbp), %rcx
leaq 0x320f6(%rip), %rax     # 0x8c790
imulq $0x50, %rcx, %rcx
addq %rcx, %rax
movslq -0x38(%rbp), %rcx
mulss (%rax,%rcx,4), %xmm0
cvtss2sd %xmm0, %xmm0
addsd -0x20(%rbp), %xmm0
movsd %xmm0, -0x20(%rbp)
movslq -0x34(%rbp), %rcx
leaq 0x320cd(%rip), %rax     # 0x8c790
imulq $0x50, %rcx, %rcx
addq %rcx, %rax
movslq -0x38(%rbp), %rcx
movss (%rax,%rcx,4), %xmm0
cvtss2sd %xmm0, %xmm1
movsd -0x48(%rbp), %xmm0
ucomisd %xmm1, %xmm0
jbe 0x5a706
movslq -0x34(%rbp), %rcx
leaq 0x320a3(%rip), %rax     # 0x8c790
imulq $0x50, %rcx, %rcx
addq %rcx, %rax
movslq -0x38(%rbp), %rcx
movss (%rax,%rcx,4), %xmm0
cvtss2sd %xmm0, %xmm0
movsd %xmm0, -0x48(%rbp)
jmp 0x5a708
movl -0x38(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x38(%rbp)
jmp 0x5a5d2
jmp 0x5a718
movl -0x34(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x34(%rbp)
jmp 0x5a5c1
cvtsi2sdl -0x3c(%rbp), %xmm1
movsd -0x8(%rbp), %xmm0
divsd %xmm1, %xmm0
movsd %xmm0, -0x8(%rbp)
movsd -0x10(%rbp), %xmm1
cvtsi2sdl -0x3c(%rbp), %xmm0
divsd %xmm0, %xmm1
movsd -0x8(%rbp), %xmm0
movsd -0x8(%rbp), %xmm2
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
callq 0x1220
movsd %xmm0, -0x10(%rbp)
movsd -0x48(%rbp), %xmm0
movsd -0x18(%rbp), %xmm1
movsd 0x21cb4(%rip), %xmm2    # 0x7c440
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x18(%rbp)
movsd -0x48(%rbp), %xmm0
mulsd -0x48(%rbp), %xmm0
movsd -0x20(%rbp), %xmm1
movsd 0x21c90(%rip), %xmm2    # 0x7c440
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x20(%rbp)
cvtsi2sdl -0x3c(%rbp), %xmm1
movsd -0x18(%rbp), %xmm0
divsd %xmm1, %xmm0
movsd %xmm0, -0x18(%rbp)
movsd -0x20(%rbp), %xmm1
cvtsi2sdl -0x3c(%rbp), %xmm0
divsd %xmm0, %xmm1
movsd -0x18(%rbp), %xmm0
movsd -0x18(%rbp), %xmm2
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
callq 0x1220
movsd %xmm0, -0x20(%rbp)
movsd -0x8(%rbp), %xmm0
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
divsd -0x10(%rbp), %xmm0
cvtsd2ss %xmm0, %xmm0
movss %xmm0, -0x4c(%rbp)
movsd -0x48(%rbp), %xmm0
subsd -0x18(%rbp), %xmm0
divsd -0x20(%rbp), %xmm0
cvtsd2ss %xmm0, %xmm0
movss %xmm0, -0x50(%rbp)
movl $0x0, -0x34(%rbp)
cmpl $0x14, -0x34(%rbp)
jge 0x5a98b
movl $0x0, -0x38(%rbp)
cmpl $0x14, -0x38(%rbp)
jge 0x5a913
leaq 0x318c4(%rip), %rax     # 0x8c140
movq (%rax), %rax
movslq -0x34(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movslq -0x38(%rbp), %rcx
movss (%rax,%rcx,4), %xmm0
cvtss2sd %xmm0, %xmm0
subsd -0x8(%rbp), %xmm0
divsd -0x10(%rbp), %xmm0
cvtsd2ss %xmm0, %xmm0
movslq -0x34(%rbp), %rcx
leaq 0x42363(%rip), %rax     # 0x9cc10
imulq $0x54, %rcx, %rcx
addq %rcx, %rax
movslq -0x38(%rbp), %rcx
movss %xmm0, (%rax,%rcx,4)
movslq -0x34(%rbp), %rcx
leaq 0x31ec8(%rip), %rax     # 0x8c790
imulq $0x50, %rcx, %rcx
addq %rcx, %rax
movslq -0x38(%rbp), %rcx
movss (%rax,%rcx,4), %xmm0
cvtss2sd %xmm0, %xmm0
subsd -0x18(%rbp), %xmm0
divsd -0x20(%rbp), %xmm0
cvtsd2ss %xmm0, %xmm0
movslq -0x34(%rbp), %rcx
leaq 0x42a0b(%rip), %rax     # 0x9d300
imulq $0x54, %rcx, %rcx
addq %rcx, %rax
movslq -0x38(%rbp), %rcx
movss %xmm0, (%rax,%rcx,4)
movl -0x38(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x38(%rbp)
jmp 0x5a86b
movss -0x4c(%rbp), %xmm0
movslq -0x34(%rbp), %rcx
leaq 0x422ed(%rip), %rax     # 0x9cc10
imulq $0x54, %rcx, %rcx
addq %rcx, %rax
movss %xmm0, 0x50(%rax)
movss -0x4c(%rbp), %xmm0
movslq -0x34(%rbp), %rcx
leaq 0x422d1(%rip), %rax     # 0x9cc10
movss %xmm0, 0x690(%rax,%rcx,4)
movss -0x50(%rbp), %xmm0
movslq -0x34(%rbp), %rcx
leaq 0x429a8(%rip), %rax     # 0x9d300
imulq $0x54, %rcx, %rcx
addq %rcx, %rax
movss %xmm0, 0x50(%rax)
movss -0x50(%rbp), %xmm0
movslq -0x34(%rbp), %rcx
leaq 0x4298c(%rip), %rax     # 0x9d300
movss %xmm0, 0x690(%rax,%rcx,4)
movl -0x34(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x34(%rbp)
jmp 0x5a85a
leaq 0x2735c(%rip), %rdi     # 0x81cee
movb $0x0, %al
callq 0x10d0
movsd -0x28(%rbp), %xmm0
movsd -0x30(%rbp), %xmm1
leaq 0x27373(%rip), %rdi     # 0x81d1d
movb $0x2, %al
callq 0x10d0
movsd -0x18(%rbp), %xmm0
movsd -0x20(%rbp), %xmm1
leaq 0x27373(%rip), %rdi     # 0x81d35
movb $0x2, %al
callq 0x10d0
movsd -0x8(%rbp), %xmm0
movsd -0x10(%rbp), %xmm1
leaq 0x27372(%rip), %rdi     # 0x81d4c
movb $0x2, %al
callq 0x10d0
addq $0x50, %rsp
popq %rbp
retq
nopw (%rax,%rax)