LLM4Binary/decompile-bench-makefile huggingface dataset:
// File: /workspace/github_c_makefile_star2_O3/ugobas[P]tnm/dof_tnm.c
// Function: Axe_prot_main
char Axe_prot_main(struct bond *bond, int r_last)
{
  char *name=bond->atom->name, *previous_name=bond->previous->atom->name;
  if(strncmp(name,"N ", 2)==0){
    if((strncmp(previous_name,"C ", 2)==0)&&
       (bond->atom->res==(bond->previous->atom->res+1))){
      return('o'); // omega 
    }else{return('q');}         // Pseudo-axe after disordered loop
  }else if(strncmp(name,"CA", 2)==0){
    if((strncmp(previous_name,"N ", 2)==0)&&
       (bond->atom->res==bond->previous->atom->res)){
      return('f'); // phi
    }else{return('q');}
  }else if(strncmp(name,"C ", 2)==0){
    if((strncmp(previous_name,"CA", 2)==0)&&
       (bond->atom->res==bond->previous->atom->res)){
      if(bond->atom->res==r_last)return(NONE);
      return('p'); // psi
    }else{return('q');}
  }else{
    return('s');
  }
  return(NONE);
}

// Assembly:
movq (%rdi), %rcx
movq 0x28(%rdi), %rax
movq (%rax), %rdx
movzbl 0x18(%rcx), %edi
movb $0x73, %al
cmpl $0x43, %edi
je 0x13e77
cmpl $0x4e, %edi
jne 0x13e76
cmpb $0x20, 0x19(%rcx)
jne 0x13e76
cmpb $0x43, 0x18(%rdx)
jne 0x13e93
cmpb $0x20, 0x19(%rdx)
jne 0x13e93
movl 0x64(%rdx), %edx
incl %edx
movb $0x6f, %al
cmpl %edx, 0x64(%rcx)
jne 0x13e93
retq
cmpb $0x41, 0x19(%rcx)
jne 0x13e96
cmpb $0x4e, 0x18(%rdx)
jne 0x13e93
cmpb $0x20, 0x19(%rdx)
jne 0x13e93
movl 0x64(%rcx), %ecx
movb $0x66, %al
cmpl 0x64(%rdx), %ecx
je 0x13e76
movb $0x71, %al
retq
cmpb $0x20, 0x19(%rcx)
jne 0x13e76
movb $0x71, %al
cmpb $0x43, 0x18(%rdx)
jne 0x13e76
cmpb $0x41, 0x19(%rdx)
jne 0x13e76
movl 0x64(%rcx), %ecx
cmpl 0x64(%rdx), %ecx
jne 0x13e76
cmpl %esi, %ecx
movzbl 0x54969(%rip), %ecx     # 0x68824
movl $0x70, %eax
cmovel %ecx, %eax
retq
// File: /workspace/github_c_makefile_star2_O3/ugobas[P]tnm/dof_tnm.c
// Function: Axe_nuc_main
char Axe_nuc_main(struct bond *bond, int r_last)
{
  char *name=bond->atom->name, *previous_name=bond->previous->atom->name;
  if(strncmp(name,"O5'", 3)==0){
    if(strncmp(previous_name,"P ", 3)==0)return('1');
  }else if(strncmp(name,"C5'", 3)==0){
    if(strncmp(previous_name,"O5'", 3)==0)return('2');
  }else if(strncmp(name,"C4'", 3)==0){
    if(strncmp(previous_name,"C5'", 3)==0)return('3');
  }else if(strncmp(name,"C3'", 3)==0){
    if(strncmp(previous_name,"C4'", 3)==0)return('4');
  }else if(strncmp(name,"O3'", 3)==0){
    if(strncmp(previous_name,"C3'", 3)==0)return('5');
  }else if(strncmp(name,"P  ", 3)==0){
    if(strncmp(previous_name,"O3'", 3)==0)return('6');
    /*}else if(strncmp(name,"C1'", 3)==0){
      if(strncmp(previous_name,"N", 1)==0)return('s');*/ // side chain!
  }else{
    return(NONE);
  }
  return(NONE);
  //return('q'); // Disordered main axis
}

// Assembly:
movq (%rdi), %rax
movq 0x28(%rdi), %rcx
movq (%rcx), %rcx
movzbl 0x18(%rax), %edx
cmpl $0x43, %edx
je 0x13f3a
cmpl $0x50, %edx
je 0x13f16
cmpl $0x4f, %edx
jne 0x13fae
cmpb $0x35, 0x19(%rax)
jne 0x13f72
cmpb $0x27, 0x1a(%rax)
jne 0x13f72
cmpb $0x50, 0x18(%rcx)
jne 0x13fae
cmpb $0x20, 0x19(%rcx)
jne 0x13fae
cmpb $0x0, 0x1a(%rcx)
jne 0x13fae
movb $0x31, %al
retq
cmpb $0x20, 0x19(%rax)
jne 0x13fae
cmpb $0x20, 0x1a(%rax)
jne 0x13fae
cmpb $0x4f, 0x18(%rcx)
jne 0x13fae
cmpb $0x33, 0x19(%rcx)
jne 0x13fae
movb $0x36, %al
jmp 0x13fa8
cmpb $0x35, 0x19(%rax)
jne 0x13f56
cmpb $0x27, 0x1a(%rax)
jne 0x13f56
cmpb $0x4f, 0x18(%rcx)
jne 0x13fae
cmpb $0x35, 0x19(%rcx)
jne 0x13fae
movb $0x32, %al
jmp 0x13fa8
cmpb $0x34, 0x19(%rax)
jne 0x13f8e
cmpb $0x27, 0x1a(%rax)
jne 0x13f8e
cmpb $0x43, 0x18(%rcx)
jne 0x13fae
cmpb $0x35, 0x19(%rcx)
jne 0x13fae
movb $0x33, %al
jmp 0x13fa8
cmpb $0x33, 0x19(%rax)
jne 0x13fae
cmpb $0x27, 0x1a(%rax)
jne 0x13fae
cmpb $0x43, 0x18(%rcx)
jne 0x13fae
cmpb $0x33, 0x19(%rcx)
jne 0x13fae
movb $0x35, %al
jmp 0x13fa8
cmpb $0x33, 0x19(%rax)
jne 0x13fae
cmpb $0x27, 0x1a(%rax)
jne 0x13fae
cmpb $0x43, 0x18(%rcx)
jne 0x13fae
cmpb $0x34, 0x19(%rcx)
jne 0x13fae
movb $0x34, %al
cmpb $0x27, 0x1a(%rcx)
je 0x13fb4
movb 0x54870(%rip), %al      # 0x68824
retq
// File: /workspace/github_c_makefile_star2_O3/ugobas[P]tnm/dof_tnm.c
// Function: Count_interactions
int Count_interactions(atom *atoms, int natoms,
		       struct interaction *Int_list, int N_int)
{
  int i;
  for(i=0; i<natoms; i++)atoms[i].N_int=0;
  for(i=0; i<N_int; i++){
    atoms[Int_list[i].i1].N_int++;
    atoms[Int_list[i].i2].N_int++;
  }
  return i;
}

