octave: liboctave/numeric/lo-mappers.h annotate

annotate liboctave/numeric/lo-mappers.h @ 15409:c898dbe14e1d

build: Update imported module list from gnulib. * bootstrap.conf: Update imported module list from gnulib. * libgnu/Makefile.am: Auto-updated by gnulib-tool. * liboctave/link-deps.mk: Add FLOORF_LIBM to list of GNULIB_DEPS. * lo-mappers.h (xfloor): Change float version to forward to xfloor in lo-mappers.cc instead of being an inline function. * lo-mappers.c (xfloor): New function with float signature which calls gnulib::floorf. (xround) : Fix typo and call gnulib::roundf for float version of function. * file-ops.cc: Include "canonicalize.h" for access to gnulib::canonicalize_file_name.

author	Rik <rik@octave.org>
date	Tue, 18 Sep 2012 10:15:05 -0700
parents	648dabbb4c6b
children	ecf5be238b4a

rev	line source
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	1 /*
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	2
14138 72c96de7a403 maint: update copyright notices for 2012 John W. Eaton <jwe@octave.org> parents: 13756 diff changeset	3 Copyright (C) 1996-2012 John W. Eaton
10521 4d1fc073fbb7 add some missing copyright stmts Jaroslav Hajek <highegg@gmail.com> parents: 10436 diff changeset	4 Copyright (C) 2010 VZLU Prague
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	5
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	6 This file is part of Octave.
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	7
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	8 Octave is free software; you can redistribute it and/or modify it
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	9 under the terms of the GNU General Public License as published by the
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6108 diff changeset	10 Free Software Foundation; either version 3 of the License, or (at your
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6108 diff changeset	11 option) any later version.
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	12
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	13 Octave is distributed in the hope that it will be useful, but WITHOUT
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	14 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	15 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	16 for more details.
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	17
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	18 You should have received a copy of the GNU General Public License
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6108 diff changeset	19 along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6108 diff changeset	20 <http://www.gnu.org/licenses/>.
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	21
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	22 */
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	23
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	24 #if !defined (octave_liboctave_mappers_h)
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	25 #define octave_liboctave_mappers_h 1
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	26
11211 2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	27 #include <limits>
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	28
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	29 #include "oct-cmplx.h"
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	30 #include "lo-math.h"
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	31
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11525 diff changeset	32 // Double Precision
11237 110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11213 diff changeset	33 extern OCTAVE_API double xtrunc (double x);
13737 30414ff19d5e Use copysignf module from gnulib. Michael Goffioul <michael.goffioul@gmail.com> parents: 13720 diff changeset	34 extern OCTAVE_API double xcopysign (double x, double y);
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	35 inline double xceil (double x) { return ceil (x); }
11525 ab231f944252 avoid exposing gnulib:: in header files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	36 extern OCTAVE_API double xfloor (double x);
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	37 inline double arg (double x) { return atan2 (0.0, x); }
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	38 inline double conj (double x) { return x; }
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	39 inline double fix (double x) { return xtrunc (x); }
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	40 inline double imag (double) { return 0.0; }
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	41 inline double real (double x) { return x; }
11237 110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11213 diff changeset	42 extern OCTAVE_API double xround (double x);
7636 99c410f7f0b0 implement mapper function for banker's rounding Jaroslav Hajek <highegg@gmail.com> parents: 7017 diff changeset	43 extern OCTAVE_API double xroundb (double x);
6108 143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 5904 diff changeset	44 extern OCTAVE_API double signum (double x);
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11525 diff changeset	45 extern OCTAVE_API double xlog2 (double x);
12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11525 diff changeset	46 extern OCTAVE_API Complex xlog2 (const Complex& x);
7740 39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7636 diff changeset	47 extern OCTAVE_API double xlog2 (double x, int& exp);
39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7636 diff changeset	48 extern OCTAVE_API Complex xlog2 (const Complex& x, int& exp);
6108 143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 5904 diff changeset	49 extern OCTAVE_API double xexp2 (double x);
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	50
7922 935be827eaf8 error for NaN values in & and \| expressions John W. Eaton <jwe@octave.org> parents: 7789 diff changeset	51 // These are used by the BOOL_OP macros in mx-op-defs.h.
935be827eaf8 error for NaN values in & and \| expressions John W. Eaton <jwe@octave.org> parents: 7789 diff changeset	52 inline bool xisnan (bool) { return false; }
935be827eaf8 error for NaN values in & and \| expressions John W. Eaton <jwe@octave.org> parents: 7789 diff changeset	53 inline bool xisnan (char) { return false; }
935be827eaf8 error for NaN values in & and \| expressions John W. Eaton <jwe@octave.org> parents: 7789 diff changeset	54
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	55 #if defined (HAVE_CMATH_ISNAN)
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	56 inline bool xisnan (double x)
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	57 { return std::isnan (x); }
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	58 #else
6108 143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 5904 diff changeset	59 extern OCTAVE_API bool xisnan (double x);
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	60 #endif
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	61 #if defined (HAVE_CMATH_ISFINITE)
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	62 inline bool xfinite (double x)
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	63 { return std::isfinite (x); }
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	64 #else
6108 143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 5904 diff changeset	65 extern OCTAVE_API bool xfinite (double x);
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	66 #endif
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	67 #if defined (HAVE_CMATH_ISINF)
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	68 inline bool xisinf (double x)
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	69 { return std::isinf (x); }
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	70 #else
6108 143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 5904 diff changeset	71 extern OCTAVE_API bool xisinf (double x);
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	72 #endif
3248 68259f410026 [project @ 1999-07-13 03:34:54 by jwe] jwe parents: 3156 diff changeset	73
6108 143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 5904 diff changeset	74 extern OCTAVE_API bool octave_is_NA (double x);
143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 5904 diff changeset	75 extern OCTAVE_API bool octave_is_NaN_or_NA (double x) GCC_ATTR_DEPRECATED;
4025 cfb762dc9259 [project @ 2002-08-09 06:32:15 by jwe] jwe parents: 3777 diff changeset	76
10268 9a16a61ed43d new optimizations for accumarray Jaroslav Hajek <highegg@gmail.com> parents: 10158 diff changeset	77 // Generic xmin, xmax definitions
9a16a61ed43d new optimizations for accumarray Jaroslav Hajek <highegg@gmail.com> parents: 10158 diff changeset	78 template <class T>
9a16a61ed43d new optimizations for accumarray Jaroslav Hajek <highegg@gmail.com> parents: 10158 diff changeset	79 inline T xmin (T x, T y)
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	80 {
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	81 return x <= y ? x : y;
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	82 }
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	83
10268 9a16a61ed43d new optimizations for accumarray Jaroslav Hajek <highegg@gmail.com> parents: 10158 diff changeset	84 template <class T>
9a16a61ed43d new optimizations for accumarray Jaroslav Hajek <highegg@gmail.com> parents: 10158 diff changeset	85 inline T xmax (T x, T y)
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	86 {
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	87 return x >= y ? x : y;
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	88 }
10268 9a16a61ed43d new optimizations for accumarray Jaroslav Hajek <highegg@gmail.com> parents: 10158 diff changeset	89
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	90 // This form is favorable. GCC will translate (x <= y ? x : y) without a
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	91 // jump, hence the only conditional jump involved will be the first
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	92 // (xisnan), infrequent and hence friendly to branch prediction.
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	93 inline double
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	94 xmin (double x, double y)
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	95 {
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	96 return xisnan (y) ? x : (x <= y ? x : y);
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	97 }
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	98
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	99 inline double
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	100 xmax (double x, double y)
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	101 {
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	102 return xisnan (y) ? x : (x >= y ? x : y);
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	103 }
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	104
6108 143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 5904 diff changeset	105 extern OCTAVE_API Complex acos (const Complex& x);
143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 5904 diff changeset	106 extern OCTAVE_API Complex acosh (const Complex& x);
143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 5904 diff changeset	107 extern OCTAVE_API Complex asin (const Complex& x);
143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 5904 diff changeset	108 extern OCTAVE_API Complex asinh (const Complex& x);
143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 5904 diff changeset	109 extern OCTAVE_API Complex atan (const Complex& x);
143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 5904 diff changeset	110 extern OCTAVE_API Complex atanh (const Complex& x);
3777 b4f260ddd748 [project @ 2001-02-06 15:57:23 by jwe] jwe parents: 3776 diff changeset	111
6108 143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 5904 diff changeset	112 extern OCTAVE_API bool octave_is_NA (const Complex& x);
143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 5904 diff changeset	113 extern OCTAVE_API bool octave_is_NaN_or_NA (const Complex& x);
4469 53ee020af847 [project @ 2003-07-26 03:45:10 by jwe] jwe parents: 4065 diff changeset	114
6108 143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 5904 diff changeset	115 extern OCTAVE_API Complex xmin (const Complex& x, const Complex& y);
143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 5904 diff changeset	116 extern OCTAVE_API Complex xmax (const Complex& x, const Complex& y);
3248 68259f410026 [project @ 1999-07-13 03:34:54 by jwe] jwe parents: 3156 diff changeset	117
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11525 diff changeset	118 // Single Precision
11237 110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11213 diff changeset	119 extern OCTAVE_API float xtrunc (float x);
13737 30414ff19d5e Use copysignf module from gnulib. Michael Goffioul <michael.goffioul@gmail.com> parents: 13720 diff changeset	120 extern OCTAVE_API float xcopysign (float x, float y);
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	121 inline float xceil (float x) { return ceilf (x); }
15409 c898dbe14e1d build: Update imported module list from gnulib. Rik <rik@octave.org> parents: 15271 diff changeset	122 extern OCTAVE_API float xfloor (float x);
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	123 inline float arg (float x) { return atan2f (0.0f, x); }
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	124 inline float conj (float x) { return x; }
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	125 inline float fix (float x) { return xtrunc (x); }
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	126 inline float imag (float) { return 0.0f; }
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	127 inline float real (float x) { return x; }
11237 110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11213 diff changeset	128 extern OCTAVE_API float xround (float x);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	129 extern OCTAVE_API float xroundb (float x);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	130 extern OCTAVE_API float signum (float x);
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11525 diff changeset	131 extern OCTAVE_API float xlog2 (float x);
12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11525 diff changeset	132 extern OCTAVE_API FloatComplex xlog2 (const FloatComplex& x);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	133 extern OCTAVE_API float xlog2 (float x, int& exp);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	134 extern OCTAVE_API FloatComplex xlog2 (const FloatComplex& x, int& exp);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	135 extern OCTAVE_API float xexp2 (float x);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	136
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	137 #if defined (HAVE_CMATH_ISNANF)
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	138 inline bool xisnan (float x)
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	139 { return std::isnan (x); }
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	140 #else
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	141 extern OCTAVE_API bool xisnan (float x);
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	142 #endif
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	143 #if defined (HAVE_CMATH_ISFINITEF)
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	144 inline bool xfinite (float x)
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	145 { return std::isfinite (x); }
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	146 #else
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	147 extern OCTAVE_API bool xfinite (float x);
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	148 #endif
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	149 #if defined (HAVE_CMATH_ISINFF)
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	150 inline bool xisinf (float x)
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	151 { return std::isinf (x); }
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	152 #else
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	153 extern OCTAVE_API bool xisinf (float x);
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	154 #endif
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	155
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	156 extern OCTAVE_API bool octave_is_NA (float x);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	157 extern OCTAVE_API bool octave_is_NaN_or_NA (float x) GCC_ATTR_DEPRECATED;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	158
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	159 inline float
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	160 xmin (float x, float y)
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	161 {
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	162 return xisnan (y) ? x : (x <= y ? x : y);
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	163 }
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	164
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	165 inline float
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	166 xmax (float x, float y)
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	167 {
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	168 return xisnan (y) ? x : (x >= y ? x : y);
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	169 }
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	170
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	171 extern OCTAVE_API FloatComplex acos (const FloatComplex& x);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	172 extern OCTAVE_API FloatComplex acosh (const FloatComplex& x);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	173 extern OCTAVE_API FloatComplex asin (const FloatComplex& x);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	174 extern OCTAVE_API FloatComplex asinh (const FloatComplex& x);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	175 extern OCTAVE_API FloatComplex atan (const FloatComplex& x);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	176 extern OCTAVE_API FloatComplex atanh (const FloatComplex& x);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	177
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	178 extern OCTAVE_API bool octave_is_NA (const FloatComplex& x);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	179 extern OCTAVE_API bool octave_is_NaN_or_NA (const FloatComplex& x);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	180
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	181 extern OCTAVE_API FloatComplex xmin (const FloatComplex& x, const FloatComplex& y);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	182 extern OCTAVE_API FloatComplex xmax (const FloatComplex& x, const FloatComplex& y);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	183
9812 f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 8998 diff changeset	184 // These map reals to Complex.
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 8998 diff changeset	185
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 8998 diff changeset	186 extern OCTAVE_API Complex rc_acos (double);
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 8998 diff changeset	187 extern OCTAVE_API FloatComplex rc_acos (float);
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 8998 diff changeset	188 extern OCTAVE_API Complex rc_acosh (double);
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 8998 diff changeset	189 extern OCTAVE_API FloatComplex rc_acosh (float);
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 8998 diff changeset	190 extern OCTAVE_API Complex rc_asin (double);
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 8998 diff changeset	191 extern OCTAVE_API FloatComplex rc_asin (float);
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 8998 diff changeset	192 extern OCTAVE_API Complex rc_atanh (double);
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 8998 diff changeset	193 extern OCTAVE_API FloatComplex rc_atanh (float);
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 8998 diff changeset	194 extern OCTAVE_API Complex rc_log (double);
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 8998 diff changeset	195 extern OCTAVE_API FloatComplex rc_log (float);
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 8998 diff changeset	196 extern OCTAVE_API Complex rc_log2 (double);
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 8998 diff changeset	197 extern OCTAVE_API FloatComplex rc_log2 (float);
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 8998 diff changeset	198 extern OCTAVE_API Complex rc_log10 (double);
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 8998 diff changeset	199 extern OCTAVE_API FloatComplex rc_log10 (float);
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 8998 diff changeset	200 extern OCTAVE_API Complex rc_sqrt (double);
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 8998 diff changeset	201 extern OCTAVE_API FloatComplex rc_sqrt (float);
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 8998 diff changeset	202
11010 9478b216752e simplify more array tests Jaroslav Hajek <highegg@gmail.com> parents: 10521 diff changeset	203 // Some useful tests, that are commonly repeated.
9478b216752e simplify more array tests Jaroslav Hajek <highegg@gmail.com> parents: 10521 diff changeset	204 // Test for a finite integer.
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	205 inline bool
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	206 xisinteger (double x)
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	207 {
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	208 return xfinite (x) && x == xround (x);
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	209 }
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	210
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	211 inline bool
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	212 xisinteger (float x)
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	213 {
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	214 return xfinite (x) && x == xround (x);
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	215 }
11010 9478b216752e simplify more array tests Jaroslav Hajek <highegg@gmail.com> parents: 10521 diff changeset	216
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11525 diff changeset	217 // Test for negative sign.
11010 9478b216752e simplify more array tests Jaroslav Hajek <highegg@gmail.com> parents: 10521 diff changeset	218 extern OCTAVE_API bool xnegative_sign (double x);
9478b216752e simplify more array tests Jaroslav Hajek <highegg@gmail.com> parents: 10521 diff changeset	219 extern OCTAVE_API bool xnegative_sign (float x);
9478b216752e simplify more array tests Jaroslav Hajek <highegg@gmail.com> parents: 10521 diff changeset	220
13756 6dfebfa334cb allow negative data log plots with OpenGL+FLTK graphics (bug #34232) John W. Eaton <jwe@octave.org> parents: 13737 diff changeset	221 // Test for positive sign.
6dfebfa334cb allow negative data log plots with OpenGL+FLTK graphics (bug #34232) John W. Eaton <jwe@octave.org> parents: 13737 diff changeset	222 inline bool xpositive_sign (double x) { return ! xnegative_sign (x); }
6dfebfa334cb allow negative data log plots with OpenGL+FLTK graphics (bug #34232) John W. Eaton <jwe@octave.org> parents: 13737 diff changeset	223 inline bool xpositive_sign (float x) { return ! xnegative_sign (x); }
6dfebfa334cb allow negative data log plots with OpenGL+FLTK graphics (bug #34232) John W. Eaton <jwe@octave.org> parents: 13737 diff changeset	224
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	225 // Some old rounding functions.
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	226
11211 2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	227 extern OCTAVE_API octave_idx_type NINTbig (double x);
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	228 extern OCTAVE_API octave_idx_type NINTbig (float x);
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	229
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	230 extern OCTAVE_API int NINT (double x);
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	231 extern OCTAVE_API int NINT (float x);
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	232
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	233 template <typename T>
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	234 T
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	235 X_NINT (T x)
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	236 {
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	237 return (xisinf (x) \|\| xisnan (x)) ? x : xfloor (x + 0.5);
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	238 }
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	239
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	240 inline OCTAVE_API double D_NINT (double x) { return X_NINT (x); }
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	241 inline OCTAVE_API float F_NINT (float x) { return X_NINT (x); }
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	242
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	243 // Template functions can have either float or double arguments.
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	244
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	245 template <typename T>
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	246 bool
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	247 xisnan (const std::complex<T>& x)
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	248 {
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	249 return (xisnan (real (x)) \|\| xisnan (imag (x)));
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	250 }
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	251
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	252 template <typename T>
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	253 bool
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	254 xfinite (const std::complex<T>& x)
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	255 {
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	256 return (xfinite (real (x)) && xfinite (imag (x)));
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	257 }
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	258
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	259 template <typename T>
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	260 bool
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	261 xisinf (const std::complex<T>& x)
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	262 {
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	263 return (xisinf (real (x)) \|\| xisinf (imag (x)));
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	264 }
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	265
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	266 template <typename T>
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	267 std::complex<T>
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	268 fix (const std::complex<T>& x)
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	269 {
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	270 return std::complex<T> (fix (real (x)), fix (imag (x)));
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	271 }
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	272
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	273 template <typename T>
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	274 std::complex<T>
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	275 ceil (const std::complex<T>& x)
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	276 {
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	277 return std::complex<T> (xceil (real (x)), xceil (imag (x)));
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	278 }
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	279
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	280 template <typename T>
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	281 std::complex<T>
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	282 floor (const std::complex<T>& x)
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	283 {
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	284 return std::complex<T> (xfloor (real (x)), xfloor (imag (x)));
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	285 }
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	286
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	287 template <typename T>
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	288 std::complex<T>
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	289 xround (const std::complex<T>& x)
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	290 {
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	291 return std::complex<T> (xround (real (x)), xround (imag (x)));
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	292 }
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	293
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	294 template <typename T>
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	295 std::complex<T>
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	296 xroundb (const std::complex<T>& x)
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	297 {
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	298 return std::complex<T> (xroundb (real (x)), xroundb (imag (x)));
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	299 }
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	300
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	301 template <typename T>
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	302 std::complex<T>
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	303 signum (const std::complex<T>& x)
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	304 {
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	305 T tmp = abs (x);
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	306
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	307 return tmp == 0 ? 0.0 : x / tmp;
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	308 }
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	309
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	310 template <typename T>
11211 2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	311 T
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	312 xmod (T x, T y)
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	313 {
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	314 T retval;
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	315
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	316 if (y == 0)
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	317 retval = x;
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	318 else
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	319 {
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	320 T q = x / y;
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	321
11361 140c8fbe8b0b use correct types in mod template function John W. Eaton <jwe@octave.org> parents: 11237 diff changeset	322 T n = xfloor (q);
11211 2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	323
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	324 if (X_NINT (y) != y)
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	325 {
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	326 if (X_NINT (q) == q)
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	327 n = q;
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	328 else
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	329 {
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	330 if (x >= -1 && x <= 1)
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	331 {
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	332 if (std::abs (q - X_NINT (q))
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	333 < std::numeric_limits<T>::epsilon ())
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	334 n = X_NINT (q);
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	335 }
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	336 else
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	337 {
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	338 if (std::abs ((q - X_NINT (q))/ X_NINT (q))
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	339 < std::numeric_limits<T>::epsilon ())
11361 140c8fbe8b0b use correct types in mod template function John W. Eaton <jwe@octave.org> parents: 11237 diff changeset	340 n = X_NINT (q);
11211 2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	341 }
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	342 }
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	343 }
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	344
11397 57868a21170e lo-mappers.h (xrem, xmod): prevent use of extra precision John W. Eaton <jwe@octave.org> parents: 11361 diff changeset	345 // Prevent use of extra precision.
57868a21170e lo-mappers.h (xrem, xmod): prevent use of extra precision John W. Eaton <jwe@octave.org> parents: 11361 diff changeset	346 volatile T tmp = y * n;
57868a21170e lo-mappers.h (xrem, xmod): prevent use of extra precision John W. Eaton <jwe@octave.org> parents: 11361 diff changeset	347
57868a21170e lo-mappers.h (xrem, xmod): prevent use of extra precision John W. Eaton <jwe@octave.org> parents: 11361 diff changeset	348 retval = x - tmp;
11211 2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	349 }
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	350
11213 009d16b010fa lo-mappers.cc (xmod, xrem): don't copy sign if result is zero John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	351 if (x != y && y != 0 && retval != 0)
11211 2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	352 retval = xcopysign (retval, y);
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	353
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	354 return retval;
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	355 }
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	356
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	357 template <typename T>
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	358 T
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	359 xrem (T x, T y)
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	360 {
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	361 T retval;
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	362
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	363 if (y == 0)
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	364 retval = x;
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	365 else
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	366 {
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	367 T q = x / y;
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	368
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	369 T n = xtrunc (q);
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	370
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	371 if (X_NINT (y) != y)
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	372 {
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	373 if (X_NINT (q) == q)
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	374 n = q;
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	375 else
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	376 {
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	377 if (x >= -1 && x <= 1)
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	378 {
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	379 if (std::abs (q - X_NINT (q))
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	380 < std::numeric_limits<T>::epsilon ())
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	381 n = X_NINT (q);
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	382 }
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	383 else
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	384 {
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	385 if (std::abs ((q - X_NINT (q))/ X_NINT (q))
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	386 < std::numeric_limits<T>::epsilon ())
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	387 n = X_NINT (q);
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	388 }
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	389 }
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	390 }
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	391
11397 57868a21170e lo-mappers.h (xrem, xmod): prevent use of extra precision John W. Eaton <jwe@octave.org> parents: 11361 diff changeset	392 // Prevent use of extra precision.
57868a21170e lo-mappers.h (xrem, xmod): prevent use of extra precision John W. Eaton <jwe@octave.org> parents: 11361 diff changeset	393 volatile T tmp = y * n;
57868a21170e lo-mappers.h (xrem, xmod): prevent use of extra precision John W. Eaton <jwe@octave.org> parents: 11361 diff changeset	394
57868a21170e lo-mappers.h (xrem, xmod): prevent use of extra precision John W. Eaton <jwe@octave.org> parents: 11361 diff changeset	395 retval = x - tmp;
11211 2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	396 }
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	397
11213 009d16b010fa lo-mappers.cc (xmod, xrem): don't copy sign if result is zero John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	398 if (x != y && y != 0 && retval != 0)
11211 2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	399 retval = xcopysign (retval, x);
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	400
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	401 return retval;
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	402 }
11010 9478b216752e simplify more array tests Jaroslav Hajek <highegg@gmail.com> parents: 10521 diff changeset	403
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	404 #endif

Mercurial > octave

annotate liboctave/numeric/lo-mappers.h @ 15409:c898dbe14e1d