def check_convergence(lines):
    """
    Returns all the geometry convergence results
    :param lines: lines of a bagel output file
    """
    convergence_start = '  *** Geometry optimization started ***'
    for i, line in enumerate(lines):
        if line[:39] == convergence_start:
            convergence_list = []
            for line in lines[i+3:]:
                line = line.strip()
                if not line:
                    break
                convergence_list.append(line)
            return convergence_list
    raise Exception('Could not find convergence')