// Assembly:
movl %ecx, %eax
testl %esi, %esi
jle 0x1437e
movl %esi, %ecx
addl $0x3, %esi
decq %rcx
movq %rcx, %xmm0
pshufd $0x44, %xmm0, %xmm0     # xmm0 = xmm0[0,1,0,1]
shrl $0x2, %esi
imulq $0x1c0, %rsi, %rcx      # imm = 0x1C0
movdqa 0x3cd81(%rip), %xmm1    # 0x51010
movdqa 0x3cd89(%rip), %xmm2    # 0x51020
xorl %esi, %esi
movdqa 0x3cd8f(%rip), %xmm3    # 0x51030
pxor %xmm3, %xmm0
pcmpeqd %xmm4, %xmm4
movdqa 0x41ecf(%rip), %xmm5    # 0x56180
movdqa %xmm2, %xmm6
pxor %xmm3, %xmm6
movdqa %xmm6, %xmm7
pcmpgtd %xmm0, %xmm7
pcmpeqd %xmm0, %xmm6
pshufd $0xf5, %xmm6, %xmm8     # xmm8 = xmm6[1,1,3,3]
pand %xmm7, %xmm8
pshufd $0xf5, %xmm7, %xmm6     # xmm6 = xmm7[1,1,3,3]
por %xmm8, %xmm6
movd %xmm6, %r8d
notl %r8d
testb $0x1, %r8b
je 0x142f0
movl $0x0, 0x60(%rdi,%rsi)
pxor %xmm4, %xmm6
pextrw $0x4, %xmm6, %r8d
testb $0x1, %r8b
je 0x1430b
movl $0x0, 0xd0(%rdi,%rsi)
movdqa %xmm1, %xmm6
pxor %xmm3, %xmm6
movdqa %xmm6, %xmm7
pcmpgtd %xmm0, %xmm7
pcmpeqd %xmm0, %xmm6
pshufd $0xf5, %xmm6, %xmm8     # xmm8 = xmm6[1,1,3,3]
pand %xmm7, %xmm8
pshufd $0xf5, %xmm7, %xmm6     # xmm6 = xmm7[1,1,3,3]
por %xmm8, %xmm6
pxor %xmm4, %xmm6
pextrw $0x0, %xmm6, %r8d
testb $0x1, %r8b
je 0x1434f
movl $0x0, 0x140(%rdi,%rsi)
pextrw $0x4, %xmm6, %r8d
testb $0x1, %r8b
je 0x14366
movl $0x0, 0x1b0(%rdi,%rsi)
paddq %xmm5, %xmm2
paddq %xmm5, %xmm1
addq $0x1c0, %rsi            # imm = 0x1C0
cmpq %rsi, %rcx
jne 0x142b1
testl %eax, %eax
jle 0x143b5
movl %eax, %ecx
imulq $0xc8, %rcx, %rcx
xorl %esi, %esi
movslq (%rdx,%rsi), %r8
imulq $0x70, %r8, %r8
incl 0x60(%rdi,%r8)
movslq 0x4(%rdx,%rsi), %r8
imulq $0x70, %r8, %r8
incl 0x60(%rdi,%r8)
addq $0xc8, %rsi
cmpq %rsi, %rcx
jne 0x1438d
retq
xorl %eax, %eax
retq
// File: /workspace/github_c_makefile_star2_O3/ugobas[P]tnm/dof_tnm.c
// Function: Select_axes_int
int Select_axes_int(int *skipped, struct axe *axes, int *N_axes, int start,
		    atom *atoms, int natoms, int MIN_INT)
{
  int n_skipped=0;
  int n=start, last=*N_axes-1, last_a, a, i;
  for(a=start; a<*N_axes; a++){
    struct axe *axe=axes+a;
    if(a<last){
      struct axe *next=axe+1;
      if(((next->type=='l')||(axe->type=='l')||(axe->type=='a'))&&
	 ((a+1)<last)){
	next++;
      }
      last_a=next->first_atom;
    }else{
      last_a=natoms;
    }
    int num_int=0;
    for(i=axe->first_atom; i< last_a; i++)num_int+=atoms[i].N_int;
    if((num_int >= MIN_INT)||(axe->rigid)){ //
      if(a!=n)axes[n]=*axe;
      n++;
    }else{
      printf("%c%d ", axe->type, a);
      skipped[n_skipped]=a;
      n_skipped++;
    }
  }
  *N_axes=n;
  return(n_skipped);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x28, %rsp
movl %r9d, 0xc(%rsp)
movl %ecx, %ebx
movq %rdx, %r14
movq %rdi, 0x10(%rsp)
movslq (%rdx), %rax
xorl %r13d, %r13d
cmpl %ecx, %eax
jle 0x144cb
movq %rsi, %r12
leaq -0x1(%rax), %rbp
movslq %ebx, %r15
addq $-0x2, %rax
movq %rax, 0x18(%rsp)
addq $0x60, %r8
movq %r8, 0x20(%rsp)
xorl %r13d, %r13d
imulq $0x158, %r15, %rsi      # imm = 0x158
addq %r12, %rsi
movl 0xc(%rsp), %ecx
cmpq %rbp, %r15
jge 0x14445
leaq 0x158(%rsi), %rax
cmpb $0x6c, 0x158(%rsi)
je 0x14432
movzbl (%rsi), %ecx
cmpl $0x6c, %ecx
je 0x14432
cmpl $0x61, %ecx
jne 0x14442
leaq 0x2b0(%rsi), %rcx
cmpq 0x18(%rsp), %r15
cmovlq %rcx, %rax
movl 0x10(%rax), %ecx
movslq 0x10(%rsi), %rdi
xorl %eax, %eax
cmpl %ecx, %edi
jge 0x1446b
movslq %ecx, %rcx
imulq $0x70, %rdi, %rdx
addq 0x20(%rsp), %rdx
subq %rdi, %rcx
xorl %eax, %eax
addl (%rdx), %eax
addq $0x70, %rdx
decq %rcx
jne 0x14460
cmpl 0x60(%rsp), %eax
jge 0x14477
cmpl $0x0, 0x4(%rsi)
je 0x144a6
cmpl %r15d, %ebx
je 0x14493
movslq %ebx, %rax
imulq $0x158, %rax, %rdi      # imm = 0x158
addq %r12, %rdi
movl $0x158, %edx            # imm = 0x158
callq 0x21b0
incl %ebx
incq %r15
movslq (%r14), %rax
cmpq %rax, %r15
jl 0x14402
jmp 0x144cb
movsbl (%rsi), %esi
leaq 0x423ed(%rip), %rdi     # 0x5689d
movl %r15d, %edx
xorl %eax, %eax
callq 0x20f0
movslq %r13d, %r13
movq 0x10(%rsp), %rax
movl %r15d, (%rax,%r13,4)
incl %r13d
jmp 0x14495
movl %ebx, (%r14)
movl %r13d, %eax
addq $0x28, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_c_makefile_star2_O3/ugobas[P]tnm/dof_tnm.c
// Function: Last_axe
int Last_axe(int *aa, struct axe *axes, int ini, int n, int iatom)
{
  int a=*aa; struct axe *axe=axes+(*aa);
  while(a<n){
    if(axe->first_atom>iatom)break;
    a++; axe++;
  }
  if(a==0)return(-1);
  if(a==n){a--; axe--;}
  if(axe->first_atom>iatom && (axe-1)->first_atom<=iatom){a--; axe--;}
  if(a<ini){*aa=a;}else{*aa=ini;}
  if((axe->first_atom<=iatom))return(a); //&&(axe->last_atom>=i)
  printf("WARNING, last_axe<0: atom %d axe %d n= %d first %d last %d\n",
	 iatom, a, n, axe->first_atom, axe->last_atom); //exit(8);
  return(-1);
}

// Assembly:
pushq %rbx
movl %edx, %eax
movslq (%rdi), %rdx
imulq $0x158, %rdx, %r9       # imm = 0x158
addq %rsi, %r9
cmpl %ecx, %edx
jge 0x14509
cmpl %r8d, 0x10(%r9)
jg 0x14509
incl %edx
addq $0x158, %r9             # imm = 0x158
cmpl %edx, %ecx
jne 0x144f4
movl %ecx, %edx
movl $0xffffffff, %ebx       # imm = 0xFFFFFFFF
testl %edx, %edx
je 0x1456c
leaq -0x158(%r9), %rsi
xorl %r10d, %r10d
cmpl %ecx, %edx
cmovneq %r9, %rsi
sete %r10b
subl %r10d, %edx
cmpl %r8d, 0x10(%rsi)
jle 0x14541
cmpl %r8d, -0x148(%rsi)
jg 0x14541
addq $-0x158, %rsi           # imm = 0xFEA8
decl %edx
cmpl %eax, %edx
cmovll %edx, %eax
movl %eax, (%rdi)
movl 0x10(%rsi), %eax
cmpl %r8d, %eax
jle 0x1456a
movl 0x14(%rsi), %r9d
leaq 0x42348(%rip), %rdi     # 0x568a3
movl %r8d, %esi
movl %eax, %r8d
xorl %eax, %eax
callq 0x20f0
jmp 0x1456c
movl %edx, %ebx
movl %ebx, %eax
popq %rbx
retq
// File: /workspace/github_c_makefile_star2_O3/ugobas[P]tnm/buildup.c
// Function: Trajectory
void Trajectory(atom *atoms, int natoms, struct axe *axes, int naxe,
		float **Tors_mode, int imode, float fact, int nmove,
		struct bond *bonds)
{
  int iter, i, iatom, j;

  Set_bonds_measure(bonds, natoms, atoms);
  double d_phi[naxe];
  for(i=0; i<naxe; i++)d_phi[i]=Tors_mode[imode][i]*fact;

  atom *atommove=malloc(natoms*sizeof(atom));
  for(iter=0; iter<nmove; iter++){
    Build_up(bonds, natoms, d_phi, naxe);
    double r2=0; float x;
    for(i=0; i<natoms; i++){
      struct bond *bond=bonds+i;
      if(bond->i_atom<0)break;
      iatom=bond->i_atom;
      for(j=0; j<3; j++){
	atommove[i].r[j]=bonds[i].r[j];
	x=atommove[i].r[j]-atoms[iatom].r[j];
	r2+=x*x;
      }
      // printf("%3d %s %.4f\n", i, bond->name, x);
    }
    printf("mode %2d step %3d d2=%.4f\n", imode, iter, r2);
  }
  free(atommove);
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x18, %rsp
movss %xmm0, -0x38(%rbp)
movl %r9d, -0x30(%rbp)
movq %r8, %r14
movl %ecx, %ebx
movq %rdi, %r15
movq 0x18(%rbp), %r12
movq %r12, %rdi
movl %esi, -0x2c(%rbp)
movq %r15, %rdx
callq 0x14cf3
movss -0x38(%rbp), %xmm1
movl %ebx, %eax
movq %rsp, %rdx
leaq 0xf(,%rax,8), %rcx
andq $-0x10, %rcx
subq %rcx, %rdx
movq %rdx, %rsp
testl %ebx, %ebx
jle 0x14c3c
movslq -0x30(%rbp), %rcx
movq (%r14,%rcx,8), %rcx
xorl %esi, %esi
movss (%rcx,%rsi,4), %xmm0
mulss %xmm1, %xmm0
cvtss2sd %xmm0, %xmm0
movsd %xmm0, (%rdx,%rsi,8)
incq %rsi
cmpq %rsi, %rax
jne 0x14c22
cmpl $0x0, 0x10(%rbp)
jle 0x14ce4
movl -0x2c(%rbp), %r14d
leaq 0x8(%r12), %rax
movq %rax, -0x38(%rbp)
xorl %ebx, %ebx
movq %rdx, -0x40(%rbp)
movl -0x2c(%rbp), %r13d
movq %r12, %rdi
movl %r13d, %esi
callq 0x15083
xorps %xmm0, %xmm0
testl %r13d, %r13d
jle 0x14cc2
movq -0x38(%rbp), %rax
xorl %ecx, %ecx
leaq (%rcx,%rcx,2), %rdx
shlq $0x6, %rdx
movslq 0x20(%r12,%rdx), %rdx
testq %rdx, %rdx
js 0x14cc2
imulq $0x70, %rdx, %rdx
addq %r15, %rdx
xorl %esi, %esi
movsd (%rax,%rsi,8), %xmm1
subsd (%rdx,%rsi,8), %xmm1
cvtsd2ss %xmm1, %xmm1
mulss %xmm1, %xmm1
cvtss2sd %xmm1, %xmm1
addsd %xmm1, %xmm0
incq %rsi
cmpq $0x3, %rsi
jne 0x14c91
incq %rcx
addq $0xc0, %rax
cmpq %r14, %rcx
jne 0x14c76
leaq 0x41cd3(%rip), %rdi     # 0x5699c
movl -0x30(%rbp), %esi
movl %ebx, %edx
movb $0x1, %al
callq 0x20f0
incl %ebx
cmpl 0x10(%rbp), %ebx
movq -0x40(%rbp), %rdx
jne 0x14c5d
leaq -0x28(%rbp), %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_c_makefile_star2_O3/ugobas[P]tnm/buildup.c
// Function: Set_bonds_topology
struct bond *Set_bonds_topology(int N_atoms, atom *atoms, struct residue *seq)
{
  struct bond *bonds=malloc((N_atoms+1)*sizeof(struct bond));
  BONDTHR1_2=BONDTHR1*BONDTHR1;

  if(DEBUG)printf("Setting bond topology\n");
  int i;
  char chain=atoms[0].chain;
  atom *atom=atoms;
  struct bond *bond=bonds, *previous;
  for(i=0; i<N_atoms; i++){
    //ALIGN
    if(atom->ali<0)continue; // for 2 str, not aligned atoms are omitted
    atom->main=Is_atom_main(atom->name, seq->type);
    bond->i_atom=i;
    bond->atom=atom;
    for(int j=0; j<3; j++)bond->r[j]=atom->r[j];
    bond->terminal=1;
    bond->d_phi=0; bond->len=-1; bond->angle=-1; bond->torsion=-1;
    // The first atom of a chain is joined to the third one
    // of the first chain (n=3)
    if(atom->chain!=chain){
      chain=atom->chain;
      previous=bonds+2;
    }else if(i){
      previous=Find_previous(bond, i, N_atoms, seq+bond->atom->res);
      if(previous==NULL){
	printf("ERROR, atom %d %s res %d has no previous atom\n",
	       i, atom->name, atom->res); exit(8);
      }
    }else{previous=NULL;}
    bond->previous=previous;
    if(previous)previous->terminal=0;
    bond++; atom++;
  }
  bonds[i].atom=NULL; // Last bond
  printf("%d atoms, %d bonds for building the molecule\n", N_atoms, i);
  return(bonds);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x18, %rsp
movq %rdx, %r15
movq %rsi, %r12
movl %edi, %ebx
movslq %edi, %r14
movq %r14, %rax
shlq $0x6, %rax
leaq (%rax,%rax,2), %rdi
addq $0xc0, %rdi
callq 0x21e0
movq %rax, %r13
movss 0x531ac(%rip), %xmm0    # 0x6882c
mulss %xmm0, %xmm0
movss %xmm0, 0x63bec(%rip)    # 0x79278
testl %r14d, %r14d
jle 0x1576f
movb 0x6c(%r12), %r14b
movq %r13, %rax
addq $0x180, %rax            # imm = 0x180
movq %rax, 0x10(%rsp)
xorl %ebp, %ebp
movq %r13, 0x8(%rsp)
cmpl $0x0, 0x54(%r12)
js 0x15750
leaq 0x18(%r12), %rdi
movswl 0x18(%r15), %esi
callq 0x157a6
movl %eax, 0x48(%r12)
movl %ebp, 0x20(%r13)
movq %r12, (%r13)
movq 0x10(%r12), %rax
movq %rax, 0x18(%r13)
movups (%r12), %xmm0
movups %xmm0, 0x8(%r13)
movl $0x1, 0x30(%r13)
movw $0x0, 0x70(%r13)
movq $-0x1, 0x34(%r13)
movl $0xffffffff, 0x3c(%r13) # imm = 0xFFFFFFFF
movl 0x6c(%r12), %eax
movsbl %r14b, %ecx
cmpl %ecx, %eax
jne 0x15732
testl %ebp, %ebp
je 0x1575c
movslq 0x64(%r12), %rcx
shlq $0x5, %rcx
addq %r15, %rcx
movq %r13, %rdi
movl %ebp, %esi
callq 0x1584f
jmp 0x1573a
movl %eax, %r14d
movq 0x10(%rsp), %rax
movq %rax, 0x28(%r13)
movl $0x0, 0x30(%rax)
addq $0xc0, %r13
addq $0x70, %r12
incl %ebp
cmpl %ebp, %ebx
jne 0x156af
jmp 0x15766
movq $0x0, 0x28(%r13)
jmp 0x15745
movl %ebx, %edx
movq 0x8(%rsp), %r13
jmp 0x15771
xorl %edx, %edx
movl %edx, %eax
leaq (%rax,%rax,2), %rax
shlq $0x6, %rax
movq $0x0, (%r13,%rax)
leaq 0x4125b(%rip), %rdi     # 0x569e6
movl %ebx, %esi
xorl %eax, %eax
callq 0x20f0
movq %r13, %rax
addq $0x18, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_c_makefile_star2_O3/ugobas[P]tnm/interactions_tnm.c
// Function: Interactions
int Interactions(struct interaction *Int_list, char *atom_type, float thr,
		 atom *atoms, int N_atoms, int N_res)
{
  int N_int=0, i1=0, i2;
  atom *atom1=atoms, *atom2;
  float thr2=thr*thr;
  //float t_max2=t_max+22; t_max2*=t_max2;
  float d, d2;
  int res1, res2;

  for(res1=0; res1<N_res; res1++){

    if(Find_atom(atoms, &i1, res1, N_atoms, atom_type)< 0)continue;
    atom1=atoms+i1; i2=i1+1;

    for(res2=res1+1; res2<N_res; res2++){
      if(Find_atom(atoms, &i2, res2, N_atoms, atom_type)< 0)continue;
      atom2=atoms+i2;

      // Interatomic distance
      d=atom2->r[0]-atom1->r[0]; d2 =d*d; if(d2>thr2)continue;
      d=atom2->r[1]-atom1->r[1]; d2+=d*d; if(d2>thr2)continue;
      d=atom2->r[2]-atom1->r[2]; d2+=d*d; if(d2>thr2)continue;

      // store interaction
      Int_list[N_int].i1=i1; Int_list[N_int].i2=i2;
      Int_list[N_int].r0=sqrt(d2); Int_list[N_int].sec_der=1;
      N_int++;
      if(N_int >= INT_MAX){
	printf("ERROR, too many interactions %d\n", N_int); exit(8);
      }

    }
  }
  return(N_int);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x28, %rsp
movl %ecx, 0x8(%rsp)
movq %rsi, 0x18(%rsp)
movq %rdi, 0x20(%rsp)
movl $0x0, 0x4(%rsp)
testl %r8d, %r8d
jle 0x1717c
movl %r8d, %ebx
movq %rdx, %r14
mulss %xmm0, %xmm0
xorl %r12d, %r12d
leaq 0x4(%rsp), %rbp
xorl %r13d, %r13d
movss %xmm0, 0x14(%rsp)
movq %r14, %rdi
movq %rbp, %rsi
movl %r12d, %edx
movl 0x8(%rsp), %ecx
movq 0x18(%rsp), %r8
callq 0x1867d
testl %eax, %eax
js 0x1716e
movslq 0x4(%rsp), %r15
leal 0x1(%r15), %eax
movl %eax, 0xc(%rsp)
incl %r12d
cmpl %ebx, %r12d
jge 0x17171
imulq $0x70, %r15, %rbp
addq %r14, %rbp
movl %r12d, 0x10(%rsp)
movq %r14, %rdi
leaq 0xc(%rsp), %rsi
movl %r12d, %edx
movl 0x8(%rsp), %ecx
movq 0x18(%rsp), %r8
callq 0x1867d
testl %eax, %eax
js 0x17156
movslq 0xc(%rsp), %rax
imulq $0x70, %rax, %rcx
movsd (%r14,%rcx), %xmm0
subsd (%rbp), %xmm0
xorps %xmm2, %xmm2
cvtsd2ss %xmm0, %xmm2
mulss %xmm2, %xmm2
movss 0x14(%rsp), %xmm0
ucomiss %xmm0, %xmm2
ja 0x17156
addq %r14, %rcx
movsd 0x8(%rcx), %xmm1
subsd 0x8(%rbp), %xmm1
cvtsd2ss %xmm1, %xmm1
mulss %xmm1, %xmm1
addss %xmm2, %xmm1
ucomiss %xmm0, %xmm1
ja 0x17156
movsd 0x10(%rcx), %xmm2
subsd 0x10(%rbp), %xmm2
cvtsd2ss %xmm2, %xmm2
mulss %xmm2, %xmm2
addss %xmm1, %xmm2
ucomiss %xmm0, %xmm2
ja 0x17156
movslq %r13d, %r13
imulq $0xc8, %r13, %rcx
movq 0x20(%rsp), %rdx
movl %r15d, (%rdx,%rcx)
movl %eax, 0x4(%rdx,%rcx)
xorps %xmm1, %xmm1
cvtss2sd %xmm2, %xmm1
sqrtsd %xmm1, %xmm1
movsd %xmm1, 0x20(%rdx,%rcx)
movabsq $0x3ff0000000000000, %rax # imm = 0x3FF0000000000000
movq %rax, 0x28(%rdx,%rcx)
incl %r13d
leaq 0x62187(%rip), %rax     # 0x792d8
cmpl (%rax), %r13d
jge 0x17191
incl %r12d
cmpl %r12d, %ebx
jne 0x1708c
movl 0x10(%rsp), %r12d
leaq 0x4(%rsp), %rbp
jmp 0x17171
incl %r12d
cmpl %ebx, %r12d
jne 0x17048
jmp 0x1717f
xorl %r13d, %r13d
movl %r13d, %eax
addq $0x28, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
leaq 0x3fdeb(%rip), %rdi     # 0x56f83
movl %r13d, %esi
xorl %eax, %eax
callq 0x20f0
movl $0x8, %edi
callq 0x2280
// File: /workspace/github_c_makefile_star2_O3/ugobas[P]tnm/interactions_tnm.c
// Function: Purge_interactions_res_pair
void Purge_interactions_res_pair(int *N_int, struct interaction *Int_list,
				 int nres)
{
  // Maximum one interaction per residue pair
  int Nnew=0, res1, res2, k; 
  struct interaction *res12;
  struct interaction *ptr, *ini=Int_list, *end=Int_list+*N_int;
  struct interaction *tmp=malloc(*N_int*sizeof(struct interaction));

  for(res1=0; res1<nres; res1++){
    for(res2=res1+1; res2<nres; res2++){
      // Find strongest res1-res2 contact
      // Find the next value of res2
      int res2min=nres; res12=NULL;
      for(ptr=ini; ptr<end; ptr++){
	if(ptr->res1 > res1)break;
	if(ptr->res1 == res1){
	  if(ptr->res2 == res2){
	    if((res12==NULL)||(ptr->sec_der>res12->sec_der)||
	       ((ptr->sec_der==res12->sec_der)&&(ptr->r0 < res12->r0))){
	      res12=ptr;
	    }
	  }else if((ptr->res2 < res2min)&&(ptr->res2 > res2)){
	    res2min=ptr->res2;
	  }
	}
      } // End pair

      if(res12){tmp[Nnew]=*res12; Nnew++;} // store

      // Prepare next value
      if(res2min==nres){ini=ptr; break;} // end res2
      else{res2=res2min-1;}

    } // end res2
  } // end res1

  for(k=0; k<Nnew; k++)Int_list[k]=tmp[k];
  *N_int=Nnew;
  free(tmp);

}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x28, %rsp
movl %edx, %r14d
movq %rsi, %r15
movq %rdi, %rbx
movslq (%rdi), %rax
imulq $0xc8, %rax, %r12
movq %r12, %rdi
callq 0x21e0
testl %r14d, %r14d
jle 0x182df
movq %rbx, 0x10(%rsp)
addq %r15, %r12
xorl %r13d, %r13d
xorl %ecx, %ecx
movq %r15, 0x20(%rsp)
movq %r15, %r8
movq %rax, 0x18(%rsp)
movl %r14d, 0xc(%rsp)
movl %ecx, %r15d
leal 0x1(%r15), %edi
movl %edi, 0x8(%rsp)
cmpl %r14d, %edi
jge 0x182a4
cmpq %r12, %r8
jae 0x182a4
movl %r14d, %ebp
movq %r8, %rbx
xorl %esi, %esi
cmpl %r15d, 0x8(%rbx)
jg 0x18262
jne 0x18256
movl 0xc(%rbx), %ecx
cmpl %edi, %ecx
jne 0x1824a
testq %rsi, %rsi
je 0x18245
movsd 0x28(%rbx), %xmm0
ucomisd 0x28(%rsi), %xmm0
ja 0x18245
jne 0x18256
jp 0x18256
movsd 0x20(%rsi), %xmm0
ucomisd 0x20(%rbx), %xmm0
jbe 0x18256
movq %rbx, %rsi
jmp 0x18256
cmpl %ebp, %ecx
movl %ebp, %edx
cmovll %ecx, %edx
cmpl %edi, %ecx
cmovgl %edx, %ebp
addq $0xc8, %rbx
cmpq %r12, %rbx
jb 0x18215
testq %rsi, %rsi
je 0x18297
movslq %r13d, %r13
imulq $0xc8, %r13, %rdi
addq %rax, %rdi
movl $0xc8, %edx
movq %r13, %r14
movq %r8, %r13
callq 0x21b0
movq %r13, %r8
movq %r14, %r13
movl 0xc(%rsp), %r14d
movq 0x18(%rsp), %rax
incl %r13d
movl %ebp, %edi
cmpl %r14d, %ebp
jne 0x181fb
jmp 0x182a7
movq %r8, %rbx
movq %rbx, %r8
movl 0x8(%rsp), %ecx
cmpl %r14d, %ecx
jne 0x181f0
testl %r13d, %r13d
jle 0x182d8
movl %r13d, %ecx
imulq $0xc8, %rcx, %rdx
movq 0x20(%rsp), %rdi
movq %rax, %rsi
callq 0x21b0
movq 0x18(%rsp), %rax
movq 0x10(%rsp), %rbx
jmp 0x182e2
xorl %r13d, %r13d
movl %r13d, (%rbx)
movq %rax, %rdi
addq $0x28, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0x2050
// File: /workspace/github_c_makefile_star2_O3/ugobas[P]tnm/interactions_tnm.c
// Function: Find_atom
int Find_atom(atom *atoms, int *i1, int res, int N_atoms, char *atom_type){
  int i; atom *atom;

  atom=atoms+*i1;
  for(i=*i1; i<N_atoms; i++){
    if(atom->res > res)break;
    if(atom->res < res)goto next1;
    if(strncmp(atom->name, atom_type, 2)==0){*i1=i; return(0);}
  next1: atom++;
  }
  // Not found? Look for CA atom
  atom=atoms+*i1;
  for(i=*i1; i<N_atoms; i++){
    if(atom->res > res)break;
    if(atom->res < res)goto next2;
    if(strncmp(atom->name, "CA", 2)==0){*i1=i; return(0);}
  next2: atom++;
  }
  return(-1);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x18, %rsp
movl (%rsi), %r14d
movl $0xffffffff, %eax       # imm = 0xFFFFFFFF
cmpl %ecx, %r14d
jge 0x18733
movq %r8, %r13
movl %ecx, %ebp
movl %edx, %r15d
movq %rsi, (%rsp)
movslq %r14d, %rax
imulq $0x70, %rax, %rax
movq %rdi, 0x10(%rsp)
movq %rax, 0x8(%rsp)
leaq (%rax,%rdi), %r12
addq $0x64, %r12
movl %r14d, %ebx
cmpl %r15d, (%r12)
jg 0x186ec
jl 0x186e2
leaq -0x4c(%r12), %rdi
movl $0x2, %edx
movq %r13, %rsi
callq 0x2070
testl %eax, %eax
je 0x1872b
incl %ebx
addq $0x70, %r12
cmpl %ebx, %ebp
jne 0x186c4
cmpl %ebp, %r14d
jge 0x18721
movq 0x10(%rsp), %rax
movq 0x8(%rsp), %rcx
addq %rcx, %rax
addq $0x18, %rax
cmpl %r15d, 0x4c(%rax)
jg 0x18721
jl 0x18715
cmpb $0x43, (%rax)
jne 0x18715
cmpb $0x41, 0x1(%rax)
je 0x18728
incl %r14d
addq $0x70, %rax
cmpl %r14d, %ebp
jne 0x18702
movl $0xffffffff, %eax       # imm = 0xFFFFFFFF
jmp 0x18733
movl %r14d, %ebx
movq (%rsp), %rax
movl %ebx, (%rax)
xorl %eax, %eax
addq $0x18, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_c_makefile_star2_O3/ugobas[P]tnm/interactions_tnm.c
// Function: Covalent
int Covalent(int res1, int res2, atom *atom1, atom *atom2){
  if(res2!=(res1+1))return(0);
  if((strncmp(atom2->name, "N ", 2)==0)&&
       ((strncmp(atom1->name, "O ", 2)==0)
	||(strncmp(atom1->name, "CA", 2)==0)
	||(strncmp(atom1->name, "C ", 2)==0)
	))return(1);
  if((strncmp(atom2->name, "CA", 2)==0)&&
     (strncmp(atom1->name, "C ", 2)==0)
     //||(strncmp(atom1->name, "O ", 2)==0)
     )return(1);
  return(0);
}

// Assembly:
incl %edi
xorl %eax, %eax
cmpl %esi, %edi
jne 0x18787
movzbl 0x18(%rcx), %eax
cmpl $0x43, %eax
je 0x1876e
cmpl $0x4e, %eax
jne 0x18785
cmpb $0x20, 0x19(%rcx)
jne 0x18785
movzbl 0x18(%rdx), %eax
cmpl $0x43, %eax
je 0x18788
cmpl $0x4f, %eax
je 0x1877a
jmp 0x18785
cmpb $0x41, 0x19(%rcx)
jne 0x18785
cmpb $0x43, 0x18(%rdx)
jne 0x18785
movl $0x1, %eax
cmpb $0x20, 0x19(%rdx)
je 0x18787
xorl %eax, %eax
retq
movl $0x1, %eax
cmpb $0x41, 0x19(%rdx)
jne 0x1877f
jmp 0x18787
// File: /workspace/github_c_makefile_star2_O3/ugobas[P]tnm/interactions_tnm.c
// Function: Heavy_atoms_interactors
void Heavy_atoms_interactors(struct interaction *Int_list, int N_int,
			     atom *atoms)
{
  int k, h=0;
  for(k=0; k<N_int; k++){
    h+=Heavy_atom(&(Int_list[k].i1), atoms);
    h+=Heavy_atom(&(Int_list[k].i2), atoms);
    if(1){
      Int_list[k].r0=
	sqrt(Distance_square(atoms[Int_list[k].i1].r,atoms[Int_list[k].i2].r));
    }

  }
  printf("%d hydrogen atoms changed to heavy atoms\n", h);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
pushq %rax
testl %esi, %esi
jle 0x18817
movq %rdx, %r14
movq %rdi, %r15
movl %esi, %r13d
xorl %ebx, %ebx
movq %r15, %rdi
movq %r14, %rsi
callq 0x18837
movl %eax, %ebp
leaq 0x4(%r15), %rdi
movq %r14, %rsi
callq 0x18837
movl %eax, %r12d
movslq (%r15), %rax
imulq $0x70, %rax, %rdi
addq %r14, %rdi
movslq 0x4(%r15), %rax
imulq $0x70, %rax, %rsi
addq %r14, %rsi
callq 0x141ea
cvtss2sd %xmm0, %xmm0
ucomisd 0x38987(%rip), %xmm0    # 0x51178
jb 0x187f9
sqrtsd %xmm0, %xmm0
jmp 0x187fe
callq 0x22a0
addl %ebp, %ebx
addl %r12d, %ebx
movsd %xmm0, 0x20(%r15)
addq $0xc8, %r15
decq %r13
jne 0x187af
jmp 0x18819
xorl %ebx, %ebx
leaq 0x3e784(%rip), %rdi     # 0x56fa4
movl %ebx, %esi
xorl %eax, %eax
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0x20f0
// File: /workspace/github_c_makefile_star2_O3/ugobas[P]tnm/interactions_tnm.c
// Function: Heavy_atom
int Heavy_atom(int *iatom, atom *atoms)
{
  if(atoms[*iatom].name[0]!='H')return(0); // Not hydrogen atoms
  atom *atom1=atoms+*iatom;
  int i, imin=*iatom, res=atom1->res;
  double dmin=100, *r=atom1->r;
  char branch=atom1->name[1];
  if((branch=='N')||(branch=='1')||(branch=='2')||(branch=='3'))branch=' ';
  for(i=*iatom-1; i>=0; i--){
    atom1--; if(atom1->res!=res)break;
    if((atom1->name[0]!='H')&&(atom1->name[1]==branch)){
      double x=r[0]-atom1->r[0], d=x*x;
      x=r[1]-atom1->r[1]; d+=x*x;
      x=r[2]-atom1->r[2]; d+=x*x;
      if(d<dmin){dmin=d; imin=i;}
    }
  }

  if(*iatom==imin){
    printf("WARNING, atom %s %d, no heavy atom found\n",
	   atoms[*iatom].name, *iatom); return(0);
  }else if(dmin > 4){
    printf("WARNING, atom %s %d, heavy atom %s %d has d=%.3f\n",
	   atoms[*iatom].name, *iatom, atoms[imin].name, imin, dmin);
  }
  *iatom=imin; return(1);
}

// Assembly:
pushq %rbp
pushq %rbx
pushq %rax
movl (%rdi), %edx
movslq %edx, %rax
imulq $0x70, %rax, %rax
xorl %ebp, %ebp
cmpb $0x48, 0x18(%rsi,%rax)
jne 0x1895d
testl %edx, %edx
jle 0x1890a
movq %rdi, %rbx
leaq (%rsi,%rax), %rcx
movsbl 0x19(%rcx), %r9d
leal -0x31(%r9), %edi
cmpb $0x3, %dil
movl $0x20, %r10d
movl %r9d, %edi
cmovbl %r10d, %edi
movl 0x64(%rcx), %r8d
cmpl $0x4e, %r9d
cmovel %r10d, %edi
leal 0x1(%rdx), %r9d
leaq (%rax,%rsi), %r10
addq $-0x70, %r10
movsd 0x3e647(%rip), %xmm0    # 0x56ee0
movl %edx, %ebp
cmpl %r8d, 0x64(%r10)
jne 0x18906
cmpb $0x48, 0x18(%r10)
je 0x188f9
movsbl 0x19(%r10), %r11d
cmpl %r11d, %edi
jne 0x188f9
movsd (%rcx), %xmm1
movsd 0x8(%rcx), %xmm2
subsd (%r10), %xmm1
mulsd %xmm1, %xmm1
subsd 0x8(%r10), %xmm2
mulsd %xmm2, %xmm2
addsd %xmm1, %xmm2
movsd 0x10(%rcx), %xmm1
subsd 0x10(%r10), %xmm1
mulsd %xmm1, %xmm1
addsd %xmm2, %xmm1
ucomisd %xmm1, %xmm0
leal -0x2(%r9), %r11d
cmoval %r11d, %ebp
minsd %xmm0, %xmm1
movapd %xmm1, %xmm0
decl %r9d
addq $-0x70, %r10
cmpl $0x1, %r9d
ja 0x1889b
cmpl %ebp, %edx
jne 0x18923
addq %rax, %rsi
addq $0x18, %rsi
leaq 0x3e6b6(%rip), %rdi     # 0x56fce
xorl %ebp, %ebp
xorl %eax, %eax
callq 0x20f0
jmp 0x1895d
ucomisd 0x3e5bd(%rip), %xmm0    # 0x56ee8
jbe 0x18956
addq %rsi, %rax
addq $0x18, %rax
movslq %ebp, %rcx
imulq $0x70, %rcx, %rcx
addq %rsi, %rcx
addq $0x18, %rcx
leaq 0x3e6af(%rip), %rdi     # 0x56ff8
movq %rax, %rsi
movl %ebp, %r8d
movb $0x1, %al
callq 0x20f0
movl %ebp, (%rbx)
movl $0x1, %ebp
movl %ebp, %eax
addq $0x8, %rsp
popq %rbx
popq %rbp
retq
// File: /workspace/github_c_makefile_star2_O3/ugobas[P]tnm/interactions_tnm.c
// Function: Force_constant
float Force_constant(float r)
{
  if(r <= C1_THR){
    if(POW){return(K0*pow(r, -EXP_HESSIAN));}
    else{return(K0*exp(-r*EXP_HESSIAN));}
  }else if(r < C_THR){
    return(K1*(C_THR-r));
  }else{
    return(0);
  }
}

// Assembly:
pushq %rax
leaq 0x5017c(%rip), %rax     # 0x69070
movss (%rax), %xmm1
ucomiss %xmm0, %xmm1
jae 0x18f33
leaq 0x50164(%rip), %rax     # 0x69068
movss (%rax), %xmm2
xorps %xmm1, %xmm1
ucomiss %xmm0, %xmm2
jbe 0x18f94
leaq 0x5016d(%rip), %rax     # 0x69088
subss %xmm0, %xmm2
xorps %xmm0, %xmm0
cvtss2sd %xmm2, %xmm0
mulsd (%rax), %xmm0
xorps %xmm1, %xmm1
cvtsd2ss %xmm0, %xmm1
jmp 0x18f94
leaq 0x4f6b2(%rip), %rax     # 0x685ec
cmpl $0x0, (%rax)
leaq 0x5013c(%rip), %rax     # 0x69080
movsd (%rax), %xmm1
movsd %xmm1, (%rsp)
leaq 0x50118(%rip), %rax     # 0x6906c
movss (%rax), %xmm1
je 0x18f70
cvtss2sd %xmm0, %xmm0
xorps 0x380cb(%rip), %xmm1    # 0x51030
cvtss2sd %xmm1, %xmm1
callq 0x20a0
jmp 0x18f87
xorps 0x380b9(%rip), %xmm0    # 0x51030
mulss %xmm0, %xmm1
xorps %xmm0, %xmm0
cvtss2sd %xmm1, %xmm0
callq 0x22b0
movsd (%rsp), %xmm1
mulsd %xmm0, %xmm1
cvtsd2ss %xmm1, %xmm1
movaps %xmm1, %xmm0
popq %rax
retq
// File: /workspace/github_c_makefile_star2_O3/ugobas[P]tnm/interactions_tnm.c
// Function: Select_pairs
void Select_pairs(int *selected, struct interaction **res12, int n12)
{
  int k, kopt=0, nsel=0; float w=0;
  //selected: 0=not yet -1=discarded 1=selected
  for(k=0; k<n12; k++)selected[k]=0;
  while(kopt>=0){
    kopt=-1;
    for(k=0; k<n12; k++){
      if(selected[k])continue;
      if(res12[k]->sec_der==1){
	if((kopt<0)||(res12[k]->r0 < w)){
	  kopt=k; w=res12[k]->r0; 
	}
      }else{
	if((kopt<0)||(res12[k]->sec_der > w)){
	  kopt=k; w=res12[k]->sec_der;
	}
      } 
    }
    if(kopt<0)return;
    selected[kopt]=1; nsel++;
    int i1=res12[kopt]->i1, i2=res12[kopt]->i2;
    for(k=0; k<n12; k++){
      if(selected[k])continue;
      if((res12[k]->i1==i1)||(res12[k]->i2==i2))selected[k]=-1;
    }
  }
  return;
}

// Assembly:
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
pushq %rax
movl %edx, %ebx
movq %rsi, %r14
movq %rdi, %r15
movl %edx, %r12d
testl %edx, %edx
jle 0x194e8
leaq (,%r12,4), %rdx
movq %r15, %rdi
xorl %esi, %esi
callq 0x2110
testl %ebx, %ebx
jle 0x195a7
xorps %xmm0, %xmm0
movl $0xffffffff, %ecx       # imm = 0xFFFFFFFF
xorl %eax, %eax
movsd 0x37c46(%rip), %xmm1    # 0x51148
cmpl $0x0, (%r15,%rax,4)
jne 0x1954c
movq (%r14,%rax,8), %rdx
movsd 0x28(%rdx), %xmm2
ucomisd %xmm1, %xmm2
jne 0x19532
jp 0x19532
movsd 0x20(%rdx), %xmm2
testl %ecx, %ecx
js 0x19543
xorps %xmm3, %xmm3
cvtss2sd %xmm0, %xmm3
ucomisd %xmm2, %xmm3
jbe 0x1954c
jmp 0x19543
testl %ecx, %ecx
js 0x19543
xorps %xmm3, %xmm3
cvtss2sd %xmm0, %xmm3
ucomisd %xmm3, %xmm2
jbe 0x1954c
xorps %xmm0, %xmm0
cvtsd2ss %xmm2, %xmm0
movl %eax, %ecx
incq %rax
cmpq %r12, %rax
jne 0x19502
testl %ecx, %ecx
js 0x195a7
movl %ecx, %eax
movl $0x1, (%r15,%rax,4)
testl %ebx, %ebx
jle 0x195a7
movq (%r14,%rax,8), %rcx
movl (%rcx), %eax
movl 0x4(%rcx), %ecx
xorl %edx, %edx
cmpl $0x0, (%r15,%rdx,4)
jne 0x1958d
movq (%r14,%rdx,8), %rsi
cmpl %eax, (%rsi)
je 0x19585
cmpl %ecx, 0x4(%rsi)
jne 0x1958d
movl $0xffffffff, (%r15,%rdx,4) # imm = 0xFFFFFFFF
incq %rdx
cmpq %rdx, %r12
jne 0x19571
movl $0x0, %eax
movl $0xffffffff, %ecx       # imm = 0xFFFFFFFF
testl %ebx, %ebx
jg 0x19502
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r14
popq %r15
retq
// File: /workspace/github_c_makefile_star2_O3/ugobas[P]tnm/interactions_tnm.c
// Function: Purge_interactions_res_pair_new
void Purge_interactions_res_pair_new(int *N_int, struct interaction *Int_list,
				     int nres)
{
  // Maximum one interaction per residue pair
  int Nnew=0, res1, res2, k; 
  struct interaction *ini=Int_list, *end=Int_list+*N_int;
  struct interaction *tmp=malloc(*N_int*sizeof(struct interaction));

  for(res1=0; res1<nres; res1++){
    for(res2=res1+1; res2<nres; res2++){
      // Find strongest res1-res2 contact
      // Find the next value of res2
      int new1;
      struct interaction *res12=
	Check_pair(&new1, res1, &res2, nres, &ini, end);

      if(res12){tmp[Nnew]=*res12; Nnew++;} // store
      if(new1)break;

    } // end res2
  } // end res1

  for(k=0; k<Nnew; k++)Int_list[k]=tmp[k];
  *N_int=Nnew;
  free(tmp);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x28, %rsp
movl %edx, %ebp
movq %rsi, %r13
movq %rdi, %rbx
movq %rsi, 0x20(%rsp)
movslq (%rdi), %rax
imulq $0xc8, %rax, %r12
movq %r12, %rdi
callq 0x21e0
movq %rax, %r15
testl %ebp, %ebp
jle 0x19686
movq %rbx, 0x10(%rsp)
movq %r13, 0x18(%rsp)
addq %r13, %r12
xorl %r13d, %r13d
xorl %ebx, %ebx
movl %ebx, %r14d
leal 0x1(%r14), %ebx
movl %ebx, 0x8(%rsp)
cmpl %ebp, %ebx
jge 0x1965f
leaq 0xc(%rsp), %rdi
movl %r14d, %esi
leaq 0x8(%rsp), %rdx
movl %ebp, %ecx
leaq 0x20(%rsp), %r8
movq %r12, %r9
callq 0x196a2
testq %rax, %rax
je 0x1964a
movslq %r13d, %r13
imulq $0xc8, %r13, %rdi
addq %r15, %rdi
movl $0xc8, %edx
movq %rax, %rsi
callq 0x21b0
incl %r13d
cmpl $0x0, 0xc(%rsp)
jne 0x1965f
movl 0x8(%rsp), %eax
incl %eax
movl %eax, 0x8(%rsp)
cmpl %ebp, %eax
jl 0x1960c
cmpl %ebp, %ebx
jne 0x195fd
testl %r13d, %r13d
jle 0x1967f
movl %r13d, %eax
imulq $0xc8, %rax, %rdx
movq 0x18(%rsp), %rdi
movq %r15, %rsi
callq 0x21b0
movq 0x10(%rsp), %rbx
jmp 0x19689
xorl %r13d, %r13d
movl %r13d, (%rbx)
movq %r15, %rdi
addq $0x28, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0x2050
// File: /workspace/github_c_makefile_star2_O3/ugobas[P]tnm/interactions_tnm.c
// Function: Check_pair
struct interaction *Check_pair(int *new1, int res1, int *res2, int nres,
			       struct interaction **ini,
			       struct interaction *end)
{
  struct interaction *res12=NULL, *ptr;
  int res2min=nres; *new1=0;
  for(ptr=*ini; ptr<end; ptr++){
    if(ptr->res1 > res1)break;
    if(ptr->res1 == res1){
      if(ptr->res2 == *res2){
	if((res12==NULL)||(ptr->sec_der>res12->sec_der)||
	   ((ptr->sec_der==res12->sec_der)&&(ptr->r0 < res12->r0))){
	  res12=ptr;
	}
      }else if((ptr->res2 < res2min)&&(ptr->res2 > *res2)){
	res2min=ptr->res2;
      }
    }
  } // End pair
  // Prepare next value
  if(res2min==nres){*ini=ptr; *new1=1;} // end res2
  else{*res2=res2min-1;}
  return(res12);
}

// Assembly:
pushq %rbp
pushq %r14
pushq %rbx
movl $0x0, (%rdi)
movq (%r8), %r10
xorl %eax, %eax
cmpq %r9, %r10
jae 0x19714
movl %ecx, %r11d
cmpl %esi, 0x8(%r10)
jg 0x1970f
jne 0x19703
movl 0xc(%r10), %ebx
movl (%rdx), %ebp
cmpl %ebp, %ebx
jne 0x196f3
testq %rax, %rax
je 0x196ee
movsd 0x28(%r10), %xmm0
ucomisd 0x28(%rax), %xmm0
ja 0x196ee
jne 0x19703
jp 0x19703
movsd 0x20(%rax), %xmm0
ucomisd 0x20(%r10), %xmm0
jbe 0x19703
movq %r10, %rax
jmp 0x19703
cmpl %r11d, %ebx
movl %r11d, %r14d
cmovll %ebx, %r14d
cmpl %ebp, %ebx
cmovgl %r14d, %r11d
addq $0xc8, %r10
cmpq %r9, %r10
jb 0x196b9
cmpl %ecx, %r11d
jne 0x1971f
movq %r10, (%r8)
movl $0x1, (%rdi)
jmp 0x19725
decl %r11d
movl %r11d, (%rdx)
popq %rbx
popq %r14
popq %rbp
retq
// File: /workspace/github_c_makefile_star2_O3/ugobas[P]tnm/interactions_tnm.c
// Function: Compute_sec_der_old
void Compute_sec_der_old(struct interaction *Int_list, int N_int, float expo)
{
  double t=0; int n;
    for(n=0; n<N_int; n++){
      if(Int_list[n].r0<=0){
	printf("ERROR, d2= %.2f\n", Int_list[n].r0); exit(8);
      }
      if(expo)Int_list[n].sec_der*=pow(Int_list[n].r0, -expo);
      t+=Int_list[n].sec_der;
    }
    t/=N_int;
    for(n=0; n<N_int; n++)Int_list[n].sec_der/=t;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
subq $0x20, %rsp
testl %esi, %esi
jle 0x19829
movaps %xmm0, %xmm2
movl %esi, %ebp
movq %rdi, %rbx
movaps 0x378bf(%rip), %xmm0    # 0x51030
xorps %xmm2, %xmm0
cvtss2sd %xmm0, %xmm0
movsd %xmm0, (%rsp)
movl %esi, %r14d
imulq $0xc8, %r14, %r15
xorl %r12d, %r12d
xorpd %xmm1, %xmm1
movaps %xmm2, 0x10(%rsp)
movsd 0x20(%rbx,%r12), %xmm0
xorpd %xmm3, %xmm3
ucomisd %xmm0, %xmm3
jae 0x19836
ucomiss 0x37969(%rip), %xmm2    # 0x51118
jne 0x197b3
jnp 0x197de
movsd %xmm1, 0x8(%rsp)
movsd (%rsp), %xmm1
callq 0x20a0
movaps 0x10(%rsp), %xmm2
movsd 0x8(%rsp), %xmm1
mulsd 0x28(%rbx,%r12), %xmm0
movsd %xmm0, 0x28(%rbx,%r12)
jmp 0x197e5
movsd 0x28(%rbx,%r12), %xmm0
addsd %xmm0, %xmm1
addq $0xc8, %r12
cmpq %r12, %r15
jne 0x19793
testl %ebp, %ebp
jle 0x19829
xorps %xmm0, %xmm0
cvtsi2sd %ebp, %xmm0
divsd %xmm0, %xmm1
imulq $0xc8, %r14, %rax
xorl %ecx, %ecx
movsd 0x28(%rbx,%rcx), %xmm0
divsd %xmm1, %xmm0
movsd %xmm0, 0x28(%rbx,%rcx)
addq $0xc8, %rcx
cmpq %rcx, %rax
jne 0x1980d
addq $0x20, %rsp
popq %rbx
popq %r12
popq %r14
popq %r15
popq %rbp
retq
leaq 0x3d973(%rip), %rdi     # 0x571b0
movb $0x1, %al
callq 0x20f0
movl $0x8, %edi
callq 0x2280
nop
// File: /workspace/github_c_makefile_star2_O3/ugobas[P]tnm/anharmonic_tnm.c
// Function: Compute_pow_ia
double **Compute_pow_ia(int npara, int nsam){
  double **pow_ia=Allocate_mat2_d(nsam, npara);
  for(int a=0; a<npara; a++){
    for(int i=1; i<nsam; i++)pow_ia[i][a]=pow(i,a);
  }
  return(pow_ia);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
pushq %rax
movl %esi, %ebx
movl %edi, %ebp
movl %esi, %edi
movl %ebp, %esi
callq 0x49804
movq %rax, %r14
testl %ebp, %ebp
jle 0x1b465
movl %ebp, %r15d
movl %ebx, %r12d
xorl %r13d, %r13d
cmpl $0x2, %ebx
jl 0x1b45d
xorps %xmm0, %xmm0
cvtsi2sd %r13d, %xmm0
movsd %xmm0, (%rsp)
movl $0x1, %ebp
xorps %xmm0, %xmm0
cvtsi2sd %ebp, %xmm0
movsd (%rsp), %xmm1
callq 0x20a0
movq (%r14,%rbp,8), %rax
movsd %xmm0, (%rax,%r13,8)
incq %rbp
cmpq %rbp, %r12
jne 0x1b43a
incq %r13
cmpq %r15, %r13
jne 0x1b423
movq %r14, %rax
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_c_makefile_star2_O3/ugobas[P]tnm/anharmonic_tnm.c
// Function: Get_fit_para
void Get_fit_para(double *para, double *y, double **Corr_inv,
		  int npara, int nsam, float scale)
{
  // p_a*d^a = sum_b C^(-1)_ab sum_i i^b*y_i = C^(-1) Y
  double YX[npara]; int a, b, i;
  for(a=0; a<npara; a++){
    YX[a]=0; for(i=1; i<nsam; i++)YX[a]+=pow_ia[i][a]*y[i]; //pow(i,a)
  }
  for(a=0; a<npara; a++){
    para[a]=0;
    for(b=0; b<npara; b++)para[a]+=Corr_inv[a][b]*YX[b];
    para[a]*=pow(scale, -a);
  }
  if(0){
    printf("Y for fit: "); for(a=1; a<npara; a++)printf(" %.2g",y[a]);
    printf("\n");
    printf("YX for fit: "); for(a=0; a<npara; a++)printf(" %.2g",YX[a]);
    printf("\n");
    for(a=0; a<npara; a++){if(isnan(YX[a]))exit(8);}
  }
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
pushq %rax
movl %ecx, %r15d
movq %rsp, %r12
leaq 0xf(,%r15,8), %rax
andq $-0x10, %rax
subq %rax, %r12
movq %r12, %rsp
testl %ecx, %ecx
jle 0x1c045
movq %rdx, %rbx
movq %rdi, %r14
movq 0x5d31a(%rip), %rax     # 0x792b0
movl %r8d, %edx
xorl %edi, %edi
movq $0x0, (%r12,%rdi,8)
cmpl $0x2, %r8d
jl 0x1bfd5
xorpd %xmm1, %xmm1
movl $0x1, %r9d
movq (%rax,%r9,8), %r10
movsd (%r10,%rdi,8), %xmm2
mulsd (%rsi,%r9,8), %xmm2
addsd %xmm2, %xmm1
incq %r9
cmpq %r9, %rdx
jne 0x1bfb3
movsd %xmm1, (%r12,%rdi,8)
incq %rdi
cmpq %r15, %rdi
jne 0x1bf9b
testl %ecx, %ecx
jle 0x1c045
cvtss2sd %xmm0, %xmm0
movsd %xmm0, -0x30(%rbp)
xorl %r13d, %r13d
movq $0x0, (%r14,%r13,8)
movq (%rbx,%r13,8), %rax
xorps %xmm0, %xmm0
xorl %ecx, %ecx
movsd (%rax,%rcx,8), %xmm1
mulsd (%r12,%rcx,8), %xmm1
addsd %xmm1, %xmm0
movsd %xmm0, (%r14,%r13,8)
incq %rcx
cmpq %rcx, %r15
jne 0x1bffe
movl %r13d, %eax
negl %eax
xorps %xmm1, %xmm1
cvtsi2sd %eax, %xmm1
movsd -0x30(%rbp), %xmm0
callq 0x20a0
mulsd (%r14,%r13,8), %xmm0
movsd %xmm0, (%r14,%r13,8)
incq %r13
cmpq %r15, %r13
jne 0x1bfed
leaq -0x28(%rbp), %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_c_makefile_star2_O3/ugobas[P]tnm/anharmonic_tnm.c
// Function: Native_interactions
int **Native_interactions(struct interaction *Int_list, int N_int,
			  atom *atoms, int nres)
{
  int **matrix=malloc(nres*sizeof (int *)); int i1, i2;
  for(i1=0; i1<nres; i1++){
    matrix[i1]=malloc(nres*sizeof(int));
    for(i2=0; i2<nres; i2++)matrix[i1][i2]=-1;
  }
  struct interaction *Int=Int_list;
  for(int n=0; n<N_int; n++){
    i1=atoms[Int->i1].res;
    i2=atoms[Int->i2].res;
    if((i1<0)||(i1>=nres)||(i2<0)||(i2>=nres)){
      printf("ERROR in Native_interactions i1=%d i2= %d nres=%d n= %d\n",
	     i1, i2, nres, n); exit(8);
    }
    matrix[i1][i2]=n;
    matrix[i2][i1]=n;
    Int++;
  }
  printf("\n%d native interactions stored\n", N_int);
  return(matrix);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x18, %rsp
movq %rdx, 0x10(%rsp)
movl %esi, 0x8(%rsp)
movq %rdi, %r12
movl %ecx, 0xc(%rsp)
movslq %ecx, %rbx
leaq (,%rbx,8), %rdi
callq 0x21e0
movq %rax, %r14
testl %ebx, %ebx
jle 0x1c149
leaq (,%rbx,4), %r13
movl %ebx, %ebx
leaq (,%rbx,4), %rbp
xorl %r15d, %r15d
movq %r13, %rdi
callq 0x21e0
movq %rax, (%r14,%r15,8)
movq %rax, %rdi
movl $0xff, %esi
movq %rbp, %rdx
callq 0x2110
incq %r15
cmpq %r15, %rbx
jne 0x1c125
movl 0x8(%rsp), %r9d
testl %r9d, %r9d
movl 0xc(%rsp), %ecx
movq 0x10(%rsp), %rdi
jle 0x1c1a9
xorl %r8d, %r8d
movslq (%r12), %rax
imulq $0x70, %rax, %rax
movslq 0x64(%rdi,%rax), %rsi
movslq 0x4(%r12), %rax
imulq $0x70, %rax, %rax
movl 0x64(%rdi,%rax), %edx
testq %rsi, %rsi
js 0x1c1cc
cmpl %ecx, %esi
jge 0x1c1cc
testl %edx, %edx
js 0x1c1cc
cmpl %ecx, %edx
jge 0x1c1cc
movq (%r14,%rsi,8), %rax
movl %r8d, (%rax,%rdx,4)
movq (%r14,%rdx,8), %rax
movl %r8d, (%rax,%rsi,4)
addq $0xc8, %r12
incl %r8d
cmpl %r8d, %r9d
jne 0x1c15f
leaq 0x3b721(%rip), %rdi     # 0x578d1
movl %r9d, %esi
xorl %eax, %eax
callq 0x20f0
movq %r14, %rax
addq $0x18, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
leaq 0x3b6c5(%rip), %rdi     # 0x57898
xorl %eax, %eax
callq 0x20f0
movl $0x8, %edi
callq 0x2280
// File: /workspace/github_c_makefile_star2_O3/ugobas[P]tnm/anharmonic_tnm.c
// Function: Find_CA
int Find_CA(atom *atoms, int *i1, int res, int natoms){
  atom *atom=atoms+*i1; int i=*i1;
  while(atom->res <= res){
    if((atom->res == res)&&(strncmp(atom->name, "CA", 2)==0)){
      *i1=i; return(0);
    }
    atom++; i++; if(i>=natoms)break;
  }
  if(*i1)(*i1)--;
  if(atoms[*i1].res==res)return(0);
  return(-1);
}

// Assembly:
movl (%rsi), %r8d
movslq %r8d, %rax
imulq $0x70, %rax, %r9
addq %rdi, %r9
addq $0x18, %r9
movl %r8d, %r10d
cmpl %edx, 0x4c(%r9)
jg 0x1c219
jne 0x1c20d
cmpb $0x43, (%r9)
jne 0x1c20d
cmpb $0x41, 0x1(%r9)
je 0x1c237
incl %r10d
addq $0x70, %r9
cmpl %ecx, %r10d
jl 0x1c1f8
testl %r8d, %r8d
je 0x1c225
decq %rax
movl %eax, (%rsi)
jmp 0x1c227
xorl %eax, %eax
imulq $0x70, %rax, %rcx
xorl %eax, %eax
cmpl %edx, 0x64(%rdi,%rcx)
setne %al
negl %eax
retq
movl %r10d, (%rsi)
xorl %eax, %eax
retq
// File: /workspace/github_c_makefile_star2_O3/ugobas[P]tnm/anharmonic_tnm.c
// Function: Initialize_energy
void Initialize_energy(struct interaction *Int_list, int N_int,
		       float *r, atom *atoms, int natoms, int nres)
{
  INI_ENERGY=1;

  E_repulsion=Int_list[0].A4;
  C22=C_THR*C_THR;
  lambda_2=2*lambda;
  lambda_4=4*lambda;
  
  // Auxiliary matrices for energy computation: interaction list
  nat_res1_res2=Native_interactions(Int_list, N_int, atoms, nres);
  // YY matrix with the index of the interaction data (in Int_list)
  INT_MAX=150*nres;
  Int_list_all=malloc(INT_MAX*sizeof(struct interaction));
  N_int_all=
    Interactions_CA(Int_list_all, d_CA_high, atoms, natoms, nres, r);
  N_int_ini=N_int_all;
  Int_list_ini=malloc(N_int_ini*sizeof(struct interaction));
  for(int i=0; i<N_int_ini; i++)Int_list_ini[i]=Int_list_all[i];
  // YY List of interactions with d(CA-CA)<22
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
movl %r9d, %ebx
movl %r8d, %ebp
movq %rcx, %r15
movq %rdx, %r14
movl $0x1, 0x5cdfe(%rip)     # 0x7928c
cvtsd2ss 0x48(%rdi), %xmm0
leaq 0x4cc3e(%rip), %rax     # 0x690d8
movss %xmm0, (%rax)
leaq 0x4cbc3(%rip), %rax     # 0x69068
movss (%rax), %xmm0
mulss %xmm0, %xmm0
movss %xmm0, 0x5ce27(%rip)    # 0x792dc
leaq 0x4cc18(%rip), %rax     # 0x690d4
movss (%rax), %xmm0
movaps %xmm0, %xmm1
addss %xmm0, %xmm1
cvtss2sd %xmm1, %xmm1
movsd %xmm1, 0x5ce0d(%rip)    # 0x792e0
mulss 0x3a9cd(%rip), %xmm0    # 0x56ea8
cvtss2sd %xmm0, %xmm0
movsd %xmm0, 0x5ce01(%rip)    # 0x792e8
movq %rcx, %rdx
movl %r9d, %ecx
callq 0x1c0db
movq %rax, 0x5cdd7(%rip)     # 0x792d0
imull $0x96, %ebx, %eax
movl %eax, 0x5cdd3(%rip)     # 0x792d8
cltq
imulq $0xc8, %rax, %rdi
callq 0x21e0
movq %rax, 0x5cda6(%rip)     # 0x792c0
movss 0x4c96e(%rip), %xmm0    # 0x68e90
movq %rax, %rdi
movq %r15, %rsi
movl %ebp, %edx
movl %ebx, %ecx
movq %r14, %r8
callq 0x1b477
movl %eax, %ebp
movl %eax, 0x5cd80(%rip)     # 0x792bc
movl %eax, 0x5cd76(%rip)     # 0x792b8
movslq %eax, %r14
imulq $0xc8, %r14, %rdi
callq 0x21e0
movq %rax, 0x5cd70(%rip)     # 0x792c8
testl %r14d, %r14d
jle 0x1c591
movq %rax, %rbx
movq 0x5cd59(%rip), %r14     # 0x792c0
movl %ebp, %eax
imulq $0xc8, %rax, %r15
xorl %r12d, %r12d
leaq (%rbx,%r12), %rdi
leaq (%r14,%r12), %rsi
movl $0xc8, %edx
callq 0x21b0
addq $0xc8, %r12
cmpq %r12, %r15
jne 0x1c573
popq %rbx
popq %r12
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_c_makefile_star2_O3/ugobas[P]tnm/anharmonic_tnm.c
// Function: Dist_thr
int Dist_thr(struct interaction *Int_cov, atom *atoms){
  if(Int_cov->i1>=0 && Int_cov->i2>=0){   
    double *r1=atoms[Int_cov->i1].r, *r2=atoms[Int_cov->i2].r;
    double d2=0, d; for(int j=0; j<3; j++){d=r1[j]-r2[j]; d*=d; d2+=d;}
    Int_cov->r0=sqrt(d2);
    if(Int_cov->r0 < C_THR){return(1);}
    else{return(0);}
  }
  return(0);
}

// Assembly:
movslq (%rdi), %rcx
xorl %eax, %eax
testq %rcx, %rcx
js 0x1d580
pushq %rbx
movq %rdi, %rbx
movslq 0x4(%rdi), %rdx
testq %rdx, %rdx
js 0x1d57f
imulq $0x70, %rcx, %rax
addq %rsi, %rax
imulq $0x70, %rdx, %rcx
addq %rsi, %rcx
xorpd %xmm0, %xmm0
xorl %edx, %edx
movsd (%rax,%rdx,8), %xmm1
subsd (%rcx,%rdx,8), %xmm1
mulsd %xmm1, %xmm1
addsd %xmm1, %xmm0
incq %rdx
cmpq $0x3, %rdx
jne 0x1d533
xorpd %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jb 0x1d55e
sqrtsd %xmm0, %xmm0
jmp 0x1d563
callq 0x22a0
movsd %xmm0, 0x20(%rbx)
leaq 0x4baf9(%rip), %rax     # 0x69068
xorps %xmm1, %xmm1
cvtss2sd (%rax), %xmm1
xorl %eax, %eax
ucomisd %xmm0, %xmm1
seta %al
popq %rbx
retq