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