Changed routine to construct stiffness matrix for a system to representing the diffusion tensor at the nodes instead of at the elements