V🇺
какого порядка число параметров и степень производной у этих более сложных функций?
Size: a a a
2019 November 04
FO
максимум сколько видел - третья производная у функции. три параметра, по которым берется производная
V🇺
Ну тогда forward AD - твой выбор
V🇺
например https://autodiff.github.io/, если плюсы не страшно
FO
а я смогу потом с этого какой-то код производных получить?
V🇺
смотря что ты имеешь в виду под этим
FO
мне бы каждый раз в коде не вызывать диффиренциирование, так как число точек порядка нескольких миллионов для простейшего расчета
V🇺
мне бы каждый раз в коде не вызывать диффиренциирование, так как число точек порядка нескольких миллионов для простейшего расчета
оно в compile time дифференциирует, условно говоря
FO
Вот в
R я могу сделать так...> f <- expression(a*x^(4/3) * (1-exp(x^(-8/3))))
> D(f,"x");
a * (x^((4/3) - 1) * (4/3)) * (1 - exp(x^(-8/3))) - a * x^(4/3) *
(exp(x^(-8/3)) * (x^((-8/3) - 1) * (-8/3)))
FO
Пару взмахов sed и код для Fortran готов
V🇺
Вот в
R я могу сделать так...> f <- expression(a*x^(4/3) * (1-exp(x^(-8/3))))
> D(f,"x");
a * (x^((4/3) - 1) * (4/3)) * (1 - exp(x^(-8/3))) - a * x^(4/3) *
(exp(x^(-8/3)) * (x^((-8/3) - 1) * (-8/3)))
это символьная производная, и считать ее влоб дороже чем AD
FO
это символьная производная, и считать ее влоб дороже чем AD
а потом я могу её записать кодом
V🇺
а потом я могу её записать кодом
да, и это будет дороже - см. количество дублированных подвыражений
FO
да, и это будет дороже - см. количество дублированных подвыражений
а компилятор меня тут не спасет?
V🇺
не скажу за фортран, но clang с gcc не спасут
https://godbolt.org/z/He4lv2
https://godbolt.org/z/He4lv2
V🇺
не скажу за фортран, но clang с gcc не спасут
https://godbolt.org/z/He4lv2
https://godbolt.org/z/He4lv2
можно там же проверить и фортран, но я в него не умею
FO
не скажу за фортран, но clang с gcc не спасут
https://godbolt.org/z/He4lv2
https://godbolt.org/z/He4lv2
увы, но без разницы.. а icc смог
V🇺
icc молодец, но дорогой зараза