octave: libinterp/corefcn/xdiv.cc annotate

annotate libinterp/corefcn/xdiv.cc @ 21200:fcac5dbbf9ed

maint: Indent #ifdef blocks in libinterp. * builtins.h, Cell.cc, __contourc__.cc, __dispatch__.cc, __dsearchn__.cc, __ichol__.cc, __ilu__.cc, __lin_interpn__.cc, __pchip_deriv__.cc, __qp__.cc, balance.cc, besselj.cc, betainc.cc, bitfcns.cc, bsxfun.cc, c-file-ptr-stream.cc, c-file-ptr-stream.h, cellfun.cc, colloc.cc, comment-list.cc, conv2.cc, daspk.cc, dasrt.cc, dassl.cc, data.cc, debug.cc, defaults.cc, defaults.in.h, defun-dld.h, defun.cc, defun.h, det.cc, dirfns.cc, display.cc, dlmread.cc, dot.cc, dynamic-ld.cc, eig.cc, ellipj.cc, error.cc, errwarn.cc, event-queue.cc, fft.cc, fft2.cc, fftn.cc, file-io.cc, filter.cc, find.cc, gammainc.cc, gcd.cc, getgrent.cc, getpwent.cc, getrusage.cc, givens.cc, gl-render.cc, gl2ps-print.cc, graphics.cc, graphics.in.h, gripes.cc, hash.cc, help.cc, hess.cc, hex2num.cc, input.cc, inv.cc, jit-ir.cc, jit-typeinfo.cc, jit-util.cc, jit-util.h, kron.cc, load-path.cc, load-save.cc, lookup.cc, ls-ascii-helper.cc, ls-hdf5.cc, ls-mat-ascii.cc, ls-mat4.cc, ls-mat5.cc, ls-oct-binary.cc, ls-oct-text.cc, ls-oct-text.h, ls-utils.cc, ls-utils.h, lsode.cc, lu.cc, luinc.cc, mappers.cc, matrix_type.cc, max.cc, mex.h, mexproto.h, mgorth.cc, nproc.cc, oct-errno.in.cc, oct-fstrm.cc, oct-hdf5-types.cc, oct-hdf5.h, oct-hist.cc, oct-iostrm.cc, oct-lvalue.cc, oct-map.cc, oct-prcstrm.cc, oct-procbuf.cc, oct-stream.cc, oct-strstrm.cc, octave-link.cc, ordschur.cc, pager.cc, pinv.cc, pr-output.cc, procstream.cc, profiler.cc, psi.cc, pt-jit.cc, quad.cc, quadcc.cc, qz.cc, rand.cc, rcond.cc, regexp.cc, schur.cc, sighandlers.cc, sparse-xdiv.cc, sparse-xpow.cc, sparse.cc, spparms.cc, sqrtm.cc, str2double.cc, strfind.cc, strfns.cc, sub2ind.cc, svd.cc, sylvester.cc, symtab.cc, syscalls.cc, sysdep.cc, sysdep.h, time.cc, toplev.cc, tril.cc, tsearch.cc, txt-eng-ft.cc, txt-eng.cc, typecast.cc, urlwrite.cc, utils.cc, variables.cc, xdiv.cc, xnorm.cc, xpow.cc, zfstream.cc, __delaunayn__.cc, __eigs__.cc, __fltk_uigetfile__.cc, __glpk__.cc, __init_fltk__.cc, __init_gnuplot__.cc, __magick_read__.cc, __osmesa_print__.cc, __voronoi__.cc, amd.cc, audiodevinfo.cc, audioread.cc, ccolamd.cc, chol.cc, colamd.cc, convhulln.cc, dmperm.cc, fftw.cc, oct-qhull.h, qr.cc, symbfact.cc, symrcm.cc, oct-conf.in.cc, ov-base-diag.cc, ov-base-int.cc, ov-base-mat.cc, ov-base-scalar.cc, ov-base-sparse.cc, ov-base.cc, ov-bool-mat.cc, ov-bool-sparse.cc, ov-bool.cc, ov-builtin.cc, ov-cell.cc, ov-ch-mat.cc, ov-class.cc, ov-classdef.cc, ov-colon.cc, ov-complex.cc, ov-cs-list.cc, ov-cx-diag.cc, ov-cx-mat.cc, ov-cx-sparse.cc, ov-dld-fcn.cc, ov-fcn-handle.cc, ov-fcn-inline.cc, ov-fcn.cc, ov-float.cc, ov-flt-complex.cc, ov-flt-cx-diag.cc, ov-flt-cx-mat.cc, ov-flt-re-diag.cc, ov-flt-re-mat.cc, ov-int16.cc, ov-int32.cc, ov-int64.cc, ov-int8.cc, ov-java.cc, ov-lazy-idx.cc, ov-mex-fcn.cc, ov-null-mat.cc, ov-oncleanup.cc, ov-perm.cc, ov-range.cc, ov-re-diag.cc, ov-re-mat.cc, ov-re-sparse.cc, ov-scalar.cc, ov-str-mat.cc, ov-struct.cc, ov-typeinfo.cc, ov-uint16.cc, ov-uint32.cc, ov-uint64.cc, ov-uint8.cc, ov-usr-fcn.cc, ov.cc, ovl.cc, octave.cc, op-b-b.cc, op-b-bm.cc, op-b-sbm.cc, op-bm-b.cc, op-bm-bm.cc, op-bm-sbm.cc, op-cdm-cdm.cc, op-cell.cc, op-chm.cc, op-class.cc, op-cm-cm.cc, op-cm-cs.cc, op-cm-m.cc, op-cm-s.cc, op-cm-scm.cc, op-cm-sm.cc, op-cs-cm.cc, op-cs-cs.cc, op-cs-m.cc, op-cs-s.cc, op-cs-scm.cc, op-cs-sm.cc, op-dm-dm.cc, op-dm-scm.cc, op-dm-sm.cc, op-dm-template.cc, op-dms-template.cc, op-double-conv.cc, op-fcdm-fcdm.cc, op-fcdm-fdm.cc, op-fcm-fcm.cc, op-fcm-fcs.cc, op-fcm-fm.cc, op-fcm-fs.cc, op-fcn.cc, op-fcs-fcm.cc, op-fcs-fcs.cc, op-fcs-fm.cc, op-fcs-fs.cc, op-fdm-fdm.cc, op-float-conv.cc, op-fm-fcm.cc, op-fm-fcs.cc, op-fm-fm.cc, op-fm-fs.cc, op-fs-fcm.cc, op-fs-fcs.cc, op-fs-fm.cc, op-fs-fs.cc, op-i16-i16.cc, op-i32-i32.cc, op-i64-i64.cc, op-i8-i8.cc, op-int-concat.cc, op-int-conv.cc, op-m-cm.cc, op-m-cs.cc, op-m-m.cc, op-m-s.cc, op-m-scm.cc, op-m-sm.cc, op-pm-pm.cc, op-pm-scm.cc, op-pm-sm.cc, op-pm-template.cc, op-range.cc, op-s-cm.cc, op-s-cs.cc, op-s-m.cc, op-s-s.cc, op-s-scm.cc, op-s-sm.cc, op-sbm-b.cc, op-sbm-bm.cc, op-sbm-sbm.cc, op-scm-cm.cc, op-scm-cs.cc, op-scm-m.cc, op-scm-s.cc, op-scm-scm.cc, op-scm-sm.cc, op-sm-cm.cc, op-sm-cs.cc, op-sm-m.cc, op-sm-s.cc, op-sm-scm.cc, op-sm-sm.cc, op-str-m.cc, op-str-s.cc, op-str-str.cc, op-struct.cc, op-ui16-ui16.cc, op-ui32-ui32.cc, op-ui64-ui64.cc, op-ui8-ui8.cc, pt-arg-list.cc, pt-array-list.cc, pt-assign.cc, pt-binop.cc, pt-bp.cc, pt-cbinop.cc, pt-cell.cc, pt-check.cc, pt-classdef.cc, pt-cmd.cc, pt-colon.cc, pt-colon.h, pt-const.cc, pt-decl.cc, pt-eval.cc, pt-except.cc, pt-exp.cc, pt-fcn-handle.cc, pt-funcall.cc, pt-id.cc, pt-idx.cc, pt-jump.cc, pt-loop.cc, pt-mat.cc, pt-misc.cc, pt-pr-code.cc, pt-select.cc, pt-stmt.cc, pt-unop.cc, pt.cc, token.cc, Array-jit.cc, Array-os.cc, Array-sym.cc, Array-tc.cc, version.cc: Indent #ifdef blocks in libinterp.

