seqLib.py 文件源码

def getComplexRepr(sequenza, K):
    """
    Crea la rappresentazione ottenuta dalla somma delle varie rappresentazioni 
    in una unica matrice. K e' la lista delle dimensioni dei k-mers.
    """
    m_init=getMatrice(sequenza, K[0])
    dim=2**K[-1]
    #out=np.zeros((dim, dim))
    for ki in K[1:]:
        m_init=espandiMatrice(m_init)
        # divido gli elementi per due per diminuire l'impatto dei termini 
        # piu' corti e quindi piu' comuni
        m_init /= 2.0
        temp=getMatrice(sequenza, ki)
        m_init=np.add(m_init,temp)
    return m_init