Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Metric intersection #382

Open
xgarnaud opened this issue Dec 10, 2021 · 0 comments
Open

Metric intersection #382

xgarnaud opened this issue Dec 10, 2021 · 0 comments

Comments

@xgarnaud
Copy link

Dear all,

Thank you very much for this library.

I was trying to understand how metric are computed and handled, and I think there might be an issue with metric intersection.
If I understand correctly when I intersect two metrics M_0$ and M_1 to get M, for any vector v I should have
v^T M_0 v < v^T M v
and
v^T M_1 v < v^T M v
or at least numerically something very close to this.

I tried to check this numerically. Usually this test passes, but when I use e.g.

M0 << 0.686582, 0.310726, -0.440491, 0.310726, 0.765968, 0.284695, -0.440491,
      0.284695, 0.657285;
M1 << 0.914328, -0.689839, -0.270308, -0.689839, 1.29997, 0.048584, -0.270308,
      0.048584, 1.15847;

I find that v^T M_1 v > v^T M v for the last two eigenvectors of M1:
v^T M_1 v = 1.14424 and v^T M v = 1.06664
and
v^T M_1 v = 1.88177 and v^T M v = 1.58967

The code I used is here:
https://gist.github.com/xgarnaud/0b6c5e54ad3398319292fb95b7c7eaf3

Have you ever seen such a behavior (it seems when it happens M_0 has a small eigenvalue, but it can work fine even if it has one)? Do you have an idea where it could comme from?
Best regards

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant