|
| sys | examples.uc.read_mesh1d.plot_mesh = sys.argv[1] |
| |
| h5py | examples.uc.read_mesh1d.f = h5py.File(plot_mesh) |
| |
| np | examples.uc.read_mesh1d.coordinates = np.asarray(f["domain/coordinates"]) |
| |
| np | examples.uc.read_mesh1d.subelement_connectivity = np.asarray(f["domain/connectivity"]) |
| |
| np | examples.uc.read_mesh1d.element_connectivity = np.asarray(f["domain/element_connectivity"]) |
| |
| h5py | examples.uc.read_mesh1d.patchNames = f["domain"].attrs["patchNames"] |
| |
| h5py | examples.uc.read_mesh1d.patchIndexes = f["domain"].attrs["patchIndexes"] |
| |
| h5py | examples.uc.read_mesh1d.startingElementIndexes = f["domain"].attrs["startingElementIndexes"] |
| |
| np | examples.uc.read_mesh1d.numcoords = coordinates.shape[0] |
| |
| np | examples.uc.read_mesh1d.numsubelems = subelement_connectivity.shape[0] |
| |
| np | examples.uc.read_mesh1d.numelems = element_connectivity.shape[0] |
| |
| len | examples.uc.read_mesh1d.numpatches = len(patchIndexes) |
| |
| np | examples.uc.read_mesh1d.x = coordinates[:,0] |
| |
| np | examples.uc.read_mesh1d.y = coordinates[:,1] |
| |
| np | examples.uc.read_mesh1d.gp_subelements = coordinates[subelement_connectivity] |
| |
| np | examples.uc.read_mesh1d.gp_elements = coordinates[element_connectivity] |
| |
| np | examples.uc.read_mesh1d.Np = gp_elements.shape[1] |
| |
| np | examples.uc.read_mesh1d.patchMarkers = np.zeros(x.shape) |
| |
| np | examples.uc.read_mesh1d.element_patches = patchMarkers[element_connectivity] |
| |
| np | examples.uc.read_mesh1d.elem_start = gp_elements[:,0,0] |
| |
| np | examples.uc.read_mesh1d.elem_end = gp_elements[:,-1,0] |
| |
| np | examples.uc.read_mesh1d.elems = gp_elements[:,:,0] |
| |
| np | examples.uc.read_mesh1d.elem_edges = np.vstack((elem_start, elem_end)).T |
| |
| plt | examples.uc.read_mesh1d.cmap = plt.cm.get_cmap('rainbow',len(patchIndexes)) |
| |
| | examples.uc.read_mesh1d.color |
| |
| | examples.uc.read_mesh1d.label |
| |
| | examples.uc.read_mesh1d.s |
| |
1.9.6