Hello, I am using mixed finite element method to solve a PDE. Theoretically，
the following two equations are obtained, in which there are two unknown quantities u, p,
A * u+B * p=b1;
C * p+D * u=b2;
How does freefem++ implement this solution process?
Is p=C^ (-1) * (b2-D*u) brought into the first equation group? If matrix C is irreversible?