def test_orbital_radial(): r = np.linspace(0, 10, 1000) f = np.exp(- r) o = sisl.SphericalOrbital(2, (r, f)) sisl.plot(o) fig = plt.figure() sisl.plot(o, axes=fig.gca()) plt.close('all')
def test_orbital_harmonics(): r = np.linspace(0, 10, 1000) f = np.exp(-r) o = sisl.SphericalOrbital(2, (r, f)) sisl.plot(o, harmonics=True) plt.close('all')