forge: main/fixed/src/fixedMatrix.cc annotate

annotate main/fixed/src/fixedMatrix.cc @ 2914:e2142eefb50d octave-forge

Add missing != to matrix classes

author	adb014
date	Mon, 15 Jan 2007 15:10:24 +0000
parents	1728eb2120ef
children	73fa4496fb07

rev	line source
2383 d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1 /*
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	2
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	3 Copyright (C) 2003 Motorola Inc
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	4 Copyright (C) 2003 David Bateman
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	5
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	6 This program is free software; you can redistribute it and/or modify it
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	7 under the terms of the GNU General Public License as published by the
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	8 Free Software Foundation; either version 2, or (at your option) any
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	9 later version.
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	10
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	11 This program is distributed in the hope that it will be useful, but WITHOUT
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	12 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	13 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	14 for more details.
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	15
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	16 You should have received a copy of the GNU General Public License
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	17 along with this program; see the file COPYING. If not, write to the Free
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	18 Software Foundation, 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	19
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	20 In addition to the terms of the GPL, you are permitted to link
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	21 this program with any Open Source program, as defined by the
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	22 Open Source Initiative (www.opensource.org)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	23
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	24 */
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	25
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	26 #include <iostream>
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	27
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	28 #include <octave/config.h>
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	29 #include <octave/lo-error.h>
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	30 #include <octave/lo-utils.h>
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	31 #include <octave/lo-error.h>
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	32 #include <octave/error.h>
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	33 #include <octave/dMatrix.h>
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	34 #include <octave/gripes.h>
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	35 #include <octave/ops.h>
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	36 #include <octave/quit.h>
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	37
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	38 #include "fixedColVector.h"
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	39 #include "fixedRowVector.h"
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	40 #include "fixedMatrix.h"
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	41 #include "fixedCMatrix.h"
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	42
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	43 // Fixed Point Matrix class.
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	44
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	45 FixedMatrix::FixedMatrix (const MArray2<int> &is, const MArray2<int> &ds)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	46 : MArray2<FixedPoint> (is.rows(), is.cols())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	47 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	48 if ((rows() != ds.rows()) \|\| (cols() != ds.cols())) {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	49 (*current_liboctave_error_handler) ("matrix size mismatch");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	50 return;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	51 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	52
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	53 for (int j = 0; j < cols (); j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	54 for (int i = 0; i < rows (); i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	55 elem (i, j) = FixedPoint((unsigned int)is(i,j), (unsigned int)ds(i,j));
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	56 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	57
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	58
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	59 FixedMatrix::FixedMatrix (const Matrix &is, const Matrix &ds)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	60 : MArray2<FixedPoint> (is.rows(), is.cols())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	61 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	62 if ((rows() != ds.rows()) \|\| (cols() != ds.cols())) {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	63 (*current_liboctave_error_handler) ("matrix size mismatch");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	64 return;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	65 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	66
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	67 for (int j = 0; j < cols (); j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	68 for (int i = 0; i < rows (); i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	69 elem (i, j) = FixedPoint((unsigned int)is(i,j), (unsigned int)ds(i,j));
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	70 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	71
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	72 FixedMatrix::FixedMatrix (const MArray2<int> &is, const MArray2<int> &ds,
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	73 const FixedMatrix& a)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	74 : MArray2<FixedPoint> (a.rows(), a.cols())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	75 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	76 if ((rows() != is.rows()) \|\| (cols() != is.cols()) \|\| (rows() != ds.rows())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	77 \|\| (cols() != ds.cols())) {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	78 (*current_liboctave_error_handler) ("matrix size mismatch");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	79 return;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	80 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	81
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	82 for (int j = 0; j < cols (); j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	83 for (int i = 0; i < rows (); i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	84 elem (i, j) = FixedPoint((unsigned int)is(i,j), (unsigned int)ds(i,j),
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	85 a.elem (i, j));
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	86 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	87
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	88 FixedMatrix::FixedMatrix (const Matrix &is, const Matrix &ds,
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	89 const FixedMatrix& a)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	90 : MArray2<FixedPoint> (a.rows(), a.cols())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	91 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	92 if ((rows() != is.rows()) \|\| (cols() != is.cols()) \|\| (rows() != ds.rows())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	93 \|\| (cols() != ds.cols())) {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	94 (*current_liboctave_error_handler) ("matrix size mismatch");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	95 return;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	96 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	97
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	98 for (int j = 0; j < cols (); j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	99 for (int i = 0; i < rows (); i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	100 elem (i, j) = FixedPoint((unsigned int)is(i,j), (unsigned int)ds(i,j),
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	101 a.elem (i, j));
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	102 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	103
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	104 FixedMatrix::FixedMatrix (unsigned int is, unsigned int ds,
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	105 const FixedMatrix& a)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	106 : MArray2<FixedPoint> (a.rows(), a.cols())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	107 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	108 for (int j = 0; j < cols (); j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	109 for (int i = 0; i < rows (); i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	110 elem (i, j) = FixedPoint(is, ds, a.elem (i, j));
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	111 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	112
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	113 FixedMatrix::FixedMatrix (const MArray2<int> &is, const MArray2<int> &ds,
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	114 const Matrix& a)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	115 : MArray2<FixedPoint> (a.rows(), a.cols())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	116 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	117 if ((rows() != is.rows()) \|\| (cols() != is.cols()) \|\| (rows() != ds.rows())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	118 \|\| (cols() != ds.cols())) {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	119 (*current_liboctave_error_handler) ("matrix size mismatch");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	120 return;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	121 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	122
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	123 for (int j = 0; j < cols (); j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	124 for (int i = 0; i < rows (); i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	125 elem (i, j) = FixedPoint((unsigned int)is(i,j), (unsigned int)ds(i,j),
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	126 a.elem (i, j));
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	127 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	128
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	129 FixedMatrix::FixedMatrix (const Matrix &is, const Matrix &ds,
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	130 const Matrix& a)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	131 : MArray2<FixedPoint> (a.rows(), a.cols())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	132 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	133 if ((rows() != is.rows()) \|\| (cols() != is.cols()) \|\| (rows() != ds.rows())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	134 \|\| (cols() != ds.cols())) {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	135 (*current_liboctave_error_handler) ("matrix size mismatch");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	136 return;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	137 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	138
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	139 for (int j = 0; j < cols (); j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	140 for (int i = 0; i < rows (); i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	141 elem (i, j) = FixedPoint((unsigned int)is(i,j), (unsigned int)ds(i,j),
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	142 a.elem (i, j));
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	143 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	144
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	145 FixedMatrix::FixedMatrix (unsigned int is, unsigned int ds, const Matrix& a)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	146 : MArray2<FixedPoint> (a.rows(), a.cols())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	147 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	148 for (int j = 0; j < cols (); j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	149 for (int i = 0; i < rows (); i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	150 elem (i, j) = FixedPoint(is, ds, a.elem (i, j));
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	151 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	152
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	153 FixedMatrix::FixedMatrix (unsigned int is, unsigned int ds, const Matrix& a,
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	154 const Matrix& b)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	155 : MArray2<FixedPoint> (a.rows(), a.cols())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	156 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	157 if ((rows() != b.rows()) \|\| (cols() != b.cols())) {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	158 (*current_liboctave_error_handler) ("matrix size mismatch");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	159 return;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	160 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	161
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	162 for (int j = 0; j < cols (); j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	163 for (int i = 0; i < rows (); i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	164 elem (i, j) = FixedPoint(is, ds, (unsigned int)a.elem (i, j),
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	165 (unsigned int)b.elem (i,j));
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	166 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	167
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	168 FixedMatrix::FixedMatrix (const MArray2<int> &is, const MArray2<int> &ds,
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	169 const Matrix& a, const Matrix& b)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	170 : MArray2<FixedPoint> (a.rows(), a.cols())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	171 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	172 if ((rows() != b.rows()) \|\| (cols() != b.cols()) \|\| (rows() != is.rows())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	173 \|\| (cols() != is.cols()) \|\| (rows() != ds.rows())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	174 \|\| (cols() != ds.cols())) {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	175 (*current_liboctave_error_handler) ("matrix size mismatch");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	176 return;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	177 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	178
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	179 for (int j = 0; j < cols (); j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	180 for (int i = 0; i < rows (); i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	181 elem (i, j) = FixedPoint((unsigned int)is(i,j), (unsigned int)ds(i,j),
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	182 (unsigned int)a.elem (i, j),
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	183 (unsigned int)b.elem (i,j));
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	184 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	185
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	186 FixedMatrix::FixedMatrix (const Matrix &is, const Matrix &ds, const Matrix& a,
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	187 const Matrix& b)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	188 : MArray2<FixedPoint> (a.rows(), a.cols())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	189 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	190 if ((rows() != b.rows()) \|\| (cols() != b.cols()) \|\| (rows() != is.rows())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	191 \|\| (cols() != is.cols()) \|\| (rows() != ds.rows())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	192 \|\| (cols() != ds.cols())) {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	193 (*current_liboctave_error_handler) ("matrix size mismatch");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	194 return;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	195 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	196
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	197 for (int j = 0; j < cols (); j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	198 for (int i = 0; i < rows (); i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	199 elem (i, j) = FixedPoint((unsigned int)is(i,j), (unsigned int)ds(i,j),
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	200 (unsigned int)a.elem (i, j),
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	201 (unsigned int)b.elem (i,j));
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	202 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	203
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	204 FixedMatrix::FixedMatrix (const MArray2<int> &a)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	205 : MArray2<FixedPoint> (a.rows(), a.cols())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	206 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	207
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	208 for (int j = 0; j < cols (); j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	209 for (int i = 0; i < rows (); i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	210 elem (i, j) = FixedPoint(a.elem (i, j));
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	211 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	212
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	213 FixedMatrix::FixedMatrix (const Matrix &a)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	214 : MArray2<FixedPoint> (a.rows(), a.cols())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	215 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	216
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	217 for (int j = 0; j < cols (); j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	218 for (int i = 0; i < rows (); i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	219 elem (i, j) = FixedPoint(a.elem (i, j));
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	220 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	221
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	222 FixedMatrix::FixedMatrix (const FixedRowVector& rv)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	223 : MArray2<FixedPoint> (1, rv.length (), FixedPoint())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	224 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	225 for (int i = 0; i < rv.length (); i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	226 elem (0, i) = rv.elem (i);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	227 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	228
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	229 FixedMatrix::FixedMatrix (const FixedColumnVector& cv)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	230 : MArray2<FixedPoint> (cv.length (), 1, FixedPoint())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	231 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	232 for (int i = 0; i < cv.length (); i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	233 elem (i, 0) = cv.elem (i);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	234 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	235
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	236 #define GET_FIXED_PROP(METHOD) \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	237 Matrix \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	238 FixedMatrix:: METHOD (void) const \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	239 { \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	240 int nr = rows(); \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	241 int nc = cols(); \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	242 Matrix retval(nr,nc); \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	243 for (int i = 0; i < nr; i++) \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	244 for (int j = 0; j < nc; j++) \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	245 retval(i,j) = (double) elem(i,j) . METHOD (); \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	246 return retval; \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	247 } \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	248
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	249 GET_FIXED_PROP(sign);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	250 GET_FIXED_PROP(signbit);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	251 GET_FIXED_PROP(getdecsize);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	252 GET_FIXED_PROP(getintsize);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	253 GET_FIXED_PROP(getnumber);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	254 GET_FIXED_PROP(fixedpoint);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	255
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	256 #undef GET_FIXED_PROP
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	257
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	258 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	259 FixedMatrix::chdecsize (const double n)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	260 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	261 int nr = rows();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	262 int nc = cols();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	263 FixedMatrix retval(nr,nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	264
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	265 for (int i = 0; i < nr; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	266 for (int j = 0; j < nc; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	267 retval(i,j) = FixedPoint(elem(i,j).getintsize(), (int)n, elem(i,j));
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	268
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	269 return retval;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	270 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	271
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	272 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	273 FixedMatrix::chdecsize (const Matrix &n)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	274 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	275 int nr = rows();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	276 int nc = cols();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	277
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	278 if ((nr != n.rows()) \|\| (nc != n.cols())) {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	279 (*current_liboctave_error_handler) ("matrix size mismatch in chdecsize");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	280 return FixedMatrix();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	281 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	282
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	283 FixedMatrix retval(nr,nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	284
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	285 for (int i = 0; i < nr; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	286 for (int j = 0; j < nc; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	287 retval(i,j) = FixedPoint(elem(i,j).getintsize(), (int)n(i,j), elem(i,j));
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	288
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	289 return retval;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	290 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	291
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	292 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	293 FixedMatrix::chintsize (const double n)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	294 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	295 int nr = rows();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	296 int nc = cols();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	297 FixedMatrix retval(nr,nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	298
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	299 for (int i = 0; i < nr; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	300 for (int j = 0; j < nc; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	301 retval(i,j) = FixedPoint((int)n, elem(i,j).getdecsize(), elem(i,j));
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	302
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	303 return retval;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	304 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	305
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	306 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	307 FixedMatrix::chintsize (const Matrix &n)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	308 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	309 int nr = rows();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	310 int nc = cols();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	311
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	312 if ((nr != n.rows()) \|\| (nc != n.cols())) {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	313 (*current_liboctave_error_handler) ("matrix size mismatch in chintsize");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	314 return FixedMatrix();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	315 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	316
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	317 FixedMatrix retval(nr,nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	318
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	319 for (int i = 0; i < nr; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	320 for (int j = 0; j < nc; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	321 retval(i,j) = FixedPoint((int)n(i,j), elem(i,j).getdecsize(), elem(i,j));
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	322
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	323 return retval;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	324 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	325
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	326 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	327 FixedMatrix::incdecsize (const double n) {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	328 return chdecsize(n + getdecsize());
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	329 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	330
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	331 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	332 FixedMatrix::incdecsize (const Matrix &n) {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	333 if ((n.rows() != rows()) \|\| (n.cols() != cols())) {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	334 (*current_liboctave_error_handler) ("matrix size mismatch in chintsize");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	335 return FixedMatrix();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	336 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	337 return chdecsize(n + getdecsize());
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	338 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	339
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	340 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	341 FixedMatrix::incdecsize () {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	342 return chdecsize(1 + getdecsize());
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	343 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	344
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	345 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	346 FixedMatrix::incintsize (const double n) {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	347 return chintsize(n + getintsize());
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	348 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	349
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	350 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	351 FixedMatrix::incintsize (const Matrix &n) {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	352 if ((n.rows() != rows()) \|\| (n.cols() != cols())) {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	353 (*current_liboctave_error_handler) ("matrix size mismatch in chintsize");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	354 return FixedMatrix();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	355 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	356 return chintsize(n + getintsize());
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	357 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	358
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	359 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	360 FixedMatrix::incintsize () {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	361 return chintsize(1 + getintsize());
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	362 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	363
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	364 bool
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	365 FixedMatrix::operator == (const FixedMatrix& a) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	366 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	367 if (rows () != a.rows () \|\| cols () != a.cols ())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	368 return false;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	369
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	370 for (int i = 0; i < rows(); i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	371 for (int j = 0; j < cols(); j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	372 if (elem(i,j) != a.elem(i,j))
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	373 return false;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	374 return true;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	375 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	376
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	377 bool
2914 e2142eefb50d Add missing != to matrix classes adb014 parents: 2394 diff changeset	378 FixedMatrix::operator != (const FixedMatrix& a) const
e2142eefb50d Add missing != to matrix classes adb014 parents: 2394 diff changeset	379 {
e2142eefb50d Add missing != to matrix classes adb014 parents: 2394 diff changeset	380 return !(*this == a);
e2142eefb50d Add missing != to matrix classes adb014 parents: 2394 diff changeset	381 }
e2142eefb50d Add missing != to matrix classes adb014 parents: 2394 diff changeset	382
e2142eefb50d Add missing != to matrix classes adb014 parents: 2394 diff changeset	383 bool
2383 d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	384 FixedMatrix::is_symmetric (void) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	385 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	386 if (is_square () && rows () > 0)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	387 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	388 for (int i = 0; i < rows (); i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	389 for (int j = i+1; j < cols (); j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	390 if (elem (i, j) != elem (j, i))
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	391 return false;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	392
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	393 return true;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	394 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	395
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	396 return false;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	397 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	398
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	399 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	400 FixedMatrix::concat (const FixedMatrix& rb, const Array<int>& ra_idx)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	401 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	402 if (rb.numel() > 0)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	403 insert (rb, ra_idx(0), ra_idx(1));
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	404 return *this;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	405 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	406
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	407 FixedComplexMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	408 FixedMatrix::concat (const FixedComplexMatrix& rb, const Array<int>& ra_idx)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	409 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	410 FixedComplexMatrix retval (*this);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	411 if (rb.numel() > 0)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	412 retval.insert (rb, ra_idx(0), ra_idx(1));
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	413 return retval;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	414 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	415
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	416 FixedMatrix&
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	417 FixedMatrix::insert (const FixedMatrix& a, int r, int c)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	418 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	419 Array2<FixedPoint>::insert (a, r, c);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	420 return *this;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	421 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	422
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	423 FixedMatrix&
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	424 FixedMatrix::insert (const FixedRowVector& a, int r, int c)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	425 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	426 int a_len = a.length ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	427
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	428 if (r < 0 \|\| r >= rows () \|\| c < 0 \|\| c + a_len > cols ())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	429 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	430 (*current_liboctave_error_handler) ("range error for insert");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	431 return *this;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	432 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	433
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	434 if (a_len > 0)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	435 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	436 make_unique ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	437
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	438 for (int i = 0; i < a_len; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	439 xelem (r, c+i) = a.elem (i);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	440 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	441
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	442 return *this;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	443 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	444
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	445 FixedMatrix&
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	446 FixedMatrix::insert (const FixedColumnVector& a, int r, int c)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	447 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	448 int a_len = a.length ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	449
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	450 if (r < 0 \|\| r + a_len > rows () \|\| c < 0 \|\| c >= cols ())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	451 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	452 (*current_liboctave_error_handler) ("range error for insert");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	453 return *this;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	454 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	455
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	456 if (a_len > 0)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	457 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	458 make_unique ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	459
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	460 for (int i = 0; i < a_len; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	461 xelem (r+i, c) = a.elem (i);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	462 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	463
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	464 return *this;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	465 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	466
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	467 FixedMatrix&
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	468 FixedMatrix::fill (FixedPoint val)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	469 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	470 int nr = rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	471 int nc = cols ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	472
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	473 if (nr > 0 && nc > 0)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	474 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	475 make_unique ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	476
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	477 for (int j = 0; j < nc; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	478 for (int i = 0; i < nr; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	479 xelem (i, j) = val;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	480 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	481
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	482 return *this;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	483 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	484
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	485 FixedMatrix&
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	486 FixedMatrix::fill (FixedPoint val, int r1, int c1, int r2, int c2)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	487 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	488 int nr = rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	489 int nc = cols ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	490
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	491 if (r1 < 0 \|\| r2 < 0 \|\| c1 < 0 \|\| c2 < 0
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	492 \|\| r1 >= nr \|\| r2 >= nr \|\| c1 >= nc \|\| c2 >= nc)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	493 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	494 (*current_liboctave_error_handler) ("range error for fill");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	495 return *this;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	496 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	497
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	498 if (r1 > r2) { int tmp = r1; r1 = r2; r2 = tmp; }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	499 if (c1 > c2) { int tmp = c1; c1 = c2; c2 = tmp; }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	500
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	501 if (r2 >= r1 && c2 >= c1)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	502 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	503 make_unique ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	504
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	505 for (int j = c1; j <= c2; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	506 for (int i = r1; i <= r2; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	507 xelem (i, j) = val;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	508 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	509
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	510 return *this;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	511 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	512
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	513 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	514 FixedMatrix::append (const FixedMatrix& a) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	515 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	516 int nr = rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	517 int nc = cols ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	518 if (nr != a.rows ())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	519 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	520 (*current_liboctave_error_handler) ("row dimension mismatch for append");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	521 return FixedMatrix ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	522 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	523
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	524 int nc_insert = nc;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	525 FixedMatrix retval (nr, nc + a.cols ());
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	526 retval.insert (*this, 0, 0);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	527 retval.insert (a, 0, nc_insert);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	528 return retval;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	529 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	530
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	531 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	532 FixedMatrix::append (const FixedRowVector& a) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	533 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	534 int nr = rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	535 int nc = cols ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	536 if (nr != 1)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	537 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	538 (*current_liboctave_error_handler) ("row dimension mismatch for append");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	539 return FixedMatrix ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	540 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	541
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	542 int nc_insert = nc;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	543 FixedMatrix retval (nr, nc + a.length ());
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	544 retval.insert (*this, 0, 0);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	545 retval.insert (a, 0, nc_insert);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	546 return retval;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	547 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	548
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	549 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	550 FixedMatrix::append (const FixedColumnVector& a) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	551 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	552 int nr = rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	553 int nc = cols ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	554 if (nr != a.length ())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	555 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	556 (*current_liboctave_error_handler) ("row dimension mismatch for append");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	557 return FixedMatrix ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	558 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	559
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	560 int nc_insert = nc;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	561 FixedMatrix retval (nr, nc + 1);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	562 retval.insert (*this, 0, 0);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	563 retval.insert (a, 0, nc_insert);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	564 return retval;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	565 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	566
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	567 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	568 FixedMatrix::stack (const FixedMatrix& a) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	569 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	570 int nr = rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	571 int nc = cols ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	572 if (nc != a.cols ())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	573 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	574 (*current_liboctave_error_handler)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	575 ("column dimension mismatch for stack");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	576 return FixedMatrix ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	577 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	578
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	579 int nr_insert = nr;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	580 FixedMatrix retval (nr + a.rows (), nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	581 retval.insert (*this, 0, 0);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	582 retval.insert (a, nr_insert, 0);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	583 return retval;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	584 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	585
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	586 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	587 FixedMatrix::stack (const FixedRowVector& a) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	588 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	589 int nr = rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	590 int nc = cols ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	591 if (nc != a.length ())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	592 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	593 (*current_liboctave_error_handler)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	594 ("column dimension mismatch for stack");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	595 return FixedMatrix ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	596 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	597
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	598 int nr_insert = nr;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	599 FixedMatrix retval (nr + 1, nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	600 retval.insert (*this, 0, 0);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	601 retval.insert (a, nr_insert, 0);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	602 return retval;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	603 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	604
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	605 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	606 FixedMatrix::stack (const FixedColumnVector& a) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	607 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	608 int nr = rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	609 int nc = cols ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	610 if (nc != 1)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	611 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	612 (*current_liboctave_error_handler)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	613 ("column dimension mismatch for stack");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	614 return FixedMatrix ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	615 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	616
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	617 int nr_insert = nr;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	618 FixedMatrix retval (nr + a.length (), nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	619 retval.insert (*this, 0, 0);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	620 retval.insert (a, nr_insert, 0);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	621 return retval;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	622 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	623
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	624 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	625 FixedMatrix::extract (int r1, int c1, int r2, int c2) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	626 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	627 if (r1 > r2) { int tmp = r1; r1 = r2; r2 = tmp; }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	628 if (c1 > c2) { int tmp = c1; c1 = c2; c2 = tmp; }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	629
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	630 int new_r = r2 - r1 + 1;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	631 int new_c = c2 - c1 + 1;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	632
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	633 FixedMatrix result (new_r, new_c);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	634
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	635 for (int j = 0; j < new_c; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	636 for (int i = 0; i < new_r; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	637 result.xelem (i, j) = elem (r1+i, c1+j);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	638
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	639 return result;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	640 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	641
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	642 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	643 FixedMatrix::extract_n (int r1, int c1, int nr, int nc) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	644 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	645 FixedMatrix result (nr, nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	646
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	647 for (int j = 0; j < nc; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	648 for (int i = 0; i < nr; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	649 result.xelem (i, j) = elem (r1+i, c1+j);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	650
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	651 return result;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	652 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	653
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	654 // extract row or column i.
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	655
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	656 FixedRowVector
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	657 FixedMatrix::row (int i) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	658 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	659 int nc = cols ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	660 if (i < 0 \|\| i >= rows ())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	661 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	662 (*current_liboctave_error_handler) ("invalid row selection");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	663 return FixedRowVector ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	664 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	665
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	666 FixedRowVector retval (nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	667 for (int j = 0; j < nc; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	668 retval.xelem (j) = elem (i, j);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	669
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	670 return retval;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	671 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	672
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	673 FixedRowVector
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	674 FixedMatrix::row (char *s) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	675 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	676 if (! s)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	677 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	678 (*current_liboctave_error_handler) ("invalid row selection");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	679 return FixedRowVector ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	680 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	681
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	682 char c = *s;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	683 if (c == 'f' \|\| c == 'F')
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	684 return row (0);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	685 else if (c == 'l' \|\| c == 'L')
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	686 return row (rows () - 1);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	687 else
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	688 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	689 (*current_liboctave_error_handler) ("invalid row selection");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	690 return FixedRowVector ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	691 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	692 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	693
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	694 FixedColumnVector
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	695 FixedMatrix::column (int i) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	696 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	697 int nr = rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	698 if (i < 0 \|\| i >= cols ())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	699 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	700 (*current_liboctave_error_handler) ("invalid column selection");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	701 return FixedColumnVector ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	702 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	703
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	704 FixedColumnVector retval (nr);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	705 for (int j = 0; j < nr; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	706 retval.xelem (j) = elem (j, i);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	707
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	708 return retval;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	709 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	710
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	711 FixedColumnVector
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	712 FixedMatrix::column (char *s) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	713 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	714 if (! s)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	715 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	716 (*current_liboctave_error_handler) ("invalid column selection");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	717 return FixedColumnVector ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	718 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	719
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	720 char c = *s;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	721 if (c == 'f' \|\| c == 'F')
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	722 return column (0);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	723 else if (c == 'l' \|\| c == 'L')
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	724 return column (cols () - 1);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	725 else
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	726 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	727 (*current_liboctave_error_handler) ("invalid column selection");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	728 return FixedColumnVector ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	729 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	730 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	731
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	732 // unary operations
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	733
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	734 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	735 FixedMatrix::operator ! (void) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	736 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	737 int nr = rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	738 int nc = cols ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	739
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	740 FixedMatrix b (nr, nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	741
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	742 for (int j = 0; j < nc; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	743 for (int i = 0; i < nr; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	744 b.elem (i, j) = ! elem (i, j) ;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	745
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	746 return b;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	747 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	748
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	749 // column vector by row vector -> matrix operations
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	750
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	751 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	752 operator * (const FixedColumnVector& v, const FixedRowVector& a)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	753 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	754 FixedMatrix retval;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	755
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	756 int len = v.length ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	757
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	758 if (len != 0)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	759 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	760 int a_len = a.length ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	761
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	762 retval.resize (len, a_len);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	763
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	764 for (int i = 0; i < len; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	765 for (int j = 0; j < a_len; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	766 retval.elem(j,i) = v.elem(i) * a.elem(j);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	767 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	768
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	769 return retval;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	770 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	771
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	772 // other operations.
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	773
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	774 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	775 FixedMatrix::map (f_f_Mapper f) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	776 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	777 FixedMatrix b (*this);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	778 return b.apply (f);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	779 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	780
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	781 FixedMatrix&
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	782 FixedMatrix::apply (f_f_Mapper f)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	783 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	784 FixedPoint *d = fortran_vec (); // Ensures only one reference to my privates!
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	785
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	786 for (int i = 0; i < length (); i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	787 d[i] = f (d[i]);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	788
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	789 return *this;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	790 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	791
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	792 boolMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	793 FixedMatrix::all (int dim) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	794 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	795 #define ROW_EXPR \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	796 if (elem (i, j) .fixedpoint () == 0.0) \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	797 { \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	798 retval.elem (i, 0) = false; \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	799 break; \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	800 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	801 #define COL_EXPR \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	802 if (elem (i, j) .fixedpoint () == 0.0) \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	803 { \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	804 retval.elem (0, j) = false; \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	805 break; \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	806 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	807
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	808 MX_BASE_REDUCTION_OP (boolMatrix, ROW_EXPR, COL_EXPR, true, true);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	809
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	810 #undef ROW_EXPR
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	811 #undef COL_EXPR
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	812 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	813
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	814 boolMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	815 FixedMatrix::any (int dim) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	816 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	817 #define ROW_EXPR \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	818 if (elem (i, j) .fixedpoint () != 0.0) \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	819 { \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	820 retval.elem (i, 0) = true; \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	821 break; \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	822 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	823 #define COL_EXPR \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	824 if (elem (i, j) .fixedpoint () != 0.0) \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	825 { \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	826 retval.elem (0, j) = true; \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	827 break; \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	828 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	829
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	830 MX_BASE_REDUCTION_OP (boolMatrix, ROW_EXPR, COL_EXPR, false, false);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	831
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	832 #undef ROW_EXPR
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	833 #undef COL_EXPR
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	834 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	835
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	836 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	837 FixedMatrix::cumprod (int dim) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	838 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	839 MX_CUMULATIVE_OP (FixedMatrix, FixedPoint, *=);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	840 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	841
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	842 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	843 FixedMatrix::cumsum (int dim) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	844 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	845 MX_CUMULATIVE_OP (FixedMatrix, FixedPoint, +=);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	846 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	847
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	848 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	849 FixedMatrix::prod (int dim) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	850 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	851 FixedPoint one(1,0,1,0);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	852 MX_REDUCTION_OP (FixedMatrix, *=, one, one);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	853 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	854
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	855 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	856 FixedMatrix::sum (int dim) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	857 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	858 FixedPoint zero;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	859 MX_REDUCTION_OP (FixedMatrix, +=, zero, zero);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	860 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	861
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	862 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	863 FixedMatrix::sumsq (int dim) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	864 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	865 FixedPoint zero;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	866 #define ROW_EXPR \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	867 FixedPoint d = elem (i, j); \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	868 retval.elem (i, 0) += d * d
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	869
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	870 #define COL_EXPR \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	871 FixedPoint d = elem (i, j); \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	872 retval.elem (0, j) += d * d
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	873
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	874 MX_BASE_REDUCTION_OP (FixedMatrix, ROW_EXPR, COL_EXPR, zero, zero);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	875
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	876 #undef ROW_EXPR
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	877 #undef COL_EXPR
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	878 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	879
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	880 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	881 FixedMatrix::abs (void) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	882 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	883 int nr = rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	884 int nc = cols ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	885
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	886 FixedMatrix retval (nr, nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	887
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	888 for (int j = 0; j < nc; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	889 for (int i = 0; i < nr; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	890 retval (i, j) = ::abs(elem (i, j));
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	891
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	892 return retval;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	893 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	894
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	895 FixedColumnVector
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	896 FixedMatrix::diag (void) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	897 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	898 return diag (0);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	899 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	900
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	901 FixedColumnVector
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	902 FixedMatrix::diag (int k) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	903 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	904 int nnr = rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	905 int nnc = cols ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	906 if (k > 0)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	907 nnc -= k;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	908 else if (k < 0)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	909 nnr += k;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	910
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	911 FixedColumnVector d;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	912
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	913 if (nnr > 0 && nnc > 0)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	914 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	915 int ndiag = (nnr < nnc) ? nnr : nnc;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	916
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	917 d.resize (ndiag);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	918
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	919 if (k > 0)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	920 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	921 for (int i = 0; i < ndiag; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	922 d.elem (i) = elem (i, i+k);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	923 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	924 else if ( k < 0)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	925 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	926 for (int i = 0; i < ndiag; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	927 d.elem (i) = elem (i-k, i);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	928 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	929 else
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	930 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	931 for (int i = 0; i < ndiag; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	932 d.elem (i) = elem (i, i);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	933 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	934 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	935 else
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	936 std::cerr << "diag: requested diagonal out of range\n";
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	937
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	938 return d;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	939 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	940
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	941 FixedColumnVector
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	942 FixedMatrix::row_min (void) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	943 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	944 Array<int> index;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	945 return row_min (index);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	946 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	947
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	948 FixedColumnVector
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	949 FixedMatrix::row_min (Array<int>& index) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	950 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	951 FixedColumnVector result;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	952
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	953 int nr = rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	954 int nc = cols ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	955
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	956 if (nr > 0 && nc > 0)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	957 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	958 result.resize (nr);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	959 index.resize (nr);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	960
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	961 for (int i = 0; i < nr; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	962 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	963 int idx_j = 0;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	964
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	965 FixedPoint tmp_min = elem (i, idx_j);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	966
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	967 for (int j = 1; j < nc; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	968 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	969 FixedPoint tmp = elem (i, j);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	970 if (tmp < tmp_min)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	971 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	972 idx_j = j;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	973 tmp_min = tmp;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	974 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	975 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	976
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	977 result.elem (i) = tmp_min;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	978 index.elem (i) = idx_j;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	979 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	980 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	981
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	982 return result;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	983 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	984
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	985 FixedColumnVector
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	986 FixedMatrix::row_max (void) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	987 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	988 Array<int> index;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	989 return row_max (index);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	990 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	991
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	992 FixedColumnVector
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	993 FixedMatrix::row_max (Array<int>& index) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	994 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	995 FixedColumnVector result;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	996
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	997 int nr = rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	998 int nc = cols ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	999
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1000 if (nr > 0 && nc > 0)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1001 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1002 result.resize (nr);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1003 index.resize (nr);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1004
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1005 for (int i = 0; i < nr; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1006 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1007 int idx_j = 0;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1008
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1009 FixedPoint tmp_max = elem (i, idx_j);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1010
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1011 for (int j = 1; j < nc; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1012 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1013 FixedPoint tmp = elem (i, j);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1014 if (tmp > tmp_max)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1015 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1016 idx_j = j;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1017 tmp_max = tmp;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1018 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1019 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1020
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1021 result.elem (i) = tmp_max;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1022 index.elem (i) = idx_j;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1023 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1024 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1025
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1026 return result;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1027 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1028
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1029 FixedRowVector
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1030 FixedMatrix::column_min (void) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1031 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1032 Array<int> index;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1033 return column_min (index);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1034 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1035
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1036 FixedRowVector
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1037 FixedMatrix::column_min (Array<int>& index) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1038 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1039 FixedRowVector result;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1040
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1041 int nr = rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1042 int nc = cols ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1043
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1044 if (nr > 0 && nc > 0)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1045 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1046 result.resize (nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1047 index.resize (nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1048
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1049 for (int j = 0; j < nc; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1050 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1051 int idx_i = 0;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1052
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1053 FixedPoint tmp_min = elem (idx_i, j);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1054
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1055 for (int i = 1; i < nr; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1056 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1057 FixedPoint tmp = elem (i, j);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1058 if (tmp < tmp_min)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1059 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1060 idx_i = i;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1061 tmp_min = tmp;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1062 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1063 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1064
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1065 result.elem (j) = tmp_min;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1066 index.elem (j) = idx_i;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1067 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1068 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1069
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1070 return result;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1071 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1072
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1073 FixedRowVector
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1074 FixedMatrix::column_max (void) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1075 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1076 Array<int> index;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1077 return column_max (index);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1078 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1079
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1080 FixedRowVector
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1081 FixedMatrix::column_max (Array<int>& index) const
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1082 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1083 FixedRowVector result;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1084
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1085 int nr = rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1086 int nc = cols ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1087
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1088 if (nr > 0 && nc > 0)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1089 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1090 result.resize (nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1091 index.resize (nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1092
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1093 for (int j = 0; j < nc; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1094 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1095 int idx_i = 0;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1096
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1097 FixedPoint tmp_max = elem (idx_i, j);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1098
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1099 for (int i = 1; i < nr; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1100 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1101 FixedPoint tmp = elem (i, j);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1102
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1103 if (tmp > tmp_max)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1104 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1105 idx_i = i;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1106 tmp_max = tmp;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1107 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1108 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1109
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1110 result.elem (j) = tmp_max;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1111 index.elem (j) = idx_i;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1112 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1113 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1114
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1115 return result;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1116 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1117
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1118 #define DO_FIXED_MAT_FUNC(FUNC) \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1119 FixedMatrix FUNC (const FixedMatrix& x) \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1120 { \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1121 FixedMatrix retval ( x.rows(), x.cols()); \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1122 for (int j = 0; j < x.cols(); j++) \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1123 for (int i = 0; i < x.rows(); i++) \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1124 retval(i,j) = FUNC ( x (i,j) ); \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1125 return retval; \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1126 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1127
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1128 DO_FIXED_MAT_FUNC(real);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1129 DO_FIXED_MAT_FUNC(imag);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1130 DO_FIXED_MAT_FUNC(conj);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1131 DO_FIXED_MAT_FUNC(abs);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1132 DO_FIXED_MAT_FUNC(cos);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1133 DO_FIXED_MAT_FUNC(cosh);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1134 DO_FIXED_MAT_FUNC(sin);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1135 DO_FIXED_MAT_FUNC(sinh);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1136 DO_FIXED_MAT_FUNC(tan);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1137 DO_FIXED_MAT_FUNC(tanh);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1138 DO_FIXED_MAT_FUNC(sqrt);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1139 DO_FIXED_MAT_FUNC(exp);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1140 DO_FIXED_MAT_FUNC(log);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1141 DO_FIXED_MAT_FUNC(log10);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1142 DO_FIXED_MAT_FUNC(round);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1143 DO_FIXED_MAT_FUNC(rint);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1144 DO_FIXED_MAT_FUNC(floor);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1145 DO_FIXED_MAT_FUNC(ceil);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1146
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1147 FixedMatrix elem_pow (const FixedMatrix &a, const FixedMatrix &b)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1148 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1149 FixedMatrix retval;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1150 int a_nr = a.rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1151 int a_nc = a.cols ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1152 int b_nr = b.rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1153 int b_nc = b.cols ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1154
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1155 if (a_nr == 1 && a_nc == 1)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1156 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1157 retval.resize(b_nr,b_nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1158 FixedPoint ad = a(0,0);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1159 for (int j = 0; j < b_nc; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1160 for (int i = 0; i < b_nr; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1161 retval(i,j) = pow(ad, b(i,j));
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1162 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1163 else if (b_nr == 1 && b_nc == 1)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1164 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1165 retval.resize(a_nr,a_nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1166 FixedPoint bd = b(0,0);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1167 for (int j = 0; j < a_nc; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1168 for (int i = 0; i < a_nr; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1169 retval(i,j) = pow(a(i,j), bd);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1170 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1171 else if ((a_nr == b_nr) && (a_nc == b_nc))
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1172 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1173 retval.resize(a_nr,a_nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1174 for (int j = 0; j < a_nc; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1175 for (int i = 0; i < a_nr; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1176 retval(i,j) = pow(a(i,j), b(i,j));
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1177 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1178 else
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1179 gripe_nonconformant ("operator .^", a_nr, a_nc, a_nr, a_nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1180
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1181 return retval;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1182 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1183
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1184 FixedMatrix elem_pow (const FixedMatrix &a, const FixedPoint &b)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1185 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1186 return elem_pow (a, FixedMatrix(1, 1, b));
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1187 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1188
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1189 FixedMatrix elem_pow (const FixedPoint &a, const FixedMatrix &b)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1190 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1191 return elem_pow (FixedMatrix(1, 1, a), b);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1192 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1193
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1194 FixedMatrix pow (const FixedMatrix& a, int b)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1195 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1196 FixedMatrix retval;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1197 int nr = a.rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1198 int nc = a.cols ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1199
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1200 if (nr == 0 \|\| nc == 0 \|\| nr != nc)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1201 (*current_liboctave_error_handler)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1202 ("for A^x, A must be square and x scalar");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1203 else if (b == 0)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1204 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1205 retval = a;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1206 for (int j = 0; j < nc; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1207 for (int i = 0; i < nr; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1208 retval(i,j) =
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1209 FixedPoint(a(i,j).getintsize(), a(i,j).getdecsize(), 1, 0);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1210 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1211 else
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1212 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1213 if (b < 0 )
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1214 (*current_liboctave_error_handler)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1215 ("can not treat matrix inversion");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1216
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1217 FixedMatrix atmp (a);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1218 retval = atmp;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1219 b--;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1220 while (b > 0)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1221 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1222 if (b & 1)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1223 retval = retval * atmp;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1224
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1225 b >>= 1;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1226
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1227 if (b > 0)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1228 atmp = atmp * atmp;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1229 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1230 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1231 return retval;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1232 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1233
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1234 FixedMatrix pow (const FixedMatrix &a, const double b)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1235 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1236 int bi = (int)b;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1237 if ((double)bi != b) {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1238 (*current_liboctave_error_handler)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1239 ("can only treat integer powers of a matrix");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1240 return FixedMatrix();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1241 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1242
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1243 return pow(a, bi);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1244 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1245
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1246 FixedMatrix pow (const FixedMatrix &a, const FixedPoint &b)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1247 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1248 return pow(a, b.fixedpoint());
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1249 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1250
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1251 FixedMatrix pow (const FixedMatrix &a, const FixedMatrix &b)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1252 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1253 if (a.rows() == 0 \|\| a.rows() == 0 \|\| a.rows() != a.rows() \|\|
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1254 b.rows() != 1 \|\| b.cols() != 1) {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1255 (*current_liboctave_error_handler)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1256 ("for A^x, A must be square and x scalar");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1257 return FixedMatrix();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1258 } else
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1259 return pow(a, b(0,0).fixedpoint());
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1260 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1261
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1262 FixedMatrix atan2 (const FixedMatrix &x, const FixedMatrix &y)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1263 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1264 FixedMatrix retval;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1265 if ((x.rows() == y.rows()) && (x.cols() == y.cols())) {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1266 retval.resize(x.rows(),x.cols());
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1267 for (int j = 0; j < x.cols(); j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1268 for (int i = 0; i < x.rows(); i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1269 retval(i,j) = atan2 ( x(i,j), y(i,j));
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1270 } else
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1271 (*current_liboctave_error_handler) ("matrix size mismatch");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1272 return retval;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1273 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1274
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1275 std::ostream&
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1276 operator << (std::ostream& os, const FixedMatrix& a)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1277 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1278 for (int i = 0; i < a.rows (); i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1279 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1280 for (int j = 0; j < a.cols (); j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1281 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1282 os << " " << a.elem(i,j);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1283 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1284 os << "\n";
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1285 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1286 return os;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1287 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1288
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1289 std::istream&
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1290 operator >> (std::istream& is, FixedMatrix& a)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1291 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1292 int nr = a.rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1293 int nc = a.cols ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1294
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1295 if (nr < 1 \|\| nc < 1)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1296 is.clear (std::ios::badbit);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1297 else
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1298 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1299 FixedPoint tmp;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1300 for (int i = 0; i < nr; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1301 for (int j = 0; j < nc; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1302 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1303 is >> tmp;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1304 if (is)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1305 a.elem (i, j) = tmp;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1306 else
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1307 goto done;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1308 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1309 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1310
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1311 done:
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1312
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1313 return is;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1314 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1315
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1316 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1317 operator * (const FixedMatrix& a, const FixedMatrix& b)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1318 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1319 FixedMatrix retval;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1320
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1321 int a_nr = a.rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1322 int a_nc = a.cols ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1323
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1324 int b_nr = b.rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1325 int b_nc = b.cols ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1326
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1327 if (a_nc != b_nr)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1328 gripe_nonconformant ("operator *", a_nr, a_nc, b_nr, b_nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1329 else
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1330 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1331 retval.resize (a_nr, b_nc, FixedPoint());
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1332 if (a_nr != 0 && a_nc != 0 && b_nc != 0)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1333 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1334 for (int j = 0; j < b_nr; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1335 for (int i = 0; i < b_nc; i++) {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1336 FixedPoint tmp = b.elem(j,i);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1337 for (int k = 0; k < a_nr; k++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1338 retval.elem (k,i) += a.elem(k,j) * tmp;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1339 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1340 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1341 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1342
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1343 return retval;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1344 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1345
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1346 // XXX FIXME XXX -- it would be nice to share code among the min/max
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1347 // functions below.
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1348
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1349 #define EMPTY_RETURN_CHECK(T) \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1350 if (nr == 0 \|\| nc == 0) \
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1351 return T (nr, nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1352
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1353 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1354 min (FixedPoint d, const FixedMatrix& m)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1355 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1356 int nr = m.rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1357 int nc = m.columns ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1358
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1359 EMPTY_RETURN_CHECK (FixedMatrix);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1360
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1361 FixedMatrix result (nr, nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1362
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1363 for (int j = 0; j < nc; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1364 for (int i = 0; i < nr; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1365 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1366 OCTAVE_QUIT;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1367 result (i, j) = m(i,j) < d ? m(i,j) : d;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1368 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1369
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1370 return result;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1371 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1372
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1373 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1374 min (const FixedMatrix& m, FixedPoint d)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1375 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1376 int nr = m.rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1377 int nc = m.columns ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1378
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1379 EMPTY_RETURN_CHECK (FixedMatrix);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1380
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1381 FixedMatrix result (nr, nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1382
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1383 for (int j = 0; j < nc; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1384 for (int i = 0; i < nr; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1385 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1386 OCTAVE_QUIT;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1387 result (i, j) = m(i,j) < d ? m(i,j) : d;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1388 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1389
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1390 return result;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1391 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1392
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1393 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1394 min (const FixedMatrix& a, const FixedMatrix& b)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1395 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1396 int nr = a.rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1397 int nc = a.columns ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1398
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1399 if (nr != b.rows () \|\| nc != b.columns ())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1400 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1401 (*current_liboctave_error_handler)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1402 ("two-arg min expecting args of same size");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1403 return FixedMatrix ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1404 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1405
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1406 EMPTY_RETURN_CHECK (FixedMatrix);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1407
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1408 FixedMatrix result (nr, nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1409
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1410 for (int j = 0; j < nc; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1411 for (int i = 0; i < nr; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1412 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1413 OCTAVE_QUIT;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1414 result (i, j) = a(i,j) < b(i,j) ? a(i,j) : b(i,j);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1415 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1416
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1417 return result;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1418 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1419
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1420 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1421 max (FixedPoint d, const FixedMatrix& m)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1422 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1423 int nr = m.rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1424 int nc = m.columns ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1425
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1426 EMPTY_RETURN_CHECK (FixedMatrix);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1427
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1428 FixedMatrix result (nr, nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1429
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1430 for (int j = 0; j < nc; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1431 for (int i = 0; i < nr; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1432 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1433 OCTAVE_QUIT;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1434 result (i, j) = m(i,j) > d ? m(i,j) : d;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1435 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1436
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1437 return result;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1438 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1439
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1440 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1441 max (const FixedMatrix& m, FixedPoint d)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1442 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1443 int nr = m.rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1444 int nc = m.columns ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1445
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1446 EMPTY_RETURN_CHECK (FixedMatrix);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1447
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1448 FixedMatrix result (nr, nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1449
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1450 for (int j = 0; j < nc; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1451 for (int i = 0; i < nr; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1452 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1453 OCTAVE_QUIT;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1454 result (i, j) = m(i,j) > d ? m(i,j) : d;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1455 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1456
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1457 return result;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1458 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1459
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1460 FixedMatrix
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1461 max (const FixedMatrix& a, const FixedMatrix& b)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1462 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1463 int nr = a.rows ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1464 int nc = a.columns ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1465
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1466 if (nr != b.rows () \|\| nc != b.columns ())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1467 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1468 (*current_liboctave_error_handler)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1469 ("two-arg max expecting args of same size");
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1470 return FixedMatrix ();
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1471 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1472
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1473 EMPTY_RETURN_CHECK (FixedMatrix);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1474
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1475 FixedMatrix result (nr, nc);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1476
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1477 for (int j = 0; j < nc; j++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1478 for (int i = 0; i < nr; i++)
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1479 {
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1480 OCTAVE_QUIT;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1481 result (i, j) = a(i,j) > b(i,j) ? a(i,j) : b(i,j);
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1482 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1483
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1484 return result;
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1485 }
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1486
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1487 MS_CMP_OPS(FixedMatrix, , FixedPoint, )
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1488 MS_BOOL_OPS(FixedMatrix, FixedPoint, FixedPoint())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1489
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1490 SM_CMP_OPS(FixedPoint, , FixedMatrix, )
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1491 SM_BOOL_OPS(FixedPoint, FixedMatrix, FixedPoint())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1492
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1493 MM_CMP_OPS(FixedMatrix, , FixedMatrix, )
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1494 MM_BOOL_OPS(FixedMatrix, FixedMatrix, FixedPoint())
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1495
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1496 /*
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1497 ;;; Local Variables: ***
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1498 ;;; mode: C++ ***
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1499 ;;; End: ***
d1930f8af7ca Changed the directory structure of fixed to match the package system hauberg parents: diff changeset	1500 */

Mercurial > forge

annotate main/fixed/src/fixedMatrix.cc @ 2914:e2142eefb50d octave-forge