标签: fortran lapack blas
我想执行一系列矩阵乘法,如果我可以避免它,我宁愿不分配新的矩阵。例如,我想计算C = C * B(所以我将C传递为A和C)
A
C
我是否可以这样做(由于dgemm已将A和B标记为in变量?
dgemm
B
in
答案 0 :(得分:3)
我所使用的BLAS实现都没有一个可以就地工作的dgemm。
答案 1 :(得分:2)
不,你永远不能这样做,它与BLAS实现无关。参见
Is it safe to pass GEMV the same output- as input vector to achieve a destructive matrix application?