biotite 0.41.1__cp310-cp310-win_amd64.whl → 1.0.0__cp310-cp310-win_amd64.whl

biotite/sequence/alphabet.py

@@ -4,58 +4,62 @@
 
 __name__ = "biotite.sequence"
 __author__ = "Patrick Kunzmann"
-__all__ = ["Alphabet", "LetterAlphabet", "AlphabetMapper", "AlphabetError",
-           "common_alphabet"]
+__all__ = [
+    "Alphabet",
+    "LetterAlphabet",
+    "AlphabetMapper",
+    "AlphabetError",
+    "common_alphabet",
+]
 
-import copy
-from numbers import Integral
 import string
+from numbers import Integral
 import numpy as np
-from .codec import encode_chars, decode_to_chars, map_sequence_code
+from biotite.sequence.codec import decode_to_chars, encode_chars, map_sequence_code
 
 
 class Alphabet(object):
     """
     This class defines the allowed symbols for a :class:`Sequence` and
     handles the encoding/decoding between symbols and symbol codes.
-    
+
     An :class:`Alphabet` is created with the list of symbols, that can
     be used in this context.
     In most cases a symbol will be simply a letter, hence a string of
     length 1. But in principle every hashable Python object can serve
     as symbol.
-    
+
     The encoding of a symbol into a symbol code is
     done in the following way: Find the first index in the symbol list,
     where the list element equals the symbol. This index is the
     symbol code. If the symbol is not found in the list, an
     :class:`AlphabetError` is raised.
-    
+
     Internally, a dictionary is used for encoding, with symbols as keys
     and symbol codes as values. Therefore, every symbol must be
     hashable. For decoding the symbol list is indexed with the symbol
     code.
-    
+
     If an alphabet *1* contains the same symbols and the same
     symbol-code-mappings like another alphabet *2*, but alphabet *1*
     introduces also new symbols, then alphabet *1* *extends* alphabet
     *2*.
     Per definition, every alphabet also extends itself.
-    
+
     Objects of this class are immutable.
-    
+
     Parameters
     ----------
     symbols : iterable object
         The symbols, that are allowed in this alphabet. The
         corresponding code for a symbol, is the index of that symbol
         in this list.
-    
+
     Examples
     --------
     Create an Alphabet containing DNA letters and encode/decode a
     letter/code:
-    
+
     >>> alph = Alphabet(["A","C","G","T"])
     >>> print(alph.encode("G"))
     2
@@ -66,9 +70,9 @@ class Alphabet(object):
     ... except Exception as e:
     ...    print(e)
     Symbol 'foo' is not in the alphabet
-    
+
     Create an Alphabet of arbitrary objects:
-    
+
     >>> alph = Alphabet(["foo", 42, (1,2,3), 5, 3.141])
     >>> print(alph.encode((1,2,3)))
     2
@@ -77,53 +81,53 @@ class Alphabet(object):
 
     On the subject of alphabet extension:
     An alphabet always extends itself.
-    
+
     >>> Alphabet(["A","C","G","T"]).extends(Alphabet(["A","C","G","T"]))
     True
 
     An alphabet extends an alphabet when it contains additional symbols...
-    
+
     >>> Alphabet(["A","C","G","T","U"]).extends(Alphabet(["A","C","G","T"]))
     True
-    
+
     ...but not vice versa
-    
+
     >>> Alphabet(["A","C","G","T"]).extends(Alphabet(["A","C","G","T","U"]))
     False
-    
+
     Two alphabets with same symbols but different symbol-code-mappings
-    
-    >>> Alphabet(["A","C","G","T"]).extends(Alphabet(["A","C","T","G"]))    
+
+    >>> Alphabet(["A","C","G","T"]).extends(Alphabet(["A","C","T","G"]))
     False
     """
-    
+
     def __init__(self, symbols):
         if len(symbols) == 0:
             raise ValueError("Symbol list is empty")
-        self._symbols = copy.deepcopy(list(symbols))
+        self._symbols = tuple(symbols)
         self._symbol_dict = {}
         for i, symbol in enumerate(symbols):
             self._symbol_dict[symbol] = i
 
     def __repr__(self):
         """Represent Alphabet as a string for debugging."""
-        return f'Alphabet({self._symbols})'
+        return f"Alphabet({self._symbols})"
 
     def get_symbols(self):
         """
         Get the symbols in the alphabet.
-        
+
         Returns
         -------
-        symbols : list
-            Copy of the internal list of symbols.
+        symbols : tuple
+            The symbols.
         """
