Python functions to wrap Cubit mesh metrics

[hovey]

Currently, we have the minimum scaled Jacobian (MSJ) from Cubit accessible through Python

• see cubit_inp_to_minsj_csv.py

We need the following (two, and possibly four) scripts to be created:

• cubit_inp_to_skew_csv.py page 97 of 107 of [2],
• cubit_inp_to_aspect_ratio_csv.py (aka Maximum edge ratio) page 87 of 107 of [2]
• cubit_inp_to_warp_csv.py # This seems like a quadrilateral metric, not a hexahedral metric.
• cubit_inp_to_min_max_angle_csv.py # This seems like a quadrilateral metric, not a hexahedral metric.

Reference:

• [1] Li 2022 Table 2 (page 14) uses these same metrics, and we may wish to benchmark our results with their results.
• [2] Stimpson C, Ernst C, Knupp P, Pébay P, Thompson D. The verdict library reference manual. Sandia National Laboratories Technical Report. 2007 May 8;9(6), pdf
  • Cubit makes use of the Verdict library for mesh quality metrics.
  • There is also a Kitware version.
• [3] The quad_quality.pdf document may also be useful.

[calehma]

Skew script completed. Will update repository with aspect ratio next.

[calehma]

Completed the second of the hexahedral quality metric scripts. Do you want to go ahead with the quadrahedral metric scripts as well?

[hovey]

I have crossed-out the two quad-related metrics. We currently have no plan to use quads in our meshes. The Li 2022 meshes used both hex elements and quad elements, which explains why warp and min/max angle were also included in their work.

So, please cancel the request for the quad-related items. I will mark this issue as completed.