We have investigated and compared the ability of numerical and Gaussian-type basis sets combined with density functional theory (DFT) to accurately describe the geometries, binding energies, and electronic properties of aluminum clusters, Al12XHn (X = Al, Si; n = 0, 1, 2). DFT results are compared against high-level benchmark calculations and experimental data where available. Properties compared include geometries, binding energies, ionization potentials, electron affinities, and HOMO-LUMO gaps. Generally, the PBE functional with the double numerical basis set with polarization (DNP) performs very well against experiment and the analytical basis sets for considerably less computational expense.