def test_break_tri3_to_tri6_sierra(self):
fixture_setup()
pm = MPI.COMM_WORLD
p_size = parallel_machine_size(pm)
print p_size, "++++++++++++++++++++"
for i in range(100):
print "-------------------"
if p_size <= 3:
n = 12
nx = n
ny = n
createEdgeSets = True
fixture = QuadFixture_3(pm, nx, ny, createEdgeSets)
isCommited = False
eMesh = PerceptMesh(fixture.meta_data, fixture.bulk_data, isCommited)
scalarDimension = 0
proc_rank_field = eMesh.add_field("proc_rank", eMesh.element_rank(), scalarDimension)
eMesh.add_field("proc_rank_ede", eMesh.edge_rank(), scalarDimension)
breaker = Refiner(eMesh, TRI3_TRI6_1, proc_rank_field)
eMesh.commit()
fixture.generate_mesh()
eMesh.print_info("tri mesh tri6")
eMesh.save_as("quad_fixture_tri3_tri6_0.e")
breaker.doBreak()
eMesh.print_info("tri mesh enriched")
eMesh.save_as("quad_fixture_tri6_tri6_0.e")
eMesh.save_as("quad_fixture_tri3_tri6_1.e")
def test_break_tri3_to_tri6_sierra(self):
fixture_setup()
pm = MPI.COMM_WORLD
p_size = parallel_machine_size(pm)
print p_size, "++++++++++++++++++++"
for i in range(100):
print "-------------------"
if p_size <= 3:
n = 12
nx = n
ny = n
createEdgeSets = True
fixture = QuadFixture_3(pm, nx, ny, createEdgeSets)
isCommited = False
eMesh = PerceptMesh(fixture.meta_data, fixture.bulk_data,
isCommited)
scalarDimension = 0
proc_rank_field = eMesh.add_field("proc_rank",
eMesh.element_rank(),
scalarDimension)
eMesh.add_field("proc_rank_ede", eMesh.edge_rank(),
scalarDimension)
breaker = Refiner(eMesh, TRI3_TRI6_1, proc_rank_field)
eMesh.commit()
fixture.generate_mesh()
eMesh.print_info("tri mesh tri6")
eMesh.save_as("quad_fixture_tri3_tri6_0.e")
breaker.doBreak()
eMesh.print_info("tri mesh enriched")
eMesh.save_as("quad_fixture_tri6_tri6_0.e")
eMesh.save_as("quad_fixture_tri3_tri6_1.e")
def test_break_tri3_to_tri6_to_tri6_sierra(self):
fixture_setup()
pm = MPI.COMM_WORLD
p_size = parallel_machine_size(pm)
if p_size <= 3:
eMesh = PerceptMesh(2)
eMesh.open("quad_fixture_tri6_tri6_0.e")
scalarDimension = 0
proc_rank_field = eMesh.add_field("proc_rank", eMesh.element_rank(), scalarDimension)
eMesh.add_field("proc_rank_ede", eMesh.edge_rank(), scalarDimension)
breaker = Refiner(eMesh, TRI6_TRI6_4, proc_rank_field)
eMesh.commit()
eMesh.print_info("tri mesh tri6")
eMesh.save_as("quad_fixture_tri6_tri6_0.e")
breaker.doBreak()
eMesh.print_info("tri mesh refined")
eMesh.save_as("quad_fixture_tri6_tri6_1.e")
def test_break_tri3_to_tri6_to_tri6_sierra(self):
fixture_setup()
pm = MPI.COMM_WORLD
p_size = parallel_machine_size(pm)
if p_size <= 3:
eMesh = PerceptMesh(2)
eMesh.open("quad_fixture_tri6_tri6_0.e")
scalarDimension = 0
proc_rank_field = eMesh.add_field("proc_rank", eMesh.element_rank(), scalarDimension)
eMesh.add_field("proc_rank_ede", eMesh.edge_rank(), scalarDimension)
breaker = Refiner(eMesh, TRI6_TRI6_4, proc_rank_field)
eMesh.commit()
eMesh.print_info("tri mesh tri6")
eMesh.save_as("quad_fixture_tri6_tri6_0.e")
breaker.doBreak()
eMesh.print_info("tri mesh refined")
eMesh.save_as("quad_fixture_tri6_tri6_1.e")
def test_break_tri_to_tri_sierra_1_test(self):
pm = MPI.COMM_WORLD
p_size = parallel_machine_size(pm)
if p_size <= 2:
n = 2
nx = n
ny = n
createEdgeSets = False
fixture = QuadFixture_4(pm,nx,ny,createEdgeSets)
isCommited = False
eMesh = PerceptMesh(fixture.meta_data, fixture.bulk_data, isCommited)
scalarDimension = 0
proc_rank_field = eMesh.add_field("proc_rank", eMesh.element_rank(), scalarDimension)
eMesh.add_field("proc_rank_edge", eMesh.edge_rank(), scalarDimension)
breaker = Refiner(eMesh, TRI3_TRI3_4, proc_rank_field)
eMesh.commit()
fixture.generate_mesh()
breaker.setRemoveOldElements(False)
breaker.doBreak()
eMesh.dump_elements_compact()
def test_break_tri_to_tri_sierra_1_test(self):
pm = MPI.COMM_WORLD
p_size = parallel_machine_size(pm)
if p_size <= 2:
n = 2
nx = n
ny = n
createEdgeSets = False
fixture = QuadFixture_4(pm,nx,ny,createEdgeSets)
isCommited = False
eMesh = PerceptMesh(fixture.meta_data, fixture.bulk_data, isCommited)
scalarDimension = 0
proc_rank_field = eMesh.add_field("proc_rank", eMesh.element_rank(), scalarDimension)
eMesh.add_field("proc_rank_edge", eMesh.edge_rank(), scalarDimension)
breaker = Refiner(eMesh, TRI3_TRI3_4, proc_rank_field)
eMesh.commit()
fixture.generate_mesh()
breaker.setRemoveOldElements(False)
breaker.doBreak()
eMesh.dump_elements_compact()
def test_break_tri_to_tri_sierra_1(self):
fixture_setup()
pm = MPI.COMM_WORLD
p_size = parallel_machine_size(pm)
if p_size <= 3:
n = 12
nx = n
ny = n
createEdgeSets = True
fixture = QuadFixture_3(pm, nx, ny, createEdgeSets)
isCommited = False
eMesh = PerceptMesh(fixture.meta_data, fixture.bulk_data, isCommited)
scalarDimension = 0
proc_rank_field = eMesh.add_field("proc_rank", stk::topology::ELEMENT_RANK, scalarDimension)
eMesh.add_field("proc_rank_ede", eMesh.edge_rank(), scalarDimension)
breaker = Refiner(eMesh, TRI3_TRI3_4, proc_rank_field)
eMesh.commit()
fixture.generate_mesh()
eMesh.print_info("tri mesh")
eMesh.save_as("quad_fixture_tri3_0.e")
breaker.doBreak()
eMesh.print_info("tri mesh refined")
eMesh.save_as("quad_fixture_tri3_1.e")
def test_perceptMesh_walk_nodes(self):
self.fixture_setup()
p_size = self.pm.size
p_rank = self.pm.rank
if p_size <= 2:
n = 12
nx = n
ny = n
sidesets_on = 1
fixture = QuadFixture_4(self.pm, nx, ny, sidesets_on)
fixture.meta_data.commit()
fixture.generate_mesh()
eMesh = PerceptMesh(fixture.meta_data, fixture.bulk_data)
eMesh.print_info("quad fixture", 2)
metaData = eMesh.get_fem_meta_data()
parts = metaData.get_parts()
nparts = len(parts)
print "Number of parts = ", nparts
surface_id = 2
surface_name = "surface_" + str(surface_id)
part = eMesh.get_part(surface_name)
in_surface_selector = Selector(part)
bulkData = eMesh.get_bulk_data()
coordField = eMesh.get_coordinates_field()
if eMesh.get_spatial_dim() == 2:
buckets_arg = eMesh.edge_rank()
else:
buckets_arg = eMesh.face_rank
buckets = bulkData.buckets(buckets_arg)
sum = 0.0
def test_perceptMesh_walk_nodes(self):
self.fixture_setup()
p_size = self.pm.size
p_rank = self.pm.rank
if p_size <= 2:
n = 12
nx = n
ny = n
sidesets_on = 1
fixture = QuadFixture_4(self.pm, nx, ny, sidesets_on)
fixture.meta_data.commit()
fixture.generate_mesh()
eMesh = PerceptMesh(fixture.meta_data, fixture.bulk_data)
eMesh.print_info("quad fixture", 2)
metaData = eMesh.get_fem_meta_data()
parts = metaData.get_parts()
nparts = len(parts)
print "Number of parts = ", nparts
surface_id = 2
surface_name = "surface_" + str(surface_id)
part = eMesh.get_part(surface_name)
in_surface_selector = Selector(part)
bulkData = eMesh.get_bulk_data()
coordField = eMesh.get_coordinates_field()
if eMesh.get_spatial_dim() == 2:
buckets_arg = eMesh.edge_rank()
else:
buckets_arg = eMesh.face_rank
buckets = bulkData.buckets(buckets_arg)
sum = 0.0
