I am a little confused with the meaning of the block_size property of a mixed element in FEniCSx. Consider the following FEniCS code: from dolfin import * mesh = UnitSquareMesh(1, 1) FE = FiniteElement("Lagrange", mesh.ufl_cell(), 1) ME = MixedElement([FE, FE]) V_ME = FunctionSpace(mesh, ME) V …

Please note that you are comparing two different block sizes. Note that DOLFINx also has a DofMap block size, which tells you if dofs are blocked in the cell_dofs, i.e. print(V_ME.dofmap.bs) 1 print(V_ME.dofmap.cell_dofs(0)) array([0, 1, 2, 3, 4, 5], dtype=int32) Note that for V, you get: pr…

Block_size of a mixed element in DOLFINx

dolfinx

violetus December 29, 2023, 9:32am 2

Consider Getting cell dofs for vector valued function space in dolfinx.

Topic		Replies	Views
Fenicsx method, which is equaivalent to dofmap().dofs() in fenics	6	902	August 6, 2023
Dofmap of fem.VectorFunctionSpace vs. fem.FunctionSpace(mesh, MixedElement) dolfinx	4	700	December 20, 2023
Inconsistent results with sub MixedElement functions	1	371	July 7, 2021
Vector quadrature_element in FEniCSx 0.7.0 General	2	272	October 12, 2023
Initial Conditions and Function update on time in mixed Element dolfinx	2	287	December 5, 2023

Block_size of a mixed element in DOLFINx

Related Topics