DB
Ну в целом вопрос с подвохом) мне в целом конкретный тип не важен, если весь остальной код после умеет работать с этим типом и что самое главное без потери производительности
Size: a a a
2021 July 06
АО
Что например, если вы захотите использовать эту библиотеку для построения какой-нибудь экономической модели? А для решения этой модели понадобится вычислять производные. И как уже было замечено выше, явная аннотация результата сломает AD.
АО
И по простому сделать уже не получится.
DB
Идея в том, чтобы математическая функция писалась в максимально шаблонном виде и умела работать с самым широким спектром чисел(f32, f64, big, complex) а не просто выдавала ошибку, когда в результате получается что то отличное от big
A
function foo(a, b, c, d)
local r
if some_condition
r = a/b
else
r = c/d
end
return r
end
Функция выше нестабильна по типам, а функция ниже стабильна
function foo(a, b, c, d)
local r::BigFloat
if some_condition
r = a/b
else
r = c/d
end
return r
end
local r
if some_condition
r = a/b
else
r = c/d
end
return r
end
Функция выше нестабильна по типам, а функция ниже стабильна
function foo(a, b, c, d)
local r::BigFloat
if some_condition
r = a/b
else
r = c/d
end
return r
end
АО
Вообще, если захочется сделать такое обобщение, то лучше взять все типы, вычислить через promote общий тип и к нему всё приводить.
АО
Тогда и стабильность будет и ограничений меньше.
DB
Первая функция вполне может быть type stable, все зависит от входных данных
DB
А вот вторая функция не работает для всех входных данных
DB
Хотя безусловно является typestable конечно, тут смысла спорить нет
2021 July 07
DB
Как вариант да, promote хорошее решение
DB
Стандартный паттерн для таких случаев это
foo(a::T, b::T) where T = …
foo(a, b) = foo(promote(a,b)…)
foo(a::T, b::T) where T = …
foo(a, b) = foo(promote(a,b)…)
DB
Конечно, бывают случаи когда можно и нужно проаннотировать типы
DB
Но идея джулии, чтобы делать это как можно меньше без явной необходимости
АО
Собственно говоря весь multiple dispatch :-)
DB
Мы в нашей лабе написали достаточно большой пакет для байесовских методов для обработки сигналов, там очень много математики…очень, но ни в одной строчке кода мы нигде явно не указываем аннотации типов как раз для того чтобы поддерживать разные типы и автодифференцирования «из коробки». Производительность при этом не страдает
DB
(Пользуясь случаем, кстати, приглашаю посмотреть на наш пакет на JuliaCon2021:)
АО
А этот пакет есть в открытом доступе?
АО