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
Example #2
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