// File: /workspace/github_c_makefile_star2_O0/ugobas[P]tnm/mutation.c
// Function: Pred_mut_str
float Pred_mut_str(float C_size, float C_stab, float C_dist,
		   float *DE, float *MSD_mut, float *MSD_cross,
		   float *A_mut, float *A_nomut, float *A_cross, float *offset,
		   float *pred_mut_3, float *pred_mut_2, float *pred_mut_tot_2,
		   int Na, int *n_cont_mut, float *Force, float *Force_coeff,
		   float *pred_mut_Cart, int N_Cart, 
		   float *pred_mut_Tors, int N_axes,
		   char *AAwt, int *Posmut, char *AAmut, int Nmut,
		   int **clist, int **cnum, struct interaction *Int_list,
		   struct residue *seq, int Nres, atom *atoms, int *resatom,
		   struct Reference Ref_kin,
		   int *kin_atom, struct Normal_Mode NM,
		   float *sigma2, float **B3, float *B_pred, 
		   float *str_diff_3, float *str_diff_2)
{
  // Compute force
  *n_cont_mut=
    Mutation_force(Force, N_Cart, C_size, C_stab, C_dist, 
		   AAwt, Posmut, AAmut, Nmut, clist, cnum, Int_list,
		   seq, Nres, atoms, Ref_kin, kin_atom);

  // Compute deformation mut_Cart and mut_Tors
  Compute_mut_def(pred_mut_Cart, N_Cart, pred_mut_Tors, N_axes,
		  DE, Force_coeff, Force, NM, sigma2, Ref_kin.mass_sum);
      
  // Deformation of resatom mut_3, mut_abs2 and MSD_mut
  float pred_mut_ave[3], *str=pred_mut_3; int i, j;
  for(j=0; j<3; j++)pred_mut_ave[j]=0;
  for(i=0; i<Na; i++){
    float *x=pred_mut_Cart+3*kin_atom[resatom[i]];
    for(j=0; j<3; j++){*str=*x; pred_mut_ave[j]+=*x; x++; str++;}
  }
  // Subtract center of mass
  for(j=0; j<3; j++)pred_mut_ave[j]/=Na;
  *MSD_mut=0; str=pred_mut_3;
  for(i=0; i<Na; i++){
    double d2=0;
    for(j=0; j<3; j++){*str-=pred_mut_ave[j]; d2+=(*str)*(*str); str++;}
    *MSD_mut+=d2;
    pred_mut_2[i]=d2;
  }
  *MSD_mut/=Na;

  //for(i=0; i<20; i++){printf("%.3g ",pred_mut_2[i]);} printf("\n");
  //printf("Predicted MSD: %.4g\n", *MSD_mut);

  *MSD_cross=0;
  // No mutant structure is present for computing observed differences
  if((str_diff_3==NULL)||(str_diff_2==NULL))return(0);


  // pred_mut_cross[i]= 2 pred_mut[i] X B3[i]
  // B3_0i = sum_a s_a va_i (B_EVEC==0)
  // Does it make sense? va_i and -va_i are equivalent,
  // the dot product should be zero... 
  float pred_mut_cross[Na];
  for(i=0; i<Na; i++){pred_mut_cross[i]=0;}
  for(int a=0; a<amax; a++){ // amax=1 (B_EVEC==0) or 3 (B_EVEC==1)
    float *bb=B3[a]; str=pred_mut_3;
    // Scalar product pred_mut_3*B3
    for(i=0; i<Na; i++){
      double d_cross=0;
      for(j=0; j<3; j++){d_cross+=(*str)*(*bb); str++; bb++;}
      pred_mut_cross[i]+=d_cross*d_cross;
    }
  }
  (*MSD_cross)=0;
  for(i=0; i<Na; i++){
    pred_mut_cross[i]=4*sqrt(pred_mut_cross[i]); // Note the factor 4!
    (*MSD_cross)+=pred_mut_cross[i];
  }
  *MSD_cross/=Na;

  // Fit observed deformation profile squared str_diff_2[i] 
  float score;
  int Npar=3;
  float D_out[Npar], Y_pred[Na], *X[Na];
  for(i=0; i<Na; i++){
    X[i]=malloc(Npar*sizeof(float));
    X[i][0]=pred_mut_2[i];
    X[i][1]=B_pred[i];
    X[i][2]=pred_mut_cross[i];
  }

  struct ridge_fit fit; fit.A=malloc(Npar*sizeof(float)); char type;
  if(ini_fit){type='C';} // Specific-heat fit
  else{type='L'; fit.Lambda=Lambda_RRR;}
  score=Ridge_regression(&fit, Y_pred, D_out, "ridge_regression",
			 X, str_diff_2, Na, Npar, type);
  printf("RRR3: A_mut= %.3g A_nomut= %.3g A_cross= %.3g\n",
	 fit.A[0], fit.A[1], fit.A[2]);

  if(fit.A[0]<=0 || fit.A[1]<=0){ //*A_mut<=0 *A_nomut<0
    // Fit without A_cross
    fit.A[2]=0; Npar=2;
    score=Ridge_regression(&fit, Y_pred, D_out, "ridge_regression",
			   X, str_diff_2, Na, Npar, type);
    printf("RRR2: A_mut= %.3g A_nomut= %.3g A_cross= %.3g\n",
	   fit.A[0], fit.A[1], fit.A[2]);
  }

  Lambda_RRR=fit.Lambda;
  *A_mut=fit.A[0];
  *A_nomut=fit.A[1];
  *A_cross=fit.A[2];
  *offset=0;
  for(i=0; i<Na; i++){
    (*offset)+=str_diff_2[i]-Y_pred[i];
    pred_mut_tot_2[i]=Y_pred[i];
    free(X[i]);
  }
  (*offset)/=Na;
  free(fit.A);

  ini_fit=0;

  if(*A_nomut<=0){
    // Fit only with the predicted effect of mutation
    score=Fit_1(str_diff_2, pred_mut_2, Na, A_mut,offset);
    *A_nomut=0; *A_cross=0;
    printf("Fit1: A_mut= %.3g A_nomut= %.3g A_cross= %.3g\n",
	   *A_mut, *A_nomut, *A_cross);
  }else if(*A_mut<=0){
    // Fit only with thermal fluctuations
    score=Fit_1(str_diff_2, B_pred, Na, A_nomut,offset);
    *A_mut=0; *A_cross=0;
    printf("Fit1: A_mut= %.3g A_nomut= %.3g A_cross= %.3g\n",
	   *A_mut, *A_nomut, *A_cross);
  }

  return(score);
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x158, %rsp            # imm = 0x158
movq 0x1f0(%rbp), %rax
movq 0x1e8(%rbp), %rax
movq 0x1e0(%rbp), %rax
movq 0x1d8(%rbp), %rax
movq 0x1d0(%rbp), %rax
leaq 0x100(%rbp), %rax
movq %rax, -0x150(%rbp)
movq 0xf8(%rbp), %rax
leaq 0xc8(%rbp), %rax
movq %rax, -0x148(%rbp)
movq 0xc0(%rbp), %r10
movq 0xb8(%rbp), %r10
movl 0xb0(%rbp), %r10d
movq 0xa8(%rbp), %r10
movq 0xa0(%rbp), %r10
movq 0x98(%rbp), %r10
movq 0x90(%rbp), %r10
movl 0x88(%rbp), %r10d
movq 0x80(%rbp), %r10
movq 0x78(%rbp), %r10
movq 0x70(%rbp), %r10
movl 0x68(%rbp), %r10d
movq 0x60(%rbp), %r10
movl 0x58(%rbp), %r10d
movq 0x50(%rbp), %r10
movq 0x48(%rbp), %r10
movq 0x40(%rbp), %r10
movq 0x38(%rbp), %r10
movl 0x30(%rbp), %r10d
movq 0x28(%rbp), %r10
movq 0x20(%rbp), %r10
movq 0x18(%rbp), %r10
movq 0x10(%rbp), %r10
movss %xmm0, -0x30(%rbp)
movss %xmm1, -0x34(%rbp)
movss %xmm2, -0x38(%rbp)
movq %rdi, -0x40(%rbp)
movq %rsi, -0x48(%rbp)
movq %rdx, -0x50(%rbp)
movq %rcx, -0x58(%rbp)
movq %r8, -0x60(%rbp)
movq %r9, -0x68(%rbp)
movq 0x40(%rbp), %rdi
movl 0x58(%rbp), %esi
movss -0x30(%rbp), %xmm0
movss -0x34(%rbp), %xmm1
movss -0x38(%rbp), %xmm2
movq 0x70(%rbp), %rdx
movq 0x78(%rbp), %rcx
movq 0x80(%rbp), %r8
movl 0x88(%rbp), %r9d
movq 0x90(%rbp), %r10
movq 0x98(%rbp), %r11
movq 0xa0(%rbp), %rbx
movq 0xa8(%rbp), %r14
movl 0xb0(%rbp), %r15d
movq 0xb8(%rbp), %r12
movq 0xf8(%rbp), %r13
subq $0x70, %rsp
movups (%rax), %xmm3
movups 0x10(%rax), %xmm4
movups 0x20(%rax), %xmm5
movq %rsp, %rax
movq %rax, -0x158(%rbp)
movups %xmm5, 0x50(%rax)
movups %xmm4, 0x40(%rax)
movups %xmm3, 0x30(%rax)
movq %r13, 0x60(%rax)
movq %r12, 0x28(%rax)
movl %r15d, 0x20(%rax)
movq %r14, 0x18(%rax)
movq %rbx, 0x10(%rax)
movq %r11, 0x8(%rax)
movq %r10, (%rax)
callq 0x59810
movq -0x150(%rbp), %rsi
addq $0x70, %rsp
movl %eax, %edx
movq -0x148(%rbp), %rax
movq 0x38(%rbp), %rcx
movl %edx, (%rcx)
movq 0x50(%rbp), %rcx
movq %rcx, -0x138(%rbp)
movl 0x58(%rbp), %ecx
movl %ecx, -0x13c(%rbp)
movq 0x60(%rbp), %rdx
movl 0x68(%rbp), %ecx
movl %ecx, -0x130(%rbp)
movq -0x40(%rbp), %r8
movq 0x48(%rbp), %r9
movq 0x40(%rbp), %r10
movq 0x1d0(%rbp), %r11
movss 0x8(%rax), %xmm0
subq $0xe0, %rsp
movq %rsp, %rax
movq %rax, %rdi
addq $0x8, %rdi
movl $0x1a, %ecx
rep  movsq (%rsi), %es:(%rdi)
movl -0x13c(%rbp), %esi
movq -0x138(%rbp), %rdi
movl -0x130(%rbp), %ecx
movq %r11, 0xd8(%rax)
movq %r10, (%rax)
callq 0x59eb0
addq $0xe0, %rsp
movq 0x18(%rbp), %rax
movq %rax, -0x80(%rbp)
movl $0x0, -0x88(%rbp)
cmpl $0x3, -0x88(%rbp)
jge 0x58e6b
movslq -0x88(%rbp), %rax
xorps %xmm0, %xmm0
movss %xmm0, -0x74(%rbp,%rax,4)
movl -0x88(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x88(%rbp)
jmp 0x58e41
movl $0x0, -0x84(%rbp)
movl -0x84(%rbp), %eax
cmpl 0x30(%rbp), %eax
jge 0x58f3f
movq 0x50(%rbp), %rax
movq 0xf8(%rbp), %rcx
movq 0xc0(%rbp), %rdx
movslq -0x84(%rbp), %rsi
movslq (%rdx,%rsi,4), %rdx
imull $0x3, (%rcx,%rdx,4), %ecx
movslq %ecx, %rcx
shlq $0x2, %rcx
addq %rcx, %rax
movq %rax, -0x90(%rbp)
movl $0x0, -0x88(%rbp)
cmpl $0x3, -0x88(%rbp)
jge 0x58f29
movq -0x90(%rbp), %rax
movss (%rax), %xmm0
movq -0x80(%rbp), %rax
movss %xmm0, (%rax)
movq -0x90(%rbp), %rax
movss (%rax), %xmm0
movslq -0x88(%rbp), %rax
addss -0x74(%rbp,%rax,4), %xmm0
movss %xmm0, -0x74(%rbp,%rax,4)
movq -0x90(%rbp), %rax
addq $0x4, %rax
movq %rax, -0x90(%rbp)
movq -0x80(%rbp), %rax
addq $0x4, %rax
movq %rax, -0x80(%rbp)
movl -0x88(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x88(%rbp)
jmp 0x58ec0
jmp 0x58f2b
movl -0x84(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x84(%rbp)
jmp 0x58e75
movl $0x0, -0x88(%rbp)
cmpl $0x3, -0x88(%rbp)
jge 0x58f7f
cvtsi2ssl 0x30(%rbp), %xmm1
movslq -0x88(%rbp), %rax
movss -0x74(%rbp,%rax,4), %xmm0
divss %xmm1, %xmm0
movss %xmm0, -0x74(%rbp,%rax,4)
movl -0x88(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x88(%rbp)
jmp 0x58f49
movq -0x48(%rbp), %rax
xorps %xmm0, %xmm0
movss %xmm0, (%rax)
movq 0x18(%rbp), %rax
movq %rax, -0x80(%rbp)
movl $0x0, -0x84(%rbp)
movl -0x84(%rbp), %eax
cmpl 0x30(%rbp), %eax
jge 0x59077
xorps %xmm0, %xmm0
movsd %xmm0, -0x98(%rbp)
movl $0x0, -0x88(%rbp)
cmpl $0x3, -0x88(%rbp)
jge 0x59027
movslq -0x88(%rbp), %rax
movss -0x74(%rbp,%rax,4), %xmm1
movq -0x80(%rbp), %rax
movss (%rax), %xmm0
subss %xmm1, %xmm0
movss %xmm0, (%rax)
movq -0x80(%rbp), %rax
movss (%rax), %xmm0
movq -0x80(%rbp), %rax
mulss (%rax), %xmm0
cvtss2sd %xmm0, %xmm0
addsd -0x98(%rbp), %xmm0
movsd %xmm0, -0x98(%rbp)
movq -0x80(%rbp), %rax
addq $0x4, %rax
movq %rax, -0x80(%rbp)
movl -0x88(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x88(%rbp)
jmp 0x58fc0
movsd -0x98(%rbp), %xmm1
movq -0x48(%rbp), %rax
movss (%rax), %xmm0
cvtss2sd %xmm0, %xmm0
addsd %xmm1, %xmm0
cvtsd2ss %xmm0, %xmm0
movss %xmm0, (%rax)
movsd -0x98(%rbp), %xmm0
cvtsd2ss %xmm0, %xmm0
movq 0x20(%rbp), %rax
movslq -0x84(%rbp), %rcx
movss %xmm0, (%rax,%rcx,4)
movl -0x84(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x84(%rbp)
jmp 0x58f9c
cvtsi2ssl 0x30(%rbp), %xmm1
movq -0x48(%rbp), %rax
movss (%rax), %xmm0
divss %xmm1, %xmm0
movss %xmm0, (%rax)
movq -0x50(%rbp), %rax
xorps %xmm0, %xmm0
movss %xmm0, (%rax)
cmpq $0x0, 0x1e8(%rbp)
je 0x590ab
cmpq $0x0, 0x1f0(%rbp)
jne 0x590b8
xorps %xmm0, %xmm0
movss %xmm0, -0x2c(%rbp)
jmp 0x597ef
movl 0x30(%rbp), %eax
movq %rsp, %rcx
movq %rcx, -0xa0(%rbp)
leaq 0xf(,%rax,4), %rdx
andq $-0x10, %rdx
movq %rsp, %rcx
subq %rdx, %rcx
movq %rcx, -0x160(%rbp)
movq %rcx, %rsp
movq %rax, -0xa8(%rbp)
movl $0x0, -0x84(%rbp)
movl -0x84(%rbp), %eax
cmpl 0x30(%rbp), %eax
jge 0x59124
movq -0x160(%rbp), %rax
movslq -0x84(%rbp), %rcx
xorps %xmm0, %xmm0
movss %xmm0, (%rax,%rcx,4)
movl -0x84(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x84(%rbp)
jmp 0x590f2
movl $0x0, -0xac(%rbp)
movl -0xac(%rbp), %eax
cmpl 0x32dee(%rip), %eax     # 0x8bf28
jge 0x59250
movq 0x1d8(%rbp), %rax
movslq -0xac(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movq %rax, -0xb8(%rbp)
movq 0x18(%rbp), %rax
movq %rax, -0x80(%rbp)
movl $0x0, -0x84(%rbp)
movl -0x84(%rbp), %eax
cmpl 0x30(%rbp), %eax
jge 0x5923a
xorps %xmm0, %xmm0
movsd %xmm0, -0xc0(%rbp)
movl $0x0, -0x88(%rbp)
cmpl $0x3, -0x88(%rbp)
jge 0x591ee
movq -0x80(%rbp), %rax
movss (%rax), %xmm0
movq -0xb8(%rbp), %rax
mulss (%rax), %xmm0
cvtss2sd %xmm0, %xmm0
addsd -0xc0(%rbp), %xmm0
movsd %xmm0, -0xc0(%rbp)
movq -0x80(%rbp), %rax
addq $0x4, %rax
movq %rax, -0x80(%rbp)
movq -0xb8(%rbp), %rax
addq $0x4, %rax
movq %rax, -0xb8(%rbp)
movl -0x88(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x88(%rbp)
jmp 0x5918f
movq -0x160(%rbp), %rax
movsd -0xc0(%rbp), %xmm0
movsd -0xc0(%rbp), %xmm2
movslq -0x84(%rbp), %rcx
movss (%rax,%rcx,4), %xmm1
cvtss2sd %xmm1, %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
cvtsd2ss %xmm0, %xmm0
movss %xmm0, (%rax,%rcx,4)
movl -0x84(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x84(%rbp)
jmp 0x5916b
jmp 0x5923c
movl -0xac(%rbp), %eax
addl $0x1, %eax
movl %eax, -0xac(%rbp)
jmp 0x5912e
movq -0x50(%rbp), %rax
xorps %xmm0, %xmm0
movss %xmm0, (%rax)
movl $0x0, -0x84(%rbp)
movl -0x84(%rbp), %eax
cmpl 0x30(%rbp), %eax
jge 0x592db
movq -0x160(%rbp), %rax
movslq -0x84(%rbp), %rcx
movss (%rax,%rcx,4), %xmm0
cvtss2sd %xmm0, %xmm0
callq 0x1220
movq -0x160(%rbp), %rax
movaps %xmm0, %xmm1
movsd 0x22d8a(%rip), %xmm0    # 0x7c028
mulsd %xmm1, %xmm0
cvtsd2ss %xmm0, %xmm0
movslq -0x84(%rbp), %rcx
movss %xmm0, (%rax,%rcx,4)
movslq -0x84(%rbp), %rcx
movss (%rax,%rcx,4), %xmm0
movq -0x50(%rbp), %rax
addss (%rax), %xmm0
movss %xmm0, (%rax)
movl -0x84(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x84(%rbp)
jmp 0x59265
cvtsi2ssl 0x30(%rbp), %xmm1
movq -0x50(%rbp), %rax
movss (%rax), %xmm0
divss %xmm1, %xmm0
movss %xmm0, (%rax)
movl $0x3, -0xc8(%rbp)
movl -0xc8(%rbp), %eax
leaq 0xf(,%rax,4), %rdx
andq $-0x10, %rdx
movq %rsp, %rcx
subq %rdx, %rcx
movq %rcx, -0x178(%rbp)
movq %rcx, %rsp
movq %rax, -0xd0(%rbp)
movl 0x30(%rbp), %eax
leaq 0xf(,%rax,4), %rdx
andq $-0x10, %rdx
movq %rsp, %rcx
subq %rdx, %rcx
movq %rcx, -0x170(%rbp)
movq %rcx, %rsp
movq %rax, -0xd8(%rbp)
movl 0x30(%rbp), %eax
leaq 0xf(,%rax,8), %rdx
andq $-0x10, %rdx
movq %rsp, %rcx
subq %rdx, %rcx
movq %rcx, -0x168(%rbp)
movq %rcx, %rsp
movq %rax, -0xe0(%rbp)
movl $0x0, -0x84(%rbp)
movl -0x84(%rbp), %eax
cmpl 0x30(%rbp), %eax
jge 0x59426
movslq -0xc8(%rbp), %rdi
shlq $0x2, %rdi
callq 0x1180
movq -0x160(%rbp), %rcx
movq %rax, %rsi
movq -0x168(%rbp), %rax
movslq -0x84(%rbp), %rdx
movq %rsi, (%rax,%rdx,8)
movq 0x20(%rbp), %rdx
movslq -0x84(%rbp), %rsi
movss (%rdx,%rsi,4), %xmm0
movslq -0x84(%rbp), %rdx
movq (%rax,%rdx,8), %rdx
movss %xmm0, (%rdx)
movq 0x1e0(%rbp), %rdx
movslq -0x84(%rbp), %rsi
movss (%rdx,%rsi,4), %xmm0
movslq -0x84(%rbp), %rdx
movq (%rax,%rdx,8), %rdx
movss %xmm0, 0x4(%rdx)
movslq -0x84(%rbp), %rdx
movss (%rcx,%rdx,4), %xmm0
movslq -0x84(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movss %xmm0, 0x8(%rax)
movl -0x84(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x84(%rbp)
jmp 0x59379
movslq -0xc8(%rbp), %rdi
shlq $0x2, %rdi
callq 0x1180
movq %rax, -0x128(%rbp)
cmpl $0x0, 0x32ae0(%rip)     # 0x8bf24
je 0x5944f
movb $0x43, -0x129(%rbp)
jmp 0x59466
movb $0x4c, -0x129(%rbp)
movss 0x430f2(%rip), %xmm0    # 0x9c550
movss %xmm0, -0x11c(%rbp)
movq -0x168(%rbp), %r8
movq -0x178(%rbp), %rdx
movq -0x170(%rbp), %rsi
movq 0x1f0(%rbp), %r9
movl 0x30(%rbp), %r10d
movl -0xc8(%rbp), %eax
leaq -0x128(%rbp), %rdi
leaq 0x2850c(%rip), %rcx     # 0x819a6
subq $0x20, %rsp
movl %r10d, (%rsp)
movl %eax, 0x8(%rsp)
movsbl -0x129(%rbp), %eax
movl %eax, 0x10(%rsp)
callq 0x3d040
addq $0x20, %rsp
movss %xmm0, -0xc4(%rbp)
movq -0x128(%rbp), %rax
movss (%rax), %xmm0
cvtss2sd %xmm0, %xmm0
movq -0x128(%rbp), %rax
movss 0x4(%rax), %xmm1
cvtss2sd %xmm1, %xmm1
movq -0x128(%rbp), %rax
movss 0x8(%rax), %xmm2
cvtss2sd %xmm2, %xmm2
leaq 0x284bf(%rip), %rdi     # 0x819b7
movb $0x3, %al
callq 0x10d0
movq -0x128(%rbp), %rax
xorps %xmm0, %xmm0
ucomiss (%rax), %xmm0
jae 0x59522
movq -0x128(%rbp), %rax
xorps %xmm0, %xmm0
ucomiss 0x4(%rax), %xmm0
jb 0x595d4
movq -0x168(%rbp), %r8
movq -0x178(%rbp), %rdx
movq -0x170(%rbp), %rsi
movq -0x128(%rbp), %rax
xorps %xmm0, %xmm0
movss %xmm0, 0x8(%rax)
movl $0x2, -0xc8(%rbp)
movq 0x1f0(%rbp), %r9
movl 0x30(%rbp), %r10d
movl -0xc8(%rbp), %eax
leaq -0x128(%rbp), %rdi
leaq 0x28437(%rip), %rcx     # 0x819a6
subq $0x20, %rsp
movl %r10d, (%rsp)
movl %eax, 0x8(%rsp)
movsbl -0x129(%rbp), %eax
movl %eax, 0x10(%rsp)
callq 0x3d040
addq $0x20, %rsp
movss %xmm0, -0xc4(%rbp)
movq -0x128(%rbp), %rax
movss (%rax), %xmm0
cvtss2sd %xmm0, %xmm0
movq -0x128(%rbp), %rax
movss 0x4(%rax), %xmm1
cvtss2sd %xmm1, %xmm1
movq -0x128(%rbp), %rax
movss 0x8(%rax), %xmm2
cvtss2sd %xmm2, %xmm2
leaq 0x28419(%rip), %rdi     # 0x819e6
movb $0x3, %al
callq 0x10d0
movss -0x11c(%rbp), %xmm0
movss %xmm0, 0x42f6c(%rip)    # 0x9c550
movq -0x128(%rbp), %rax
movss (%rax), %xmm0
movq -0x58(%rbp), %rax
movss %xmm0, (%rax)
movq -0x128(%rbp), %rax
movss 0x4(%rax), %xmm0
movq -0x60(%rbp), %rax
movss %xmm0, (%rax)
movq -0x128(%rbp), %rax
movss 0x8(%rax), %xmm0
movq -0x68(%rbp), %rax
movss %xmm0, (%rax)
movq 0x10(%rbp), %rax
xorps %xmm0, %xmm0
movss %xmm0, (%rax)
movl $0x0, -0x84(%rbp)
movl -0x84(%rbp), %eax
cmpl 0x30(%rbp), %eax
jge 0x596b8
movq -0x168(%rbp), %rax
movq -0x170(%rbp), %rcx
movq 0x1f0(%rbp), %rdx
movslq -0x84(%rbp), %rsi
movss (%rdx,%rsi,4), %xmm0
movslq -0x84(%rbp), %rdx
subss (%rcx,%rdx,4), %xmm0
movq 0x10(%rbp), %rdx
addss (%rdx), %xmm0
movss %xmm0, (%rdx)
movslq -0x84(%rbp), %rdx
movss (%rcx,%rdx,4), %xmm0
movq 0x28(%rbp), %rcx
movslq -0x84(%rbp), %rdx
movss %xmm0, (%rcx,%rdx,4)
movslq -0x84(%rbp), %rcx
movq (%rax,%rcx,8), %rdi
callq 0x1040
movl -0x84(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x84(%rbp)
jmp 0x59634
cvtsi2ssl 0x30(%rbp), %xmm1
movq 0x10(%rbp), %rax
movss (%rax), %xmm0
divss %xmm1, %xmm0
movss %xmm0, (%rax)
movq -0x128(%rbp), %rdi
callq 0x1040
movl $0x0, 0x32841(%rip)     # 0x8bf24
movq -0x60(%rbp), %rax
xorps %xmm0, %xmm0
ucomiss (%rax), %xmm0
jb 0x5975c
movq 0x1f0(%rbp), %rdi
movq 0x20(%rbp), %rsi
movl 0x30(%rbp), %edx
movq -0x58(%rbp), %rcx
movq 0x10(%rbp), %r8
callq 0x5cb70
movss %xmm0, -0xc4(%rbp)
movq -0x60(%rbp), %rax
xorps %xmm0, %xmm0
movss %xmm0, (%rax)
movq -0x68(%rbp), %rax
xorps %xmm0, %xmm0
movss %xmm0, (%rax)
movq -0x58(%rbp), %rax
movss (%rax), %xmm0
cvtss2sd %xmm0, %xmm0
movq -0x60(%rbp), %rax
movss (%rax), %xmm1
cvtss2sd %xmm1, %xmm1
movq -0x68(%rbp), %rax
movss (%rax), %xmm2
cvtss2sd %xmm2, %xmm2
leaq 0x282c2(%rip), %rdi     # 0x81a15
movb $0x3, %al
callq 0x10d0
jmp 0x597d8
movq -0x58(%rbp), %rax
xorps %xmm0, %xmm0
ucomiss (%rax), %xmm0
jb 0x597d6
movq 0x1f0(%rbp), %rdi
movq 0x1e0(%rbp), %rsi
movl 0x30(%rbp), %edx
movq -0x60(%rbp), %rcx
movq 0x10(%rbp), %r8
callq 0x5cb70
movss %xmm0, -0xc4(%rbp)
movq -0x58(%rbp), %rax
xorps %xmm0, %xmm0
movss %xmm0, (%rax)
movq -0x68(%rbp), %rax
xorps %xmm0, %xmm0
movss %xmm0, (%rax)
movq -0x58(%rbp), %rax
movss (%rax), %xmm0
cvtss2sd %xmm0, %xmm0
movq -0x60(%rbp), %rax
movss (%rax), %xmm1
cvtss2sd %xmm1, %xmm1
movq -0x68(%rbp), %rax
movss (%rax), %xmm2
cvtss2sd %xmm2, %xmm2
leaq 0x28246(%rip), %rdi     # 0x81a15
movb $0x3, %al
callq 0x10d0
jmp 0x597d8
movss -0xc4(%rbp), %xmm0
movss %xmm0, -0x2c(%rbp)
movq -0xa0(%rbp), %rax
movq %rax, %rsp
movss -0x2c(%rbp), %xmm0
leaq -0x28(%rbp), %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
nopw %cs:(%rax,%rax)