sfepy.fem.domain module

Computational domain, consisting of the mesh and regions.

class sfepy.fem.domain.Domain(name, mesh, verbose=False)

Domain is divided into groups, whose purpose is to have homogeneous data shapes.

clear_surface_groups()

Remove surface group data.

create_region(name, select, flags=None, check_parents=True, functions=None, add_to_regions=True)

Region factory constructor. Append the new region to self.regions list.

create_regions(region_defs, functions=None)

create_surface_group(region)

Create a new surface group corresponding to region if it does not exist yet.

Notes

Surface groups define surface facet connectivity that is needed for sfepy.fem.mappings.SurfaceMapping.

fix_element_orientation()

Ensure element nodes ordering giving positive element volume.

The groups with elements of lower dimension than the space dimension are skipped.

get_cell_offsets()

get_conns()

Return the element connectivity groups of the underlying mesh.

get_diameter()

Return the diameter of the domain.

Notes

The diameter corresponds to the Friedrichs constant.

get_element_diameters(ig, cells, vg, mode, square=True)

get_facets(force_faces=False)

Return edge and face descriptions.

get_mesh_bounding_box()

Return the bounding box of the underlying mesh.

Returns:

bbox : ndarray (2, dim)

The bounding box with min. values in the first row and max. values in the second row.

get_mesh_coors(actual=False)

Return the coordinates of the underlying mesh vertices.

has_faces()

iter_groups(igs=None)

refine()

Uniformly refine the domain mesh.

Returns:

domain : Domain instance

The new domain with the refined mesh.

Notes

Works only for meshes with single element type! Does not preserve node groups!

reset_regions()

Reset the list of regions associated with the domain.

save_regions(filename, region_names=None)

Save regions as individual meshes.

Parameters:

filename : str

The output filename.

region_names : list, optional

If given, only the listed regions are saved.

save_regions_as_groups(filename, region_names=None)

Save regions in a single mesh but mark them by using different element/node group numbers.

If regions overlap, the result is undetermined, with exception of the whole domain region, which is marked by group id 0.

Region masks are also saved as scalar point data for output formats that support this.

Parameters:

filename : str

The output filename.

region_names : list, optional

If given, only the listed regions are saved.

setup_facets(create_edges=True, create_faces=True, verbose=False)

Setup the edges and faces (in 3D) of domain elements.

setup_groups()

surface_faces()

sfepy.fem.domain.region_leaf(domain, regions, rdef, functions)

Create/setup a region instance according to rdef.

sfepy.fem.domain.region_op(level, op, item1, item2)