Source code for show_mesh_info

#!/usr/bin/env python
"""
Print various information about a mesh.
"""
from argparse import RawDescriptionHelpFormatter, ArgumentParser

from sfepy.base.base import output
from sfepy.discrete.fem import Mesh

helps = {
    'filename' :
    'mesh file name',
}

[docs]def main(): parser = ArgumentParser(description=__doc__.rstrip(), formatter_class=RawDescriptionHelpFormatter) parser.add_argument('filename', help=helps['filename']) options = parser.parse_args() mesh = Mesh.from_file(options.filename) output(mesh.cmesh) output('element types:', mesh.descs) output('nodal BCs:', mesh.nodal_bcs) bbox = mesh.get_bounding_box() output('bounding box: %s' % ', '.join('%s: [%s, %s]' % (name, bbox[0, ii], bbox[1, ii]) for ii, name in enumerate('xyz'[:mesh.dim]))) output('centre:', mesh.coors.mean(0))
if __name__ == '__main__': main()