Geometric statistics

meshioplusplus.compute_stats(mesh) reports geometric statistics of a mesh — the geometric complement to the topological info verb. It computes the bounding box + extents, centroid, per-cell-type counts, total surface area (of 2D cells plus the boundary of 3D cells), total signed & unsigned volume of 3D cells, and the count of inverted (negative signed-volume) cells — without modifying the mesh. It is a mesh operation (like quality metrics), not a file format, and uses only standard C++/numpy, so it runs under every mesh backend.

import meshioplusplus

mesh = meshioplusplus.read("part.vtu")
s = meshioplusplus.compute_stats(mesh)

print(s["bbox_min"], s["bbox_max"], s["extent"])
print(s["centroid"])
print(s["cell_type_counts"])          # {"tetra": 1234, ...}
print(s["total_area"], s["signed_volume"], s["unsigned_volume"])
print(s["num_inverted"])              # cells with negative volume

compute_stats returns a dict with keys num_points, num_cells, bbox_min/bbox_max/extent/centroid (3-tuples), cell_type_counts ({type: count}), total_area, signed_volume, unsigned_volume, and num_inverted.

What is measured

• Bounding box / extent / centroid — over all points (a parallel reduction).
• Per-cell-type counts — the number of cells of each type.
• Total area — the area of 2D cells (triangles, quads) plus the boundary surface area of 3D cells (computed via surface extraction).
• Signed & unsigned volume — of 3D cells (tetra / hex / wedge / pyramid), via the divergence theorem over each cell's outward-wound faces. signed sums the volumes; unsigned sums their magnitudes.
• Inverted cells — 3D cells with negative signed volume (wrong orientation).

CLI

meshioplusplus stats mesh.vtu             # human-readable table
meshioplusplus stats mesh.vtu --json      # machine-readable JSON

See the CLI reference.

Other languages

• C API — mio_stats(mesh, &report) fills a mio_stats_report struct (the scalar measures; per-cell-type counts are not carried across the flat ABI). See the C API reference.
• Fortran — report = mesh%stats() returns a mio_stats_report derived type. See the Fortran reference.
• WebAssembly / JavaScript — stats(mesh) returns an object with numPoints, bboxMin/…, cellTypeCounts, totalArea, signedVolume, etc. See the WebAssembly reference.