Sorry for the delay. I followed the directions from this post: ode-for-cell-annotation/3762
What I did was run Resolvi in unsupervised mode and perform leiden clustering on the X_resolVI embeddings to get clusters.
Then I used the X_resolVI embeddings along with the cluster labels as input to scVIVA. From there I used the differential_expression module scVIVA uses.