First time getting started with Fenics (and using a finite element analysis package in general), could I have some advice or help on a problem I'm working on?

[dokken]

@fea

The vertex_to_dof_map_vectorized and perturb_mesh certainly look more advanced than what I’m capable of understanding at the moment, but thankfully it looks like the only thing I need to mess around with are the user-defined expressions to see if I can’t implement a way to randomly generate y-deltas like you probably saw in my code.

I’ll

Another alternative is to generate the grid by hand, as sketched out in:

L-shape mesh generation with right element

This is way easier in DOLFINx, as you can work directly with numpy arrays from mpi4py import MPI import dolfinx import basix.ufl import ufl import numpy geometrical_dim = 2 c_el = basix.ufl.element("Lagrange", "triangle", 1, shape=(geometrical_dim,)) ufl_domain = ufl.Mesh(c_el) if MPI.COMM_WORLD.rank == 0: nodes = numpy.array( [ [0.0, 0.0], [1.0, 0.0], [0.0, 1.0], [1.0, 1.0], [0.0, 2.0], [1.0, 2.0], …