-        return copy.deepcopy(self._symbols)
-    
+        return self._symbols
+
     def extends(self, alphabet):
         """
         Check, if this alphabet extends another alphabet.
-        
+
         Parameters
         ----------
         alphabet : Alphabet
@@ -139,23 +143,22 @@ class Alphabet(object):
         elif len(alphabet) > len(self):
             return False
         else:
-            return alphabet.get_symbols() \
-                == self.get_symbols()[:len(alphabet)]
-    
+            return alphabet.get_symbols() == self.get_symbols()[: len(alphabet)]
+
     def encode(self, symbol):
         """
         Use the alphabet to encode a symbol.
-        
+
         Parameters
         ----------
         symbol : object
             The object to encode into a symbol code.
-        
+
         Returns
         -------
         code : int
             The symbol code of `symbol`.
-        
+
         Raises
         ------
         AlphabetError
@@ -164,24 +167,22 @@ class Alphabet(object):
         try:
             return self._symbol_dict[symbol]
         except KeyError:
-            raise AlphabetError(
-                f"Symbol {repr(symbol)} is not in the alphabet"
-            )
-    
+            raise AlphabetError(f"Symbol {repr(symbol)} is not in the alphabet")
+
     def decode(self, code):
         """
         Use the alphabet to decode a symbol code.
-        
+
         Parameters
         ----------
         code : int
             The symbol code to be decoded.
-        
+
         Returns
         -------
         symbol : object
             The symbol corresponding to `code`.
-        
+
         Raises
         ------
         AlphabetError
@@ -190,41 +191,41 @@ class Alphabet(object):
         if code < 0 or code >= len(self._symbols):
             raise AlphabetError(f"'{code:d}' is not a valid code")
         return self._symbols[code]
-    
+
     def encode_multiple(self, symbols, dtype=np.int64):
         """
         Encode a list of symbols.
-        
+
         Parameters
         ----------
         symbols : array-like
             The symbols to encode.
         dtype : dtype, optional
             The dtype of the output ndarray. (Default: `int64`)
-            
+
         Returns
         -------
         code : ndarray
             The sequence code.
         """
         return np.array([self.encode(e) for e in symbols], dtype=dtype)
-    
+
     def decode_multiple(self, code):
         """
         Decode a sequence code into a list of symbols.
-        
+
         Parameters
         ----------
         code : ndarray
             The sequence code to decode.
-        
+
         Returns
         -------
         symbols : list
             The decoded list of symbols.
         """
         return [self.decode(c) for c in code]
-    
+
     def is_letter_alphabet(self):
         """
         Check whether the symbols in this alphabet are single printable
@@ -238,30 +239,33 @@ class Alphabet(object):
             have length 1 and are printable.
         """
         for symbol in self:
-            if not isinstance(symbol, (str, bytes)) \
-                or len(symbol) > 1:
-                    return False
+            if not isinstance(symbol, (str, bytes)) or len(symbol) > 1:
+                return False
             if isinstance(symbol, str):
                 symbol = symbol.encode("ASCII")
-            if symbol not in LetterAlphabet.PRINATBLES:
+            if symbol not in LetterAlphabet.PRINTABLES:
                 return False
         return True
-    
+
     def __str__(self):
         return str(self.get_symbols())
-    
+
     def __len__(self):
         return len(self.get_symbols())
-    
+
     def __iter__(self):
         return self.get_symbols().__iter__()
-    
+
     def __contains__(self, symbol):
         return symbol in self.get_symbols()
-    
+
     def __hash__(self):
-        return hash(tuple(self._symbols))
-    
+        symbols = self.get_symbols()
+        if isinstance(symbols, tuple):
+            return hash(symbols)
+        else:
+            return hash(tuple(symbols))
+
     def __eq__(self, item):
         if item is self:
             return True
@@ -291,9 +295,10 @@ class LetterAlphabet(Alphabet):
         corresponding code for a symbol, is the index of that symbol
         in this list.
     """
-    
-    PRINATBLES = (string.digits + string.ascii_letters + string.punctuation) \
-                 .encode("ASCII")
+
+    PRINTABLES = (string.digits + string.ascii_letters + string.punctuation).encode(
+        "ASCII"
+    )
 
     def __init__(self, symbols):
         if len(symbols) == 0:
