octave: liboctave/array/dRowVector.cc annotate

annotate liboctave/array/dRowVector.cc @ 33636:3ec6fcce7715 default tip @

gui: Avoid using HAVE_QSCINTILLA in more header files. * gui-settings.h, settings-dialog.h: Don't include QScintilla header. Forward-declare QSciLexer class instead if necessary. Declare all member functions unconditionally. * gui-settings.cc (gui_settings::get_valid_lexer_styles, gui_settings::read_lexer_settings), settings-dialog.cc (settings_dialog::update_lexer, settings_dialog::get_lexer_settings, settings_dialog::write_lexer_settings): Define functions unconditionally. * gui-preferences-ed.h: Don't include QScintilla header. Remove definition of local variable os_eol_mode from header. * gui-preferences-ed.cc (os_eol_mode): Move definition of local variable here.

author	Markus Mützel <markus.muetzel@gmx.de>
date	Tue, 28 May 2024 14:54:58 +0200
parents	f53ac65ffba6
children

rev	line source
27923 bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	1 ////////////////////////////////////////////////////////////////////////
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	2 //
32632 2e484f9f1f18 maint: update Octave Project Developers copyright for the new year John W. Eaton <jwe@octave.org> parents: 31927 diff changeset	3 // Copyright (C) 1994-2024 The Octave Project Developers
27923 bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	4 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	5 // See the file COPYRIGHT.md in the top-level directory of this
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	6 // distribution or <https://octave.org/copyright/>.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	7 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	8 // This file is part of Octave.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	9 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	10 // Octave is free software: you can redistribute it and/or modify it
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	11 // under the terms of the GNU General Public License as published by
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	12 // the Free Software Foundation, either version 3 of the License, or
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	13 // (at your option) any later version.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	14 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	15 // Octave is distributed in the hope that it will be useful, but
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	16 // WITHOUT ANY WARRANTY; without even the implied warranty of
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	17 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	18 // GNU General Public License for more details.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	19 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	20 // You should have received a copy of the GNU General Public License
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	21 // along with Octave; see the file COPYING. If not, see
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	22 // <https://www.gnu.org/licenses/>.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	23 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	24 ////////////////////////////////////////////////////////////////////////
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	25
21724 aba2e6293dd8 use "#if ..." consistently instead of "#ifdef" and "#ifndef" John W. Eaton <jwe@octave.org> parents: 21301 diff changeset	26 #if defined (HAVE_CONFIG_H)
21301 40de9f8f23a6 Use '#include "config.h"' rather than <config.h>. Rik <rik@octave.org> parents: 21202 diff changeset	27 # include "config.h"
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	28 #endif
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	29
25438 cb1606f78f6b prefer <istream>, <ostream>, or <iosfwd> to <iostream> where possible John W. Eaton <jwe@octave.org> parents: 25237 diff changeset	30 #include <istream>
cb1606f78f6b prefer <istream>, <ostream>, or <iosfwd> to <iostream> where possible John W. Eaton <jwe@octave.org> parents: 25237 diff changeset	31 #include <ostream>
27435 a8a5d2e8807f Produce symmetric sequences for linspace (bug #56659). Rik <rik@octave.org> parents: 27180 diff changeset	32 #include <type_traits>
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	33
4669 334a27c8f453 [project @ 2003-11-26 07:02:42 by jwe] jwe parents: 4587 diff changeset	34 #include "Array-util.h"
22322 93b3cdd36854 move most f77 function decls to separate header files John W. Eaton <jwe@octave.org> parents: 22133 diff changeset	35 #include "lo-blas-proto.h"
1368 33bb7975f866 [project @ 1995-09-06 00:12:48 by jwe] jwe parents: 1360 diff changeset	36 #include "lo-error.h"
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	37 #include "mx-base.h"
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	38 #include "mx-inlines.cc"
1650 23aa282707e8 [project @ 1995-12-20 06:53:12 by jwe] jwe parents: 1368 diff changeset	39 #include "oct-cmplx.h"
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	40
1360 7eb93d12654c [project @ 1995-09-05 21:51:54 by jwe] jwe parents: 1315 diff changeset	41 // Row Vector class.
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	42
2386 4fc9fd1424a9 [project @ 1996-10-12 18:31:34 by jwe] jwe parents: 1993 diff changeset	43 bool
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	44 RowVector::operator == (const RowVector& a) const
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	45 {
20232 a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	46 octave_idx_type len = numel ();
a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	47 if (len != a.numel ())
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	48 return 0;
9550 3d6a9aea2aea refactor binary & bool ops in liboctave Jaroslav Hajek <highegg@gmail.com> parents: 8999 diff changeset	49 return mx_inline_equal (len, data (), a.data ());
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	50 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	51
2386 4fc9fd1424a9 [project @ 1996-10-12 18:31:34 by jwe] jwe parents: 1993 diff changeset	52 bool
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	53 RowVector::operator != (const RowVector& a) const
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	54 {
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	55 return !(*this == a);
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	56 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	57
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	58 RowVector&
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	59 RowVector::insert (const RowVector& a, octave_idx_type c)
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	60 {
20232 a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	61 octave_idx_type a_len = a.numel ();
4316 236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	62
20232 a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	63 if (c < 0 \|\| c + a_len > numel ())
21136 7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21120 diff changeset	64 (*current_liboctave_error_handler) ("range error for insert");
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	65
4316 236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	66 if (a_len > 0)
236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	67 {
236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	68 make_unique ();
236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	69
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	70 for (octave_idx_type i = 0; i < a_len; i++)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	71 xelem (c+i) = a.elem (i);
4316 236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	72 }
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	73
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	74 return *this;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	75 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	76
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	77 RowVector&
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	78 RowVector::fill (double val)
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	79 {
20232 a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	80 octave_idx_type len = numel ();
4316 236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	81
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	82 if (len > 0)
4316 236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	83 {
236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	84 make_unique ();
236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	85
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	86 for (octave_idx_type i = 0; i < len; i++)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	87 xelem (i) = val;
4316 236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	88 }
236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	89
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	90 return *this;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	91 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	92
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	93 RowVector&
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	94 RowVector::fill (double val, octave_idx_type c1, octave_idx_type c2)
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	95 {
20232 a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	96 octave_idx_type len = numel ();
4316 236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	97
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	98 if (c1 < 0 \|\| c2 < 0 \|\| c1 >= len \|\| c2 >= len)
21136 7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21120 diff changeset	99 (*current_liboctave_error_handler) ("range error for fill");
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	100
17663 7975d75f933c Use std::swap in liboctave instead of temporary variable. Rik <rik@octave.org> parents: 15271 diff changeset	101 if (c1 > c2) { std::swap (c1, c2); }
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	102
4316 236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	103 if (c2 >= c1)
236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	104 {
236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	105 make_unique ();
236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	106
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	107 for (octave_idx_type i = c1; i <= c2; i++)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	108 xelem (i) = val;
4316 236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	109 }
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	110
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	111 return *this;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	112 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	113
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	114 RowVector
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	115 RowVector::append (const RowVector& a) const
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	116 {
20232 a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	117 octave_idx_type len = numel ();
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	118 octave_idx_type nc_insert = len;
20232 a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	119 RowVector retval (len + a.numel ());
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	120 retval.insert (*this, 0);
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	121 retval.insert (a, nc_insert);
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	122 return retval;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	123 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	124
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	125 ColumnVector
31771 21f9b34eb893 maint: Eliminate "(void)" in C++ function prototypes/declarations. Rik <rik@octave.org> parents: 31706 diff changeset	126 RowVector::transpose () const
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	127 {
14846 460a3c6d8bf1 maint: Use Octave coding convention for cuddled parenthis in function calls with empty argument lists. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	128 return MArray<double>::transpose ();
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	129 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	130
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	131 RowVector
1205 8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	132 real (const ComplexRowVector& a)
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	133 {
10363 a0728e81ed25 improve diag matrix interface & implementation Jaroslav Hajek <highegg@gmail.com> parents: 10314 diff changeset	134 return do_mx_unary_op<double, Complex> (a, mx_inline_real);
1205 8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	135 }
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	136
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	137 RowVector
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	138 imag (const ComplexRowVector& a)
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	139 {
10363 a0728e81ed25 improve diag matrix interface & implementation Jaroslav Hajek <highegg@gmail.com> parents: 10314 diff changeset	140 return do_mx_unary_op<double, Complex> (a, mx_inline_imag);
1205 8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	141 }
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	142
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	143 RowVector
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	144 RowVector::extract (octave_idx_type c1, octave_idx_type c2) const
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	145 {
17663 7975d75f933c Use std::swap in liboctave instead of temporary variable. Rik <rik@octave.org> parents: 15271 diff changeset	146 if (c1 > c2) { std::swap (c1, c2); }
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	147
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	148 octave_idx_type new_c = c2 - c1 + 1;
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	149
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	150 RowVector result (new_c);
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	151
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	152 for (octave_idx_type i = 0; i < new_c; i++)
4316 236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	153 result.xelem (i) = elem (c1+i);
236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	154
236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	155 return result;
236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	156 }
236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	157
236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	158 RowVector
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	159 RowVector::extract_n (octave_idx_type r1, octave_idx_type n) const
4316 236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	160 {
236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	161 RowVector result (n);
236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	162
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	163 for (octave_idx_type i = 0; i < n; i++)
4316 236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	164 result.xelem (i) = elem (r1+i);
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	165
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	166 return result;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	167 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	168
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	169 // row vector by matrix -> row vector
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	170
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	171 RowVector
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	172 operator * (const RowVector& v, const Matrix& a)
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	173 {
1947 5ab3c25a3cf9 [project @ 1996-02-14 01:50:42 by jwe] jwe parents: 1882 diff changeset	174 RowVector retval;
5ab3c25a3cf9 [project @ 1996-02-14 01:50:42 by jwe] jwe parents: 1882 diff changeset	175
22988 cd33c785e80e put to_f77_int inside octave namespace John W. Eaton <jwe@octave.org> parents: 22981 diff changeset	176 F77_INT len = octave::to_f77_int (v.numel ());
1947 5ab3c25a3cf9 [project @ 1996-02-14 01:50:42 by jwe] jwe parents: 1882 diff changeset	177
22988 cd33c785e80e put to_f77_int inside octave namespace John W. Eaton <jwe@octave.org> parents: 22981 diff changeset	178 F77_INT a_nr = octave::to_f77_int (a.rows ());
cd33c785e80e put to_f77_int inside octave namespace John W. Eaton <jwe@octave.org> parents: 22981 diff changeset	179 F77_INT a_nc = octave::to_f77_int (a.cols ());
2386 4fc9fd1424a9 [project @ 1996-10-12 18:31:34 by jwe] jwe parents: 1993 diff changeset	180
4fc9fd1424a9 [project @ 1996-10-12 18:31:34 by jwe] jwe parents: 1993 diff changeset	181 if (a_nr != len)
22327 d0562b3159c7 move more classes inside octave namespace John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	182 octave::err_nonconformant ("operator *", 1, len, a_nr, a_nc);
21120 499b851fbfae Replace pattern if/err_XXX/else/code with if/err_XXX/ code. Rik <rik@octave.org> parents: 21100 diff changeset	183
499b851fbfae Replace pattern if/err_XXX/else/code with if/err_XXX/ code. Rik <rik@octave.org> parents: 21100 diff changeset	184 if (len == 0)
499b851fbfae Replace pattern if/err_XXX/else/code with if/err_XXX/ code. Rik <rik@octave.org> parents: 21100 diff changeset	185 retval.resize (a_nc, 0.0);
1947 5ab3c25a3cf9 [project @ 1996-02-14 01:50:42 by jwe] jwe parents: 1882 diff changeset	186 else
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	187 {
21120 499b851fbfae Replace pattern if/err_XXX/else/code with if/err_XXX/ code. Rik <rik@octave.org> parents: 21100 diff changeset	188 // Transpose A to form A'x == (x'A)'
1947 5ab3c25a3cf9 [project @ 1996-02-14 01:50:42 by jwe] jwe parents: 1882 diff changeset	189
22981 4a53c121cdec fix a few mistakes in F77_INT changes John W. Eaton <jwe@octave.org> parents: 22847 diff changeset	190 F77_INT ld = a_nr;
1947 5ab3c25a3cf9 [project @ 1996-02-14 01:50:42 by jwe] jwe parents: 1882 diff changeset	191
21120 499b851fbfae Replace pattern if/err_XXX/else/code with if/err_XXX/ code. Rik <rik@octave.org> parents: 21100 diff changeset	192 retval.resize (a_nc);
32660 f53ac65ffba6 maint: New method rwdata() as clearer alternative to fortran_vec(). Rik <rik@octave.org> parents: 32632 diff changeset	193 double *y = retval.rwdata ();
1947 5ab3c25a3cf9 [project @ 1996-02-14 01:50:42 by jwe] jwe parents: 1882 diff changeset	194
21120 499b851fbfae Replace pattern if/err_XXX/else/code with if/err_XXX/ code. Rik <rik@octave.org> parents: 21100 diff changeset	195 F77_XFCN (dgemv, DGEMV, (F77_CONST_CHAR_ARG2 ("T", 1),
499b851fbfae Replace pattern if/err_XXX/else/code with if/err_XXX/ code. Rik <rik@octave.org> parents: 21100 diff changeset	196 a_nr, a_nc, 1.0, a.data (),
499b851fbfae Replace pattern if/err_XXX/else/code with if/err_XXX/ code. Rik <rik@octave.org> parents: 21100 diff changeset	197 ld, v.data (), 1, 0.0, y, 1
499b851fbfae Replace pattern if/err_XXX/else/code with if/err_XXX/ code. Rik <rik@octave.org> parents: 21100 diff changeset	198 F77_CHAR_ARG_LEN (1)));
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	199 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	200
1947 5ab3c25a3cf9 [project @ 1996-02-14 01:50:42 by jwe] jwe parents: 1882 diff changeset	201 return retval;
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	202 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	203
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	204 // other operations
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	205
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	206 double
31771 21f9b34eb893 maint: Eliminate "(void)" in C++ function prototypes/declarations. Rik <rik@octave.org> parents: 31706 diff changeset	207 RowVector::min () const
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	208 {
20232 a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	209 octave_idx_type len = numel ();
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	210 if (len == 0)
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	211 return 0;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	212
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	213 double res = elem (0);
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	214
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	215 for (octave_idx_type i = 1; i < len; i++)
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	216 if (elem (i) < res)
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	217 res = elem (i);
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	218
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	219 return res;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	220 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	221
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	222 double
31771 21f9b34eb893 maint: Eliminate "(void)" in C++ function prototypes/declarations. Rik <rik@octave.org> parents: 31706 diff changeset	223 RowVector::max () const
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	224 {
20232 a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	225 octave_idx_type len = numel ();
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	226 if (len == 0)
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	227 return 0;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	228
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	229 double res = elem (0);
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	230
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	231 for (octave_idx_type i = 1; i < len; i++)
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	232 if (elem (i) > res)
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	233 res = elem (i);
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	234
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	235 return res;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	236 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	237
3504 5eef8a2294bd [project @ 2000-02-01 10:06:51 by jwe] jwe parents: 3503 diff changeset	238 std::ostream&
5eef8a2294bd [project @ 2000-02-01 10:06:51 by jwe] jwe parents: 3503 diff changeset	239 operator << (std::ostream& os, const RowVector& a)
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	240 {
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	241 // int field_width = os.precision () + 7;
1360 7eb93d12654c [project @ 1995-09-05 21:51:54 by jwe] jwe parents: 1315 diff changeset	242
20232 a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	243 for (octave_idx_type i = 0; i < a.numel (); i++)
23807 336f89b6208b Use character literals 'c' rather than string literals "c" when possible. Rik <rik@octave.org> parents: 23461 diff changeset	244 os << ' ' /* setw (field_width) */ << a.elem (i);
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	245 return os;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	246 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	247
3504 5eef8a2294bd [project @ 2000-02-01 10:06:51 by jwe] jwe parents: 3503 diff changeset	248 std::istream&
5eef8a2294bd [project @ 2000-02-01 10:06:51 by jwe] jwe parents: 3503 diff changeset	249 operator >> (std::istream& is, RowVector& a)
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	250 {
20232 a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	251 octave_idx_type len = a.numel ();
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	252
8999 dc07bc4157b8 allow empty matrices in stream input operators Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	253 if (len > 0)
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	254 {
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	255 double tmp;
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	256 for (octave_idx_type i = 0; i < len; i++)
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	257 {
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	258 is >> tmp;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	259 if (is)
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	260 a.elem (i) = tmp;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	261 else
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	262 break;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	263 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	264 }
532 2ca256b77602 [project @ 1994-07-20 19:56:55 by jwe] jwe parents: 458 diff changeset	265 return is;
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	266 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	267
1205 8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	268 // other operations
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	269
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	270 RowVector
27435 a8a5d2e8807f Produce symmetric sequences for linspace (bug #56659). Rik <rik@octave.org> parents: 27180 diff changeset	271 linspace (double x1, double x2, octave_idx_type n_in)
1205 8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	272 {
27666 3e8faed1b7d8 Remove wrapper template class NoAlias<T> (bug #56752) Carlo de Falco <carlo.defalco@polimi.it> parents: 27435 diff changeset	273 RowVector retval;
9653 e087d7c77ff9 improve linspace in liboctave Jaroslav Hajek <highegg@gmail.com> parents: 9550 diff changeset	274
27435 a8a5d2e8807f Produce symmetric sequences for linspace (bug #56659). Rik <rik@octave.org> parents: 27180 diff changeset	275 if (n_in < 1)
20466 16b9ec39ff46 Return empty matrix when linspace called with 0 points (bug #45820) Rik <rik@octave.org> parents: 20232 diff changeset	276 return retval;
27435 a8a5d2e8807f Produce symmetric sequences for linspace (bug #56659). Rik <rik@octave.org> parents: 27180 diff changeset	277 else if (n_in == 1)
a8a5d2e8807f Produce symmetric sequences for linspace (bug #56659). Rik <rik@octave.org> parents: 27180 diff changeset	278 {
a8a5d2e8807f Produce symmetric sequences for linspace (bug #56659). Rik <rik@octave.org> parents: 27180 diff changeset	279 retval.resize (1, x2);
a8a5d2e8807f Produce symmetric sequences for linspace (bug #56659). Rik <rik@octave.org> parents: 27180 diff changeset	280 return retval;
a8a5d2e8807f Produce symmetric sequences for linspace (bug #56659). Rik <rik@octave.org> parents: 27180 diff changeset	281 }
31927 fef004fa8575 Improve Matlab compatibility of linspace for Inf endpoints. Rik <rik@octave.org> parents: 31771 diff changeset	282 else if (x1 == x2)
fef004fa8575 Improve Matlab compatibility of linspace for Inf endpoints. Rik <rik@octave.org> parents: 31771 diff changeset	283 {
fef004fa8575 Improve Matlab compatibility of linspace for Inf endpoints. Rik <rik@octave.org> parents: 31771 diff changeset	284 retval.resize (n_in, x2);
fef004fa8575 Improve Matlab compatibility of linspace for Inf endpoints. Rik <rik@octave.org> parents: 31771 diff changeset	285 return retval;
fef004fa8575 Improve Matlab compatibility of linspace for Inf endpoints. Rik <rik@octave.org> parents: 31771 diff changeset	286 }
20466 16b9ec39ff46 Return empty matrix when linspace called with 0 points (bug #45820) Rik <rik@octave.org> parents: 20232 diff changeset	287
27435 a8a5d2e8807f Produce symmetric sequences for linspace (bug #56659). Rik <rik@octave.org> parents: 27180 diff changeset	288 // Use unsigned type (guaranteed n_in > 1 at this point) so that divisions
a8a5d2e8807f Produce symmetric sequences for linspace (bug #56659). Rik <rik@octave.org> parents: 27180 diff changeset	289 // by 2 can be replaced by compiler with shift right instructions.
a8a5d2e8807f Produce symmetric sequences for linspace (bug #56659). Rik <rik@octave.org> parents: 27180 diff changeset	290 typedef std::make_unsigned<octave_idx_type>::type unsigned_octave_idx_type;
a8a5d2e8807f Produce symmetric sequences for linspace (bug #56659). Rik <rik@octave.org> parents: 27180 diff changeset	291
a8a5d2e8807f Produce symmetric sequences for linspace (bug #56659). Rik <rik@octave.org> parents: 27180 diff changeset	292 unsigned_octave_idx_type n = n_in;
a8a5d2e8807f Produce symmetric sequences for linspace (bug #56659). Rik <rik@octave.org> parents: 27180 diff changeset	293
a8a5d2e8807f Produce symmetric sequences for linspace (bug #56659). Rik <rik@octave.org> parents: 27180 diff changeset	294 // Set endpoints, rather than calculate, for maximum accuracy.
a8a5d2e8807f Produce symmetric sequences for linspace (bug #56659). Rik <rik@octave.org> parents: 27180 diff changeset	295 retval.clear (n);
27666 3e8faed1b7d8 Remove wrapper template class NoAlias<T> (bug #56752) Carlo de Falco <carlo.defalco@polimi.it> parents: 27435 diff changeset	296 retval.xelem (0) = x1;
3e8faed1b7d8 Remove wrapper template class NoAlias<T> (bug #56752) Carlo de Falco <carlo.defalco@polimi.it> parents: 27435 diff changeset	297 retval.xelem (n-1) = x2;
1205 8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	298
27435 a8a5d2e8807f Produce symmetric sequences for linspace (bug #56659). Rik <rik@octave.org> parents: 27180 diff changeset	299 // Construct linspace symmetrically from both ends.
31927 fef004fa8575 Improve Matlab compatibility of linspace for Inf endpoints. Rik <rik@octave.org> parents: 31771 diff changeset	300 bool isnan_delta = false;
27435 a8a5d2e8807f Produce symmetric sequences for linspace (bug #56659). Rik <rik@octave.org> parents: 27180 diff changeset	301 double delta = (x2 - x1) / (n - 1);
31927 fef004fa8575 Improve Matlab compatibility of linspace for Inf endpoints. Rik <rik@octave.org> parents: 31771 diff changeset	302 if (octave::math::isinf (delta))
fef004fa8575 Improve Matlab compatibility of linspace for Inf endpoints. Rik <rik@octave.org> parents: 31771 diff changeset	303 {
fef004fa8575 Improve Matlab compatibility of linspace for Inf endpoints. Rik <rik@octave.org> parents: 31771 diff changeset	304 delta = octave::numeric_limits<double>::NaN ();
fef004fa8575 Improve Matlab compatibility of linspace for Inf endpoints. Rik <rik@octave.org> parents: 31771 diff changeset	305 isnan_delta = true;
fef004fa8575 Improve Matlab compatibility of linspace for Inf endpoints. Rik <rik@octave.org> parents: 31771 diff changeset	306 }
fef004fa8575 Improve Matlab compatibility of linspace for Inf endpoints. Rik <rik@octave.org> parents: 31771 diff changeset	307
27435 a8a5d2e8807f Produce symmetric sequences for linspace (bug #56659). Rik <rik@octave.org> parents: 27180 diff changeset	308 unsigned_octave_idx_type n2 = n/2;
a8a5d2e8807f Produce symmetric sequences for linspace (bug #56659). Rik <rik@octave.org> parents: 27180 diff changeset	309 for (unsigned_octave_idx_type i = 1; i < n2; i++)
a8a5d2e8807f Produce symmetric sequences for linspace (bug #56659). Rik <rik@octave.org> parents: 27180 diff changeset	310 {
27666 3e8faed1b7d8 Remove wrapper template class NoAlias<T> (bug #56752) Carlo de Falco <carlo.defalco@polimi.it> parents: 27435 diff changeset	311 retval.xelem (i) = x1 + i*delta;
3e8faed1b7d8 Remove wrapper template class NoAlias<T> (bug #56752) Carlo de Falco <carlo.defalco@polimi.it> parents: 27435 diff changeset	312 retval.xelem (n-1-i) = x2 - i*delta;
27435 a8a5d2e8807f Produce symmetric sequences for linspace (bug #56659). Rik <rik@octave.org> parents: 27180 diff changeset	313 }
a8a5d2e8807f Produce symmetric sequences for linspace (bug #56659). Rik <rik@octave.org> parents: 27180 diff changeset	314 if (n % 2 == 1) // Middle element if number of elements is odd.
31927 fef004fa8575 Improve Matlab compatibility of linspace for Inf endpoints. Rik <rik@octave.org> parents: 31771 diff changeset	315 {
fef004fa8575 Improve Matlab compatibility of linspace for Inf endpoints. Rik <rik@octave.org> parents: 31771 diff changeset	316 if (x1 == -x2)
fef004fa8575 Improve Matlab compatibility of linspace for Inf endpoints. Rik <rik@octave.org> parents: 31771 diff changeset	317 retval.xelem (n2) = 0;
fef004fa8575 Improve Matlab compatibility of linspace for Inf endpoints. Rik <rik@octave.org> parents: 31771 diff changeset	318 else if (isnan_delta)
fef004fa8575 Improve Matlab compatibility of linspace for Inf endpoints. Rik <rik@octave.org> parents: 31771 diff changeset	319 retval.xelem (n2) = octave::numeric_limits<double>::NaN ();
fef004fa8575 Improve Matlab compatibility of linspace for Inf endpoints. Rik <rik@octave.org> parents: 31771 diff changeset	320 else
fef004fa8575 Improve Matlab compatibility of linspace for Inf endpoints. Rik <rik@octave.org> parents: 31771 diff changeset	321 retval.xelem (n2) = (x1 + x2) / 2;
fef004fa8575 Improve Matlab compatibility of linspace for Inf endpoints. Rik <rik@octave.org> parents: 31771 diff changeset	322 }
1205 8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	323
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	324 return retval;
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	325 }
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	326
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	327 // row vector by column vector -> scalar
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	328
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	329 double
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	330 operator * (const RowVector& v, const ColumnVector& a)
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	331 {
1947 5ab3c25a3cf9 [project @ 1996-02-14 01:50:42 by jwe] jwe parents: 1882 diff changeset	332 double retval = 0.0;
5ab3c25a3cf9 [project @ 1996-02-14 01:50:42 by jwe] jwe parents: 1882 diff changeset	333
22988 cd33c785e80e put to_f77_int inside octave namespace John W. Eaton <jwe@octave.org> parents: 22981 diff changeset	334 F77_INT len = octave::to_f77_int (v.numel ());
1947 5ab3c25a3cf9 [project @ 1996-02-14 01:50:42 by jwe] jwe parents: 1882 diff changeset	335
22988 cd33c785e80e put to_f77_int inside octave namespace John W. Eaton <jwe@octave.org> parents: 22981 diff changeset	336 F77_INT a_len = octave::to_f77_int (a.numel ());
2386 4fc9fd1424a9 [project @ 1996-10-12 18:31:34 by jwe] jwe parents: 1993 diff changeset	337
4fc9fd1424a9 [project @ 1996-10-12 18:31:34 by jwe] jwe parents: 1993 diff changeset	338 if (len != a_len)
22327 d0562b3159c7 move more classes inside octave namespace John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	339 octave::err_nonconformant ("operator *", len, a_len);
21120 499b851fbfae Replace pattern if/err_XXX/else/code with if/err_XXX/ code. Rik <rik@octave.org> parents: 21100 diff changeset	340
499b851fbfae Replace pattern if/err_XXX/else/code with if/err_XXX/ code. Rik <rik@octave.org> parents: 21100 diff changeset	341 if (len != 0)
5983 ae09df27153f [project @ 2006-09-12 02:15:47 by jwe] jwe parents: 5307 diff changeset	342 F77_FUNC (xddot, XDDOT) (len, v.data (), 1, a.data (), 1, retval);
1205 8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	343
1947 5ab3c25a3cf9 [project @ 1996-02-14 01:50:42 by jwe] jwe parents: 1882 diff changeset	344 return retval;
1205 8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	345 }
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	346
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	347 Complex
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	348 operator * (const RowVector& v, const ComplexColumnVector& a)
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	349 {
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	350 ComplexRowVector tmp (v);
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	351 return tmp * a;
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	352 }

Mercurial > octave

annotate liboctave/array/dRowVector.cc @ 33636:3ec6fcce7715 default tip @