return xyz,direction,speed
(xyz,direction,speed) = load_input_file("input.txt")
mesh = Mesh(mesh="../geom/conv_pipe_mesh.h5m")
# load dag geometry
load("../geom/conv_pipe_geom.h5m")
# discretise geom onto mesh
cell_fracs = discretize_geom(mesh)
mesh.vol_id = IMeshTag(1,int)
mesh.source_strength = IMeshTag(1,float)
id_list = []
for element in cell_fracs:
id_list.append(element["cell"])
mesh.vol_id[:]=id_list
# loop along path, use point in volume to determine the vol id
volumes = get_all_cells(xyz,direction,speed)
# loop along path, use point in volume to determine the vol id
# but tag each tet in the range with right source strength
dump = get_all_cells_detailed(xyz,direction,speed,mesh)
