start = math.log(10**start_exp)
end = math.log(10**end_exp)
step = (end - start) / (n * number_of_decades - 1)
uniform_time = []
for i in range(n * number_of_decades):
uniform_time.append(b**(start + step * i))
existingfiles = glob.glob(prefix + '*')
existingfiles = filter(lambda u: u[u.index('.', 5) + 1:].isdigit(),
existingfiles)
existingfiles = [int(u[u.index('.', 5) + 1:]) for u in existingfiles]
existingfiles.sort()
uniform_time = [
pypolymer.find_closest_existing_file(existingfiles, u)
for u in uniform_time
]
uniform_time = pypolymer.unique(uniform_time)
blockfiles = [prefix + str(u) for u in uniform_time]
else:
sys.stderr.write('ERROR: blockfiles_method\n')
sys.exit()
sys.stdout.write('\n\n')
sys.stdout.write('=================================================\n')
sys.stdout.write(' S(q) \n')
sys.stdout.write('=================================================\n')
sys.stdout.write('\n\nTotal number of files: ' + str(len(blockfiles)) + '\n')
if (len(blockfiles) > 2):
sys.stdout.write(' blockfiles[0] = ' + blockfiles[0] + '\n')
sys.stdout.write(' blockfiles[1] = ' + blockfiles[1] + '\n')
# base
b = math.e
number_of_decades = end_exp - start_exp
start = math.log(10 ** start_exp)
end = math.log(10 ** end_exp)
step = (end - start) / (n * number_of_decades - 1)
uniform_time = []
for i in range(n * number_of_decades):
uniform_time.append(b ** (start + step * i))
existingfiles = glob.glob(prefix + "*")
existingfiles = filter(lambda u: u[u.index(".", 5) + 1 :].isdigit(), existingfiles)
existingfiles = [int(u[u.index(".", 5) + 1 :]) for u in existingfiles]
existingfiles.sort()
uniform_time = [pypolymer.find_closest_existing_file(existingfiles, u) for u in uniform_time]
uniform_time = pypolymer.unique(uniform_time)
blockfiles = [prefix + str(u) for u in uniform_time]
else:
sys.stderr.write("ERROR: blockfiles_method\n")
sys.exit()
sys.stdout.write("\n\n")
sys.stdout.write("=================================================\n")
sys.stdout.write(" S(q) \n")
sys.stdout.write("=================================================\n")
sys.stdout.write("\n\nTotal number of files: " + str(len(blockfiles)) + "\n")
if len(blockfiles) > 2:
sys.stdout.write(" blockfiles[0] = " + blockfiles[0] + "\n")
sys.stdout.write(" blockfiles[1] = " + blockfiles[1] + "\n")