Is there a non-allocating division for two matrices? For example: A = randn(10,10); B = randn(10,10); test(A,B) = A / B @benchmark test($A, $B) Output: BenchmarkTools.Trial: memory estimate: 2.81 KiB allocs estimate: 7 -------------- minimum time: 9.443 μs (0.00% GC) median time…

I found a way without the need of a temporary storage: function rdiv_lu!(A, B) B_lu = lu!(B) rdiv!(rdiv!(A, UpperTriangular(B_lu.factors)), UnitLowerTriangular(B_lu.factors)) _apply_inverse_ipiv_cols!(B_lu, A) end The function _apply_inverse_ipiv_cols! is the equivilant of _apply_inver…

Non-allocating A / B where A and B are matrices

General Usage

Tamas_Papp March 6, 2019, 12:41pm 2

Please read the help, eg ?rdiv!, and if that is not enough provide a self-contained example of how you are trying to use it and what goes wrong.