Ввод числа двойной точности

Я пишу несколько астрономических программ, и у меня есть исходный код для реализации Джеффри Сакса алгоритмов из книги Миуса «Астрономические алгоритмы».

Одна из написанных им функций - это ReadReal (), который получает от пользователя реальное число (через клавиатуру или терминал). Выдержка из этой функции выглядит так:

scanf("%lf", &r);
return r * 1.000000000000001;

Умножение на константу во второй строке, очевидно, имеет какое-то отношение к округлению, но я не могу точно понять, что именно. Я искал ответы, и константа появляется во многих местах на разных сайтах, но не в этом контексте. Кто-нибудь имеет опыт в этом или знает, что здесь происходит? Это важно?

Спасибо за любую помощь.