MixedElement vs MixedFunctionSpace

variational formulation
11

Hi @Thomas2,

I don’t know how to thank you, I was FINALLY able to run the simulation thanks to your suggestions! For clarity, I needed to change your pseudo code a bit, this is the final version:

F = # my weak form
a, L = lhs(F), rhs(F)
A = assemble(a, keep_diagonal=True)
b = assemble(L)
for bc in bcs: bc.apply(A, b)
A.ident_zeros()
sol = PETScLUSolver()
x = Function(V)
sol.solve(A, x.vector(), b)

This solution does not need the use of MixedFunctionSpace but just of MixedElement, I also run the simulation for my system of PDEs and everything run smoothly :smiley:!