KT
Size: a a a
2020 January 30
KT
Там и числа в формулу подставить можно, и продифференцировать, и проинтегрировать, и простые ODE решить символьно
АО
Вот ещё было какое-то старое обсуждение
https://discourse.julialang.org/t/sympy-jl-vs-symengine-jl-vs-reduce-jl-vs/10381
https://discourse.julialang.org/t/sympy-jl-vs-symengine-jl-vs-reduce-jl-vs/10381
KT
Symengine мне не особо понравился, там нет нифига почти
FO
Библиотека выглядит немного outdated
Дёшево и сердито - использовать интерполяцию
Но вообще, выглядит так, что лучше пользоваться библиотеками типа
Дёшево и сердито - использовать интерполяцию
K_sigma =:(energy/rho^(4/3))
a_sigma = :(mu * sqrt($K_sigma / pi) / (3*rho^(1/3)))
differentiate(a_sigma, :rho)
Но вообще, выглядит так, что лучше пользоваться библиотеками типа
ForwardDiffСпасибо большое! Сейчас поизучаю, точно ли он те формулы выдает при дифференциировании ф-ий многих переменных :)
Но вроде все ок, когда я отдельно дифференцирую и по энергии, и по плотности..
Но вроде все ок, когда я отдельно дифференцирую и по энергии, и по плотности..
FO
Я за sympy
по некоторым причинам мне не очень нравится SymPy
KT
Да мне тоже, он здоровый и не быстрый, к тому же странно топить за Julia и при этом пользоваться питоном, но на Julia имхо ничего близкого по функционалу не предвидится и не понятно, нужно ли вообще )
FO
Да мне тоже, он здоровый и не быстрый, к тому же странно топить за Julia и при этом пользоваться питоном, но на Julia имхо ничего близкого по функционалу не предвидится и не понятно, нужно ли вообще )
ну я код с него брал (после дифференциирования) и вставлял в фортран, чтобы вычислять выражения. он давал скорость в два раза меньшую, чем код из maxima
FO
У julia поплыла точность после симплификации :(
KT
Если конечная цель вычислять выражения, то можно смотреть автоматическое дифференцирование, а не символьное
KT
Оно по идее и с big float должно работать
KT
И с интервальной арифметикой (наверное)
FO
Если конечная цель вычислять выражения, то можно смотреть автоматическое дифференцирование, а не символьное
Надо формулу отдебажить, которая сгенерировалась, но где-то косяк есть в ней
АО
Мне кажется, мы столько обсуждаем, что уже можно было бы её руками посчитать ^_^
АО
Но вообще, если нужно формулы символьно дифференцировать только для того, чтобы их потом копипастить в фортран, то джулия это имхо оверкилл.
Проще какой-нибудь онлайн математикой воспользоваться.
Проще какой-нибудь онлайн математикой воспользоваться.
FO
Мне кажется, мы столько обсуждаем, что уже можно было бы её руками посчитать ^_^
Уже считал. Для дебага не помогло :(
FO
а как можно удалить переменную?