Fixed a bug in the routine to fill the stiffness matrix for elliptic systems