def main(filename): """ Takes a spacetime file, inverts the coordinates, and outputs an inverted file. """ spacetime_data = read_3simplices_from_data_file(filename) new_spacetime = invert_time_for_spacetime(spacetime_data) print_output_file(spacetime_data[0],new_spacetime,filename)
def main(filename): """ Takes a spacetime file, inverts the coordinates, and outputs an inverted file. """ spacetime_data = read_3simplices_from_data_file(filename) new_spacetime = invert_time_for_spacetime(spacetime_data) print_output_file(spacetime_data[0], new_spacetime, filename)
def extract_2_volume_ensemble(filename_list): """ Reads in each *.3sx2p1 file in filename_list and produces a list of lists of 2-volume information as a function of proper time. Returns a vector. The first element contains header data, which is assumed to be the same for all elements of the ensemble. """ # Do the first spacetime seperately to extract header data first_spacetime = vs.read_3simplices_from_data_file(filename_list[0]) header_data = first_spacetime[0] spacetimes = [first_spacetime[1]] # We need to calculate the average of N3-31 and N3-13 for all spacetimes N3 = eval(header_data[11])/2 # Extract the remaining spacetimes for f in filename_list[1:]: temp = vs.read_3simplices_from_data_file(f) spacetimes.append(temp[1]) N3 += eval(temp[0][11])/2 # Extract the 2-simplex information sl2simplices = [vs.get_all_sl2simplices([header_data,s]) \ for s in spacetimes] # Run garbage collection del spacetimes # Average the N3s instead of summing them N3 = N3/len(filename_list) # Extract volumes as a function of time volumes = [vs.make_v_of_t(s) for s in sl2simplices] # Delete the old lists del sl2simplices return [header_data,volumes,N3]
def extract_2_volume_ensemble(filename_list): """ Reads in each *.3sx2p1 file in filename_list and produces a list of lists of 2-volume information as a function of proper time. Returns a vector. The first element contains header data, which is assumed to be the same for all elements of the ensemble. """ # Do the first spacetime seperately to extract header data first_spacetime = vs.read_3simplices_from_data_file(filename_list[0]) header_data = first_spacetime[0] spacetimes = [first_spacetime[1]] # We need to calculate the average of N3-31 and N3-13 for all spacetimes N3 = eval(header_data[11]) / 2 # Extract the remaining spacetimes for f in filename_list[1:]: temp = vs.read_3simplices_from_data_file(f) spacetimes.append(temp[1]) N3 += eval(temp[0][11]) / 2 # Extract the 2-simplex information sl2simplices = [vs.get_all_sl2simplices([header_data,s]) \ for s in spacetimes] # Run garbage collection del spacetimes # Average the N3s instead of summing them N3 = N3 / len(filename_list) # Extract volumes as a function of time volumes = [vs.make_v_of_t(s) for s in sl2simplices] # Delete the old lists del sl2simplices return [header_data, volumes, N3]
def get_raw_time_slices(filename, t_low, t_high): """ Given the filename/filepath of a *.3sx2p1 file, extract the 2-simplices for the high and low simplices we'll use to make our boundary conditions as a pair of lists (in a list). """ # Load spacetime spacetime = vs.read_3simplices_from_data_file(filename) # Assign time slices t_slice_low = vs.get_spacelike_2simplices_at_t(spacetime[1], t_low) t_slice_high = vs.get_spacelike_2simplices_at_t(spacetime[1], t_high) # spacetime[0] contains coupling constant and spacetime # size. Useful for file naming conventions. return [t_slice_low, t_slice_high, spacetime[0]]
def get_raw_time_slices(filename,t_low,t_high): """ Given the filename/filepath of a *.3sx2p1 file, extract the 2-simplices for the high and low simplices we'll use to make our boundary conditions as a pair of lists (in a list). """ # Load spacetime spacetime = vs.read_3simplices_from_data_file(filename) # Assign time slices t_slice_low = vs.get_spacelike_2simplices_at_t(spacetime[1],t_low) t_slice_high = vs.get_spacelike_2simplices_at_t(spacetime[1],t_high) # spacetime[0] contains coupling constant and spacetime # size. Useful for file naming conventions. return [t_slice_low, t_slice_high,spacetime[0]]