“…We do not explicitly build in the coupling condition by constructing a basis of the space V h , but incorporate it implicitly, by replacing V h by Vh and adding (18) as additional constraint. Since Vh , Q h and Λ h are finite dimensional, the bilinear forms a, b, c and d can be represented as matrices A h , B h , C h and D h acting on vectors of real numbers x h , u h and λ h representing the elements in Vh , Q h and Λ h , respectively, with respect to the chosen basis.…”