changed loop in hommola_cospeciation to numpy vectorization