@@ -304,7 +309,7 @@ class LetterAlphabet(Alphabet):
                 raise ValueError(f"Symbol '{symbol}' is not a single letter")
             if isinstance(symbol, str):
                 symbol = symbol.encode("ASCII")
-            if symbol not in LetterAlphabet.PRINATBLES:
+            if symbol not in LetterAlphabet.PRINTABLES:
                 raise ValueError(
                     f"Symbol {repr(symbol)} is not printable or whitespace"
                 )
@@ -312,57 +317,43 @@ class LetterAlphabet(Alphabet):
         # Direct 'astype' conversion is not allowed by numpy
         # -> frombuffer()
         self._symbols = np.frombuffer(
-            np.array(self._symbols, dtype="|S1"),
-            dtype=np.ubyte
+            np.array(self._symbols, dtype="|S1"), dtype=np.ubyte
         )
 
     def __repr__(self):
         """Represent LetterAlphabet as a string for debugging."""
-        return f'LetterAlphabet({self.get_symbols()})'
-    
+        return f"LetterAlphabet({self.get_symbols()})"
+
     def extends(self, alphabet):
         if alphabet is self:
             return True
-        elif type(alphabet) == LetterAlphabet:
+        elif isinstance(alphabet, LetterAlphabet):
             if len(alphabet._symbols) > len(self._symbols):
                 return False
-            return np.all(
-                alphabet._symbols == self._symbols[:len(alphabet._symbols)]
-            )
+            return np.all(alphabet._symbols == self._symbols[: len(alphabet._symbols)])
         else:
             return super().extends(alphabet)
 
     def get_symbols(self):
-        """
-        Get the symbols in the alphabet.
-        
-        Returns
-        -------
-        symbols : list
-            Copy of the internal list of symbols.
-        """
-        return [symbol.decode("ASCII") for symbol
-                in self._symbols_as_bytes()]
-    
+        return tuple([symbol.decode("ASCII") for symbol in self._symbols_as_bytes()])
+
     def encode(self, symbol):
         if not isinstance(symbol, (str, bytes)) or len(symbol) > 1:
             raise AlphabetError(f"Symbol '{symbol}' is not a single letter")
         indices = np.where(self._symbols == ord(symbol))[0]
         if len(indices) == 0:
-            raise AlphabetError(
-                f"Symbol {repr(symbol)} is not in the alphabet"
-            )
-        return indices[0]
-    
+            raise AlphabetError(f"Symbol {repr(symbol)} is not in the alphabet")
+        return indices[0].item()
+
     def decode(self, code, as_bytes=False):
         if code < 0 or code >= len(self._symbols):
             raise AlphabetError(f"'{code:d}' is not a valid code")
         return chr(self._symbols[code])
-    
+
     def encode_multiple(self, symbols, dtype=None):
         """
         Encode multiple symbols.
-        
+
         Parameters
         ----------
         symbols : iterable object or str or bytes
@@ -371,7 +362,7 @@ class LetterAlphabet(Alphabet):
             containing the symbols is provided, instead of e.g. a list.
         dtype : dtype, optional
             For compatibility with superclass. The value is ignored
-            
+
         Returns
         -------
         code : ndarray
@@ -382,20 +373,17 @@ class LetterAlphabet(Alphabet):
         elif isinstance(symbols, bytes):
             symbols = np.frombuffer(symbols, dtype=np.ubyte)
         elif isinstance(symbols, np.ndarray):
-            symbols = np.frombuffer(
-                symbols.astype(dtype="|S1"), dtype=np.ubyte
-            )
+            symbols = np.frombuffer(symbols.astype(dtype="|S1"), dtype=np.ubyte)
         else:
             symbols = np.frombuffer(
-                np.array(list(symbols), dtype="|S1"),
-                dtype=np.ubyte
+                np.array(list(symbols), dtype="|S1"), dtype=np.ubyte
             )
         return encode_chars(alphabet=self._symbols, symbols=symbols)
-    
+
     def decode_multiple(self, code, as_bytes=False):
         """
         Decode a sequence code into a list of symbols.
-        
+
         Parameters
         ----------
         code : ndarray, dtype=uint8
@@ -421,20 +409,19 @@ class LetterAlphabet(Alphabet):
         if not as_bytes:
             symbols = symbols.astype("U1")
         return symbols
-    
+
     def __contains__(self, symbol):
         if not isinstance(symbol, (str, bytes)):
             return False
         return ord(symbol) in self._symbols
