Кстати, недавно кто-то рекламировал пакет, что-то там kernel назывался.

Скорости не прибавит (скорее даже наоброт), но зато должен уменьшить длину кода, так как вроде бы там достаточно компактно можно было ядра прокатывать по матрицам.

Секунду, сейчас найду.

Хм нет, на самом деле это старый пакет, видимо недавно реанимировали.

https://github.com/stev47/StaticKernels.jl

Тем не менее, вдруг пригодится (сам не пользовался и гарантировать применимость к данной задаче не могу 😊 ).

Нестабильностей вроде не было. А что, лучше передавать callable struct?

Нет, не лучше.

Там просто есть особые случаи.

То есть можно написать что-то вроде

function adv_oper!(fout::stvec_t, fin::stvec_t, params::params_t, m::Int, fluxf::F) where {F}

Но по-моему это не тот случай.

" Julia will always specialize when the argument is used within the method, but not if the argument is just passed through to another function."
Да, точно не тот, тут проблем быть не должно.

Спасибо всем за обсуждение, если удастся что-то вытянуть из оптимизатора-отпишу. Маленькая ремарка. (По моим скромным и возможно ложным представлениям) джулия должна в этой задаче иметь преимущество за счет JIT, тк может оптимизировать функцию adv_oper, когда уже точно известно, какая fluxf туда пришла. Отчасти это подтверждается тем, что вариант с потоком аргументом и хардкод вариант работа ют одинаково быстро. Фортран в этом месте теряет около 30%. Но возможно, все это не так.

Чтобы наша дискуссия не была совсем уж абстрактной - решение тестовой задачи

Что-то гифка не работает(

Работает, только надо мышку навести

перемешивается

А если @benchmark показывает 1-2 аллокации, это же имеется в виду не Int или ещё что-то типа того, а нечто на куче?

Кому попадались сайты подобные "The Benchmark Game" где примеры или задачи на оптимизацию кода на разных языках есть?

Да, что-то в куче, стек всегда как 0 идёт.

интересно, откуда ) если у меня ни одного массива в коде нет

Внутренние операции могут аллоцировать.
Можно Julia запустить с трекингом аллокаций, должна показать в какой строке они происходят.