Mercurial > octave-dspies

comparison libinterp/corefcn/qz.cc @ 18678:6113e0c6920b

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
maint: Clean up extra spaces before/after parentheses.
author Rik <rik@octave.org>
date Fri, 25 Apr 2014 13:25:25 -0700
parents 175b392e91fe
children
comparison
equal deleted inserted replaced
18677:f684f7075aee 18678:6113e0c6920b
926 // Real zero. 926 // Real zero.
927 #ifdef DEBUG_EIG 927 #ifdef DEBUG_EIG
928 std::cout << " single gen eig:" << std::endl; 928 std::cout << " single gen eig:" << std::endl;
929 std::cout << " alphar(" << jj << ") = " << aa(jj,jj) 929 std::cout << " alphar(" << jj << ") = " << aa(jj,jj)
930 << std::endl; 930 << std::endl;
931 std::cout << " betar( " << jj << ") = " << bb(jj,jj) 931 std::cout << " betar(" << jj << ") = " << bb(jj,jj)
932 << std::endl; 932 << std::endl;
933 std::cout << " alphai(" << jj << ") = 0" << std::endl; 933 std::cout << " alphai(" << jj << ") = 0" << std::endl;
934 #endif 934 #endif
935 935
936 alphar(jj) = aa(jj,jj); 936 alphar(jj) = aa(jj,jj);
1275 %! 1 1 2]; 1275 %! 1 1 2];
1276 %! b = reshape (1:9,3,3); 1276 %! b = reshape (1:9,3,3);
1277 %! [aa, bb, q, z, v, w, lambda] = qz (a, b); 1277 %! [aa, bb, q, z, v, w, lambda] = qz (a, b);
1278 %! sz = length (lambda); 1278 %! sz = length (lambda);
1279 %! observed = (b * v * diag ([lambda;0])) (:, 1:sz); 1279 %! observed = (b * v * diag ([lambda;0])) (:, 1:sz);
1280 %! assert ( (a*v) (:, 1:sz), observed, norm (observed) * 1e-14); 1280 %! assert ((a*v)(:, 1:sz), observed, norm (observed) * 1e-14);
1281 %! observed = (diag ([lambda;0]) * w' * b) (1:sz, :); 1281 %! observed = (diag ([lambda;0]) * w' * b) (1:sz, :);
1282 %! assert ( (w'*a) (1:sz, :) , observed, norm (observed) * 1e-13); 1282 %! assert ((w'*a)(1:sz, :) , observed, norm (observed) * 1e-13);
1283 %! assert (q * a * z, aa, norm (aa) * 1e-14); 1283 %! assert (q * a * z, aa, norm (aa) * 1e-14);
1284 %! assert (q * b * z, bb, norm (bb) * 1e-14); 1284 %! assert (q * b * z, bb, norm (bb) * 1e-14);
1285 1285
1286 %!test 1286 %!test
1287 %! A = [0, 0, -1, 0; 1, 0, 0, 0; -1, 0, -2, -1; 0, -1, 1, 0]; 1287 %! A = [0, 0, -1, 0; 1, 0, 0, 0; -1, 0, -2, -1; 0, -1, 1, 0];