LHS and RHS of global solve implemented, ready for eigen solver