def find_weight_match(current_w, w_list): for weight in w_list: for item in pw_dict().items(): if abs(item[1] - (weight - current)) < 0.01: return item[0] return -1
return item[0]
return -1
if __name__ == '__main__':
with open('data/rosalind_full.txt') as input_data:
weights = [float(line.strip()) for line in input_data.readlines()]
# Given that len(weights) = 2n+3
n = (len(weights)-3)/2
# Initialize Variables
protein = ''
current = weights[1]
myw = [w for w in weights[2:]]
# Iteratively build the protein.
while len(protein) < n:
temp = find_weight_match(current, myw)
if temp == -1:
break
else:
protein += temp
current += pw_dict()[temp]
myw = filter(lambda w: w-current > 0, myw)
# Print and save the output.
print protein
with open('output/059_FULL.txt', 'w') as output_file:
output_file.write(protein)