-    
+
     def __len__(self):
         return len(self._symbols)
-    
+
     def _symbols_as_bytes(self):
         "Properly convert from dtype 'np.ubyte' to '|S1'"
         return np.frombuffer(self._symbols, dtype="|S1")
 
-    
 
 class AlphabetMapper(object):
     """
@@ -445,7 +432,7 @@ class AlphabetMapper(object):
     alphabet so that the symbol itself is preserved.
     This class works for single symbol codes or an entire sequence code
     likewise.
-    
+
     Parameters
     ----------
     source_alphabet, target_alphabet : Alphabet
@@ -454,7 +441,7 @@ class AlphabetMapper(object):
         The target alphabet must contain at least all symbols of the
         source alphabet, but it is not required that the shared symbols
         are in the same order.
-    
+
     Examples
     --------
 
@@ -470,56 +457,54 @@ class AlphabetMapper(object):
     >>> in_sequence = GeneralSequence(source_alph, "GCCTAT")
     >>> print(in_sequence.code)
     [2 1 1 3 0 3]
-    >>> print(in_sequence)
+    >>> print("".join(in_sequence.symbols))
     GCCTAT
     >>> out_sequence = GeneralSequence(target_alph)
     >>> out_sequence.code = mapper[in_sequence.code]
     >>> print(out_sequence.code)
     [3 4 4 0 2 0]
-    >>> print(out_sequence)
+    >>> print("".join(out_sequence.symbols))
     GCCTAT
     """
-    
+
     def __init__(self, source_alphabet, target_alphabet):
         if target_alphabet.extends(source_alphabet):
             self._necessary_mapping = False
         else:
             self._necessary_mapping = True
             self._mapper = np.zeros(
-                len(source_alphabet),
-                dtype=AlphabetMapper._dtype(len(target_alphabet))
+                len(source_alphabet), dtype=AlphabetMapper._dtype(len(target_alphabet))
             )
             for old_code in range(len(source_alphabet)):
                 symbol = source_alphabet.decode(old_code)
                 new_code = target_alphabet.encode(symbol)
                 self._mapper[old_code] = new_code
-        
+
     def __getitem__(self, code):
         if isinstance(code, Integral):
             if self._necessary_mapping:
                 return self._mapper[code]
             else:
                 return code
-        if not isinstance(code, np.ndarray) \
-           or code.dtype not in (np.uint8, np.uint16, np.uint32, np.uint64):
-                code = np.array(code, dtype=np.uint64)
+        if not isinstance(code, np.ndarray) or code.dtype not in (
+            np.uint8,
+            np.uint16,
+            np.uint32,
+            np.uint64,
+        ):
+            code = np.array(code, dtype=np.uint64)
         if self._necessary_mapping:
             mapped_code = np.empty(len(code), dtype=self._mapper.dtype)
-            map_sequence_code(
-                self._mapper,
-                code,
-                mapped_code
-            )
+            map_sequence_code(self._mapper, code, mapped_code)
             return mapped_code
         else:
             return code
 
-
     @staticmethod
     def _dtype(alphabet_size):
-        _size_uint8  = np.iinfo(np.uint8 ).max +1
-        _size_uint16 = np.iinfo(np.uint16).max +1
-        _size_uint32 = np.iinfo(np.uint32).max +1
+        _size_uint8 = np.iinfo(np.uint8).max + 1
+        _size_uint16 = np.iinfo(np.uint16).max + 1
+        _size_uint32 = np.iinfo(np.uint32).max + 1
         if alphabet_size <= _size_uint8:
             return np.uint8
         elif alphabet_size <= _size_uint16:
@@ -535,6 +520,7 @@ class AlphabetError(Exception):
     This exception is raised, when a code or a symbol is not in an
     :class:`Alphabet`.
     """
+
     pass
 
 
@@ -552,7 +538,7 @@ def common_alphabet(alphabets):
     -------
     common_alphabet : Alphabet or None
         The alphabet from `alphabets` that extends all alphabets.
-        ``None`` if no such common alphabet exists. 
+        ``None`` if no such common alphabet exists.
     """
     common_alphabet = None
     for alphabet in alphabets:
@@ -563,4 +549,4 @@ def common_alphabet(alphabets):
                 common_alphabet = alphabet
             else:
                 return None
-    return common_alphabet
+    return common_alphabet