Ускорение аппроксиматоров нелинейных функций

Задача
Построить нелинейный аппроксиматор специальной функции lgamma(x) аргумента двойной точности. Аппроксиматор должен обеспечивать высокую точность (в битах мантиссы) в интервале (-10, +10) и отрабатывать за минимальное время. Работа индивидуальная.

Кроме lgamma(x) можно взять tgamma(x), erf(x), j0(x), atan(x). Эти функции имеют разную сложность к аппроксимации; можно выбрать функцию попроще или посложнее, и обеспечить различную точность, насколько усидчивости/интереса хватит.

Замер производительности студент будет делать сам. Замер точности делается с помощью машинерии, доступной инженерам Huawei (код аппроксиматора нужно прислать по почте, обратно отсылается картинка).

Кураторы
Павел Буровский
Ведущий инженер ключевых проектов,
Huawei Technologies Co.,Ltd.
Что вы узнаете?
Аппроксиматоры нелинейных функций
Вы получите опыт построения аппроксиматоров нелинейных функций с
заданными требованиями точности и максимальной производительности.
Алгоритмы
Вы изучите некоторые факты из функционального анализа, необходимые для проекта.
Система Sollya
Вы получите первый опыт работы со специальной
системой компьютерной алгебры Sollya и внимательнее посмотрите на бинарное представление чисел с плавающей точкой стандарта IEEE-754.
Требуемые навыки
Компьютерная алгебра, язык С, знакомство с бинарным представлением
вещественных чисел IEEE-754, работа в Linux. Идеально (но не обязательно) понимание, как на принципиальном уровне работает процессорный конвейер.
Дополнительные материалы
Более подробно ознакомиться с проектом можно в презентации по ссылке