octave-nkf: src/xdiv.cc annotate

annotate src/xdiv.cc @ 14237:11949c9795a0

Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Add clf() to all demos using plot features to get reproducibility. Use 64 as input to all colormaps (jet (64)) to get reproducibility. * bicubic.m, cell2mat.m, celldisp.m, cplxpair.m, interp1.m, interp2.m, interpft.m, interpn.m, profile.m, profshow.m, convhull.m, delaunay.m, griddata.m, inpolygon.m, voronoi.m, autumn.m, bone.m, contrast.m, cool.m, copper.m, flag.m, gmap40.m, gray.m, hot.m, hsv.m, image.m, imshow.m, jet.m, ocean.m, pink.m, prism.m, rainbow.m, spring.m, summer.m, white.m, winter.m, condest.m, onenormest.m, axis.m, clabel.m, colorbar.m, comet.m, comet3.m, compass.m, contour.m, contour3.m, contourf.m, cylinder.m, daspect.m, ellipsoid.m, errorbar.m, ezcontour.m, ezcontourf.m, ezmesh.m, ezmeshc.m, ezplot.m, ezplot3.m, ezpolar.m, ezsurf.m, ezsurfc.m, feather.m, fill.m, fplot.m, grid.m, hold.m, isosurface.m, legend.m, loglog.m, loglogerr.m, pareto.m, patch.m, pbaspect.m, pcolor.m, pie.m, pie3.m, plot3.m, plotmatrix.m, plotyy.m, polar.m, quiver.m, quiver3.m, rectangle.m, refreshdata.m, ribbon.m, rose.m, scatter.m, scatter3.m, semilogx.m, semilogxerr.m, semilogy.m, semilogyerr.m, shading.m, slice.m, sombrero.m, stairs.m, stem.m, stem3.m, subplot.m, surf.m, surfc.m, surfl.m, surfnorm.m, text.m, title.m, trimesh.m, triplot.m, trisurf.m, uigetdir.m, uigetfile.m, uimenu.m, uiputfile.m, waitbar.m, xlim.m, ylim.m, zlim.m, mkpp.m, pchip.m, polyaffine.m, spline.m, bicgstab.m, cgs.m, gplot.m, pcg.m, pcr.m, treeplot.m, strtok.m, demo.m, example.m, rundemos.m, speed.m, test.m, calendar.m, datestr.m, datetick.m, weekday.m: Revamp %!demos to use Octave coding conventions on spacing, etc.

