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_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