dolfinx-tutorial
dolfinx-tutorial copied to clipboard
jorgensd
A reimplementation of the Springer book: https://github.com/hplgit/fenics-tutorial/, covering new topics as well as transitioning from dolfin to dolfinx
```python # Computing a onesided integral over an interior facet with consistent orientation # Enabled by: https://github.com/FEniCS/dolfinx/pull/2269 # Copyright 2023 Jørgen S. Dokken # SPDX: MIT import dolfinx import numpy...
Minimal example ```python # Copyright (C) 2023 Jørgen S. Dokken # # Solve a singular Poisson problem with a null-space # # SPDX-License-Identifier: MIT from dolfinx import mesh, fem, io...
Gather all data on one process, as shown in: https://fenicsproject.discourse.group/t/how-to-collect-the-global-mesh-without-writing-a-file/12445/4?u=dokken ```python # Copyright (C) 2023 Jørgen S. Dokken # # SPDX-License-Identifier: MIT # # Code to plot mesh with pyvista...
Ref: https://fenicsproject.discourse.group/t/mixed-formulation-of-poisson-problem-in-dolfinx/10149/7?u=dokken and https://fenicsproject.org/olddocs/dolfin/1.5.0/python/demo/documented/mixed-poisson/python/documentation.html
Hi @jorgensd Are there plans to add a demo for solving the poisson equation using the DG method? Thanks! Rem
Can use `ufl.reconstruct` to increase the element degree: ```python def error_L2(uh, u_ex, degree_raise=3): # Create higher order function space ufl_el = uh.function_space.ufl_element() new_el = ufl_el.reconstruct(degree=degree_raise+ufl_el.degree()) mesh = uh.function_space.mesh W =...
My operating system is Ubuntu 22.04 in WSL for Windows 11. When I run the code, there are two points that need attention. First, `pyvista.start_xvfb()` should be deleted. I tried...
Hyperelasticity demo uses F as deformation gradient and then again as the left hand side of the nonlinear variational problem. While the code works fine; not sure how fem.petsc.NonlinearProblem(F, u,...
