2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




Начать новую тему Ответить на тему
 
 [C] Гамма - функция
Сообщение10.01.2016, 18:56 
Заслуженный участник


27/07/12
1405
САФУ Архангельск
потребовалось, к сожалению, перевести свои вычисления с Wolfram M. на более быстрый С

Однако, столкнулся с проблемой.

Использую ненормализированную неполную гамма функцию, то что в математике выглядит как

Код:
Gamma [a,z]


в библиотеке gsl, которую решил использовать, вызываю аналогичную (по описанию в хелпе):

Код:
double gsl_sf_gamma_inc (double a, double x)


при попытке посчитать, ответы у них разные, причем, серьезно.

возможно, кто-то уже сталкивался с данной проблемой? в чем разница между этими реализациями?

 Профиль  
                  
 
 Re: гамма - функция
Сообщение11.01.2016, 00:17 
Заслуженный участник


25/02/11
1786
C си реализациями не сталкивался. А компилирование в самой математике не устраивает? Это порядочно ускоряет процесс. В частности, гамма-функция там компилируется. Можно подключать для этого сторонний компилятор СИ.

 Профиль  
                  
 
 Re: гамма - функция
Сообщение11.01.2016, 07:22 
Заслуженный участник


27/07/12
1405
САФУ Архангельск
по мимо скорости, там к сожалению проблема в том, что одно из разросшихся финальных выражений - подсчет интеграла численно - а в него входит эта гамма функция - она не хочет считать вообще никак(

 Профиль  
                  
 
 Re: гамма - функция
Сообщение11.01.2016, 11:07 
Заслуженный участник


25/02/11
1786
Можно сравнить с Wolfram Alpha. Например, она выдает
Код:
Gamma[2,3]=4/E^3=0.1991482734714558

В математике то же самое. Это правильный ответ — интеграл по частям берется.

 Профиль  
                  
 
 Re: [C] Гамма - функция
Сообщение25.01.2016, 23:33 
Заслуженный участник
Аватара пользователя


11/03/08
9490
Москва
А там верхняя или нижняя неполная гамма?

 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ Сообщений: 5 ] 

Модераторы: Karan, Toucan, PAV, maxal, Супермодераторы



Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group