Hello,
I am working on optimizing a code that looks at 3D deformation fields. At
one point I have to take the QR decompositions at many points. My problem
is that regardless of how I pre-allocate my variables to store the results
of the QR, they are allocated each time through the loop. I tried to use
the qrfact! method but found it suffered from the same problem. Any help
in getting around this would be great. A basic function that has the same
problem is below.
function test_qr()
tmp=rand(3,3)
a=Array(Float64,2)
b=Array(Float64,2)
for i=1:100000
a,b=qr(tmp)
end
end
