def run():
infileName = os.path.join(pkg_path, 'converters', 'panair', 'M100', 'M100.inp')
model = PanairGrid(log=None, debug=True)
model.read_panair(infileName)
p3d_name = 'M100.plt'
model.write_plot3d(p3d_name)
def run():
infilename = os.path.join(PKG_PATH, 'converters', 'panair', 'M100',
'M100.inp')
model = PanairGrid(log=None, debug=True)
model.read_panair(infilename)
p3d_name = 'M100.plt'
model.write_plot3d(p3d_name)
def test_panair_io_03(self):
"""tests the SWB model"""
log = get_logger(level='warning')
in_filename = os.path.join(TEST_PATH, 'SWB', 'SWB.inp')
#out_filename = os.path.join(TEST_PATH, 'M100', 'M100_out.inp')
model = PanairGrid(log=log, debug=False)
model.read_panair(in_filename)
(points, elements, regions, kt, cp_nrom) = model.get_points_elements_regions()
model.write_panair('junk_swb.inp')
os.remove('junk_swb.inp')
def test_panair_io_01(self):
in_filename = os.path.join(test_path, 'M100', 'M100.inp')
out_filename = os.path.join(test_path, 'M100', 'M100_out.inp')
#with open(infile_name, 'w') as f:
#f.write(lines)
model = PanairGrid(log=None, debug=False)
model.read_panair(in_filename)
#assert len(cart3d.points) == 7, 'npoints=%s' % len(cart3d.points)
#assert len(cart3d.elements) == 6, 'nelements=%s' % len(cart3d.elements)
#assert len(cart3d.regions) == 6, 'nregions=%s' % len(cart3d.regions)
#assert len(cart3d.loads) == 0, 'nloads=%s' % len(cart3d.loads)
model.write_panair(out_filename)
(points, elements, regions, kt, cp_nrom) = model.get_points_elements_regions()
os.remove(out_filename)
def test_panair_io_01(self):
log = get_logger(level='warning')
in_filename = os.path.join(test_path, 'M100', 'M100.inp')
out_filename = os.path.join(test_path, 'M100', 'M100_out.inp')
#with open(infile_name, 'w') as f:
#f.write(lines)
model = PanairGrid(log=log, debug=False)
model.read_panair(in_filename)
#assert len(cart3d.points) == 7, 'npoints=%s' % len(cart3d.points)
#assert len(cart3d.elements) == 6, 'nelements=%s' % len(cart3d.elements)
#assert len(cart3d.regions) == 6, 'nregions=%s' % len(cart3d.regions)
#assert len(cart3d.loads) == 0, 'nloads=%s' % len(cart3d.loads)
model.write_panair(out_filename)
(points, elements, regions, kt,
cp_nrom) = model.get_points_elements_regions()
os.remove(out_filename)
def test_panair_io_01(self):
"""test the M100 model"""
log = get_logger(level='warning')
in_filename = os.path.join(TEST_PATH, 'M100', 'M100.inp')
out_filename = os.path.join(TEST_PATH, 'M100', 'M100_out.inp')
#with open(infile_name, 'w') as f:
#f.write(lines)
model = PanairGrid(log=log, debug=False)
model.read_panair(in_filename)
model.write_panair(out_filename)
(points, elements, regions, kt, cp_nrom) = model.get_points_elements_regions()
model.write_panair('junk_m100.inp')
os.remove('junk_m100.inp')
model.print_options()
model.print_abutments()
model.print_grid_summary()
model.print_out_header()
os.remove(out_filename)
def load_panair_geometry(self,
panair_filename,
dirname,
name='main',
plot=True):
self.nid_map = {}
#key = self.case_keys[self.icase]
#case = self.result_cases[key]
skip_reading = self._remove_old_geometry(panair_filename)
if skip_reading:
return
model = PanairGrid(log=self.log, debug=self.debug)
self.model_type = model.model_type
model.read_panair(panair_filename)
nodes, elements, regions, kt, cp_norm = model.get_points_elements_regions(
)
#for nid, node in enumerate(nodes):
#print "node[%s] = %s" % (nid, str(node))
self.nNodes = len(nodes)
self.nElements = len(elements)
#print("nNodes = ",self.nNodes)
#print("nElements = ", self.nElements)
self.grid.Allocate(self.nElements, 1000)
#self.gridResult.SetNumberOfComponents(self.nElements)
points = vtk.vtkPoints()
points.SetNumberOfPoints(self.nNodes)
#self.gridResult.Allocate(self.nNodes, 1000)
#vectorReselt.SetNumberOfComponents(3)
#elem.SetNumberOfPoints(nNodes)
if 0:
fraction = 1. / nnodes # so you can color the nodes by ID
for nid, node in sorted(iteritems(nodes)):
points.InsertPoint(nid - 1, *point)
self.gridResult.InsertNextValue(nid * fraction)
#print str(element)
#elem = vtk.vtkVertex()
#elem.GetPointIds().SetId(0, i)
#self.aQuadGrid.InsertNextCell(elem.GetCellType(), elem.GetPointIds())
#vectorResult.InsertTuple3(0, 0.0, 0.0, 1.0)
assert len(nodes) > 0
mmax = amax(nodes, axis=0)
mmin = amin(nodes, axis=0)
dim_max = (mmax - mmin).max()
self.create_global_axes(dim_max)
for nid, node in enumerate(nodes):
points.InsertPoint(nid, *node)
assert len(elements) > 0
elem = vtkQuad()
quad_type = elem.GetCellType()
for eid, element in enumerate(elements):
(p1, p2, p3, p4) = element
elem = vtkQuad()
elem.GetPointIds().SetId(0, p1)
elem.GetPointIds().SetId(1, p2)
elem.GetPointIds().SetId(2, p3)
elem.GetPointIds().SetId(3, p4)
self.grid.InsertNextCell(quad_type, elem.GetPointIds())
self.grid.SetPoints(points)
#self.grid.GetPointData().SetScalars(self.gridResult)
#print dir(self.grid) #.SetNumberOfComponents(0)
#self.grid.GetCellData().SetNumberOfTuples(1);
#self.grid.GetCellData().SetScalars(self.gridResult)
self.grid.Modified()
if hasattr(self.grid, 'Update'):
self.grid.Update()
# loadPanairResults - regions/loads
if plot:
self.scalarBar.VisibilityOn()
self.scalarBar.Modified()
self.iSubcaseNameMap = {1: ['Panair', '']}
cases = {}
ID = 1
#print "nElements = ", nElements
loads = []
form, cases = self._fill_panair_geometry_case(cases, ID, nodes,
elements, regions, kt,
cp_norm, loads)
#if plot:
self._finish_results_io2(form, cases)
def load_panair_geometry(self, panair_filename, name='main', plot=True):
model_name = name
self.gui.nid_map = {}
#key = self.case_keys[self.icase]
#case = self.result_cases[key]
skip_reading = self.gui._remove_old_geometry(panair_filename)
if skip_reading:
return
model = PanairGrid(log=self.gui.log, debug=self.gui.debug)
self.gui.model_type = model.model_type
model.read_panair(panair_filename)
self.gui.geom_model = model
# get_wakes=True shows explicit wakes
#
# TODO: bad for results...what about just not adding it to the patches/bcs?
nodes, elements, regions, kt, cp_norm = model.get_points_elements_regions(
get_wakes=True)
self.gui.nnodes = len(nodes)
self.gui.nelements = len(elements)
#print("nnodes = ",self.nnodes)
#print("nelements = ", self.nelements)
grid = self.gui.grid
grid.Allocate(self.gui.nelements, 1000)
points = vtk.vtkPoints()
points.SetNumberOfPoints(self.gui.nnodes)
assert len(nodes) > 0
mmax = amax(nodes, axis=0)
mmin = amin(nodes, axis=0)
dim_max = (mmax - mmin).max()
self.gui.create_global_axes(dim_max)
points = numpy_to_vtk_points(nodes)
assert len(elements) > 0
elem = vtkQuad()
quad_type = elem.GetCellType()
create_vtk_cells_of_constant_element_type(grid, elements, quad_type)
grid.SetPoints(points)
grid.Modified()
# loadPanairResults - regions/loads
if plot:
self.gui.scalar_bar_actor.VisibilityOn()
self.gui.scalar_bar_actor.Modified()
self.gui.isubcase_name_map = {1: ['Panair', '']}
cases = OrderedDict()
ID = 1
loads = []
form, cases, node_ids, element_ids = self._fill_panair_geometry_case(
cases, ID, nodes, elements, regions, kt, cp_norm, loads)
self.gui.node_ids = node_ids
self.gui.element_ids = element_ids
#if plot:
self.gui._finish_results_io2(model_name, form, cases)
def load_panair_geometry(self, panair_filename, name='main', plot=True):
self.gui.nid_map = {}
#key = self.case_keys[self.icase]
#case = self.result_cases[key]
skip_reading = self.gui._remove_old_geometry(panair_filename)
if skip_reading:
return
model = PanairGrid(log=self.gui.log, debug=self.gui.debug)
self.gui.model_type = model.model_type
model.read_panair(panair_filename)
nodes, elements, regions, kt, cp_norm = model.get_points_elements_regions(
)
#for nid, node in enumerate(nodes):
#print "node[%s] = %s" % (nid, str(node))
self.gui.nnodes = len(nodes)
self.gui.nelements = len(elements)
#print("nnodes = ",self.nnodes)
#print("nelements = ", self.nelements)
grid = self.gui.grid
grid.Allocate(self.gui.nelements, 1000)
points = vtk.vtkPoints()
points.SetNumberOfPoints(self.gui.nnodes)
assert len(nodes) > 0
mmax = amax(nodes, axis=0)
mmin = amin(nodes, axis=0)
dim_max = (mmax - mmin).max()
self.gui.create_global_axes(dim_max)
for nid, node in enumerate(nodes):
points.InsertPoint(nid, *node)
assert len(elements) > 0
elem = vtkQuad()
quad_type = elem.GetCellType()
for eid, element in enumerate(elements):
(p1, p2, p3, p4) = element
elem = vtkQuad()
elem.GetPointIds().SetId(0, p1)
elem.GetPointIds().SetId(1, p2)
elem.GetPointIds().SetId(2, p3)
elem.GetPointIds().SetId(3, p4)
grid.InsertNextCell(quad_type, elem.GetPointIds())
grid.SetPoints(points)
grid.Modified()
if hasattr(grid, 'Update'):
grid.Update()
# loadPanairResults - regions/loads
if plot:
self.gui.scalarBar.VisibilityOn()
self.gui.scalarBar.Modified()
self.gui.isubcase_name_map = {1: ['Panair', '']}
cases = OrderedDict()
ID = 1
#print "nElements = ", nElements
loads = []
form, cases, node_ids, element_ids = self._fill_panair_geometry_case(
cases, ID, nodes, elements, regions, kt, cp_norm, loads)
self.gui.node_ids = node_ids
self.gui.element_ids = element_ids
#if plot:
self.gui._finish_results_io2(form, cases)
def load_panair_geometry(self, panair_filename, dirname, name='main', plot=True):
self.nid_map = {}
#key = self.case_keys[self.icase]
#case = self.result_cases[key]
skip_reading = self._remove_old_geometry(panair_filename)
if skip_reading:
return
model = PanairGrid(log=self.log, debug=self.debug)
self.model_type = model.model_type
model.read_panair(panair_filename)
nodes, elements, regions, kt, cp_norm = model.get_points_elements_regions()
#for nid, node in enumerate(nodes):
#print "node[%s] = %s" % (nid, str(node))
self.nNodes = len(nodes)
self.nElements = len(elements)
#print("nNodes = ",self.nNodes)
#print("nElements = ", self.nElements)
self.grid.Allocate(self.nElements, 1000)
#self.gridResult.SetNumberOfComponents(self.nElements)
points = vtk.vtkPoints()
points.SetNumberOfPoints(self.nNodes)
#self.gridResult.Allocate(self.nNodes, 1000)
#vectorReselt.SetNumberOfComponents(3)
#elem.SetNumberOfPoints(nNodes)
if 0:
fraction = 1. / nnodes # so you can color the nodes by ID
for nid, node in sorted(iteritems(nodes)):
points.InsertPoint(nid - 1, *point)
self.gridResult.InsertNextValue(nid * fraction)
#print str(element)
#elem = vtk.vtkVertex()
#elem.GetPointIds().SetId(0, i)
#self.aQuadGrid.InsertNextCell(elem.GetCellType(), elem.GetPointIds())
#vectorResult.InsertTuple3(0, 0.0, 0.0, 1.0)
assert len(nodes) > 0
mmax = amax(nodes, axis=0)
mmin = amin(nodes, axis=0)
dim_max = (mmax - mmin).max()
self.create_global_axes(dim_max)
for nid, node in enumerate(nodes):
points.InsertPoint(nid, *node)
assert len(elements) > 0
elem = vtkQuad()
quad_type = elem.GetCellType()
for eid, element in enumerate(elements):
(p1, p2, p3, p4) = element
elem = vtkQuad()
elem.GetPointIds().SetId(0, p1)
elem.GetPointIds().SetId(1, p2)
elem.GetPointIds().SetId(2, p3)
elem.GetPointIds().SetId(3, p4)
self.grid.InsertNextCell(quad_type, elem.GetPointIds())
self.grid.SetPoints(points)
#self.grid.GetPointData().SetScalars(self.gridResult)
#print dir(self.grid) #.SetNumberOfComponents(0)
#self.grid.GetCellData().SetNumberOfTuples(1);
#self.grid.GetCellData().SetScalars(self.gridResult)
self.grid.Modified()
if hasattr(self.grid, 'Update'):
self.grid.Update()
# loadPanairResults - regions/loads
self. turn_text_on()
if plot:
self.scalarBar.VisibilityOn()
self.scalarBar.Modified()
self.iSubcaseNameMap = {1: ['Panair', '']}
cases = {}
ID = 1
#print "nElements = ", nElements
loads = []
form, cases = self._fill_panair_geometry_case(cases, ID, nodes, elements, regions, kt, cp_norm, loads)
#if plot:
self._finish_results_io2(form, cases)