import firedrake
/home/firedrake/firedrake/lib/python3.12/site-packages/pytools/persistent_dict.py:52: RecommendedHashNotFoundWarning: Unable to import recommended hash 'siphash24.siphash13', falling back to 'hashlib.sha256'. Run 'python3 -m pip install siphash24' to install the recommended hash.
warn("Unable to import recommended hash 'siphash24.siphash13', "
import viskex
import common_firedrake as common # isort: skip
Generate meshes of the unit cube by dividing each edge of the cube in 6 segments, using either a tetrahedral or hexahedral mesh.
cube_tetra = firedrake.UnitCubeMesh(6, 6, 6, hexahedral=False)
cube_hexa = firedrake.UnitCubeMesh(6, 6, 6, hexahedral=True)
Interpolate the scalar field $x^3 + y^2 + z^4$, and store it either in a firedrake function or a UFL expression.
cube_tetra_scalar_field, cube_tetra_scalar_field_ufl = common.prepare_scalar_field_cases(
cube_tetra, lambda x: x[0]**3 + x[1]**2 + x[2]**4)
cube_hexa_scalar_field, cube_hexa_scalar_field_ufl = common.prepare_scalar_field_cases(
cube_hexa, lambda x: x[0]**3 + x[1]**2 + x[2]**4)
/home/firedrake/firedrake/src/loopy/loopy/schedule/__init__.py:2255: LoopyWarning: in kernel expression_kernel: Falling back to a slow scheduler implementation due to: cannot schedule kernels with priority dependencies between sibling loop nests (add 'v1_scheduler_fallback' to silenced_warnings kernel attribute to disable) warn_with_kernel(
/home/firedrake/firedrake/src/loopy/loopy/schedule/__init__.py:2255: LoopyWarning: in kernel expression_kernel: Falling back to a slow scheduler implementation due to: cannot schedule kernels with priority dependencies between sibling loop nests (add 'v1_scheduler_fallback' to silenced_warnings kernel attribute to disable) warn_with_kernel(
Plot the scalar field with a contour plot.
viskex.firedrake.plot_scalar_field(cube_tetra_scalar_field, "scalar")
error: XDG_RUNTIME_DIR is invalid or not set in the environment. MESA: error: ZINK: failed to choose pdev glx: failed to create drisw screen
viskex.firedrake.plot_scalar_field(cube_tetra_scalar_field_ufl, "scalar")
error: XDG_RUNTIME_DIR is invalid or not set in the environment. MESA: error: ZINK: failed to choose pdev glx: failed to create drisw screen
# PYTEST_XFAIL: Temporarily broken due to firedrake#3315
"""Expect this cell to fail.
viskex.firedrake.plot_scalar_field(cube_hexa_scalar_field, "scalar")
"""
'Expect this cell to fail.\n\nviskex.firedrake.plot_scalar_field(cube_hexa_scalar_field, "scalar")\n'
# PYTEST_XFAIL: Temporarily broken due to firedrake#3315
"""Expect this cell to fail.
viskex.firedrake.plot_scalar_field(cube_hexa_scalar_field_ufl, "scalar")
"""
'Expect this cell to fail.\n\nviskex.firedrake.plot_scalar_field(cube_hexa_scalar_field_ufl, "scalar")\n'
Warp mesh according to scalar field.
viskex.firedrake.plot_scalar_field(cube_tetra_scalar_field, "scalar", warp_factor=0.1)
error: XDG_RUNTIME_DIR is invalid or not set in the environment. MESA: error: ZINK: failed to choose pdev glx: failed to create drisw screen
# PYTEST_XFAIL: Temporarily broken due to firedrake#3315
"""Expect this cell to fail.
viskex.firedrake.plot_scalar_field(cube_hexa_scalar_field, "scalar", warp_factor=0.1)
"""
'Expect this cell to fail.\n\nviskex.firedrake.plot_scalar_field(cube_hexa_scalar_field, "scalar", warp_factor=0.1)\n'