def generate_ammonia_sites():
    'Generate the locations for the atoms in the ammonia molecule'
    (x, y) = (np.array([1.0, 0.0, 0.0]), np.array([0.0, 1.0, 0.0]))
    a = 1.4
    n = np.array([0.0, 0.0, 0.0])
    h1 = (n + (a * y))
    h2 = ((n - ((a * y) / 2.0)) + ((a * x) * (np.sqrt(3) / 2)))
    h3 = (h2 - ((a * x) * np.sqrt(3)))
    return (n, h1, h2, h3)
Generate the locations for the atoms in the ammonia molecule
def generate_ammonia_sites():
    
    (x, y) = (np.array([1.0, 0.0, 0.0]), np.array([0.0, 1.0, 0.0]))
    a = 1.4
    n = np.array([0.0, 0.0, 0.0])
    h1 = (n + (a * y))
    h2 = ((n - ((a * y) / 2.0)) + ((a * x) * (np.sqrt(3) / 2)))
    h3 = (h2 - ((a * x) * np.sqrt(3)))
    return (n, h1, h2, h3)