FEniCS/Dolfinメッシュを確認.#

PyVistaは meshio を活用して, FEniCS/Dolfin XMLフォーマットを含む,VTKがネイティブにサポートしていない多くのメッシュフォーマットを読み込みます.

import pyvista as pv
from pyvista import examples

サンプルデータリポジトリからFEniCS/Dolfinメッシュのサンプルをダウンロードしましょう.これにより,XML Dolfinメッシュがダウンロードされ,PyVistaのデータディレクトリに保存されます.

saved_file = examples.download_file("dolfin_fine.xml")
print(saved_file)
/home/runner/.cache/pyvista_3/dolfin_fine.xml

これで,XML Dolfinメッシュをローカルに保存できました.このファイル名は,PyVistaの pyvista.read() メソッドに直接渡して,PyVistaメッシュに読み込むことができます.

dolfin = pv.read(saved_file)
dolfin
UnstructuredGridInformation
N Cells5400
N Points2868
X Bounds0.000e+00, 1.000e+00
Y Bounds0.000e+00, 1.000e+00
Z Bounds0.000e+00, 0.000e+00
N Arrays0


Dolfinメッシュを加工してプロットすることができます.

qual = dolfin.compute_cell_quality()
qual.plot(show_edges=True, cpos="xy")
read dolfin

Total running time of the script: (0 minutes 0.731 seconds)

Sphinx-Galleryによるギャラリー