Having the residual stiffness 'k' in two variables k and k_local
was a constant source of confusion.  This patch eliminates
the k_local variable, and makes sure the correct values
are used.