Commit 9bcb1883 authored by MichelJuillard's avatar MichelJuillard
Browse files

fixing bug related to use of QZ decomposition in Octave

parent 3170d931
......@@ -38,7 +38,7 @@ end
if exist('OCTAVE_VERSION')
[aa,bb,qq,zz]=qz(full(a),full(b));
for j=1:p,
d(:,:,j)=qq'*d(:,:,j)*u;
d(:,:,j)=qq*d(:,:,j)*u;
end
else
[aa,bb,qq,zz]=qz(full(a),full(b),'real'); % available in Matlab version 6.0
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment