it should work now

Author: rflamary

ot/gromov/_gw.py

@@ -586,6 +586,7 @@ def fused_gromov_wasserstein2(M, C1, C2, p=None, q=None, loss_fun='square_loss',
     lin_term = nx.sum(T * M)
     log_fgw['quad_loss'] = (fgw_dist - (1 - alpha) * lin_term)
     log_fgw['lin_loss'] = lin_term * (1 - alpha)
+    gw_term = log_fgw['quad_loss'] / alpha
 
     if loss_fun == 'square_loss':
         gC1 = 2 * C1 * nx.outer(p, p) - 2 * nx.dot(T, nx.dot(C2, T.T))