Fixing #38: replaced interpolate( 0 ) with setToZero() again in inv diag
Fixing #38 (closed): replaced interpolate( 0 ) with setToZero() again in inv diag kernel type. Also extended test such that the diagonal is assemble more than once (without this fix that failed).