// File: /workspace/github_c_makefile_star2_O2/ugobas[P]tnm/simulation.c
// Function: Optimize_amplitude
float Optimize_amplitude(struct bond *bonds_min, float *f_min,
			 float *Diff_Tors, int naxes,
			 struct bond *bonds,
			 int natoms1, struct Reference Ref1,
			 float *coord_ref_2)
{
  float EPS=0.005, RMSD_min=10000, RMSD_RRR=0; *f_min=-1;
  struct bond bonds_tmp[natoms1];
  float f[3], y[3], fmin=0.01, fmax=20; int i;
  f[0]=1.3; f[1]=1; f[2]=0.7;
  printf("Optimizing amplitude f of torsion angle change:\n");
  for(i=0; i<3; i++){
    y[i]=Compute_RMSD(f[i],Diff_Tors,naxes,
		      bonds,bonds_tmp,natoms1,Ref1,coord_ref_2);
    printf("RMSD=%.5f f=%.5f\n", y[i], f[i]);
    if(f[i]==1)RMSD_RRR=y[i];
    if(y[i]<RMSD_min){
      RMSD_min=y[i]; *f_min=f[i]; Copy_bonds(bonds_min, bonds_tmp, natoms1);
    }
    y[i]=exp(-y[i]);
  }

  for(int iter=0; iter<20; iter++){
    float ff=Find_max_quad(f[0],f[1],f[2],y[0],y[1],y[2],fmin,fmax);
    float yy=Compute_RMSD(ff, Diff_Tors, naxes, bonds, bonds_tmp, natoms1,
			  Ref1, coord_ref_2);
    printf("RMSD=%.5f f=%.5f\n", yy, ff);
    if(yy<RMSD_min){
      RMSD_min=yy; *f_min=ff; Copy_bonds(bonds_min, bonds_tmp, natoms1);
    }else{
      if(iter>10)break;
    }
    yy=exp(-yy);
    if(ff<f[0]){
      if((f[0]-ff)<EPS)break;
      f[2]=f[1]; y[2]=y[1]; f[1]=f[0]; y[1]=y[0]; f[0]=ff; y[0]=yy;
    }else if(ff < f[1]){
      if(((f[1]-ff)<EPS)||((ff-f[0])<EPS))break;
      f[2]=f[1]; y[2]=y[1]; f[1]=ff; y[1]=yy;
    }else if(ff < f[2]){
      if(((f[2]-ff)<EPS)||((ff-f[1])<EPS))break;
      f[0]=f[1]; y[0]=y[1]; f[1]=ff; y[1]=yy;
    }else{
      if((ff-f[2])<EPS)break;
      f[0]=f[1]; y[0]=y[1]; f[1]=f[2]; y[1]=y[2]; f[2]=ff; y[2]=yy;
    }
  }
  printf("RMSD RRR= %.4f factor= %.4g\n", RMSD_RRR, 1.0);
  printf("Minimal RMSD= %.4f factor= %.4g\n", RMSD_min, *f_min);
  return(RMSD_min);
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0xa8, %rsp
movl %r9d, %r14d
movq %r8, -0x98(%rbp)
movl %ecx, -0x3c(%rbp)
movq %rdx, -0x90(%rbp)
movq %rdi, -0x88(%rbp)
leaq 0x10(%rbp), %r12
movq %rsi, -0x58(%rbp)
movl $0xbf800000, (%rsi)     # imm = 0xBF800000
movl %r9d, %eax
movq %rsp, %rbx
imulq $0xc0, %rax, %rax
subq %rax, %rbx
movq %rbx, %rsp
movsd 0x1ef60(%rip), %xmm0    # 0x54cd0
movsd %xmm0, -0xd0(%rbp)
movl $0x3f333333, -0xc8(%rbp) # imm = 0x3F333333
leaq 0x1ffdb(%rip), %rdi     # 0x55d64
callq 0x2090
movss 0x1b06a(%rip), %xmm0    # 0x50e00
movss %xmm0, -0x2c(%rbp)
xorps %xmm0, %xmm0
movaps %xmm0, -0xc0(%rbp)
xorl %r13d, %r13d
leaq 0x1fd2f(%rip), %r15     # 0x55ade
movss 0x15259(%rip), %xmm1    # 0x4b010
cmpq $0x3, %r13
je 0x35ebd
movss -0xd0(%rbp,%r13,4), %xmm0
movss %xmm0, -0x80(%rbp)
subq $0x30, %rsp
movups (%r12), %xmm4
movups 0x10(%r12), %xmm3
movups 0x20(%r12), %xmm2
movups %xmm2, 0x20(%rsp)
movups %xmm3, 0x10(%rsp)
movups %xmm4, (%rsp)
movq -0x90(%rbp), %rdi
movl -0x3c(%rbp), %esi
movq -0x98(%rbp), %rdx
movq %rbx, %rcx
movl %r14d, %r8d
movq 0x40(%rbp), %r9
callq 0x3613e
addq $0x30, %rsp
movaps %xmm0, -0x70(%rbp)
cvtss2sd %xmm0, %xmm0
cvtss2sd -0x80(%rbp), %xmm1
movq %r15, %rdi
movb $0x2, %al
callq 0x20f0
movss -0x80(%rbp), %xmm2
movaps -0x70(%rbp), %xmm1
movaps %xmm2, %xmm0
cmpeqss 0x151cd(%rip), %xmm0    # 0x4b010
movaps %xmm0, -0x50(%rbp)
movss -0x2c(%rbp), %xmm0
ucomiss %xmm1, %xmm0
jbe 0x35e74
movq -0x58(%rbp), %rax
movss %xmm2, (%rax)
movq -0x88(%rbp), %rdi
movq %rbx, %rsi
movl %r14d, %edx
callq 0x301ca
movaps -0x70(%rbp), %xmm1
movss %xmm1, -0x2c(%rbp)
movaps %xmm1, %xmm0
movaps -0x50(%rbp), %xmm2
andps %xmm2, %xmm0
andnps -0xc0(%rbp), %xmm2
orps %xmm0, %xmm2
movaps %xmm2, -0x50(%rbp)
xorps 0x1a9cd(%rip), %xmm1    # 0x50860
cvtss2sd %xmm1, %xmm0
callq 0x22a0
cvtsd2ss %xmm0, %xmm0
movss %xmm0, -0xa4(%rbp,%r13,4)
incq %r13
movaps -0x50(%rbp), %xmm0
movaps %xmm0, -0xc0(%rbp)
jmp 0x35daf
movss -0xa4(%rbp), %xmm3
movss -0xa0(%rbp), %xmm4
xorl %r15d, %r15d
movss 0x1ee3c(%rip), %xmm0    # 0x54d14
movss 0x1ee38(%rip), %xmm2    # 0x54d18
movss -0x9c(%rbp), %xmm5
leaq 0x1fbef(%rip), %r13     # 0x55ade
cmpl $0x14, %r15d
je 0x360f1
movss %xmm0, -0x60(%rbp)
movss %xmm1, -0x70(%rbp)
movss %xmm2, -0x5c(%rbp)
movss %xmm3, -0x38(%rbp)
movss %xmm4, -0x30(%rbp)
movss %xmm5, -0x34(%rbp)
movss 0x1b2e1(%rip), %xmm6    # 0x51200
movss 0x1aecd(%rip), %xmm7    # 0x50df4
callq 0x3bda0
movss %xmm0, -0x50(%rbp)
subq $0x30, %rsp
movups (%r12), %xmm3
movups 0x10(%r12), %xmm1
movups 0x20(%r12), %xmm2
movups %xmm2, 0x20(%rsp)
movups %xmm1, 0x10(%rsp)
movups %xmm3, (%rsp)
movq -0x90(%rbp), %rdi
movl -0x3c(%rbp), %esi
movq -0x98(%rbp), %rdx
movq %rbx, %rcx
movl %r14d, %r8d
movq 0x40(%rbp), %r9
callq 0x3613e
addq $0x30, %rsp
movaps %xmm0, -0x80(%rbp)
cvtss2sd %xmm0, %xmm0
cvtss2sd -0x50(%rbp), %xmm1
movq %r13, %rdi
movb $0x2, %al
callq 0x20f0
movss -0x2c(%rbp), %xmm0
movaps -0x80(%rbp), %xmm1
ucomiss %xmm1, %xmm0
jbe 0x35fc7
movq -0x58(%rbp), %rax
movss -0x50(%rbp), %xmm0
movss %xmm0, (%rax)
movq -0x88(%rbp), %rdi
movq %rbx, %rsi
movl %r14d, %edx
callq 0x301ca
movaps -0x80(%rbp), %xmm0
movss %xmm0, -0x2c(%rbp)
jmp 0x35fd4
movaps %xmm1, %xmm0
cmpl $0xa, %r15d
ja 0x360f1
xorps 0x1a885(%rip), %xmm0    # 0x50860
cvtss2sd %xmm0, %xmm0
callq 0x22a0
cvtsd2ss %xmm0, %xmm0
movss -0x60(%rbp), %xmm8
movss -0x50(%rbp), %xmm9
ucomiss %xmm9, %xmm8
jbe 0x3602c
movaps %xmm8, %xmm1
subss %xmm9, %xmm1
movss 0x1ed11(%rip), %xmm2    # 0x54d1c
ucomiss %xmm1, %xmm2
movss -0x30(%rbp), %xmm2
movss -0x38(%rbp), %xmm4
movss -0x70(%rbp), %xmm6
movaps %xmm8, %xmm7
jbe 0x360d9
jmp 0x360f1
movss -0x70(%rbp), %xmm3
ucomiss %xmm9, %xmm3
movss 0x1ecdf(%rip), %xmm5    # 0x54d1c
jbe 0x36078
movaps %xmm3, %xmm1
subss %xmm9, %xmm1
ucomiss %xmm1, %xmm5
ja 0x360f1
movaps %xmm9, %xmm7
movaps %xmm9, %xmm1
subss %xmm8, %xmm1
ucomiss %xmm1, %xmm5
movss -0x30(%rbp), %xmm2
movaps %xmm0, %xmm4
movss -0x38(%rbp), %xmm0
movaps %xmm3, %xmm6
movaps %xmm8, %xmm9
jbe 0x360d9
jmp 0x360f1
movss -0x5c(%rbp), %xmm7
ucomiss %xmm9, %xmm7
jbe 0x360b7
movaps %xmm7, %xmm1
subss %xmm9, %xmm1
ucomiss %xmm1, %xmm5
ja 0x360f1
movaps %xmm7, %xmm6
movaps %xmm9, %xmm7
movaps %xmm9, %xmm2
subss %xmm3, %xmm2
ucomiss %xmm2, %xmm5
movss -0x34(%rbp), %xmm2
movaps %xmm0, %xmm4
movss -0x30(%rbp), %xmm0
movaps %xmm3, %xmm9
jbe 0x360d9
jmp 0x360f1
movaps %xmm9, %xmm1
subss %xmm7, %xmm1
ucomiss %xmm1, %xmm5
movaps %xmm9, %xmm6
movaps %xmm0, %xmm2
movss -0x30(%rbp), %xmm0
movaps %xmm3, %xmm9
movss -0x34(%rbp), %xmm4
ja 0x360f1
movaps %xmm7, %xmm1
incl %r15d
movaps %xmm2, %xmm5
movaps %xmm0, %xmm3
movaps %xmm6, %xmm2
movaps %xmm9, %xmm0
jmp 0x35eef
cvtss2sd -0xc0(%rbp), %xmm0
leaq 0x1f9f0(%rip), %rdi     # 0x55af0
movsd 0x14f40(%rip), %xmm1    # 0x4b048
movb $0x2, %al
callq 0x20f0
cvtss2sd -0x2c(%rbp), %xmm0
movq -0x58(%rbp), %rax
cvtss2sd (%rax), %xmm1
leaq 0x1f9ea(%rip), %rdi     # 0x55b0d
movb $0x2, %al
callq 0x20f0
movss -0x2c(%rbp), %xmm0
leaq -0x28(%rbp), %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq