Constructing stiffness matrix for system of elliptic PDE works now. Added a unit test for FEM module.