I am aware that there exists a flowchart about how to determine the point group of a molecule with a given geometry, such as given in the diagram available in ChemLibre Texts. If we just look at the geometry visually or having it pictured in mind, it wouldn't be difficult to follow this flowchart. But I would like to automate this procedure inside my electronic-structure code.
Is there a well-established computational method to determine, e.g. how many rotational axes of symmetry a molecule has? Is it possible to perform this identification generally, that is, in a way in which I don't need to write a function or routine dedicated for each symmetry type? References are much appreciated.