2 # -*- coding: utf-8 -*-
5 from numpy
.matlib
import zeros
,isinf
6 from penalty_lookup_new
import *
8 def compute_donacc(donor_supp
, acceptor_supp
, d
, a
):
10 assert(len(donor_supp
)==len(acceptor_supp
))
11 size
= len(donor_supp
)
16 for idx
in range(size
):
17 if isinf(donor_supp
[idx
]):
18 donor
[idx
] = donor_supp
[idx
]
20 donor
[idx
] = penalty_lookup_new(d
, donor_supp
[idx
])
22 if isinf(acceptor_supp
[idx
]):
23 acceptor
[idx
] = acceptor_supp
[idx
]
25 acceptor
[idx
] = penalty_lookup_new(a
,acceptor_supp
[idx
])