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

What does "tol" in "igapfill" means? #184

Closed
furiousbean opened this issue Sep 25, 2014 · 8 comments
Closed

What does "tol" in "igapfill" means? #184

furiousbean opened this issue Sep 25, 2014 · 8 comments

Comments

@furiousbean
Copy link
Contributor

Today I'd tested a new "igapfill" method with default parameters on finite rank series with a gap (but without any noise), and found out that precision of gap filling is about 10^-3, i.e. it reconstructs only 3 digits after zero at gap measurements. Thus, there are few questions:

  1. Is it correct that "tol" parameter means squared distance between reconstructions at last two iterations of algorithm? Then maybe "tol" must mean just a distance? Or it should be clearly mentioned at "igapfill" help page?
  2. Is it enough to have just three significant digits by default?
@asl
Copy link
Owner

asl commented Sep 26, 2014

As everywhere - tol is squared distance used for converged by itself. And surely the tolerance for prediction by itself corresponds quite poorly with the true limit.

@neg99
Copy link

neg99 commented Sep 26, 2014

Поддерживаю вопросы Никиты.
Я не нашла в хелпе информации, что это расстояние в квадрате. Поэтому вопрос Никиты по-существу. (И это не стандартно, по-моему, что не в масштабе рядов.)
В статье по iossa, eps это корень из средней суммы квадратов.
Саша, у вас в iossa tol=1.e-5 (не 6 и не 3) - что это значит?.

@asl
Copy link
Owner

asl commented Sep 26, 2014

У нас с этим разброд и шатание, да. В iossa - это норма. В cadzow - квадрат нормы. igapfill я делал в соответствии с cadzow.

@neg99
Copy link

neg99 commented Sep 26, 2014

Думаю, нужно везде сделать норму.
А вот какие значения tol по умолчанию, надо выбрать. 10^-6 или ^-5 везде - чем может быть плохо? Так как есть еще maxiter. это можно ограничить.
Сейчас еще в iossa maxiter=100 (логично, так как может быть очень плохая сходимость), в cadzow и igapfill отключено по умолчанию (т.е. нули).

@eodus
Copy link
Contributor

eodus commented Sep 26, 2014

От себя добавлю, что в cadzow/igapfill норма это максимальный квадрат (max(v^2)), а в ossa --- средний (sqrt(mean(v^2))).

@asl, а какой резон в был использовании именно максимума?... Если он есть, надо везде сделать его. Или везде среднее.

@asl
Copy link
Owner

asl commented Sep 28, 2014

Я думаю, что с точки зрения сходимости в себе максимум выглядит более адекватной метрикой.

@neg99
Copy link

neg99 commented Oct 4, 2014

В iossa есть параметр norm = function(x) sqrt(mean(x^2)).
Соответственно, нужно также сделать и в остальных итеративных методах.
По умолчанию можно задать и max, не так важно.

@asl
Copy link
Owner

asl commented Oct 5, 2014

Fixed in c28d1dc

@asl asl closed this as completed Oct 5, 2014
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

4 participants