#!python
#!/usr/bin/env python
#Python script:
import sys
import numpy as np
from pybdm import BDM
from pybdm import options
options.set(raise_if_zero=False)
file = sys.argv[1]
with open(file, 'r') as file:
data = file.read()
imagelist = data.splitlines()
number_of_lines=len(imagelist[0])
s = ""
image_array = np.empty((0,number_of_lines),dtype=int)
for line in imagelist:
for char in line:
s += char + "-"
image_array=np.vstack([image_array, np.fromstring(s, dtype=int, sep="-")])
s = ""
bdm = BDM(ndim=2,nsymbols=2)
nbdm = bdm.nbdm(image_array)
print(nbdm)
#Module sys has to be imported:
import sys
import numpy as np
from pybdm import BDM
from pybdm import options
options.set(raise_if_zero=False)
file = sys.argv[1]
with open(file, 'r') as file:
data = file.read().replace('\n', '')
tape = np.fromstring(data, dtype=int, sep='-')
bdm = BDM(ndim=1, nsymbols=2)
value = bdm.nbdm(tape)
value2 = bdm.nent(tape)
if value2 == data:
value2 = 0
if (value == 0 and (np.isnan(value2) is False)
and value2 > 0.0) or (value2 < value and value2 > 0.0):
value = value2
print(value)
#!python
#!/usr/bin/env python
#Python script:
import sys
import numpy as np
from pybdm import BDM
from pybdm import options
import warnings
warnings.filterwarnings("ignore")
options.set(raise_if_zero=False)
options.set(warn_if_missing_ctm=False)
file = sys.argv[1]
with open(file, 'r') as file:
data = file.read()
Translator = {"A": 0, "C": 1, "G": 2, "T": 3}
datafile = list(map(lambda x: Translator[x], data))
datafile = np.asarray(datafile)
bdm = BDM(ndim=1, nsymbols=4)
nbdm = bdm.nbdm(datafile)
print(nbdm)