BEGIN_PROVIDER [integer, nnk]
implicit none
BEGIN_DOC
! provides nnk = total number of non zero elements in H
END_DOC
nnk=0
END_PROVIDER