author	Rik <octave@nomad.inbox5.com>
date	Fri, 20 Jan 2012 12:59:53 -0800
parents	72c96de7a403
children

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
14138 72c96de7a403 maint: update copyright notices for 2012 John W. Eaton <jwe@octave.org> parents: 12986 diff changeset	3 Copyright (C) 1993-2012 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
1192 b6360f2d4fa6 [project @ 1995-03-30 21:38:35 by jwe] jwe parents: 1009 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"
4153 6b96ce9f5743 [project @ 2002-11-06 20:38:49 by jwe] jwe parents: 3719 diff changeset	45 #include "quit.h"
1352 19c10b8657d5 [project @ 1995-09-05 08:11:57 by jwe] jwe parents: 1343 diff changeset	46
19c10b8657d5 [project @ 1995-09-05 08:11:57 by jwe] jwe parents: 1343 diff changeset	47 #include "error.h"
19c10b8657d5 [project @ 1995-09-05 08:11:57 by jwe] jwe parents: 1343 diff changeset	48 #include "xdiv.h"
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	49
3480 45742a3b1f7c [project @ 2000-01-26 06:16:41 by jwe] jwe parents: 3195 diff changeset	50 static inline bool
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	51 result_ok (octave_idx_type info)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	52 {
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	53 assert (info != -1);
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	54
3480 45742a3b1f7c [project @ 2000-01-26 06:16:41 by jwe] jwe parents: 3195 diff changeset	55 return (info != -2);
45742a3b1f7c [project @ 2000-01-26 06:16:41 by jwe] jwe parents: 3195 diff changeset	56 }
932 1c9812d0cf59 [project @ 1994-11-16 15:40:32 by jwe] jwe parents: 931 diff changeset	57
3480 45742a3b1f7c [project @ 2000-01-26 06:16:41 by jwe] jwe parents: 3195 diff changeset	58 static void
45742a3b1f7c [project @ 2000-01-26 06:16:41 by jwe] jwe parents: 3195 diff changeset	59 solve_singularity_warning (double rcond)
45742a3b1f7c [project @ 2000-01-26 06:16:41 by jwe] jwe parents: 3195 diff changeset	60 {
10630 832732c2eb80 use ID for singularity warning in matrix division Jaroslav Hajek <highegg@gmail.com> parents: 10521 diff changeset	61 warning_with_id ("Octave:singular-matrix-div",
832732c2eb80 use ID for singularity warning in matrix division Jaroslav Hajek <highegg@gmail.com> parents: 10521 diff changeset	62 "matrix singular to machine precision, rcond = %g", 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
2364 5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	65 template <class T1, class T2>
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 {
12986 f217edac2c71 fix dimension check for A'\B (bug #33997) John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	74 octave_idx_type a_nc = blas_trans == blas_no_trans ? a.cols () : a.rows ();
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	75 octave_idx_type b_nc = b.cols ();
2364 5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	76
5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	77 gripe_nonconformant ("operator \\", a_nr, a_nc, b_nr, b_nc);
5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	78 return false;
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
5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	92 template <class T1, class T2>
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
5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	104 gripe_nonconformant ("operator /", a_nr, a_nc, b_nr, b_nc);
5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	105 return false;
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	106 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	107
2364 5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	108 return true;
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	109 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	110
3195 7a5a5da64756 [project @ 1998-10-28 22:54:04 by jwe] jwe parents: 3180 diff changeset	111 #define INSTANTIATE_MX_DIV_CONFORM(T1, T2) \
7a5a5da64756 [project @ 1998-10-28 22:54:04 by jwe] jwe parents: 3180 diff changeset	112 template bool mx_div_conform (const T1&, const T2&)
7a5a5da64756 [project @ 1998-10-28 22:54:04 by jwe] jwe parents: 3180 diff changeset	113
7a5a5da64756 [project @ 1998-10-28 22:54:04 by jwe] jwe parents: 3180 diff changeset	114 INSTANTIATE_MX_DIV_CONFORM (Matrix, Matrix);
7a5a5da64756 [project @ 1998-10-28 22:54:04 by jwe] jwe parents: 3180 diff changeset	115 INSTANTIATE_MX_DIV_CONFORM (Matrix, ComplexMatrix);
7a5a5da64756 [project @ 1998-10-28 22:54:04 by jwe] jwe parents: 3180 diff changeset	116 INSTANTIATE_MX_DIV_CONFORM (ComplexMatrix, Matrix);
7a5a5da64756 [project @ 1998-10-28 22:54:04 by jwe] jwe parents: 3180 diff changeset	117 INSTANTIATE_MX_DIV_CONFORM (ComplexMatrix, ComplexMatrix);
2364 5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	118
767 42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	119 // Right division functions.
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	120 //
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	121 // op2 / op1: m cm
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	122 // +-- +---+----+
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	123 // matrix \| 1 \| 3 \|
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	124 // +---+----+
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	125 // complex_matrix \| 2 \| 4 \|
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	126 // +---+----+
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	127
767 42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	128 // -- 1 --
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	129 Matrix
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	130 xdiv (const Matrix& a, const Matrix& b, MatrixType &typ)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	131 {
2364 5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	132 if (! mx_div_conform (a, b))
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	133 return Matrix ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	134
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	135 octave_idx_type info;
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	136 double rcond = 0.0;
3480 45742a3b1f7c [project @ 2000-01-26 06:16:41 by jwe] jwe parents: 3195 diff changeset	137
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	138 Matrix result
12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	139 = b.solve (typ, a.transpose (), info, rcond,
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	140 solve_singularity_warning, true, blas_trans);
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	141
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	142 return result.transpose ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	143 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	144
767 42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	145 // -- 2 --
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	146 ComplexMatrix
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	147 xdiv (const Matrix& a, const ComplexMatrix& b, MatrixType &typ)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	148 {
2364 5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	149 if (! mx_div_conform (a, b))
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	150 return ComplexMatrix ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	151
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	152 octave_idx_type info;
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	153 double rcond = 0.0;
3480 45742a3b1f7c [project @ 2000-01-26 06:16:41 by jwe] jwe parents: 3195 diff changeset	154
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	155 ComplexMatrix result
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	156 = b.solve (typ, a.transpose (), info, rcond,
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	157 solve_singularity_warning, true, blas_trans);
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	158
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	159 return result.transpose ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	160 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	161
767 42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	162 // -- 3 --
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	163 ComplexMatrix
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	164 xdiv (const ComplexMatrix& a, const Matrix& b, MatrixType &typ)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	165 {
2364 5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	166 if (! mx_div_conform (a, b))
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	167 return ComplexMatrix ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	168
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	169 octave_idx_type info;
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	170 double rcond = 0.0;
3480 45742a3b1f7c [project @ 2000-01-26 06:16:41 by jwe] jwe parents: 3195 diff changeset	171
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	172 ComplexMatrix result
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	173 = b.solve (typ, a.transpose (), info, rcond,
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	174 solve_singularity_warning, true, blas_trans);
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	175
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	176 return result.transpose ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	177 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	178
767 42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	179 // -- 4 --
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	180 ComplexMatrix
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	181 xdiv (const ComplexMatrix& a, const ComplexMatrix& b, MatrixType &typ)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	182 {
2364 5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	183 if (! mx_div_conform (a, b))
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	184 return ComplexMatrix ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	185
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	186 octave_idx_type info;
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	187 double rcond = 0.0;
3480 45742a3b1f7c [project @ 2000-01-26 06:16:41 by jwe] jwe parents: 3195 diff changeset	188
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	189 ComplexMatrix result
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	190 = b.solve (typ, a.transpose (), info, rcond,
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	191 solve_singularity_warning, true, blas_trans);
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	192
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	193 return result.transpose ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	194 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	195
767 42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	196 // Funny element by element division operations.
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	197 //
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	198 // op2 \ op1: s cs
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	199 // +-- +---+----+
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	200 // matrix \| 1 \| 3 \|
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	201 // +---+----+
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	202 // complex_matrix \| 2 \| 4 \|
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	203 // +---+----+
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	204
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	205 Matrix
164 e2c950dd96d2 [project @ 1993-10-18 19:32:00 by jwe] jwe parents: 1 diff changeset	206 x_el_div (double a, const Matrix& b)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	207 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	208 octave_idx_type nr = b.rows ();
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	209 octave_idx_type nc = b.columns ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	210
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	211 Matrix result (nr, nc);
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	212
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	213 for (octave_idx_type j = 0; j < nc; j++)
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	214 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	215 {
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	216 octave_quit ();
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	217 result (i, j) = a / b (i, j);
4153 6b96ce9f5743 [project @ 2002-11-06 20:38:49 by jwe] jwe parents: 3719 diff changeset	218 }
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	219
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	220 return result;
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	221 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	222
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	223 ComplexMatrix
164 e2c950dd96d2 [project @ 1993-10-18 19:32:00 by jwe] jwe parents: 1 diff changeset	224 x_el_div (double a, const ComplexMatrix& b)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	225 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	226 octave_idx_type nr = b.rows ();
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	227 octave_idx_type nc = b.columns ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	228
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	229 ComplexMatrix result (nr, nc);
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	230
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	231 for (octave_idx_type j = 0; j < nc; j++)
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	232 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	233 {
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	234 octave_quit ();
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	235 result (i, j) = a / b (i, j);
4153 6b96ce9f5743 [project @ 2002-11-06 20:38:49 by jwe] jwe parents: 3719 diff changeset	236 }
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	237
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	238 return result;
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	239 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	240
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	241 ComplexMatrix
164 e2c950dd96d2 [project @ 1993-10-18 19:32:00 by jwe] jwe parents: 1 diff changeset	242 x_el_div (const Complex a, const Matrix& b)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	243 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	244 octave_idx_type nr = b.rows ();
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	245 octave_idx_type nc = b.columns ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	246
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	247 ComplexMatrix result (nr, nc);
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	248
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	249 for (octave_idx_type j = 0; j < nc; j++)
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	250 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	251 {
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	252 octave_quit ();
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	253 result (i, j) = a / b (i, j);
4153 6b96ce9f5743 [project @ 2002-11-06 20:38:49 by jwe] jwe parents: 3719 diff changeset	254 }
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	255
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	256 return result;
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	257 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	258
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	259 ComplexMatrix
164 e2c950dd96d2 [project @ 1993-10-18 19:32:00 by jwe] jwe parents: 1 diff changeset	260 x_el_div (const Complex a, const ComplexMatrix& b)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	261 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	262 octave_idx_type nr = b.rows ();
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	263 octave_idx_type nc = b.columns ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	264
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	265 ComplexMatrix result (nr, nc);
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	266
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	267 for (octave_idx_type j = 0; j < nc; j++)
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	268 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	269 {
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	270 octave_quit ();
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	271 result (i, j) = a / b (i, j);
4153 6b96ce9f5743 [project @ 2002-11-06 20:38:49 by jwe] jwe parents: 3719 diff changeset	272 }
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	273
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	274 return result;
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	275 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	276
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	277 // Funny element by element division operations.
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	278 //
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	279 // op2 \ op1: s cs
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	280 // +-- +---+----+
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	281 // N-d array \| 1 \| 3 \|
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	282 // +---+----+
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	283 // complex N-d array \| 2 \| 4 \|
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
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	286 NDArray
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	287 x_el_div (double a, const NDArray& b)
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	288 {
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	289 NDArray result (b.dims ());
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	290
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	291 for (octave_idx_type i = 0; i < b.length (); i++)
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	292 {
10142 829e69ec3110 make OCTAVE_QUIT a function Jaroslav Hajek <highegg@gmail.com> parents: 9661 diff changeset	293 octave_quit ();
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	294 result (i) = a / b (i);
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
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	297 return result;
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
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	300 ComplexNDArray
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	301 x_el_div (double a, const ComplexNDArray& b)
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	302 {
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	303 ComplexNDArray result (b.dims ());
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	304
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	305 for (octave_idx_type i = 0; i < b.length (); i++)
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	306 {
10142 829e69ec3110 make OCTAVE_QUIT a function Jaroslav Hajek <highegg@gmail.com> parents: 9661 diff changeset	307 octave_quit ();
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	308 result (i) = a / b (i);
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
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	311 return result;
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
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	314 ComplexNDArray
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	315 x_el_div (const Complex a, const NDArray& b)
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	316 {
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	317 ComplexNDArray result (b.dims ());
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	318
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	319 for (octave_idx_type i = 0; i < b.length (); i++)
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	320 {
10142 829e69ec3110 make OCTAVE_QUIT a function Jaroslav Hajek <highegg@gmail.com> parents: 9661 diff changeset	321 octave_quit ();
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	322 result (i) = a / b (i);
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
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	325 return result;
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
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	328 ComplexNDArray
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	329 x_el_div (const Complex a, const ComplexNDArray& b)
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	330 {
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	331 ComplexNDArray result (b.dims ());
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	332
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	333 for (octave_idx_type i = 0; i < b.length (); i++)
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	334 {
10142 829e69ec3110 make OCTAVE_QUIT a function Jaroslav Hajek <highegg@gmail.com> parents: 9661 diff changeset	335 octave_quit ();
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	336 result (i) = a / b (i);
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
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	339 return result;
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	340 }
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4153 diff changeset	341
767 42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	342 // Left division functions.
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	343 //
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	344 // op2 \ op1: m cm
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	345 // +-- +---+----+
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	346 // matrix \| 1 \| 3 \|
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	347 // +---+----+
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	348 // complex_matrix \| 2 \| 4 \|
42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	349 // +---+----+
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	350
767 42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	351 // -- 1 --
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	352 Matrix
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	353 xleftdiv (const Matrix& a, const Matrix& b, MatrixType &typ, 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
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	365 xleftdiv (const Matrix& a, const ComplexMatrix& b, MatrixType &typ, blas_trans_type transt)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	366 {
12986 f217edac2c71 fix dimension check for A'\B (bug #33997) John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	367 if (! mx_leftdiv_conform (a, b, transt))
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	368 return ComplexMatrix ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	369
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	370 octave_idx_type info;
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	371 double rcond = 0.0;
3480 45742a3b1f7c [project @ 2000-01-26 06:16:41 by jwe] jwe parents: 3195 diff changeset	372
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	373 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	374 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	375
767 42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	376 // -- 3 --
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	377 ComplexMatrix
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	378 xleftdiv (const ComplexMatrix& a, const Matrix& b, MatrixType &typ, blas_trans_type transt)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	379 {
12986 f217edac2c71 fix dimension check for A'\B (bug #33997) John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	380 if (! mx_leftdiv_conform (a, b, transt))
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	381 return ComplexMatrix ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	382
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	383 octave_idx_type info;
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	384 double rcond = 0.0;
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	385 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	386 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	387
767 42731861ee09 [project @ 1994-10-05 21:26:54 by jwe] jwe parents: 453 diff changeset	388 // -- 4 --
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	389 ComplexMatrix
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	390 xleftdiv (const ComplexMatrix& a, const ComplexMatrix& b, MatrixType &typ, blas_trans_type transt)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	391 {
12986 f217edac2c71 fix dimension check for A'\B (bug #33997) John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	392 if (! mx_leftdiv_conform (a, b, transt))
1800 024c75af53f1 [project @ 1996-01-29 05:55:45 by jwe] jwe parents: 1742 diff changeset	393 return ComplexMatrix ();
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	394
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	395 octave_idx_type info;
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	396 double rcond = 0.0;
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	397 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	398 }
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	399
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	400 static void
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	401 solve_singularity_warning (float rcond)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	402 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	403 warning ("matrix singular to machine precision, rcond = %g", rcond);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	404 warning ("attempting to find minimum norm solution");
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	405 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	406
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	407 INSTANTIATE_MX_LEFTDIV_CONFORM (FloatMatrix, FloatMatrix);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	408 INSTANTIATE_MX_LEFTDIV_CONFORM (FloatMatrix, FloatComplexMatrix);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	409 INSTANTIATE_MX_LEFTDIV_CONFORM (FloatComplexMatrix, FloatMatrix);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	410 INSTANTIATE_MX_LEFTDIV_CONFORM (FloatComplexMatrix, FloatComplexMatrix);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	411
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	412 INSTANTIATE_MX_DIV_CONFORM (FloatMatrix, FloatMatrix);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	413 INSTANTIATE_MX_DIV_CONFORM (FloatMatrix, FloatComplexMatrix);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	414 INSTANTIATE_MX_DIV_CONFORM (FloatComplexMatrix, FloatMatrix);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	415 INSTANTIATE_MX_DIV_CONFORM (FloatComplexMatrix, FloatComplexMatrix);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	416
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	417 // Right division functions.
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 // op2 / op1: m cm
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 // matrix \| 1 \| 3 \|
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 // complex_matrix \| 2 \| 4 \|
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
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	426 // -- 1 --
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	427 FloatMatrix
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	428 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	429 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	430 if (! mx_div_conform (a, b))
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	431 return FloatMatrix ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	432
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	433 octave_idx_type info;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	434 float rcond = 0.0;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	435
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	436 FloatMatrix result
12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	437 = b.solve (typ, a.transpose (), info, rcond,
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	438 solve_singularity_warning, true, blas_trans);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	439
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	440 return result.transpose ();
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
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	443 // -- 2 --
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	444 FloatComplexMatrix
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	445 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	446 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	447 if (! mx_div_conform (a, b))
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	448 return FloatComplexMatrix ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	449
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	450 octave_idx_type info;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	451 float rcond = 0.0;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	452
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	453 FloatComplexMatrix result
12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	454 = b.solve (typ, a.transpose (), info, rcond,
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	455 solve_singularity_warning, true, blas_trans);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	456
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	457 return result.transpose ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	458 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	459
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	460 // -- 3 --
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	461 FloatComplexMatrix
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	462 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	463 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	464 if (! mx_div_conform (a, b))
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	465 return FloatComplexMatrix ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	466
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	467 octave_idx_type info;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	468 float rcond = 0.0;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	469
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	470 FloatComplexMatrix result
12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	471 = b.solve (typ, a.transpose (), info, rcond,
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	472 solve_singularity_warning, true, blas_trans);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	473
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	474 return result.transpose ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	475 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	476
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	477 // -- 4 --
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	478 FloatComplexMatrix
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	479 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	480 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	481 if (! mx_div_conform (a, b))
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	482 return FloatComplexMatrix ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	483
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	484 octave_idx_type info;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	485 float rcond = 0.0;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	486
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	487 FloatComplexMatrix result
12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	488 = b.solve (typ, a.transpose (), info, rcond,
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	489 solve_singularity_warning, true, blas_trans);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	490
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	491 return result.transpose ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	492 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	493
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	494 // Funny element by element division operations.
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 // op2 \ op1: s cs
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 // matrix \| 1 \| 3 \|
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 // complex_matrix \| 2 \| 4 \|
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
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	503 FloatMatrix
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	504 x_el_div (float a, const FloatMatrix& b)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	505 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	506 octave_idx_type nr = b.rows ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	507 octave_idx_type nc = b.columns ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	508
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	509 FloatMatrix result (nr, nc);
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 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	512 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	513 {
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	514 octave_quit ();
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	515 result (i, j) = a / b (i, j);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	516 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	517
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	518 return result;
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
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	521 FloatComplexMatrix
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	522 x_el_div (float a, const FloatComplexMatrix& b)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	523 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	524 octave_idx_type nr = b.rows ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	525 octave_idx_type nc = b.columns ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	526
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	527 FloatComplexMatrix result (nr, nc);
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 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	530 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	531 {
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	532 octave_quit ();
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	533 result (i, j) = a / b (i, j);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	534 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	535
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	536 return result;
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
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	539 FloatComplexMatrix
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	540 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	541 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	542 octave_idx_type nr = b.rows ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	543 octave_idx_type nc = b.columns ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	544
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	545 FloatComplexMatrix result (nr, nc);
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 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	548 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	549 {
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	550 octave_quit ();
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	551 result (i, j) = a / b (i, j);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	552 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	553
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	554 return result;
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
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	557 FloatComplexMatrix
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	558 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	559 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	560 octave_idx_type nr = b.rows ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	561 octave_idx_type nc = b.columns ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	562
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	563 FloatComplexMatrix result (nr, nc);
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 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	566 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	567 {
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	568 octave_quit ();
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	569 result (i, j) = a / b (i, j);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	570 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	571
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	572 return result;
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
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	575 // Funny element by element division operations.
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 // op2 \ op1: s cs
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 // N-d array \| 1 \| 3 \|
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	580 // +---+----+
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	581 // complex N-d array \| 2 \| 4 \|
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	582 // +---+----+
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	583
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	584 FloatNDArray
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	585 x_el_div (float a, const FloatNDArray& b)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	586 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	587 FloatNDArray result (b.dims ());
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 for (octave_idx_type i = 0; i < b.length (); i++)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	590 {
10142 829e69ec3110 make OCTAVE_QUIT a function Jaroslav Hajek <highegg@gmail.com> parents: 9661 diff changeset	591 octave_quit ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	592 result (i) = a / b (i);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	593 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	594
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	595 return result;
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
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	598 FloatComplexNDArray
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	599 x_el_div (float a, const FloatComplexNDArray& b)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	600 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	601 FloatComplexNDArray result (b.dims ());
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 for (octave_idx_type i = 0; i < b.length (); i++)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	604 {
10142 829e69ec3110 make OCTAVE_QUIT a function Jaroslav Hajek <highegg@gmail.com> parents: 9661 diff changeset	605 octave_quit ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	606 result (i) = a / b (i);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	607 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	608
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	609 return result;
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
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	612 FloatComplexNDArray
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	613 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	614 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	615 FloatComplexNDArray result (b.dims ());
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 for (octave_idx_type i = 0; i < b.length (); i++)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	618 {
10142 829e69ec3110 make OCTAVE_QUIT a function Jaroslav Hajek <highegg@gmail.com> parents: 9661 diff changeset	619 octave_quit ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	620 result (i) = a / b (i);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	621 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	622
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	623 return result;
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
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	626 FloatComplexNDArray
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	627 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	628 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	629 FloatComplexNDArray result (b.dims ());
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 for (octave_idx_type i = 0; i < b.length (); i++)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	632 {
10142 829e69ec3110 make OCTAVE_QUIT a function Jaroslav Hajek <highegg@gmail.com> parents: 9661 diff changeset	633 octave_quit ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	634 result (i) = a / b (i);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	635 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	636
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	637 return result;
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
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	640 // Left division functions.
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 // op2 \ op1: m cm
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 // matrix \| 1 \| 3 \|
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 // complex_matrix \| 2 \| 4 \|
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
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	649 // -- 1 --
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	650 FloatMatrix
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	651 xleftdiv (const FloatMatrix& a, const FloatMatrix& b, MatrixType &typ, blas_trans_type transt)
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	652 {
12986 f217edac2c71 fix dimension check for A'\B (bug #33997) John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	653 if (! mx_leftdiv_conform (a, b, transt))
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	654 return FloatMatrix ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	655
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	656 octave_idx_type info;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	657 float rcond = 0.0;
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	658 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	659 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	660
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	661 // -- 2 --
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	662 FloatComplexMatrix
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	663 xleftdiv (const FloatMatrix& a, const FloatComplexMatrix& b, MatrixType &typ, blas_trans_type transt)
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	664 {
12986 f217edac2c71 fix dimension check for A'\B (bug #33997) John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	665 if (! mx_leftdiv_conform (a, b, transt))
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	666 return FloatComplexMatrix ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	667
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	668 octave_idx_type info;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	669 float rcond = 0.0;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	670
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	671 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	672 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	673
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	674 // -- 3 --
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	675 FloatComplexMatrix
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	676 xleftdiv (const FloatComplexMatrix& a, const FloatMatrix& b, MatrixType &typ, blas_trans_type transt)
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	677 {
12986 f217edac2c71 fix dimension check for A'\B (bug #33997) John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	678 if (! mx_leftdiv_conform (a, b, transt))
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	679 return FloatComplexMatrix ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	680
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	681 octave_idx_type info;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	682 float rcond = 0.0;
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	683 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	684 }
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 // -- 4 --
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	687 FloatComplexMatrix
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	688 xleftdiv (const FloatComplexMatrix& a, const FloatComplexMatrix& b, MatrixType &typ, blas_trans_type transt)
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	689 {
12986 f217edac2c71 fix dimension check for A'\B (bug #33997) John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	690 if (! mx_leftdiv_conform (a, b, transt))
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	691 return FloatComplexMatrix ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	692
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	693 octave_idx_type info;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	694 float rcond = 0.0;
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	695 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	696 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	697
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	698 // Diagonal matrix division.
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	699
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	700 template <class MT, class DMT>
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	701 MT
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	702 mdm_div_impl (const MT& a, const DMT& d)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	703 {
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	704 if (! mx_div_conform (a, d))
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	705 return MT ();
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	706
8931 92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	707 octave_idx_type m = a.rows (), n = d.rows (), l = d.length ();
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	708 MT x (m, n);
8931 92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	709 typedef typename DMT::element_type S;
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	710 typedef typename MT::element_type T;
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	711 const T *aa = a.data ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	712 const S *dd = d.data ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	713 T *xx = x.fortran_vec ();
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	714
8931 92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	715 for (octave_idx_type j = 0; j < l; j++)
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	716 {
8931 92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	717 const S del = dd[j];
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	718 if (del != S ())
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	719 for (octave_idx_type i = 0; i < m; i++)
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	720 xx[i] = aa[i] / del;
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	721 else
8931 92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	722 for (octave_idx_type i = 0; i < m; i++)
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	723 xx[i] = T ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	724 aa += m; xx += m;
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	725 }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	726
8931 92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	727 for (octave_idx_type i = lm; i < nm; i++)
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	728 xx[i] = T ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	729
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	730 return x;
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	731 }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	732
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	733 // Right division functions.
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	734 //
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	735 // op2 / op1: dm cdm
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	736 // +-- +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	737 // matrix \| 1 \| \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	738 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	739 // complex_matrix \| 2 \| 3 \|
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
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	742 // -- 1 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	743 Matrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	744 xdiv (const Matrix& a, const DiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	745 { return mdm_div_impl (a, b); }
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 // -- 2 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	748 ComplexMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	749 xdiv (const ComplexMatrix& a, const DiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	750 { return mdm_div_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	751
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	752 // -- 3 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	753 ComplexMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	754 xdiv (const ComplexMatrix& a, const ComplexDiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	755 { return mdm_div_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	756
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	757 // Right division functions, float type.
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	758 //
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	759 // op2 / op1: dm cdm
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	760 // +-- +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	761 // matrix \| 1 \| \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	762 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	763 // complex_matrix \| 2 \| 3 \|
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
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	766 // -- 1 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	767 FloatMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	768 xdiv (const FloatMatrix& a, const FloatDiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	769 { return mdm_div_impl (a, b); }
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 // -- 2 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	772 FloatComplexMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	773 xdiv (const FloatComplexMatrix& a, const FloatDiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	774 { return mdm_div_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	775
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	776 // -- 3 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	777 FloatComplexMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	778 xdiv (const FloatComplexMatrix& a, const FloatComplexDiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	779 { return mdm_div_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	780
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	781 template <class MT, class DMT>
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	782 MT
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	783 dmm_leftdiv_impl (const DMT& d, const MT& a)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	784 {
12986 f217edac2c71 fix dimension check for A'\B (bug #33997) John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	785 if (! mx_leftdiv_conform (d, a, blas_no_trans))
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	786 return MT ();
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	787
8931 92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	788 octave_idx_type m = d.cols (), n = a.cols (), k = a.rows (), l = d.length ();
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	789 MT x (m, n);
8931 92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	790 typedef typename DMT::element_type S;
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	791 typedef typename MT::element_type T;
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	792 const T *aa = a.data ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	793 const S *dd = d.data ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	794 T *xx = x.fortran_vec ();
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	795
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	796 for (octave_idx_type j = 0; j < n; j++)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	797 {
8931 92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	798 for (octave_idx_type i = 0; i < l; i++)
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	799 xx[i] = dd[i] != S () ? aa[i] / dd[i] : T ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	800 for (octave_idx_type i = l; i < m; i++)
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	801 xx[i] = T ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	802 aa += k; xx += m;
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	803 }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	804
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	805 return x;
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
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	808 // Left division functions.
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	809 //
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	810 // op2 \ op1: m cm
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	811 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	812 // diag_matrix \| 1 \| 2 \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	813 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	814 // complex_diag_matrix \| \| 3 \|
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
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	817 // -- 1 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	818 Matrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	819 xleftdiv (const DiagMatrix& a, const Matrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	820 { return dmm_leftdiv_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	821
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	822 // -- 2 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	823 ComplexMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	824 xleftdiv (const DiagMatrix& a, const ComplexMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	825 { return dmm_leftdiv_impl (a, b); }
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 // -- 3 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	828 ComplexMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	829 xleftdiv (const ComplexDiagMatrix& a, const ComplexMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	830 { return dmm_leftdiv_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	831
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	832 // Left division functions, float type.
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	833 //
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	834 // op2 \ op1: m cm
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	835 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	836 // diag_matrix \| 1 \| 2 \|
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 // complex_diag_matrix \| \| 3 \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	839 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	840
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	841 // -- 1 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	842 FloatMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	843 xleftdiv (const FloatDiagMatrix& a, const FloatMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	844 { return dmm_leftdiv_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	845
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	846 // -- 2 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	847 FloatComplexMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	848 xleftdiv (const FloatDiagMatrix& a, const FloatComplexMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	849 { return dmm_leftdiv_impl (a, b); }
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 // -- 3 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	852 FloatComplexMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	853 xleftdiv (const FloatComplexDiagMatrix& a, const FloatComplexMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	854 { return dmm_leftdiv_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	855
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	856 // Diagonal by diagonal matrix division.
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	857
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	858 template <class MT, class DMT>
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	859 MT
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	860 dmdm_div_impl (const MT& a, const DMT& d)
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 if (! mx_div_conform (a, d))
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	863 return MT ();
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	864
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	865 octave_idx_type m = a.rows (), n = d.rows (), k = d.cols ();
8931 92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	866 octave_idx_type l = std::min (m, n), lk = std::min (l, k);
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	867 MT x (m, n);
8931 92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	868 typedef typename DMT::element_type S;
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	869 typedef typename MT::element_type T;
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	870 const T *aa = a.data ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	871 const S *dd = d.data ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	872 T *xx = x.fortran_vec ();
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	873
8931 92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	874 for (octave_idx_type i = 0; i < lk; i++)
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	875 xx[i] = dd[i] != S () ? aa[i] / dd[i] : T ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	876 for (octave_idx_type i = lk; i < l; i++)
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	877 xx[i] = T ();
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	878
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	879 return x;
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	880 }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	881
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	882 // Right division functions.
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	883 //
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	884 // op2 / op1: dm cdm
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	885 // +-- +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	886 // diag_matrix \| 1 \| \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	887 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	888 // complex_diag_matrix \| 2 \| 3 \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	889 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	890
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	891 // -- 1 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	892 DiagMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	893 xdiv (const DiagMatrix& a, const DiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	894 { return dmdm_div_impl (a, b); }
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 // -- 2 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	897 ComplexDiagMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	898 xdiv (const ComplexDiagMatrix& a, const DiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	899 { return dmdm_div_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	900
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	901 // -- 3 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	902 ComplexDiagMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	903 xdiv (const ComplexDiagMatrix& a, const ComplexDiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	904 { return dmdm_div_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	905
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	906 // Right division functions, float type.
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	907 //
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	908 // op2 / op1: dm cdm
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 // diag_matrix \| 1 \| \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	911 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	912 // complex_diag_matrix \| 2 \| 3 \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	913 // +---+----+
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 // -- 1 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	916 FloatDiagMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	917 xdiv (const FloatDiagMatrix& a, const FloatDiagMatrix& 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 // -- 2 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	921 FloatComplexDiagMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	922 xdiv (const FloatComplexDiagMatrix& a, const FloatDiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	923 { return dmdm_div_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	924
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	925 // -- 3 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	926 FloatComplexDiagMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	927 xdiv (const FloatComplexDiagMatrix& a, const FloatComplexDiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	928 { return dmdm_div_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	929
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	930 template <class MT, class DMT>
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	931 MT
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	932 dmdm_leftdiv_impl (const DMT& d, const MT& a)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	933 {
12986 f217edac2c71 fix dimension check for A'\B (bug #33997) John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	934 if (! mx_leftdiv_conform (d, a, blas_no_trans))
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	935 return MT ();
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	936
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	937 octave_idx_type m = d.cols (), n = a.cols (), k = d.rows ();
8931 92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	938 octave_idx_type l = std::min (m, n), lk = std::min (l, k);
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	939 MT x (m, n);
8931 92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	940 typedef typename DMT::element_type S;
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	941 typedef typename MT::element_type T;
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	942 const T *aa = a.data ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	943 const S *dd = d.data ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	944 T *xx = x.fortran_vec ();
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	945
8931 92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	946 for (octave_idx_type i = 0; i < lk; i++)
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	947 xx[i] = dd[i] != S () ? aa[i] / dd[i] : T ();
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	948 for (octave_idx_type i = lk; i < l; i++)
92dd386f0f13 optimize diag matrix division Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	949 xx[i] = T ();
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	950
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	951 return x;
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	952 }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	953
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	954 // Left division functions.
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	955 //
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	956 // op2 \ op1: dm cdm
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	957 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	958 // diag_matrix \| 1 \| 2 \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	959 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	960 // complex_diag_matrix \| \| 3 \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	961 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	962
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	963 // -- 1 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	964 DiagMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	965 xleftdiv (const DiagMatrix& a, const DiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	966 { return dmdm_leftdiv_impl (a, b); }
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 // -- 2 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	969 ComplexDiagMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	970 xleftdiv (const DiagMatrix& a, const ComplexDiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	971 { return dmdm_leftdiv_impl (a, b); }
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 // -- 3 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	974 ComplexDiagMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	975 xleftdiv (const ComplexDiagMatrix& a, const ComplexDiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	976 { return dmdm_leftdiv_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	977
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	978 // Left division functions, float type.
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 // op2 \ op1: dm cdm
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	981 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	982 // diag_matrix \| 1 \| 2 \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	983 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	984 // complex_diag_matrix \| \| 3 \|
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	985 // +---+----+
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	986
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	987 // -- 1 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	988 FloatDiagMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	989 xleftdiv (const FloatDiagMatrix& a, const FloatDiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	990 { return dmdm_leftdiv_impl (a, b); }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	991
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	992 // -- 2 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	993 FloatComplexDiagMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	994 xleftdiv (const FloatDiagMatrix& a, const FloatComplexDiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	995 { return dmdm_leftdiv_impl (a, b); }
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 // -- 3 --
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	998 FloatComplexDiagMatrix
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	999 xleftdiv (const FloatComplexDiagMatrix& a, const FloatComplexDiagMatrix& b)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	1000 { return dmdm_leftdiv_impl (a, b); }

Mercurial > hg > octave-nkf

annotate src/xdiv.cc @ 14237:11949c9795a0