Commit c4bfc3ba by Chloe Dequeker

bugfix for atomRes interface

parent a14ec87d
...@@ -172,12 +172,18 @@ void getInterface(struct pdb_values* pdbR, struct pdb_values* pdbL, int confID){ ...@@ -172,12 +172,18 @@ void getInterface(struct pdb_values* pdbR, struct pdb_values* pdbL, int confID){
*/ */
for(i=0;i<nbCandR;i++){ for(i=0;i<nbCandR;i++){
t_candidateR[i]->isCandidate = 0; t_candidateR[i]->isCandidate = 0;
for(j=0;j<t_candidateR[i]->nbAtom;j++){
t_candidateR[i]->isAtomCandidate[j] = 0;
}
} }
/* The first nbCandL values of t_candidateL correspond to the possible candidates /* The first nbCandL values of t_candidateL correspond to the possible candidates
* for the ligand * for the ligand
*/ */
for(i=0;i<nbCandL;i++){ for(i=0;i<nbCandL;i++){
t_candidateL[i]->isCandidate = 0; t_candidateL[i]->isCandidate = 0;
for(j=0;j<t_candidateL[i]->nbAtom;j++){
t_candidateL[i]->isAtomCandidate[j] = 0;
}
} }
/* Don't worry for that, this is only O(n^2) over the remaining /* Don't worry for that, this is only O(n^2) over the remaining
...@@ -195,10 +201,6 @@ void getInterface(struct pdb_values* pdbR, struct pdb_values* pdbL, int confID){ ...@@ -195,10 +201,6 @@ void getInterface(struct pdb_values* pdbR, struct pdb_values* pdbL, int confID){
/* For each candidate residue of the ligand */ /* For each candidate residue of the ligand */
for(k=0;k<nbCandL;k++){ for(k=0;k<nbCandL;k++){
/* If both the residues we are looking at are already selected */
if(force && t_candidateL[k]->isCandidate && t_candidateR[i]->isCandidate){
continue;
}
for(l=0;l<t_candidateL[k]->nbAtom;l++){ for(l=0;l<t_candidateL[k]->nbAtom;l++){
x2 = t_candidateL[k]->x[l]; x2 = t_candidateL[k]->x[l];
y2 = t_candidateL[k]->y[l]; y2 = t_candidateL[k]->y[l];
...@@ -211,7 +213,6 @@ void getInterface(struct pdb_values* pdbR, struct pdb_values* pdbL, int confID){ ...@@ -211,7 +213,6 @@ void getInterface(struct pdb_values* pdbR, struct pdb_values* pdbL, int confID){
if(clash >= TOO_MANY_CLASHES){ if(clash >= TOO_MANY_CLASHES){
break; break;
} }
continue;
} }
if(dist <= DIST_FOR_CONTACT){ if(dist <= DIST_FOR_CONTACT){
t_candidateR[i]->isCandidate = 1; t_candidateR[i]->isCandidate = 1;
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment