Я собираю некоторый научный код в Fortran 77, и у меня есть дискуссия о том, что будет быстрее.Передать в BLAS или сделать это самостоятельно?
В принципе, у меня есть матрица MxN, назовем ее A. M больше N. Позже в коде мне нужно умножить транспонирование (A) на кучу векторов.
Мой вопрос: было бы быстрее взять A, перенести его самостоятельно и сохранить это, или когда я позвоню BLAS, просто дайте ему флаг транспонирования?
Спасибо! -Patrick
мое чувство кишки говорит мне использовать флаг транспонирования. в этом случае вы делаете много точечных продуктов с шагом одного.
В действительности, это очень сложно сказать без реально действующих кодов. Современный блас использует методы блокировки кеша, которые в лучшем случае затрудняют простой анализ.
Спасибо! У меня было ощущение, что на этот вопрос не будет прямого ответа. Это не будет слишком сложно проверить, вероятно, зависит от размера A и сколько раз я в конечном итоге умножаю A '... – Patrick
@Patrick, если вы проводите тест, вы можете поместить свои результаты в качестве ответа. – Anycorn
вместо умножения на кучу векторов, можете ли вы объединить эти векторы в матрицу? – Anycorn
Это зависит от ситуации. В основном это алгоритм, который сжимает матрицу. После сжатия вы можете быстро применить его к векторам, и именно здесь A 'умножается. В любом случае, это будет зависеть от того, сколько векторов вам дал пользователь. – Patrick