def test_break_quad_to_quad_sierra_2(self): fixture_setup() pm = MPI.COMM_WORLD p_size = parallel_machine_size(pm) doGenSideSets = True if p_size <= 3: n = 12 nx = n ny = n fixture = QuadFixture_4(pm, nx, ny, doGenSideSets) isCommited = False eMesh = PerceptMesh(fixture.meta_data, fixture.bulk_data, isCommited) eMesh.commit() fixture.generate_mesh() eMesh.save_as("quad_fixture_mbreak_0.e") eMesh.close() eMesh1 = PerceptMesh(2) eMesh1.open("quad_fixture_mbreak_0.e") scalarDimension = 0 proc_rank_field = eMesh1.add_field("proc_rank", eMesh.element_rank(), scalarDimension) breaker = Refiner(eMesh1, QUAD4_QUAD4_4_SIERRA, proc_rank_field) eMesh1.commit() i = 0 while i < 2: print "\n\n\n ================ tmp Refine Pass = "******"quad_fixture_mbreak_" + str(i) + ".e") i = i + 1
def test_break_quad_to_quad_sierra_2(self): fixture_setup() pm = MPI.COMM_WORLD p_size = parallel_machine_size(pm) doGenSideSets = True if p_size <= 3: n = 12 nx = n ny = n fixture = QuadFixture_4(pm,nx,ny,doGenSideSets) isCommited = False eMesh = PerceptMesh(fixture.meta_data, fixture.bulk_data, isCommited) eMesh.commit() fixture.generate_mesh() eMesh.save_as("quad_fixture_mbreak_0.e") eMesh.close() eMesh1 = PerceptMesh(2) eMesh1.open("quad_fixture_mbreak_0.e") scalarDimension = 0 proc_rank_field = eMesh1.add_field("proc_rank", eMesh.element_rank(), scalarDimension) breaker = Refiner(eMesh1, QUAD4_QUAD4_4_SIERRA, proc_rank_field) eMesh1.commit() i = 0 while i < 2: print "\n\n\n ================ tmp Refine Pass = "******"quad_fixture_mbreak_" + str(i) + ".e") i = i + 1
def test_heterogeneous_mesh(self): fixture_setup() pm = MPI.COMM_WORLD p_size = parallel_machine_size(MPI.COMM_WORLD) if p_size <= 1: mesh = HeterogeneousFixture(MPI.COMM_WORLD, False) #put_io_part_attribute(mesh.m_block_hex) #put_io_part_attribute(mesh.m_block_wedge) #put_io_part_attribute(mesh.m_block_tet) mesh.m_metaData.commit() mesh.populate() isCommited = True em1 = PerceptMesh(mesh.m_metaData, mesh.m_bulkData, isCommited) em1.save_as("heterogeneous_0.e") em1.close() eMesh = PerceptMesh(3) eMesh.open("heterogeneous_0.e") scalarDimension = 0 proc_rank_field = eMesh.add_field("proc_rank", eMesh.element_rank(), scalarDimension)