This possible problem was pointed out by @jaimeaxel. https://github.com/dftlibs/numgrid/blob/master/numgrid/AtomGrid.cpp#L94-L95