Incorrect PETSc.Mat block size when using enriched_element (MINI element)
youre not doing anything wrong here. the block size comes from the dofmap (check V.dofmap.bs), which dolfinx derives from the basix element. a plain vector lagrange element is built as a blocked element, so bs = gdim and you get a baij matrix. enriched_element builds a generic combined element that doesnt advertise an…
Problem outcome: unknown