1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
|
"""Tests for module gaussian"""
# Author: Theo Gnassounou <theo.gnassounou@inria.fr>
# Remi Flamary <remi.flamary@polytehnique.edu>
#
# License: MIT License
import numpy as np
import pytest
import ot
from ot.datasets import make_data_classif
def test_bures_wasserstein_mapping(nx):
ns = 50
nt = 50
Xs, ys = make_data_classif('3gauss', ns)
Xt, yt = make_data_classif('3gauss2', nt)
ms = np.mean(Xs, axis=0)[None, :]
mt = np.mean(Xt, axis=0)[None, :]
Cs = np.cov(Xs.T)
Ct = np.cov(Xt.T)
Xsb, msb, mtb, Csb, Ctb = nx.from_numpy(Xs, ms, mt, Cs, Ct)
A_log, b_log, log = ot.gaussian.bures_wasserstein_mapping(msb, mtb, Csb, Ctb, log=True)
A, b = ot.gaussian.bures_wasserstein_mapping(msb, mtb, Csb, Ctb, log=False)
Xst = nx.to_numpy(nx.dot(Xsb, A) + b)
Xst_log = nx.to_numpy(nx.dot(Xsb, A_log) + b_log)
Cst = np.cov(Xst.T)
Cst_log = np.cov(Xst_log.T)
np.testing.assert_allclose(Cst_log, Cst, rtol=1e-2, atol=1e-2)
np.testing.assert_allclose(Ct, Cst, rtol=1e-2, atol=1e-2)
@pytest.mark.parametrize("bias", [True, False])
def test_empirical_bures_wasserstein_mapping(nx, bias):
ns = 50
nt = 50
Xs, ys = make_data_classif('3gauss', ns)
Xt, yt = make_data_classif('3gauss2', nt)
if not bias:
ms = np.mean(Xs, axis=0)[None, :]
mt = np.mean(Xt, axis=0)[None, :]
Xs = Xs - ms
Xt = Xt - mt
Xsb, Xtb = nx.from_numpy(Xs, Xt)
A, b, log = ot.gaussian.empirical_bures_wasserstein_mapping(Xsb, Xtb, log=True, bias=bias)
A_log, b_log = ot.gaussian.empirical_bures_wasserstein_mapping(Xsb, Xtb, log=False, bias=bias)
Xst = nx.to_numpy(nx.dot(Xsb, A) + b)
Xst_log = nx.to_numpy(nx.dot(Xsb, A_log) + b_log)
Ct = np.cov(Xt.T)
Cst = np.cov(Xst.T)
Cst_log = np.cov(Xst_log.T)
np.testing.assert_allclose(Cst_log, Cst, rtol=1e-2, atol=1e-2)
np.testing.assert_allclose(Ct, Cst, rtol=1e-2, atol=1e-2)
def test_bures_wasserstein_distance(nx):
ms, mt = np.array([0]), np.array([10])
Cs, Ct = np.array([[1]]).astype(np.float32), np.array([[1]]).astype(np.float32)
msb, mtb, Csb, Ctb = nx.from_numpy(ms, mt, Cs, Ct)
Wb_log, log = ot.gaussian.bures_wasserstein_distance(msb, mtb, Csb, Ctb, log=True)
Wb = ot.gaussian.bures_wasserstein_distance(msb, mtb, Csb, Ctb, log=False)
np.testing.assert_allclose(nx.to_numpy(Wb_log), nx.to_numpy(Wb), rtol=1e-2, atol=1e-2)
np.testing.assert_allclose(10, nx.to_numpy(Wb), rtol=1e-2, atol=1e-2)
@pytest.mark.parametrize("bias", [True, False])
def test_empirical_bures_wasserstein_distance(nx, bias):
ns = 400
nt = 400
rng = np.random.RandomState(10)
Xs = rng.normal(0, 1, ns)[:, np.newaxis]
Xt = rng.normal(10 * bias, 1, nt)[:, np.newaxis]
Xsb, Xtb = nx.from_numpy(Xs, Xt)
Wb_log, log = ot.gaussian.empirical_bures_wasserstein_distance(Xsb, Xtb, log=True, bias=bias)
Wb = ot.gaussian.empirical_bures_wasserstein_distance(Xsb, Xtb, log=False, bias=bias)
np.testing.assert_allclose(nx.to_numpy(Wb_log), nx.to_numpy(Wb), rtol=1e-2, atol=1e-2)
np.testing.assert_allclose(10 * bias, nx.to_numpy(Wb), rtol=1e-2, atol=1e-2)
|