[WIP] Fix quantized (F)GW solvers

RELEASES.md

@@ -22,6 +22,7 @@
 - Fix mean centering in `ot.dr.fda` and `ot.dr.wda`: `np.mean(X)` returned a scalar instead of the per-feature mean, so `proj` did not center the data as documented. In `ot.dr.fda` the same pattern in the class means made the between-class scatter matrix independent of which features separate the classes, and FDA returned a non-discriminant direction (PR #840)
 - `ot.dr.fda` and `ot.dr.wda` no longer modify the input array `X` in place (PR #840)
 - Fix `UnbalancedSinkhornTransport` `transform` failing with `AttributeError: 'NoneType' object has no attribute 'array_equal'` when `fit` was called with missing parameters (PR #837, Issue #650)
+- Fix quantized (F)GW solvers that ordered OT based on clusters and not initial node ordering (PR #857, Issue #786)
 
 ## 0.9.7.post1
 

examples/gromov/plot_quantized_gromov_wasserstein.py

@@ -259,6 +259,8 @@ def draw_graph(
     list_R2,
     list_h1,
     list_h2,
+    part1,
+    part2,
     MR=None,
     alpha=1.0,
     build_OT=True,