from __future__ import division import meshio import mesh_waves as msw import numpy as np import fileinput import glob #%% points, cells, point_data, cell_data , field_data = \ meshio.read("point.msh") # nodes_array = msw.node_writer(points, point_data) print len(nodes_array) #%% # Domain elements. # nfin, els1_array = msw.ele_writer(cells, cell_data, 'quad9', 30000, 3, 1, 18, 9, 0) nini = nfin # #%% # Absorbing boundaries # nfin, els2_array = msw.boundary_writer(cells, cell_data, 100, 2, nini) # #%% # np.savetxt("2nodes.txt", nodes_array, fmt=("%d", "%d", "%d", "%d", "%.4f", "%.4f")) np.savetxt("5eles.txt", els1_array, fmt="%d") np.savetxt("6eles.txt", els2_array, fmt="%d") file_list = glob.glob("*.txt")
""" from __future__ import division import meshio import mesh_waves as msw import numpy as np import fileinput import glob #%% points, cells, point_data, cell_data , field_data = \ meshio.read("transfer.msh") # nodes_array = msw.node_writer(points, point_data) print len(nodes_array) # #%% nfin, els1_array = msw.ele_writer(cells, cell_data, 'quad9', 11000, 3, 5, 18, 9, 0) nini = nfin nfin, els2_array = msw.ele_writer(cells, cell_data, 'quad9', 12000, 3, 6, 18, 9, nini) nini = nfin nfin, els3_array = msw.face_recognition(cells, cell_data, 10000, 5000, nini) nini = nfin nfin, els4_array = msw.boundary_writer(cells, cell_data, 7000, 6, nini) # #%% # np.savetxt("2nodes.txt", nodes_array, fmt=("%d", "%d", "%d", "%d", "%.4f", "%.4f")) np.savetxt("5eles.txt", els1_array, fmt="%d") np.savetxt("6eles.txt", els2_array, fmt="%d")
from __future__ import division import meshio import mesh_waves as msw import numpy as np import fileinput import glob #%% points, cells, point_data, cell_data , field_data = \ meshio.read("estrato.msh") # nodes_array = msw.node_writer(points, point_data) print len(nodes_array) #%% # Domain elements. # nfin, els1_array = msw.ele_writer(cells, cell_data, 'quad9', 20000, 3, 1, 18, 9, 0) nini = nfin nfin, els2_array = msw.ele_writer(cells, cell_data, 'quad9', 10000, 3, 2, 18, 9, nini) nini = nfin # #%% # Absorbing boundaries # nfin, els3_array = msw.boundary_writer(cells, cell_data, 100, 3, nini) # #%% # np.savetxt("2nodes.txt", nodes_array, fmt=("%d", "%d", "%d", "%d", "%.4f", "%.4f"))