Error in Barycentre computation
Dear all I followed the 1D Wasserstein Barycentre demo on Github https://github.com/rflamary/POT/blob/master/notebooks/plot_barycenter_1D.ipy... and extended it to my use case where I have 100 pdfs/distributions with 5 bins each. I assign equal weights to each distribution and regularisation parameter value is same as the demo. So the dimensions of dataframes in my code are: A (5,100) M (5,5) weights (100) But on executing my code, I get the following error: C:\Users\Shikha\AppData\Roaming\Python\Python36\site-packages\ot\bregman.py:1049: RuntimeWarning: overflow encountered in true_divide u = (u.T * geometricBar(weights, UKv)).T / UKv Seeking help in resolving this error. Best, Shikha
Hello Shikha, This warning occurs for a too small regularization term in the sinkhorn iterations. In order to avoid it you can make the regularization term larger (at the cost of doing an approximation of transport). Or you can use the stabilized version of the algorithm with option method='sinkhorn_stabilized'. Note that this option appeared in the last 0.6 release so you might have to update POT. Best, Rémi Le 22/09/2019 à 17:59, fpm16shikhav@iima.ac.in a écrit :
-- Rémi Flamary Web: https://remi.flamary.com Tel: +33 (0)4 92 07 63 80 Laboratoire Lagrange, UMR CNRS 7293 Observatoire de la Côte d'Azur Université Côte d'Azur
Hello Shikha, This warning occurs for a too small regularization term in the sinkhorn iterations. In order to avoid it you can make the regularization term larger (at the cost of doing an approximation of transport). Or you can use the stabilized version of the algorithm with option method='sinkhorn_stabilized'. Note that this option appeared in the last 0.6 release so you might have to update POT. Best, Rémi Le 22/09/2019 à 17:59, fpm16shikhav@iima.ac.in a écrit :
-- Rémi Flamary Web: https://remi.flamary.com Tel: +33 (0)4 92 07 63 80 Laboratoire Lagrange, UMR CNRS 7293 Observatoire de la Côte d'Azur Université Côte d'Azur
participants (2)
-
fpm16shikhav＠iima.ac.in -
Rémi Flamary