author	Rik <rik@octave.org>
date	Fri, 05 Feb 2016 16:29:08 -0800
parents	538b57866b90
children	40de9f8f23a6

rev	line source
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	1 /*
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	2
19697 4197fc428c7d maint: Update copyright notices for 2015. John W. Eaton <jwe@octave.org> parents: 19409 diff changeset	3 Copyright (C) 1993-2015 John W. Eaton
8920 eb63fbe60fab update copyright notices John W. Eaton <jwe@octave.org> parents: 8366 diff changeset	4 Copyright (C) 2008 Jaroslav Hajek
11523 fd0a3ac60b0e update copyright notices John W. Eaton <jwe@octave.org> parents: 10630 diff changeset	5 Copyright (C) 2009-2010 VZLU Prague
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	6
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	7 This file is part of Octave.
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	8
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	9 Octave is free software; you can redistribute it and/or modify it
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	10 under the terms of the GNU General Public License as published by the
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5785 diff changeset	11 Free Software Foundation; either version 3 of the License, or (at your
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5785 diff changeset	12 option) any later version.
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	13
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	14 Octave is distributed in the hope that it will be useful, but WITHOUT
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	15 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	16 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	17 for more details.
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	18
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	19 You should have received a copy of the GNU General Public License
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5785 diff changeset	20 along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5785 diff changeset	21 <http://www.gnu.org/licenses/>.
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	22
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	23 */
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	24
240 a99f28f5e351 [project @ 1993-11-30 20:24:36 by jwe] jwe parents: 164 diff changeset	25 #ifdef HAVE_CONFIG_H
21200 fcac5dbbf9ed maint: Indent #ifdef blocks in libinterp. Rik <rik@octave.org> parents: 21139 diff changeset	26 # include <config.h>
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	27 #endif
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	28
1343 94bedeb289e5 [project @ 1995-09-04 00:29:21 by jwe] jwe parents: 1315 diff changeset	29 #include <cassert>
94bedeb289e5 [project @ 1995-09-04 00:29:21 by jwe] jwe parents: 1315 diff changeset	30
4669 334a27c8f453 [project @ 2003-11-26 07:02:42 by jwe] jwe parents: 4543 diff changeset	31 #include "Array-util.h"
1352 19c10b8657d5 [project @ 1995-09-05 08:11:57 by jwe] jwe parents: 1343 diff changeset	32 #include "CMatrix.h"
453 393e95f46b51 [project @ 1994-06-06 00:05:20 by jwe] jwe parents: 240 diff changeset	33 #include "dMatrix.h"
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	34 #include "CNDArray.h"
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	35 #include "dNDArray.h"
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	36 #include "fCMatrix.h"
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	37 #include "fMatrix.h"
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	38 #include "fCNDArray.h"
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	39 #include "fNDArray.h"
1651 e846e361a265 [project @ 1995-12-20 06:59:12 by jwe] jwe parents: 1488 diff changeset	40 #include "oct-cmplx.h"
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	41 #include "dDiagMatrix.h"
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	42 #include "fDiagMatrix.h"
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	43 #include "CDiagMatrix.h"
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	44 #include "fCDiagMatrix.h"
21100 e39e05d90788 Switch gripe_XXX to either err_XXX or warn_XXX naming scheme. Rik <rik@octave.org> parents: 21024 diff changeset	45 #include "lo-array-errwarn.h"
4153 6b96ce9f5743 [project @ 2002-11-06 20:38:49 by jwe] jwe parents: 3719 diff changeset	46 #include "quit.h"
1352 19c10b8657d5 [project @ 1995-09-05 08:11:57 by jwe] jwe parents: 1343 diff changeset	47
19c10b8657d5 [project @ 1995-09-05 08:11:57 by jwe] jwe parents: 1343 diff changeset	48 #include "error.h"
19c10b8657d5 [project @ 1995-09-05 08:11:57 by jwe] jwe parents: 1343 diff changeset	49 #include "xdiv.h"
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	50
3480 45742a3b1f7c [project @ 2000-01-26 06:16:41 by jwe] jwe parents: 3195 diff changeset	51 static inline bool
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	52 result_ok (octave_idx_type info)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	53 {
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	54 assert (info != -1);
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	55
3480 45742a3b1f7c [project @ 2000-01-26 06:16:41 by jwe] jwe parents: 3195 diff changeset	56 return (info != -2);
45742a3b1f7c [project @ 2000-01-26 06:16:41 by jwe] jwe parents: 3195 diff changeset	57 }
932 1c9812d0cf59 [project @ 1994-11-16 15:40:32 by jwe] jwe parents: 931 diff changeset	58
3480 45742a3b1f7c [project @ 2000-01-26 06:16:41 by jwe] jwe parents: 3195 diff changeset	59 static void
45742a3b1f7c [project @ 2000-01-26 06:16:41 by jwe] jwe parents: 3195 diff changeset	60 solve_singularity_warning (double rcond)
45742a3b1f7c [project @ 2000-01-26 06:16:41 by jwe] jwe parents: 3195 diff changeset	61 {
21112 358aa7fcbd33 Rename errwarn_singular_matrix to warn_singular_matrix. Rik <rik@octave.org> parents: 21100 diff changeset	62 warn_singular_matrix (rcond);
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	63 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	64
21139 538b57866b90 consistently use "typename" intead of "class" in template declarations John W. Eaton <jwe@octave.org> parents: 21112 diff changeset	65 template <typename T1, typename T2>
2364 5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	66 bool
12986 f217edac2c71 fix dimension check for A'\B (bug #33997) John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	67 mx_leftdiv_conform (const T1& a, const T2& b, blas_trans_type blas_trans)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	68 {
12986 f217edac2c71 fix dimension check for A'\B (bug #33997) John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	69 octave_idx_type a_nr = blas_trans == blas_no_trans ? a.rows () : a.cols ();
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	70 octave_idx_type b_nr = b.rows ();
2364 5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	71
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	72 if (a_nr != b_nr)
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	73 {
17787 175b392e91fe Use GNU style coding conventions for code in libinterp/ Rik <rik@octave.org> parents: 17744 diff changeset	74 octave_idx_type a_nc = blas_trans == blas_no_trans ? a.cols ()
175b392e91fe Use GNU style coding conventions for code in libinterp/ Rik <rik@octave.org> parents: 17744 diff changeset	75 : a.rows ();
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	76 octave_idx_type b_nc = b.cols ();
2364 5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	77
21100 e39e05d90788 Switch gripe_XXX to either err_XXX or warn_XXX naming scheme. Rik <rik@octave.org> parents: 21024 diff changeset	78 err_nonconformant ("operator \\", a_nr, a_nc, b_nr, b_nc);
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	79 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	80
2364 5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	81 return true;
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	82 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	83
3195 7a5a5da64756 [project @ 1998-10-28 22:54:04 by jwe] jwe parents: 3180 diff changeset	84 #define INSTANTIATE_MX_LEFTDIV_CONFORM(T1, T2) \
12986 f217edac2c71 fix dimension check for A'\B (bug #33997) John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	85 template bool mx_leftdiv_conform (const T1&, const T2&, blas_trans_type)
3195 7a5a5da64756 [project @ 1998-10-28 22:54:04 by jwe] jwe parents: 3180 diff changeset	86
7a5a5da64756 [project @ 1998-10-28 22:54:04 by jwe] jwe parents: 3180 diff changeset	87 INSTANTIATE_MX_LEFTDIV_CONFORM (Matrix, Matrix);
7a5a5da64756 [project @ 1998-10-28 22:54:04 by jwe] jwe parents: 3180 diff changeset	88 INSTANTIATE_MX_LEFTDIV_CONFORM (Matrix, ComplexMatrix);
7a5a5da64756 [project @ 1998-10-28 22:54:04 by jwe] jwe parents: 3180 diff changeset	89 INSTANTIATE_MX_LEFTDIV_CONFORM (ComplexMatrix, Matrix);
7a5a5da64756 [project @ 1998-10-28 22:54:04 by jwe] jwe parents: 3180 diff changeset	90 INSTANTIATE_MX_LEFTDIV_CONFORM (ComplexMatrix, ComplexMatrix);
2364 5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	91
21139 538b57866b90 consistently use "typename" intead of "class" in template declarations John W. Eaton <jwe@octave.org> parents: 21112 diff changeset	92 template <typename T1, typename T2>
2364 5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	93 bool
3195 7a5a5da64756 [project @ 1998-10-28 22:54:04 by jwe] jwe parents: 3180 diff changeset	94 mx_div_conform (const T1& a, const T2& b)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	95 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	96 octave_idx_type a_nc = a.cols ();
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	97 octave_idx_type b_nc = b.cols ();
2364 5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	98
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	99 if (a_nc != b_nc)
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	100 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	101 octave_idx_type a_nr = a.rows ();
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	102 octave_idx_type b_nr = b.rows ();
2364 5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	103
21100 e39e05d90788 Switch gripe_XXX to either err_XXX or warn_XXX naming scheme. Rik <rik@octave.org> parents: 21024 diff changeset	104 err_nonconformant ("operator /", a_nr, a_nc, b_nr, b_nc);
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	105 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	106
2364 5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	107 return true;
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	108 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	109
3195 7a5a5da64756 [project @ 1998-10-28 22:54:04 by jwe] jwe parents: 3180 diff changeset	110 #define INSTANTIATE_MX_DIV_CONFORM(T1, T2) \
7a5a5da64756 [project @ 1998-10-28 22:54:04 by jwe] jwe parents: 3180 diff changeset	111 template bool mx_div_conform (const T1&, const T2&)
7a5a5da64756 [project @ 1998-10-28 22:54:04 by jwe] jwe parents: 3180 diff changeset	112
7a5a5da64756 [project @ 1998-10-28 22:54:04 by jwe] jwe parents: 3180 diff changeset	113 INSTANTIATE_MX_DIV_CONFORM (Matrix, Matrix);
7a5a5da64756 [project @ 1998-10-28 22:54:04 by jwe] jwe parents: 3180 diff changeset	114 INSTANTIATE_MX_DIV_CONFORM (Matrix, ComplexMatrix);
7a5a5da64756 [project @ 1998-10-28 22:54:04 by jwe] jwe parents: 3180 diff changeset	115 INSTANTIATE_MX_DIV_CONFORM (ComplexMatrix, Matrix);
7a5a5da64756 [project @ 1998-10-28 22:54:04 by jwe] jwe parents: 3180 diff changeset	116 INSTANTIATE_MX_DIV_CONFORM (ComplexMatrix, ComplexMatrix);
2364 5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	117
767 42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	118 // Right division functions.
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	119 //
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	120 // op2 / op1: m cm
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	121 // +-- +---+----+
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	122 // matrix \| 1 \| 3 \|
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	123 // +---+----+
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	124 // complex_matrix \| 2 \| 4 \|
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	125 // +---+----+
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	126
767 42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	127 // -- 1 --
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	128 Matrix
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	129 xdiv (const Matrix& a, const Matrix& b, MatrixType &typ)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	130 {
2364 5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	131 if (! mx_div_conform (a, b))
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	132 return Matrix ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	133
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	134 octave_idx_type info;
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	135 double rcond = 0.0;
3480 45742a3b1f7c [project @ 2000-01-26 06:16:41 by jwe] jwe parents: 3195 diff changeset	136
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	137 Matrix result
12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	138 = b.solve (typ, a.transpose (), info, rcond,
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	139 solve_singularity_warning, true, blas_trans);
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	140
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	141 return result.transpose ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	142 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	143
767 42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	144 // -- 2 --
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	145 ComplexMatrix
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	146 xdiv (const Matrix& a, const ComplexMatrix& b, MatrixType &typ)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	147 {
2364 5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	148 if (! mx_div_conform (a, b))
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	149 return ComplexMatrix ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	150
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	151 octave_idx_type info;
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	152 double rcond = 0.0;
3480 45742a3b1f7c [project @ 2000-01-26 06:16:41 by jwe] jwe parents: 3195 diff changeset	153
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	154 ComplexMatrix result
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	155 = b.solve (typ, a.transpose (), info, rcond,
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	156 solve_singularity_warning, true, blas_trans);
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	157
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	158 return result.transpose ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	159 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	160
767 42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	161 // -- 3 --
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	162 ComplexMatrix
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	163 xdiv (const ComplexMatrix& a, const Matrix& b, MatrixType &typ)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	164 {
2364 5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	165 if (! mx_div_conform (a, b))
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	166 return ComplexMatrix ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	167
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	168 octave_idx_type info;
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	169 double rcond = 0.0;
3480 45742a3b1f7c [project @ 2000-01-26 06:16:41 by jwe] jwe parents: 3195 diff changeset	170
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	171 ComplexMatrix result
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	172 = b.solve (typ, a.transpose (), info, rcond,
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	173 solve_singularity_warning, true, blas_trans);
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	174
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	175 return result.transpose ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	176 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	177
767 42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	178 // -- 4 --
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	179 ComplexMatrix
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	180 xdiv (const ComplexMatrix& a, const ComplexMatrix& b, MatrixType &typ)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	181 {
2364 5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	182 if (! mx_div_conform (a, b))
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	183 return ComplexMatrix ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	184
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	185 octave_idx_type info;
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	186 double rcond = 0.0;
3480 45742a3b1f7c [project @ 2000-01-26 06:16:41 by jwe] jwe parents: 3195 diff changeset	187
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	188 ComplexMatrix result
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	189 = b.solve (typ, a.transpose (), info, rcond,
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	190 solve_singularity_warning, true, blas_trans);
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	191
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	192 return result.transpose ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	193 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	194
767 42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	195 // Funny element by element division operations.
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	196 //
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	197 // op2 \ op1: s cs
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	198 // +-- +---+----+
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	199 // matrix \| 1 \| 3 \|
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	200 // +---+----+
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	201 // complex_matrix \| 2 \| 4 \|
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	202 // +---+----+
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	203
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	204 Matrix
164 e2c950dd96d2 [project @ 1993-10-18 19:32:00 by jwe] jwe parents: 1 diff changeset	205 x_el_div (double a, const Matrix& b)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	206 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	207 octave_idx_type nr = b.rows ();
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	208 octave_idx_type nc = b.columns ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	209
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	210 Matrix result (nr, nc);
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	211
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	212 for (octave_idx_type j = 0; j < nc; j++)
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	213 for (octave_idx_type i = 0; i < nr; i++)
4153 6b96ce9f5743 [project @ 2002-11-06 20:38:49 by jwe] jwe parents: 3719 diff changeset	214 {
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	215 octave_quit ();
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	216 result (i, j) = a / b (i, j);
4153 6b96ce9f5743 [project @ 2002-11-06 20:38:49 by jwe] jwe parents: 3719 diff changeset	217 }
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	218
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	219 return result;
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	220 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	221
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	222 ComplexMatrix
164 e2c950dd96d2 [project @ 1993-10-18 19:32:00 by jwe] jwe parents: 1 diff changeset	223 x_el_div (double a, const ComplexMatrix& b)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	224 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	225 octave_idx_type nr = b.rows ();
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	226 octave_idx_type nc = b.columns ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	227
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	228 ComplexMatrix result (nr, nc);
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	229
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	230 for (octave_idx_type j = 0; j < nc; j++)
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	231 for (octave_idx_type i = 0; i < nr; i++)
4153 6b96ce9f5743 [project @ 2002-11-06 20:38:49 by jwe] jwe parents: 3719 diff changeset	232 {
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	233 octave_quit ();
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	234 result (i, j) = a / b (i, j);
4153 6b96ce9f5743 [project @ 2002-11-06 20:38:49 by jwe] jwe parents: 3719 diff changeset	235 }
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	236
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	237 return result;
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	238 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	239
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	240 ComplexMatrix
164 e2c950dd96d2 [project @ 1993-10-18 19:32:00 by jwe] jwe parents: 1 diff changeset	241 x_el_div (const Complex a, const Matrix& b)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	242 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	243 octave_idx_type nr = b.rows ();
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	244 octave_idx_type nc = b.columns ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	245
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	246 ComplexMatrix result (nr, nc);
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	247
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	248 for (octave_idx_type j = 0; j < nc; j++)
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	249 for (octave_idx_type i = 0; i < nr; i++)
4153 6b96ce9f5743 [project @ 2002-11-06 20:38:49 by jwe] jwe parents: 3719 diff changeset	250 {
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	251 octave_quit ();
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	252 result (i, j) = a / b (i, j);
4153 6b96ce9f5743 [project @ 2002-11-06 20:38:49 by jwe] jwe parents: 3719 diff changeset	253 }
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	254
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	255 return result;
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	256 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	257
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	258 ComplexMatrix
164 e2c950dd96d2 [project @ 1993-10-18 19:32:00 by jwe] jwe parents: 1 diff changeset	259 x_el_div (const Complex a, const ComplexMatrix& b)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	260 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	261 octave_idx_type nr = b.rows ();
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	262 octave_idx_type nc = b.columns ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	263
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	264 ComplexMatrix result (nr, nc);
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	265
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	266 for (octave_idx_type j = 0; j < nc; j++)
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	267 for (octave_idx_type i = 0; i < nr; i++)
4153 6b96ce9f5743 [project @ 2002-11-06 20:38:49 by jwe] jwe parents: 3719 diff changeset	268 {
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	269 octave_quit ();
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	270 result (i, j) = a / b (i, j);
4153 6b96ce9f5743 [project @ 2002-11-06 20:38:49 by jwe] jwe parents: 3719 diff changeset	271 }
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	272
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	273 return result;
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	274 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	275
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	276 // Funny element by element division operations.
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	277 //
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	278 // op2 \ op1: s cs
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	279 // +-- +---+----+
21014 b9d4c3615e89 maint: Capitalize D for dimension in messages and comments. Rik <rik@octave.org> parents: 20232 diff changeset	280 // N-D array \| 1 \| 3 \|
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	281 // +---+----+
21014 b9d4c3615e89 maint: Capitalize D for dimension in messages and comments. Rik <rik@octave.org> parents: 20232 diff changeset	282 // complex N-D array \| 2 \| 4 \|
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	283 // +---+----+
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	284
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	285 NDArray
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	286 x_el_div (double a, const NDArray& b)
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	287 {
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	288 NDArray result (b.dims ());
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	289
20232 a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	290 for (octave_idx_type i = 0; i < b.numel (); i++)
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	291 {
10142 829e69ec3110 make OCTAVE_QUIT a function Jaroslav Hajek <highegg@gmail.com> parents: 9661 diff changeset	292 octave_quit ();
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	293 result (i) = a / b (i);
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	294 }
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	295
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	296 return result;
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	297 }
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	298
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	299 ComplexNDArray
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	300 x_el_div (double a, const ComplexNDArray& b)
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	301 {
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	302 ComplexNDArray result (b.dims ());
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	303
20232 a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	304 for (octave_idx_type i = 0; i < b.numel (); i++)
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	305 {
10142 829e69ec3110 make OCTAVE_QUIT a function Jaroslav Hajek <highegg@gmail.com> parents: 9661 diff changeset	306 octave_quit ();
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	307 result (i) = a / b (i);
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	308 }
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	309
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	310 return result;
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	311 }
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	312
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	313 ComplexNDArray
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	314 x_el_div (const Complex a, const NDArray& b)
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	315 {
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	316 ComplexNDArray result (b.dims ());
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	317
20232 a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	318 for (octave_idx_type i = 0; i < b.numel (); i++)
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	319 {
10142 829e69ec3110 make OCTAVE_QUIT a function Jaroslav Hajek <highegg@gmail.com> parents: 9661 diff changeset	320 octave_quit ();
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	321 result (i) = a / b (i);
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	322 }
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	323
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	324 return result;
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	325 }
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	326
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	327 ComplexNDArray
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	328 x_el_div (const Complex a, const ComplexNDArray& b)
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	329 {
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	330 ComplexNDArray result (b.dims ());
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	331
20232 a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	332 for (octave_idx_type i = 0; i < b.numel (); i++)
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	333 {
10142 829e69ec3110 make OCTAVE_QUIT a function Jaroslav Hajek <highegg@gmail.com> parents: 9661 diff changeset	334 octave_quit ();
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	335 result (i) = a / b (i);
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	336 }
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	337
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	338 return result;
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	339 }
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	340
767 42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	341 // Left division functions.
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	342 //
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	343 // op2 \ op1: m cm
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	344 // +-- +---+----+
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	345 // matrix \| 1 \| 3 \|
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	346 // +---+----+
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	347 // complex_matrix \| 2 \| 4 \|
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	348 // +---+----+
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	349
767 42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	350 // -- 1 --
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	351 Matrix
17787 175b392e91fe Use GNU style coding conventions for code in libinterp/ Rik <rik@octave.org> parents: 17744 diff changeset	352 xleftdiv (const Matrix& a, const Matrix& b, MatrixType &typ,
175b392e91fe Use GNU style coding conventions for code in libinterp/ Rik <rik@octave.org> parents: 17744 diff changeset	353 blas_trans_type transt)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	354 {
12986 f217edac2c71 fix dimension check for A'\B (bug #33997) John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	355 if (! mx_leftdiv_conform (a, b, transt))
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	356 return Matrix ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	357
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	358 octave_idx_type info;
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	359 double rcond = 0.0;
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	360 return a.solve (typ, b, info, rcond, solve_singularity_warning, true, transt);
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	361 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	362
767 42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	363 // -- 2 --
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	364 ComplexMatrix
17787 175b392e91fe Use GNU style coding conventions for code in libinterp/ Rik <rik@octave.org> parents: 17744 diff changeset	365 xleftdiv (const Matrix& a, const ComplexMatrix& b, MatrixType &typ,
175b392e91fe Use GNU style coding conventions for code in libinterp/ Rik <rik@octave.org> parents: 17744 diff changeset	366 blas_trans_type transt)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	367 {
12986 f217edac2c71 fix dimension check for A'\B (bug #33997) John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	368 if (! mx_leftdiv_conform (a, b, transt))
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	369 return ComplexMatrix ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	370
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	371 octave_idx_type info;
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	372 double rcond = 0.0;
3480 45742a3b1f7c [project @ 2000-01-26 06:16:41 by jwe] jwe parents: 3195 diff changeset	373
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	374 return a.solve (typ, b, info, rcond, solve_singularity_warning, true, transt);
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	375 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	376
767 42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	377 // -- 3 --
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	378 ComplexMatrix
17787 175b392e91fe Use GNU style coding conventions for code in libinterp/ Rik <rik@octave.org> parents: 17744 diff changeset	379 xleftdiv (const ComplexMatrix& a, const Matrix& b, MatrixType &typ,
175b392e91fe Use GNU style coding conventions for code in libinterp/ Rik <rik@octave.org> parents: 17744 diff changeset	380 blas_trans_type transt)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	381 {
12986 f217edac2c71 fix dimension check for A'\B (bug #33997) John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	382 if (! mx_leftdiv_conform (a, b, transt))
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	383 return ComplexMatrix ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	384
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	385 octave_idx_type info;
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	386 double rcond = 0.0;
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	387 return a.solve (typ, b, info, rcond, solve_singularity_warning, true, transt);
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	388 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	389
767 42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	390 // -- 4 --
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	391 ComplexMatrix
17787 175b392e91fe Use GNU style coding conventions for code in libinterp/ Rik <rik@octave.org> parents: 17744 diff changeset	392 xleftdiv (const ComplexMatrix& a, const ComplexMatrix& b, MatrixType &typ,
175b392e91fe Use GNU style coding conventions for code in libinterp/ Rik <rik@octave.org> parents: 17744 diff changeset	393 blas_trans_type transt)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	394 {
12986 f217edac2c71 fix dimension check for A'\B (bug #33997) John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	395 if (! mx_leftdiv_conform (a, b, transt))
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	396 return ComplexMatrix ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	397
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	398 octave_idx_type info;
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	399 double rcond = 0.0;
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	400 return a.solve (typ, b, info, rcond, solve_singularity_warning, true, transt);
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	401 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	402
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	403 static void
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	404 solve_singularity_warning (float rcond)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	405 {
21112 358aa7fcbd33 Rename errwarn_singular_matrix to warn_singular_matrix. Rik <rik@octave.org> parents: 21100 diff changeset	406 warn_singular_matrix (rcond);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	407 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	408
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	409 INSTANTIATE_MX_LEFTDIV_CONFORM (FloatMatrix, FloatMatrix);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	410 INSTANTIATE_MX_LEFTDIV_CONFORM (FloatMatrix, FloatComplexMatrix);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	411 INSTANTIATE_MX_LEFTDIV_CONFORM (FloatComplexMatrix, FloatMatrix);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	412 INSTANTIATE_MX_LEFTDIV_CONFORM (FloatComplexMatrix, FloatComplexMatrix);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	413
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	414 INSTANTIATE_MX_DIV_CONFORM (FloatMatrix, FloatMatrix);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	415 INSTANTIATE_MX_DIV_CONFORM (FloatMatrix, FloatComplexMatrix);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	416 INSTANTIATE_MX_DIV_CONFORM (FloatComplexMatrix, FloatMatrix);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	417 INSTANTIATE_MX_DIV_CONFORM (FloatComplexMatrix, FloatComplexMatrix);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	418
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	419 // Right division functions.
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	420 //
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	421 // op2 / op1: m cm
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	422 // +-- +---+----+
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	423 // matrix \| 1 \| 3 \|
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	424 // +---+----+
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	425 // complex_matrix \| 2 \| 4 \|
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	426 // +---+----+
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	427
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	428 // -- 1 --
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	429 FloatMatrix
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	430 xdiv (const FloatMatrix& a, const FloatMatrix& b, MatrixType &typ)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	431 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	432 if (! mx_div_conform (a, b))
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	433 return FloatMatrix ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	434
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	435 octave_idx_type info;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	436 float rcond = 0.0;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	437
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	438 FloatMatrix result
12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	439 = b.solve (typ, a.transpose (), info, rcond,
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	440 solve_singularity_warning, true, blas_trans);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	441
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	442 return result.transpose ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	443 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	444
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	445 // -- 2 --
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	446 FloatComplexMatrix
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	447 xdiv (const FloatMatrix& a, const FloatComplexMatrix& b, MatrixType &typ)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	448 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	449 if (! mx_div_conform (a, b))
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	450 return FloatComplexMatrix ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	451
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	452 octave_idx_type info;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	453 float rcond = 0.0;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	454
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	455 FloatComplexMatrix result
12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	456 = b.solve (typ, a.transpose (), info, rcond,
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	457 solve_singularity_warning, true, blas_trans);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	458
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	459 return result.transpose ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	460 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	461
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	462 // -- 3 --
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	463 FloatComplexMatrix
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	464 xdiv (const FloatComplexMatrix& a, const FloatMatrix& b, MatrixType &typ)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	465 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	466 if (! mx_div_conform (a, b))
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	467 return FloatComplexMatrix ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	468
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	469 octave_idx_type info;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	470 float rcond = 0.0;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	471
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	472 FloatComplexMatrix result
12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	473 = b.solve (typ, a.transpose (), info, rcond,
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	474 solve_singularity_warning, true, blas_trans);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	475
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	476 return result.transpose ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	477 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	478
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	479 // -- 4 --
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	480 FloatComplexMatrix
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	481 xdiv (const FloatComplexMatrix& a, const FloatComplexMatrix& b, MatrixType &typ)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	482 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	483 if (! mx_div_conform (a, b))
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	484 return FloatComplexMatrix ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	485
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	486 octave_idx_type info;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	487 float rcond = 0.0;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	488
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	489 FloatComplexMatrix result
12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	490 = b.solve (typ, a.transpose (), info, rcond,
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	491 solve_singularity_warning, true, blas_trans);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	492
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	493 return result.transpose ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	494 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	495
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	496 // Funny element by element division operations.
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	497 //
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	498 // op2 \ op1: s cs
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	499 // +-- +---+----+
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	500 // matrix \| 1 \| 3 \|
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	501 // +---+----+
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	502 // complex_matrix \| 2 \| 4 \|
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	503 // +---+----+
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	504
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	505 FloatMatrix
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	506 x_el_div (float a, const FloatMatrix& b)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	507 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	508 octave_idx_type nr = b.rows ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	509 octave_idx_type nc = b.columns ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	510
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	511 FloatMatrix result (nr, nc);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	512
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	513 for (octave_idx_type j = 0; j < nc; j++)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	514 for (octave_idx_type i = 0; i < nr; i++)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	515 {
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	516 octave_quit ();
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	517 result (i, j) = a / b (i, j);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	518 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	519
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	520 return result;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	521 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	522
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	523 FloatComplexMatrix
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	524 x_el_div (float a, const FloatComplexMatrix& b)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	525 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	526 octave_idx_type nr = b.rows ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	527 octave_idx_type nc = b.columns ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	528
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	529 FloatComplexMatrix result (nr, nc);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	530
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	531 for (octave_idx_type j = 0; j < nc; j++)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	532 for (octave_idx_type i = 0; i < nr; i++)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	533 {
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	534 octave_quit ();
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	535 result (i, j) = a / b (i, j);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	536 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	537
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	538 return result;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	539 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	540
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	541 FloatComplexMatrix
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	542 x_el_div (const FloatComplex a, const FloatMatrix& b)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	543 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	544 octave_idx_type nr = b.rows ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	545 octave_idx_type nc = b.columns ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	546
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	547 FloatComplexMatrix result (nr, nc);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	548
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	549 for (octave_idx_type j = 0; j < nc; j++)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	550 for (octave_idx_type i = 0; i < nr; i++)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	551 {
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	552 octave_quit ();
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	553 result (i, j) = a / b (i, j);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	554 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	555
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	556 return result;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	557 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	558
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	559 FloatComplexMatrix
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	560 x_el_div (const FloatComplex a, const FloatComplexMatrix& b)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	561 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	562 octave_idx_type nr = b.rows ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	563 octave_idx_type nc = b.columns ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	564
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	565 FloatComplexMatrix result (nr, nc);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	566
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	567 for (octave_idx_type j = 0; j < nc; j++)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	568 for (octave_idx_type i = 0; i < nr; i++)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	569 {
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	570 octave_quit ();
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	571 result (i, j) = a / b (i, j);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	572 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	573
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	574 return result;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	575 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	576
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	577 // Funny element by element division operations.
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	578 //
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	579 // op2 \ op1: s cs
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	580 // +-- +---+----+
21014 b9d4c3615e89 maint: Capitalize D for dimension in messages and comments. Rik <rik@octave.org> parents: 20232 diff changeset	581 // N-D array \| 1 \| 3 \|
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	582 // +---+----+
21014 b9d4c3615e89 maint: Capitalize D for dimension in messages and comments. Rik <rik@octave.org> parents: 20232 diff changeset	583 // complex N-D array \| 2 \| 4 \|
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	584 // +---+----+
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	585
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	586 FloatNDArray
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	587 x_el_div (float a, const FloatNDArray& b)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	588 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	589 FloatNDArray result (b.dims ());
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	590
20232 a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	591 for (octave_idx_type i = 0; i < b.numel (); i++)
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	592 {
10142 829e69ec3110 make OCTAVE_QUIT a function Jaroslav Hajek <highegg@gmail.com> parents: 9661 diff changeset	593 octave_quit ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	594 result (i) = a / b (i);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	595 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	596
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	597 return result;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	598 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	599
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	600 FloatComplexNDArray
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	601 x_el_div (float a, const FloatComplexNDArray& b)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	602 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	603 FloatComplexNDArray result (b.dims ());
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	604
20232 a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	605 for (octave_idx_type i = 0; i < b.numel (); i++)
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	606 {
10142 829e69ec3110 make OCTAVE_QUIT a function Jaroslav Hajek <highegg@gmail.com> parents: 9661 diff changeset	607 octave_quit ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	608 result (i) = a / b (i);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	609 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	610
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	611 return result;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	612 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	613
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	614 FloatComplexNDArray
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	615 x_el_div (const FloatComplex a, const FloatNDArray& b)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	616 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	617 FloatComplexNDArray result (b.dims ());
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	618
20232 a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	619 for (octave_idx_type i = 0; i < b.numel (); i++)
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	620 {
10142 829e69ec3110 make OCTAVE_QUIT a function Jaroslav Hajek <highegg@gmail.com> parents: 9661 diff changeset	621 octave_quit ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	622 result (i) = a / b (i);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	623 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	624
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	625 return result;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	626 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	627
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	628 FloatComplexNDArray
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	629 x_el_div (const FloatComplex a, const FloatComplexNDArray& b)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	630 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	631 FloatComplexNDArray result (b.dims ());
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	632
20232 a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	633 for (octave_idx_type i = 0; i < b.numel (); i++)
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	634 {
10142 829e69ec3110 make OCTAVE_QUIT a function Jaroslav Hajek <highegg@gmail.com> parents: 9661 diff changeset	635 octave_quit ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	636 result (i) = a / b (i);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	637 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	638
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	639 return result;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	640 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	641
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	642 // Left division functions.
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	643 //
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	644 // op2 \ op1: m cm
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	645 // +-- +---+----+
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	646 // matrix \| 1 \| 3 \|
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	647 // +---+----+
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	648 // complex_matrix \| 2 \| 4 \|
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	649 // +---+----+
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	650
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	651 // -- 1 --
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	652 FloatMatrix
17787 175b392e91fe Use GNU style coding conventions for code in libinterp/ Rik <rik@octave.org> parents: 17744 diff changeset	653 xleftdiv (const FloatMatrix& a, const FloatMatrix& b, MatrixType &typ,
175b392e91fe Use GNU style coding conventions for code in libinterp/ Rik <rik@octave.org> parents: 17744 diff changeset	654 blas_trans_type transt)
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	655 {
12986 f217edac2c71 fix dimension check for A'\B (bug #33997) John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	656 if (! mx_leftdiv_conform (a, b, transt))
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	657 return FloatMatrix ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	658
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	659 octave_idx_type info;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	660 float rcond = 0.0;
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	661 return a.solve (typ, b, info, rcond, solve_singularity_warning, true, transt);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	662 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	663
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	664 // -- 2 --
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	665 FloatComplexMatrix
17787 175b392e91fe Use GNU style coding conventions for code in libinterp/ Rik <rik@octave.org> parents: 17744 diff changeset	666 xleftdiv (const FloatMatrix& a, const FloatComplexMatrix& b, MatrixType &typ,
175b392e91fe Use GNU style coding conventions for code in libinterp/ Rik <rik@octave.org> parents: 17744 diff changeset	667 blas_trans_type transt)
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	668 {
12986 f217edac2c71 fix dimension check for A'\B (bug #33997) John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	669 if (! mx_leftdiv_conform (a, b, transt))
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	670 return FloatComplexMatrix ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	671
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	672 octave_idx_type info;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	673 float rcond = 0.0;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	674
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	675 return a.solve (typ, b, info, rcond, solve_singularity_warning, true, transt);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	676 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	677
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	678 // -- 3 --
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	679 FloatComplexMatrix
17787 175b392e91fe Use GNU style coding conventions for code in libinterp/ Rik <rik@octave.org> parents: 17744 diff changeset	680 xleftdiv (const FloatComplexMatrix& a, const FloatMatrix& b, MatrixType &typ,
175b392e91fe Use GNU style coding conventions for code in libinterp/ Rik <rik@octave.org> parents: 17744 diff changeset	681 blas_trans_type transt)
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	682 {
12986 f217edac2c71 fix dimension check for A'\B (bug #33997) John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	683 if (! mx_leftdiv_conform (a, b, transt))
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	684 return FloatComplexMatrix ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	685
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	686 octave_idx_type info;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	687 float rcond = 0.0;
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	688 return a.solve (typ, b, info, rcond, solve_singularity_warning, true, transt);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	689 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	690
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	691 // -- 4 --
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	692 FloatComplexMatrix
17787 175b392e91fe Use GNU style coding conventions for code in libinterp/ Rik <rik@octave.org> parents: 17744 diff changeset	693 xleftdiv (const FloatComplexMatrix& a, const FloatComplexMatrix& b,
175b392e91fe Use GNU style coding conventions for code in libinterp/ Rik <rik@octave.org> parents: 17744 diff changeset	694 MatrixType &typ, blas_trans_type transt)
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	695 {
12986 f217edac2c71 fix dimension check for A'\B (bug #33997) John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	696 if (! mx_leftdiv_conform (a, b, transt))
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	697 return FloatComplexMatrix ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	698
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	699 octave_idx_type info;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	700 float rcond = 0.0;
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	701 return a.solve (typ, b, info, rcond, solve_singularity_warning, true, transt);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	702 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	703
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	704 // Diagonal matrix division.
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	705
21139 538b57866b90 consistently use "typename" intead of "class" in template declarations John W. Eaton <jwe@octave.org> parents: 21112 diff changeset	706 template <typename MT, typename DMT>
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	707 MT
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	708 mdm_div_impl (const MT& a, const DMT& d)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	709 {
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	710 if (! mx_div_conform (a, d))
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	711 return MT ();
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	712
18100 6a71e5030df5 Follow coding convention of defining and initializing only 1 variable per line in liboctinterp. Rik <rik@octave.org> parents: 17787 diff changeset	713 octave_idx_type m = a.rows ();
6a71e5030df5 Follow coding convention of defining and initializing only 1 variable per line in liboctinterp. Rik <rik@octave.org> parents: 17787 diff changeset	714 octave_idx_type n = d.rows ();
6a71e5030df5 Follow coding convention of defining and initializing only 1 variable per line in liboctinterp. Rik <rik@octave.org> parents: 17787 diff changeset	715 octave_idx_type l = d.length ();
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	716 MT x (m, n);
8931 92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	717 typedef typename DMT::element_type S;
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	718 typedef typename MT::element_type T;
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	719 const T *aa = a.data ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	720 const S *dd = d.data ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	721 T *xx = x.fortran_vec ();
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	722
8931 92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	723 for (octave_idx_type j = 0; j < l; j++)
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	724 {
8931 92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	725 const S del = dd[j];
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	726 if (del != S ())
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	727 for (octave_idx_type i = 0; i < m; i++)
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	728 xx[i] = aa[i] / del;
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	729 else
8931 92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	730 for (octave_idx_type i = 0; i < m; i++)
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	731 xx[i] = T ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	732 aa += m; xx += m;
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	733 }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	734
8931 92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	735 for (octave_idx_type i = lm; i < nm; i++)
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	736 xx[i] = T ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	737
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	738 return x;
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	739 }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	740
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	741 // Right division functions.
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	742 //
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	743 // op2 / op1: dm cdm
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	744 // +-- +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	745 // matrix \| 1 \| \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	746 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	747 // complex_matrix \| 2 \| 3 \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	748 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	749
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	750 // -- 1 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	751 Matrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	752 xdiv (const Matrix& a, const DiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	753 { return mdm_div_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	754
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	755 // -- 2 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	756 ComplexMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	757 xdiv (const ComplexMatrix& a, const DiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	758 { return mdm_div_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	759
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	760 // -- 3 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	761 ComplexMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	762 xdiv (const ComplexMatrix& a, const ComplexDiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	763 { return mdm_div_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	764
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	765 // Right division functions, float type.
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	766 //
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	767 // op2 / op1: dm cdm
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	768 // +-- +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	769 // matrix \| 1 \| \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	770 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	771 // complex_matrix \| 2 \| 3 \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	772 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	773
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	774 // -- 1 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	775 FloatMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	776 xdiv (const FloatMatrix& a, const FloatDiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	777 { return mdm_div_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	778
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	779 // -- 2 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	780 FloatComplexMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	781 xdiv (const FloatComplexMatrix& a, const FloatDiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	782 { return mdm_div_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	783
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	784 // -- 3 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	785 FloatComplexMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	786 xdiv (const FloatComplexMatrix& a, const FloatComplexDiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	787 { return mdm_div_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	788
21139 538b57866b90 consistently use "typename" intead of "class" in template declarations John W. Eaton <jwe@octave.org> parents: 21112 diff changeset	789 template <typename MT, typename DMT>
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	790 MT
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	791 dmm_leftdiv_impl (const DMT& d, const MT& a)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	792 {
12986 f217edac2c71 fix dimension check for A'\B (bug #33997) John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	793 if (! mx_leftdiv_conform (d, a, blas_no_trans))
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	794 return MT ();
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	795
18100 6a71e5030df5 Follow coding convention of defining and initializing only 1 variable per line in liboctinterp. Rik <rik@octave.org> parents: 17787 diff changeset	796 octave_idx_type m = d.cols ();
6a71e5030df5 Follow coding convention of defining and initializing only 1 variable per line in liboctinterp. Rik <rik@octave.org> parents: 17787 diff changeset	797 octave_idx_type n = a.cols ();
6a71e5030df5 Follow coding convention of defining and initializing only 1 variable per line in liboctinterp. Rik <rik@octave.org> parents: 17787 diff changeset	798 octave_idx_type k = a.rows ();
6a71e5030df5 Follow coding convention of defining and initializing only 1 variable per line in liboctinterp. Rik <rik@octave.org> parents: 17787 diff changeset	799 octave_idx_type l = d.length ();
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	800 MT x (m, n);
8931 92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	801 typedef typename DMT::element_type S;
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	802 typedef typename MT::element_type T;
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	803 const T *aa = a.data ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	804 const S *dd = d.data ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	805 T *xx = x.fortran_vec ();
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	806
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	807 for (octave_idx_type j = 0; j < n; j++)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	808 {
8931 92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	809 for (octave_idx_type i = 0; i < l; i++)
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	810 xx[i] = dd[i] != S () ? aa[i] / dd[i] : T ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	811 for (octave_idx_type i = l; i < m; i++)
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	812 xx[i] = T ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	813 aa += k; xx += m;
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	814 }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	815
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	816 return x;
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	817 }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	818
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	819 // Left division functions.
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	820 //
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	821 // op2 \ op1: m cm
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	822 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	823 // diag_matrix \| 1 \| 2 \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	824 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	825 // complex_diag_matrix \| \| 3 \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	826 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	827
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	828 // -- 1 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	829 Matrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	830 xleftdiv (const DiagMatrix& a, const Matrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	831 { return dmm_leftdiv_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	832
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	833 // -- 2 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	834 ComplexMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	835 xleftdiv (const DiagMatrix& a, const ComplexMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	836 { return dmm_leftdiv_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	837
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	838 // -- 3 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	839 ComplexMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	840 xleftdiv (const ComplexDiagMatrix& a, const ComplexMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	841 { return dmm_leftdiv_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	842
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	843 // Left division functions, float type.
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	844 //
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	845 // op2 \ op1: m cm
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	846 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	847 // diag_matrix \| 1 \| 2 \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	848 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	849 // complex_diag_matrix \| \| 3 \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	850 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	851
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	852 // -- 1 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	853 FloatMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	854 xleftdiv (const FloatDiagMatrix& a, const FloatMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	855 { return dmm_leftdiv_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	856
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	857 // -- 2 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	858 FloatComplexMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	859 xleftdiv (const FloatDiagMatrix& a, const FloatComplexMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	860 { return dmm_leftdiv_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	861
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	862 // -- 3 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	863 FloatComplexMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	864 xleftdiv (const FloatComplexDiagMatrix& a, const FloatComplexMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	865 { return dmm_leftdiv_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	866
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	867 // Diagonal by diagonal matrix division.
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	868
21139 538b57866b90 consistently use "typename" intead of "class" in template declarations John W. Eaton <jwe@octave.org> parents: 21112 diff changeset	869 template <typename MT, typename DMT>
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	870 MT
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	871 dmdm_div_impl (const MT& a, const DMT& d)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	872 {
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	873 if (! mx_div_conform (a, d))
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	874 return MT ();
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	875
18100 6a71e5030df5 Follow coding convention of defining and initializing only 1 variable per line in liboctinterp. Rik <rik@octave.org> parents: 17787 diff changeset	876 octave_idx_type m = a.rows ();
6a71e5030df5 Follow coding convention of defining and initializing only 1 variable per line in liboctinterp. Rik <rik@octave.org> parents: 17787 diff changeset	877 octave_idx_type n = d.rows ();
6a71e5030df5 Follow coding convention of defining and initializing only 1 variable per line in liboctinterp. Rik <rik@octave.org> parents: 17787 diff changeset	878 octave_idx_type k = d.cols ();
6a71e5030df5 Follow coding convention of defining and initializing only 1 variable per line in liboctinterp. Rik <rik@octave.org> parents: 17787 diff changeset	879 octave_idx_type l = std::min (m, n);
6a71e5030df5 Follow coding convention of defining and initializing only 1 variable per line in liboctinterp. Rik <rik@octave.org> parents: 17787 diff changeset	880 octave_idx_type lk = std::min (l, k);
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	881 MT x (m, n);
8931 92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	882 typedef typename DMT::element_type S;
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	883 typedef typename MT::element_type T;
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	884 const T *aa = a.data ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	885 const S *dd = d.data ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	886 T *xx = x.fortran_vec ();
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	887
8931 92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	888 for (octave_idx_type i = 0; i < lk; i++)
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	889 xx[i] = dd[i] != S () ? aa[i] / dd[i] : T ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	890 for (octave_idx_type i = lk; i < l; i++)
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	891 xx[i] = T ();
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	892
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	893 return x;
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	894 }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	895
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	896 // Right division functions.
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	897 //
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	898 // op2 / op1: dm cdm
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	899 // +-- +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	900 // diag_matrix \| 1 \| \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	901 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	902 // complex_diag_matrix \| 2 \| 3 \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	903 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	904
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	905 // -- 1 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	906 DiagMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	907 xdiv (const DiagMatrix& a, const DiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	908 { return dmdm_div_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	909
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	910 // -- 2 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	911 ComplexDiagMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	912 xdiv (const ComplexDiagMatrix& a, const DiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	913 { return dmdm_div_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	914
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	915 // -- 3 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	916 ComplexDiagMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	917 xdiv (const ComplexDiagMatrix& a, const ComplexDiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	918 { return dmdm_div_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	919
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	920 // Right division functions, float type.
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	921 //
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	922 // op2 / op1: dm cdm
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	923 // +-- +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	924 // diag_matrix \| 1 \| \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	925 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	926 // complex_diag_matrix \| 2 \| 3 \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	927 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	928
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	929 // -- 1 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	930 FloatDiagMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	931 xdiv (const FloatDiagMatrix& a, const FloatDiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	932 { return dmdm_div_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	933
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	934 // -- 2 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	935 FloatComplexDiagMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	936 xdiv (const FloatComplexDiagMatrix& a, const FloatDiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	937 { return dmdm_div_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	938
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	939 // -- 3 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	940 FloatComplexDiagMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	941 xdiv (const FloatComplexDiagMatrix& a, const FloatComplexDiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	942 { return dmdm_div_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	943
21139 538b57866b90 consistently use "typename" intead of "class" in template declarations John W. Eaton <jwe@octave.org> parents: 21112 diff changeset	944 template <typename MT, typename DMT>
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	945 MT
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	946 dmdm_leftdiv_impl (const DMT& d, const MT& a)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	947 {
12986 f217edac2c71 fix dimension check for A'\B (bug #33997) John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	948 if (! mx_leftdiv_conform (d, a, blas_no_trans))
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	949 return MT ();
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	950
18100 6a71e5030df5 Follow coding convention of defining and initializing only 1 variable per line in liboctinterp. Rik <rik@octave.org> parents: 17787 diff changeset	951 octave_idx_type m = d.cols ();
6a71e5030df5 Follow coding convention of defining and initializing only 1 variable per line in liboctinterp. Rik <rik@octave.org> parents: 17787 diff changeset	952 octave_idx_type n = a.cols ();
6a71e5030df5 Follow coding convention of defining and initializing only 1 variable per line in liboctinterp. Rik <rik@octave.org> parents: 17787 diff changeset	953 octave_idx_type k = d.rows ();
6a71e5030df5 Follow coding convention of defining and initializing only 1 variable per line in liboctinterp. Rik <rik@octave.org> parents: 17787 diff changeset	954 octave_idx_type l = std::min (m, n);
6a71e5030df5 Follow coding convention of defining and initializing only 1 variable per line in liboctinterp. Rik <rik@octave.org> parents: 17787 diff changeset	955 octave_idx_type lk = std::min (l, k);
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	956 MT x (m, n);
8931 92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	957 typedef typename DMT::element_type S;
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	958 typedef typename MT::element_type T;
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	959 const T *aa = a.data ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	960 const S *dd = d.data ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	961 T *xx = x.fortran_vec ();
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	962
8931 92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	963 for (octave_idx_type i = 0; i < lk; i++)
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	964 xx[i] = dd[i] != S () ? aa[i] / dd[i] : T ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	965 for (octave_idx_type i = lk; i < l; i++)
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	966 xx[i] = T ();
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	967
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	968 return x;
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	969 }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	970
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	971 // Left division functions.
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	972 //
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	973 // op2 \ op1: dm cdm
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	974 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	975 // diag_matrix \| 1 \| 2 \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	976 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	977 // complex_diag_matrix \| \| 3 \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	978 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	979
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	980 // -- 1 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	981 DiagMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	982 xleftdiv (const DiagMatrix& a, const DiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	983 { return dmdm_leftdiv_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	984
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	985 // -- 2 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	986 ComplexDiagMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	987 xleftdiv (const DiagMatrix& a, const ComplexDiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	988 { return dmdm_leftdiv_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	989
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	990 // -- 3 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	991 ComplexDiagMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	992 xleftdiv (const ComplexDiagMatrix& a, const ComplexDiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	993 { return dmdm_leftdiv_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	994
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	995 // Left division functions, float type.
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	996 //
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	997 // op2 \ op1: dm cdm
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	998 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	999 // diag_matrix \| 1 \| 2 \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	1000 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	1001 // complex_diag_matrix \| \| 3 \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	1002 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	1003
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	1004 // -- 1 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	1005 FloatDiagMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	1006 xleftdiv (const FloatDiagMatrix& a, const FloatDiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	1007 { return dmdm_leftdiv_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	1008
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	1009 // -- 2 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	1010 FloatComplexDiagMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	1011 xleftdiv (const FloatDiagMatrix& a, const FloatComplexDiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	1012 { return dmdm_leftdiv_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	1013
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	1014 // -- 3 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	1015 FloatComplexDiagMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	1016 xleftdiv (const FloatComplexDiagMatrix& a, const FloatComplexDiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	1017 { return dmdm_leftdiv_impl (a, b); }

Mercurial > octave

annotate libinterp/corefcn/xdiv.cc @ 21200:fcac5dbbf9ed