gene sets score and gsva

Hi,
In single cell data analysis, we usually score a set of genes. The score is the average expression of a set of genes subtracted with the average expression of a reference set of genes: https://scanpy.readthedocs.io/en/1.9.x/generated/scanpy.tl.score_genes.html

I feel it's function is a bit similar with gsva. Which one is better to be a representation of cell cluster gene signature ?