Skip to content
Snippets Groups Projects
Commit 37662cdf authored by Gauthier Roy's avatar Gauthier Roy
Browse files

modification of the fitness

parent 75a31e9c
No related branches found
No related tags found
No related merge requests found
No preview for this file type
...@@ -17,26 +17,27 @@ class Individu(): ...@@ -17,26 +17,27 @@ class Individu():
def evaluate(self, brin): def evaluate(self, brin):
traj = Traj3D() traj = Traj3D()
numb_ajout = 3 numb_ajout = 6
fisrt_seq = brin[0:numb_ajout] fisrt_seq = brin[0:numb_ajout]
last_seq = brin[-numb_ajout:] last_seq = brin[-numb_ajout:]
traj.compute(last_seq + brin + fisrt_seq, self.table) traj.compute(last_seq + brin + fisrt_seq, self.table)
traj_array = np.array(traj.getTraj()) traj_array = np.array(traj.getTraj())
print(traj_array)
list_distance = [] list_distance = []
for i in range(numb_ajout): begining = traj_array[0:2*numb_ajout, 0:3]
first_nuc_coordonate = traj_array[numb_ajout+i, 0:3] end = traj_array[-2*numb_ajout:, 0:3]
first_nuc_coordonate_compute = traj_array[-(numb_ajout-i), 0:3]
last_nuc_coordonate = traj_array[-(2*numb_ajout-i), 0:3]
last_nuc_coordonate_compute = traj_array[i, 0:3]
distance_first_nuc = np.linalg.norm(first_nuc_coordonate - first_nuc_coordonate_compute, ord=2) for i in range(numb_ajout):
distance_last_nuc = np.linalg.norm(last_nuc_coordonate - last_nuc_coordonate_compute, ord=2)
list_distance += [distance_first_nuc, distance_last_nuc] nuc_coordonate_beg = begining[i]
print(nuc_coordonate_beg)
nuc_coordonate_end = end[i]
distance_nuc = np.linalg.norm(nuc_coordonate_beg - nuc_coordonate_end, ord=2)
list_distance += [distance_nuc]
self.score = max(list_distance) self.score = max(list_distance)
...@@ -63,10 +64,10 @@ class Individu(): ...@@ -63,10 +64,10 @@ class Individu():
# print(individu1.table.rot_table) # print(individu1.table.rot_table)
# individu1.mutation() # individu1.mutation()
# table = RotTable() table = RotTable()
# test = Individu(table) test = Individu(table)
# test.evaluate("AAAGGATCTTCTTGAGATCCTTTTTTTCTGCGCGTAATCTGCTGCCAGTAAACGAAAAAACCGCCTGGGGAGGCGGTTTAGTCGAA") test.evaluate("AAAGGATCTTCTTGAGATCCTTTTTTTCTGCGCGTAATCTGCTGCCAGTAAACGAAAAAACCGCCTGGGGAGGCGGTTTAGTCGAA")
# print(test.score) print(test.score)
# qqun=Individu(RotTable()) # qqun=Individu(RotTable())
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment