Multi-RHS in FEM - Basic Implementation

Simple modifications - Complete implementation would require to modify
also the call of boundary, to make sure BOUNDARY_DATA contains the info
from all load cases.