As continuation of PR #36504, I am working on some functionalities related to counting (graph) homomorphisms. Currently, Sage already has has_homomorphism_to function, but it only returns one homomorphism even if multiple exist.

It seems that this function is implemented with some LP solver. Would it be possible to extend it to support returning all homomorphisms?

Edit: This one will be addressed in a new branch.