def Mol(self, sequence, popSize, initialKE, minStem):
self.sequence = sequence
dotplot = population.Checkerboard(sequence)
self.stemTable = population.FindDiagonal(sequence, dotplot, minStem)
self.infoTable = self.stemTable[:] # Copying the varible to keep original one
# print(self.infoTable,'before')
# RepairKHP operator
Op = Operators()
self.infoTable = Op.RepairKHP(self.infoTable)
# print(self.infoTable,'after')
self.molecules, self.stemPool, self.infoEnergy, self.moleculeEnergy, self.moleculeTable, self.basePairs, self.infoTable, self.moleculeShort, self.elements = population.GenerateMolecule(
sequence, len(sequence), popSize, self.infoTable)
#population.PrintInfo(molecule,stemPool,infoEnergy,moleculeEnergy)
self.PE = self.moleculeEnergy
self.PE1 = self.moleculeEnergy
for i in range(len(self.moleculeTable)):
self.KE.append(initialKE)
self.KE1.append(initialKE)
self.numHit.append(0)
self.minStruct.append(self.moleculeTable[i])
self.minPE.append(self.moleculeEnergy[i])
self.minHit.append(0)
#endfor
#end
#endclass
#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
# IMPORTANT
#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
# Need to uniquify: moleculeTable
