Mercurial > octave

comparison libinterp/corefcn/schur.cc @ 29961:7d6709900da7

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more Files affected: __betainc__.cc, __contourc__.cc, __eigs__.cc, __expint__.cc, __ftp__.cc, __gammainc__.cc, __ichol__.cc, __ilu__.cc, __magick_read__.cc, __pchip_deriv__.cc, __qp__.cc, amd.cc, balance.cc, besselj.cc, bsxfun.cc, call-stack.cc, ccolamd.cc, cellfun.cc, chol.cc, colamd.cc, colloc.cc, conv2.cc, daspk.cc, dasrt.cc, dassl.cc, data.cc, defaults.cc, dirfns.cc, display.cc, dlmread.cc, dmperm.cc, dot.cc, eig.cc, ellipj.cc, environment.cc, error.cc, event-manager.cc, fft.cc, fft2.cc, fftn.cc, file-io.cc, find.cc, gcd.cc, getgrent.cc, getpwent.cc, getrusage.cc, gsvd.cc, hash.cc, help.cc, hess.cc, hex2num.cc, input.cc, inv.cc, jsondecode.cc, jsonencode.cc, load-path.cc, load-save.cc, lookup.cc, lsode.cc, lu.cc, max.cc, mgorth.cc, oct-hist.cc, ordqz.cc, ordschur.cc, pager.cc, pr-output.cc, psi.cc, qr.cc, quad.cc, quadcc.cc, qz.cc, rand.cc, regexp.cc, schur.cc, settings.cc, sighandlers.cc, sparse.cc, spparms.cc, sqrtm.cc, stream-euler.cc, strfind.cc, strfns.cc, sub2ind.cc, svd.cc, symbfact.cc, symtab.cc, syscalls.cc, sysdep.cc, time.cc, toplev.cc, tril.cc, typecast.cc, urlwrite.cc, utils.cc, variables.cc, __delaunayn__.cc, __fltk_uigetfile__.cc, __glpk__.cc, __init_gnuplot__.cc, __ode15__.cc, __voronoi__.cc, audiodevinfo.cc, audioread.cc, convhulln.cc, fftw.cc, gzip.cc, ov-cell.cc, ov-class.cc, ov-classdef.cc, ov-fcn-handle.cc, ov-struct.cc, ov-typeinfo.cc, ov-usr-fcn.cc, octave.cc, lex.ll, oct-parse.yy, profiler.cc, andpt-eval.cc.
author John W. Eaton <jwe@octave.org>
date Sat, 14 Aug 2021 22:48:52 -0400
parents 32c3a5805893
children 1e277c6b6626
comparison
equal deleted inserted replaced
29960:939bef0b66e0 29961:7d6709900da7
182 { 182 {
183 FloatMatrix tmp = arg.float_matrix_value (); 183 FloatMatrix tmp = arg.float_matrix_value ();
184 184
185 if (nargout <= 1) 185 if (nargout <= 1)
186 { 186 {
187 octave::math::schur<FloatMatrix> result (tmp, ord, false); 187 math::schur<FloatMatrix> result (tmp, ord, false);
188 retval = ovl (result.schur_matrix ()); 188 retval = ovl (result.schur_matrix ());
189 } 189 }
190 else 190 else
191 { 191 {
192 octave::math::schur<FloatMatrix> result (tmp, ord, true); 192 math::schur<FloatMatrix> result (tmp, ord, true);
193 retval = ovl (result.unitary_matrix (), 193 retval = ovl (result.unitary_matrix (),
194 result.schur_matrix ()); 194 result.schur_matrix ());
195 } 195 }
196 } 196 }
197 else 197 else
198 { 198 {
199 FloatComplexMatrix ctmp = arg.float_complex_matrix_value (); 199 FloatComplexMatrix ctmp = arg.float_complex_matrix_value ();
200 200
201 if (nargout <= 1) 201 if (nargout <= 1)
202 { 202 {
203 octave::math::schur<FloatComplexMatrix> result (ctmp, ord, false); 203 math::schur<FloatComplexMatrix> result (ctmp, ord, false);
204 retval = ovl (mark_upper_triangular (result.schur_matrix ())); 204 retval = ovl (mark_upper_triangular (result.schur_matrix ()));
205 } 205 }
206 else 206 else
207 { 207 {
208 octave::math::schur<FloatComplexMatrix> result (ctmp, ord, true); 208 math::schur<FloatComplexMatrix> result (ctmp, ord, true);
209 retval = ovl (result.unitary_matrix (), 209 retval = ovl (result.unitary_matrix (),
210 mark_upper_triangular (result.schur_matrix ())); 210 mark_upper_triangular (result.schur_matrix ()));
211 } 211 }
212 } 212 }
213 } 213 }
217 { 217 {
218 Matrix tmp = arg.matrix_value (); 218 Matrix tmp = arg.matrix_value ();
219 219
220 if (nargout <= 1) 220 if (nargout <= 1)
221 { 221 {
222 octave::math::schur<Matrix> result (tmp, ord, false); 222 math::schur<Matrix> result (tmp, ord, false);
223 retval = ovl (result.schur_matrix ()); 223 retval = ovl (result.schur_matrix ());
224 } 224 }
225 else 225 else
226 { 226 {
227 octave::math::schur<Matrix> result (tmp, ord, true); 227 math::schur<Matrix> result (tmp, ord, true);
228 retval = ovl (result.unitary_matrix (), 228 retval = ovl (result.unitary_matrix (),
229 result.schur_matrix ()); 229 result.schur_matrix ());
230 } 230 }
231 } 231 }
232 else 232 else
233 { 233 {
234 ComplexMatrix ctmp = arg.complex_matrix_value (); 234 ComplexMatrix ctmp = arg.complex_matrix_value ();
235 235
236 if (nargout <= 1) 236 if (nargout <= 1)
237 { 237 {
238 octave::math::schur<ComplexMatrix> result (ctmp, ord, false); 238 math::schur<ComplexMatrix> result (ctmp, ord, false);
239 retval = ovl (mark_upper_triangular (result.schur_matrix ())); 239 retval = ovl (mark_upper_triangular (result.schur_matrix ()));
240 } 240 }
241 else 241 else
242 { 242 {
243 octave::math::schur<ComplexMatrix> result (ctmp, ord, true); 243 math::schur<ComplexMatrix> result (ctmp, ord, true);
244 retval = ovl (result.unitary_matrix (), 244 retval = ovl (result.unitary_matrix (),
245 mark_upper_triangular (result.schur_matrix ())); 245 mark_upper_triangular (result.schur_matrix ()));
246 } 246 }
247 } 247 }
248 } 248 }
304 if (args(0).is_single_type () || args(1).is_single_type ()) 304 if (args(0).is_single_type () || args(1).is_single_type ())
305 { 305 {
306 FloatMatrix u = args(0).float_matrix_value (); 306 FloatMatrix u = args(0).float_matrix_value ();
307 FloatMatrix t = args(1).float_matrix_value (); 307 FloatMatrix t = args(1).float_matrix_value ();
308 308
309 octave::math::schur<FloatComplexMatrix> cs 309 math::schur<FloatComplexMatrix> cs
310 = octave::math::rsf2csf<FloatComplexMatrix, FloatMatrix> (t, u); 310 = math::rsf2csf<FloatComplexMatrix, FloatMatrix> (t, u);
311 311
312 return ovl (cs.unitary_matrix (), cs.schur_matrix ()); 312 return ovl (cs.unitary_matrix (), cs.schur_matrix ());
313 } 313 }
314 else 314 else
315 { 315 {
316 Matrix u = args(0).matrix_value (); 316 Matrix u = args(0).matrix_value ();
317 Matrix t = args(1).matrix_value (); 317 Matrix t = args(1).matrix_value ();
318 318
319 octave::math::schur<ComplexMatrix> cs 319 math::schur<ComplexMatrix> cs
320 = octave::math::rsf2csf<ComplexMatrix, Matrix> (t, u); 320 = math::rsf2csf<ComplexMatrix, Matrix> (t, u);
321 321
322 return ovl (cs.unitary_matrix (), cs.schur_matrix ()); 322 return ovl (cs.unitary_matrix (), cs.schur_matrix ());
323 } 323 }
324 } 324 }
325 325