Thanks for the suggestion to use the solvers in scipy which might be an easier approach during the initial development phase on a workstation. Unfortunately I am getting the following error:
Traceback (most recent call last):
File "/home/andy/work/fenicsx/modanal/main_sp.py", line 143, in <module>
K = fem.assemble_matrix(kc, bcs=[bc])
File "/usr/lib/python3.10/functools.py", line 889, in wrapper
return dispatch(args[0].__class__)(*args, **kw)
File "/usr/lib/petscdir/petsc3.15/x86_64-linux-gnu-complex/lib/python3/dist-packages/dolfinx/fem/assemble.py", line 273, in _assemble_matrix_form
A: la.MatrixCSRMetaClass = create_matrix(a)
File "/usr/lib/petscdir/petsc3.15/x86_64-linux-gnu-complex/lib/python3/dist-packages/dolfinx/fem/assemble.py", line 97, in create_matrix
return la.matrix_csr(sp, dtype=a.dtype)
File "/usr/lib/petscdir/petsc3.15/x86_64-linux-gnu-complex/lib/python3/dist-packages/dolfinx/la.py", line 81, in matrix_csr
return matrixcls(sp)
File "/usr/lib/petscdir/petsc3.15/x86_64-linux-gnu-complex/lib/python3/dist-packages/dolfinx/la.py", line 54, in __init__
super().__init__(sp)
RuntimeError: Block size not yet supported
which has been seen before with the work around to use petsc. I am also initially looking at elasticity problems.