+ fixed some index bugs in the weights calculation
[qpalma.git] / python / compute_donacc.py
index b198877..bad0872 100644 (file)
@@ -8,13 +8,12 @@ from penalty_lookup_new import *
 def compute_donacc(donor_supp, acceptor_supp, d, a):
 
    assert(len(donor_supp)==len(acceptor_supp))
+   size = len(donor_supp)
   
-   #donor = zeros((len(donor_supp),1))
-   #acceptor= zeros((len(acceptor_supp),1))
-   donor = [0.0]*len(donor_supp)
-   acceptor= [0.0]*len(acceptor_supp)
+   donor    = [0.0]*size
+   acceptor = [0.0]*size
 
-   for idx in range(len(donor_supp)):
+   for idx in range(size):
      if isinf(donor_supp[idx]):
        donor[idx] = donor_supp[idx]
      else: