octave-nkf: src/xdiv.cc annotate

annotate src/xdiv.cc @ 11542:695141f1c05c ss-3-3-55

snapshot 3.3.55

author	John W. Eaton <jwe@octave.org>
date	Sat, 15 Jan 2011 04:53:04 -0500
parents	fd0a3ac60b0e
children	12df7854fa7c

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
11523 fd0a3ac60b0e update copyright notices John W. Eaton <jwe@octave.org> parents: 10630 diff changeset	3 Copyright (C) 1993-2011 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
3195 7a5a5da64756 [project @ 1998-10-28 22:54:04 by jwe] jwe parents: 3180 diff changeset	67 mx_leftdiv_conform (const T1& a, const T2& b)
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	68 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	69 octave_idx_type a_nr = a.rows ();
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 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4669 diff changeset	74 octave_idx_type a_nc = a.cols ();
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) \
7a5a5da64756 [project @ 1998-10-28 22:54:04 by jwe] jwe parents: 3180 diff changeset	85 template bool mx_leftdiv_conform (const T1&, const T2&)
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
5785 6b9cec830d72 [project @ 2006-05-03 19:32:46 by dbateman] dbateman parents: 5307 diff changeset	138 Matrix result
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	139 = b.solve (typ, a.transpose (), info, rcond,
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
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	156 = b.solve (typ, a.transpose (), info, rcond,
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
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	173 = b.solve (typ, a.transpose (), info, rcond,
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
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	190 = b.solve (typ, a.transpose (), info, rcond,
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 {
2364 5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	355 if (! mx_leftdiv_conform (a, b))
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 {
2364 5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	367 if (! mx_leftdiv_conform (a, b))
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 {
2364 5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	380 if (! mx_leftdiv_conform (a, b))
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 {
2364 5eb0af0730d6 [project @ 1996-10-11 22:57:03 by jwe] jwe parents: 2305 diff changeset	392 if (! mx_leftdiv_conform (a, b))
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
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	436 FloatMatrix result
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	437 = b.solve (typ, a.transpose (), info, rcond,
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
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	453 FloatComplexMatrix result
afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	454 = b.solve (typ, a.transpose (), info, rcond,
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
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	470 FloatComplexMatrix result
afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	471 = b.solve (typ, a.transpose (), info, rcond,
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
9661 afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	487 FloatComplexMatrix result
afcf852256d2 optimize / and '\ for triangular matrices Jaroslav Hajek <highegg@gmail.com> parents: 8931 diff changeset	488 = b.solve (typ, a.transpose (), info, rcond,
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 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	653 if (! mx_leftdiv_conform (a, b))
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 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	665 if (! mx_leftdiv_conform (a, b))
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 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	678 if (! mx_leftdiv_conform (a, b))
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 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	690 if (! mx_leftdiv_conform (a, b))
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 {
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	785 if (! mx_leftdiv_conform (d, a))
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 {
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	934 if (! mx_leftdiv_conform (d, a))
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 @ 11542:695141f1c05c ss-3-3-55