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")
np.savetxt("7eles.txt", els3_array, fmt="%d")
np.savetxt("8eles.txt", els4_array, fmt="%d")
file_list = glob.glob("*.txt")
with open('transfer.inp', 'w') as file:
input_lines = fileinput.input(file_list)
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")
with open('point.inp', 'w') as file:
input_lines = fileinput.input(file_list)
file.writelines(input_lines)