Multisource Geodesic Distance not Working

[paul0noah]

Hi,

I wanted to use the multisource functions for geodesic distance computation (basically i wanted to compute full geodesic distance matrix), i.e. i invoked with:

solver = pp3d.MeshHeatMethodDistanceSolver(V,F)
dist = solver.compute_distance_multisource(np.arange(V.shape[0]).tolist())  

From that, i would expect output dist beeing |V|x|V| dimensional but it just returns a |V|x1 array.

My fix for now is, but i guess this is not the intended use for multisource:

solver = pp3d.MeshHeatMethodDistanceSolver(V,F)
dist = np.zeros((V.shape[0], V.shape[0]))
for i in range(V.shape[0]):
        dist[i, :] = solver.compute_distance(i)

Note, I also tried the one-off version which yields same issue:

dist = pp3d.compute_distance_multisource(V, F, np.arange(V.shape[0]).tolist())  

Thanks for looking into it and thanks for the great work! 🙏