Beamlet formation in negative ion sources for fusion applications: A computational grid system comparison