octave-nkf: liboctave/lo-mappers.cc annotate

annotate liboctave/lo-mappers.cc @ 11525:ab231f944252

avoid exposing gnulib:: in header files

author	John W. Eaton <jwe@octave.org>
date	Fri, 14 Jan 2011 06:20:06 -0500
parents	fd0a3ac60b0e
children	12df7854fa7c

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
11523 fd0a3ac60b0e update copyright notices John W. Eaton <jwe@octave.org> parents: 11237 diff changeset	3 Copyright (C) 1996-2011 John W. Eaton
10521 4d1fc073fbb7 add some missing copyright stmts Jaroslav Hajek <highegg@gmail.com> parents: 10459 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: 6596 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: 6596 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: 6596 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: 6596 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 #ifdef HAVE_CONFIG_H
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	25 #include <config.h>
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	26 #endif
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	27
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	28 #include <cfloat>
3268 fdc7dd08cd85 [project @ 1999-09-10 05:16:49 by jwe] jwe parents: 3248 diff changeset	29
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	30 #include "lo-error.h"
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	31 #include "lo-ieee.h"
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	32 #include "lo-mappers.h"
7231 2eb392d058bb [project @ 2007-11-30 18:53:29 by jwe] jwe parents: 7017 diff changeset	33 #include "lo-math.h"
3121 923049908004 [project @ 1997-12-01 00:01:00 by jwe] jwe parents: 3092 diff changeset	34 #include "lo-specfun.h"
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	35 #include "lo-utils.h"
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	36 #include "oct-cmplx.h"
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	37
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	38 #include "f77-fcn.h"
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	39
3248 68259f410026 [project @ 1999-07-13 03:34:54 by jwe] jwe parents: 3206 diff changeset	40 // double -> double mappers.
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	41
11237 110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	42 // Both xtrunc and xround belong here so we can keep gnulib:: out of
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	43 // lo-mappers.h.
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	44
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	45 double
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	46 xtrunc (double x)
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	47 {
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	48 return gnulib::trunc (x);
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	49 }
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	50
11525 ab231f944252 avoid exposing gnulib:: in header files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	51 double xfloor (double x)
ab231f944252 avoid exposing gnulib:: in header files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	52 {
ab231f944252 avoid exposing gnulib:: in header files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	53 return gnulib::floor (x);
ab231f944252 avoid exposing gnulib:: in header files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	54 }
ab231f944252 avoid exposing gnulib:: in header files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	55
11237 110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	56 double
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	57 xround (double x)
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	58 {
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	59 return gnulib::round (x);
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	60 }
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	61
7636 99c410f7f0b0 implement mapper function for banker's rounding Jaroslav Hajek <highegg@gmail.com> parents: 7231 diff changeset	62 double
99c410f7f0b0 implement mapper function for banker's rounding Jaroslav Hajek <highegg@gmail.com> parents: 7231 diff changeset	63 xroundb (double x)
99c410f7f0b0 implement mapper function for banker's rounding Jaroslav Hajek <highegg@gmail.com> parents: 7231 diff changeset	64 {
99c410f7f0b0 implement mapper function for banker's rounding Jaroslav Hajek <highegg@gmail.com> parents: 7231 diff changeset	65 double t = xround (x);
99c410f7f0b0 implement mapper function for banker's rounding Jaroslav Hajek <highegg@gmail.com> parents: 7231 diff changeset	66
99c410f7f0b0 implement mapper function for banker's rounding Jaroslav Hajek <highegg@gmail.com> parents: 7231 diff changeset	67 if (fabs (x - t) == 0.5)
99c410f7f0b0 implement mapper function for banker's rounding Jaroslav Hajek <highegg@gmail.com> parents: 7231 diff changeset	68 t = 2 * xtrunc (0.5 * t);
99c410f7f0b0 implement mapper function for banker's rounding Jaroslav Hajek <highegg@gmail.com> parents: 7231 diff changeset	69
99c410f7f0b0 implement mapper function for banker's rounding Jaroslav Hajek <highegg@gmail.com> parents: 7231 diff changeset	70 return t;
99c410f7f0b0 implement mapper function for banker's rounding Jaroslav Hajek <highegg@gmail.com> parents: 7231 diff changeset	71 }
99c410f7f0b0 implement mapper function for banker's rounding Jaroslav Hajek <highegg@gmail.com> parents: 7231 diff changeset	72
99c410f7f0b0 implement mapper function for banker's rounding Jaroslav Hajek <highegg@gmail.com> parents: 7231 diff changeset	73 double
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	74 signum (double x)
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	75 {
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	76 double tmp = 0.0;
4984 98e4b7da360a [project @ 2004-09-10 14:02:41 by jwe] jwe parents: 4968 diff changeset	77
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	78 if (x < 0.0)
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	79 tmp = -1.0;
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	80 else if (x > 0.0)
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	81 tmp = 1.0;
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	82
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	83 return xisnan (x) ? octave_NaN : tmp;
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	84 }
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	85
5904 80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	86 double
80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	87 xlog2 (double x)
80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	88 {
80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	89 #if defined (HAVE_LOG2)
80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	90 return log2 (x);
80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	91 #else
80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	92 #if defined (M_LN2)
80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	93 static double ln2 = M_LN2;
80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	94 #else
7740 39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	95 static double ln2 = log (2);
5904 80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	96 #endif
80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	97
80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	98 return log (x) / ln2;
80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	99 #endif
80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	100 }
80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	101
7740 39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	102 Complex
39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	103 xlog2 (const Complex& x)
39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	104 {
39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	105 #if defined (M_LN2)
39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	106 static double ln2 = M_LN2;
39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	107 #else
39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	108 static double ln2 = log (2);
39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	109 #endif
39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	110
39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	111 return std::log (x) / ln2;
39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	112 }
39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	113
5904 80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	114 double
80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	115 xexp2 (double x)
80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	116 {
80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	117 #if defined (HAVE_EXP2)
80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	118 return exp2 (x);
80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	119 #else
80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	120 #if defined (M_LN2)
80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	121 static double ln2 = M_LN2;
80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	122 #else
7740 39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	123 static double ln2 = log (2);
5904 80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	124 #endif
80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	125
80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	126 return exp (x * ln2);
80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	127 #endif
80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	128 }
80d3933fb8b6 [project @ 2006-07-26 17:19:10 by jwe] jwe parents: 5775 diff changeset	129
7740 39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	130 double
39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	131 xlog2 (double x, int& exp)
39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	132 {
39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	133 return frexp (x, &exp);
39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	134 }
39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	135
39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	136 Complex
39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	137 xlog2 (const Complex& x, int& exp)
39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	138 {
39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	139 double ax = std::abs (x);
39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	140 double lax = xlog2 (ax, exp);
9319 0d9178575dd7 fix log2 with 2 outargs, loosen tests to meet IEEE Jaroslav Hajek <highegg@gmail.com> parents: 8998 diff changeset	141 return (ax != lax) ? (x / ax) * lax : x;
7740 39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	142 }
39930366b709 implement builtin log2 Jaroslav Hajek <highegg@gmail.com> parents: 7707 diff changeset	143
3248 68259f410026 [project @ 1999-07-13 03:34:54 by jwe] jwe parents: 3206 diff changeset	144 // double -> bool mappers.
68259f410026 [project @ 1999-07-13 03:34:54 by jwe] jwe parents: 3206 diff changeset	145
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	146 #if ! defined(HAVE_CMATH_ISNAN)
3248 68259f410026 [project @ 1999-07-13 03:34:54 by jwe] jwe parents: 3206 diff changeset	147 bool
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	148 xisnan (double x)
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	149 {
4074 a67f76924703 [project @ 2002-09-27 21:06:02 by jwe] jwe parents: 4065 diff changeset	150 return lo_ieee_isnan (x);
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	151 }
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	152 #endif
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	153
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	154 #if ! defined(HAVE_CMATH_ISFINITE)
3248 68259f410026 [project @ 1999-07-13 03:34:54 by jwe] jwe parents: 3206 diff changeset	155 bool
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	156 xfinite (double x)
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	157 {
4074 a67f76924703 [project @ 2002-09-27 21:06:02 by jwe] jwe parents: 4065 diff changeset	158 return lo_ieee_finite (x);
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	159 }
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	160 #endif
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	161
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	162 #if ! defined(HAVE_CMATH_ISINF)
3248 68259f410026 [project @ 1999-07-13 03:34:54 by jwe] jwe parents: 3206 diff changeset	163 bool
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	164 xisinf (double x)
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	165 {
4074 a67f76924703 [project @ 2002-09-27 21:06:02 by jwe] jwe parents: 4065 diff changeset	166 return lo_ieee_isinf (x);
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	167 }
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	168 #endif
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	169
4025 cfb762dc9259 [project @ 2002-08-09 06:32:15 by jwe] jwe parents: 3873 diff changeset	170 bool
cfb762dc9259 [project @ 2002-08-09 06:32:15 by jwe] jwe parents: 3873 diff changeset	171 octave_is_NA (double x)
cfb762dc9259 [project @ 2002-08-09 06:32:15 by jwe] jwe parents: 3873 diff changeset	172 {
cfb762dc9259 [project @ 2002-08-09 06:32:15 by jwe] jwe parents: 3873 diff changeset	173 return lo_ieee_is_NA (x);
cfb762dc9259 [project @ 2002-08-09 06:32:15 by jwe] jwe parents: 3873 diff changeset	174 }
cfb762dc9259 [project @ 2002-08-09 06:32:15 by jwe] jwe parents: 3873 diff changeset	175
cfb762dc9259 [project @ 2002-08-09 06:32:15 by jwe] jwe parents: 3873 diff changeset	176 bool
cfb762dc9259 [project @ 2002-08-09 06:32:15 by jwe] jwe parents: 3873 diff changeset	177 octave_is_NaN_or_NA (double x)
cfb762dc9259 [project @ 2002-08-09 06:32:15 by jwe] jwe parents: 3873 diff changeset	178 {
5389 25c8956d2204 [project @ 2005-06-15 03:45:46 by jwe] jwe parents: 5307 diff changeset	179 return lo_ieee_isnan (x);
4025 cfb762dc9259 [project @ 2002-08-09 06:32:15 by jwe] jwe parents: 3873 diff changeset	180 }
cfb762dc9259 [project @ 2002-08-09 06:32:15 by jwe] jwe parents: 3873 diff changeset	181
3248 68259f410026 [project @ 1999-07-13 03:34:54 by jwe] jwe parents: 3206 diff changeset	182 // (double, double) -> double mappers.
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	183
3248 68259f410026 [project @ 1999-07-13 03:34:54 by jwe] jwe parents: 3206 diff changeset	184 // complex -> complex mappers.
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	185
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	186 Complex
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	187 acos (const Complex& x)
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	188 {
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	189 static Complex i (0, 1);
3056 b08a007edd15 [project @ 1997-06-16 02:12:15 by jwe] jwe parents: 3049 diff changeset	190
6414 250eebbd7ccc [project @ 2007-03-15 16:04:52 by jwe] jwe parents: 5904 diff changeset	191 return -i * (log (x + i * (sqrt (1.0 - x*x))));
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	192 }
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	193
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	194 Complex
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	195 acosh (const Complex& x)
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	196 {
3056 b08a007edd15 [project @ 1997-06-16 02:12:15 by jwe] jwe parents: 3049 diff changeset	197 return log (x + sqrt (x*x - 1.0));
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	198 }
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	199
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	200 Complex
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	201 asin (const Complex& x)
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	202 {
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	203 static Complex i (0, 1);
3056 b08a007edd15 [project @ 1997-06-16 02:12:15 by jwe] jwe parents: 3049 diff changeset	204
b08a007edd15 [project @ 1997-06-16 02:12:15 by jwe] jwe parents: 3049 diff changeset	205 return -i * log (ix + sqrt (1.0 - xx));
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	206 }
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	207
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	208 Complex
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	209 asinh (const Complex& x)
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	210 {
3056 b08a007edd15 [project @ 1997-06-16 02:12:15 by jwe] jwe parents: 3049 diff changeset	211 return log (x + sqrt (x*x + 1.0));
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	212 }
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	213
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	214 Complex
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	215 atan (const Complex& x)
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	216 {
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	217 static Complex i (0, 1);
3056 b08a007edd15 [project @ 1997-06-16 02:12:15 by jwe] jwe parents: 3049 diff changeset	218
b08a007edd15 [project @ 1997-06-16 02:12:15 by jwe] jwe parents: 3049 diff changeset	219 return i * log ((i + x) / (i - x)) / 2.0;
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	220 }
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	221
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	222 Complex
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	223 atanh (const Complex& x)
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	224 {
3092 4bb976b250bf [project @ 1997-10-17 07:48:48 by jwe] jwe parents: 3069 diff changeset	225 return log ((1.0 + x) / (1.0 - x)) / 2.0;
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	226 }
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	227
3248 68259f410026 [project @ 1999-07-13 03:34:54 by jwe] jwe parents: 3206 diff changeset	228 // complex -> bool mappers.
68259f410026 [project @ 1999-07-13 03:34:54 by jwe] jwe parents: 3206 diff changeset	229
68259f410026 [project @ 1999-07-13 03:34:54 by jwe] jwe parents: 3206 diff changeset	230 bool
4025 cfb762dc9259 [project @ 2002-08-09 06:32:15 by jwe] jwe parents: 3873 diff changeset	231 octave_is_NA (const Complex& x)
cfb762dc9259 [project @ 2002-08-09 06:32:15 by jwe] jwe parents: 3873 diff changeset	232 {
cfb762dc9259 [project @ 2002-08-09 06:32:15 by jwe] jwe parents: 3873 diff changeset	233 return (octave_is_NA (real (x)) \|\| octave_is_NA (imag (x)));
cfb762dc9259 [project @ 2002-08-09 06:32:15 by jwe] jwe parents: 3873 diff changeset	234 }
cfb762dc9259 [project @ 2002-08-09 06:32:15 by jwe] jwe parents: 3873 diff changeset	235
cfb762dc9259 [project @ 2002-08-09 06:32:15 by jwe] jwe parents: 3873 diff changeset	236 bool
cfb762dc9259 [project @ 2002-08-09 06:32:15 by jwe] jwe parents: 3873 diff changeset	237 octave_is_NaN_or_NA (const Complex& x)
cfb762dc9259 [project @ 2002-08-09 06:32:15 by jwe] jwe parents: 3873 diff changeset	238 {
5389 25c8956d2204 [project @ 2005-06-15 03:45:46 by jwe] jwe parents: 5307 diff changeset	239 return (xisnan (real (x)) \|\| xisnan (imag (x)));
4025 cfb762dc9259 [project @ 2002-08-09 06:32:15 by jwe] jwe parents: 3873 diff changeset	240 }
cfb762dc9259 [project @ 2002-08-09 06:32:15 by jwe] jwe parents: 3873 diff changeset	241
3248 68259f410026 [project @ 1999-07-13 03:34:54 by jwe] jwe parents: 3206 diff changeset	242 // (complex, complex) -> complex mappers.
68259f410026 [project @ 1999-07-13 03:34:54 by jwe] jwe parents: 3206 diff changeset	243
5775 ace8d8d26933 [project @ 2006-04-24 19:13:06 by jwe] jwe parents: 5389 diff changeset	244 // FIXME -- need to handle NA too?
4025 cfb762dc9259 [project @ 2002-08-09 06:32:15 by jwe] jwe parents: 3873 diff changeset	245
3248 68259f410026 [project @ 1999-07-13 03:34:54 by jwe] jwe parents: 3206 diff changeset	246 Complex
68259f410026 [project @ 1999-07-13 03:34:54 by jwe] jwe parents: 3206 diff changeset	247 xmin (const Complex& x, const Complex& y)
68259f410026 [project @ 1999-07-13 03:34:54 by jwe] jwe parents: 3206 diff changeset	248 {
3757 574711ce9070 [project @ 2000-12-15 18:19:15 by jwe] jwe parents: 3268 diff changeset	249 return abs (x) <= abs (y) ? x : (xisnan (x) ? x : y);
3248 68259f410026 [project @ 1999-07-13 03:34:54 by jwe] jwe parents: 3206 diff changeset	250 }
68259f410026 [project @ 1999-07-13 03:34:54 by jwe] jwe parents: 3206 diff changeset	251
68259f410026 [project @ 1999-07-13 03:34:54 by jwe] jwe parents: 3206 diff changeset	252 Complex
68259f410026 [project @ 1999-07-13 03:34:54 by jwe] jwe parents: 3206 diff changeset	253 xmax (const Complex& x, const Complex& y)
68259f410026 [project @ 1999-07-13 03:34:54 by jwe] jwe parents: 3206 diff changeset	254 {
3757 574711ce9070 [project @ 2000-12-15 18:19:15 by jwe] jwe parents: 3268 diff changeset	255 return abs (x) >= abs (y) ? x : (xisnan (x) ? x : y);
3248 68259f410026 [project @ 1999-07-13 03:34:54 by jwe] jwe parents: 3206 diff changeset	256 }
68259f410026 [project @ 1999-07-13 03:34:54 by jwe] jwe parents: 3206 diff changeset	257
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	258
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	259 // float -> float mappers.
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	260
11237 110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	261 // Both xtrunc and xround belong here so we can keep gnulib:: out of
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	262 // lo-mappers.h.
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	263
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	264 float
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	265 xtrunc (float x)
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	266 {
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	267 return gnulib::truncf (x);
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	268 }
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	269
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	270 float
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	271 xround (float x)
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	272 {
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	273 return gnulib::round (x);
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	274 }
110e570e5f8d keep gnulib out of lo-mappers.h John W. Eaton <jwe@octave.org> parents: 11212 diff changeset	275
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	276 float
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	277 xroundb (float x)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	278 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	279 float t = xround (x);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	280
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	281 if (fabsf (x - t) == 0.5)
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	282 t = 2 * xtrunc (0.5 * t);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	283
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	284 return t;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	285 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	286
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	287 float
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	288 signum (float x)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	289 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	290 float tmp = 0.0;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	291
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	292 if (x < 0.0)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	293 tmp = -1.0;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	294 else if (x > 0.0)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	295 tmp = 1.0;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	296
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	297 return xisnan (x) ? octave_Float_NaN : tmp;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	298 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	299
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	300 float
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	301 xlog2 (float x)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	302 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	303 #if defined (HAVE_LOG2)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	304 return log2 (x);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	305 #else
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	306 #if defined (M_LN2)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	307 static float ln2 = M_LN2;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	308 #else
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	309 static float ln2 = log2 (2);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	310 #endif
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	311
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	312 return log (x) / ln2;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	313 #endif
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	314 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	315
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	316 FloatComplex
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	317 xlog2 (const FloatComplex& x)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	318 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	319 #if defined (M_LN2)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	320 static float ln2 = M_LN2;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	321 #else
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	322 static float ln2 = log (2);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	323 #endif
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	324
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	325 return std::log (x) / ln2;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	326 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	327
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	328 float
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	329 xexp2 (float x)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	330 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	331 #if defined (HAVE_EXP2)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	332 return exp2 (x);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	333 #else
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	334 #if defined (M_LN2)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	335 static float ln2 = M_LN2;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	336 #else
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	337 static float ln2 = log2 (2);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	338 #endif
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	339
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	340 return exp (x * ln2);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	341 #endif
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	342 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	343
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	344 float
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	345 xlog2 (float x, int& exp)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	346 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	347 return frexpf (x, &exp);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	348 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	349
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	350 FloatComplex
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	351 xlog2 (const FloatComplex& x, int& exp)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	352 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	353 float ax = std::abs (x);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	354 float lax = xlog2 (ax, exp);
9319 0d9178575dd7 fix log2 with 2 outargs, loosen tests to meet IEEE Jaroslav Hajek <highegg@gmail.com> parents: 8998 diff changeset	355 return (ax != lax) ? (x / ax) * lax : x;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	356 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	357
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	358 // float -> bool mappers.
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	359
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	360 #if ! defined(HAVE_CMATH_ISNANF)
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	361 bool
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	362 xisnan (float x)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	363 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	364 return lo_ieee_isnan (x);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	365 }
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	366 #endif
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	367
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	368 #if ! defined(HAVE_CMATH_ISFINITEF)
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	369 bool
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	370 xfinite (float x)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	371 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	372 return lo_ieee_finite (x);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	373 }
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	374 #endif
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	375
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	376 #if ! defined(HAVE_CMATH_ISINFF)
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	377 bool
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	378 xisinf (float x)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	379 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	380 return lo_ieee_isinf (x);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	381 }
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	382 #endif
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	383
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	384 bool
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	385 octave_is_NA (float x)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	386 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	387 return lo_ieee_is_NA (x);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	388 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	389
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	390 bool
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	391 octave_is_NaN_or_NA (float x)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	392 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	393 return lo_ieee_isnan (x);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	394 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	395
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	396 // (float, float) -> float mappers.
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	397
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	398 // complex -> complex mappers.
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	399
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	400 FloatComplex
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	401 acos (const FloatComplex& x)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	402 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	403 static FloatComplex i (0, 1);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	404
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	405 return -i * (log (x + i * (sqrt (static_cast<float>(1.0) - x*x))));
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	406 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	407
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	408 FloatComplex
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	409 acosh (const FloatComplex& x)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	410 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	411 return log (x + sqrt (x*x - static_cast<float>(1.0)));
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	412 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	413
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	414 FloatComplex
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	415 asin (const FloatComplex& x)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	416 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	417 static FloatComplex i (0, 1);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	418
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	419 return -i * log (ix + sqrt (static_cast<float>(1.0) - xx));
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	420 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	421
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	422 FloatComplex
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	423 asinh (const FloatComplex& x)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	424 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	425 return log (x + sqrt (x*x + static_cast<float>(1.0)));
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	426 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	427
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	428 FloatComplex
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	429 atan (const FloatComplex& x)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	430 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	431 static FloatComplex i (0, 1);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	432
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	433 return i * log ((i + x) / (i - x)) / static_cast<float>(2.0);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	434 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	435
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	436 FloatComplex
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	437 atanh (const FloatComplex& x)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	438 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	439 return log ((static_cast<float>(1.0) + x) / (static_cast<float>(1.0) - x)) / static_cast<float>(2.0);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	440 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	441
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	442 // complex -> bool mappers.
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	443
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	444 bool
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	445 octave_is_NA (const FloatComplex& x)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	446 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	447 return (octave_is_NA (real (x)) \|\| octave_is_NA (imag (x)));
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	448 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	449
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	450 bool
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	451 octave_is_NaN_or_NA (const FloatComplex& x)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	452 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	453 return (xisnan (real (x)) \|\| xisnan (imag (x)));
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	454 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	455
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	456 // (complex, complex) -> complex mappers.
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	457
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	458 // FIXME -- need to handle NA too?
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	459
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	460 FloatComplex
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	461 xmin (const FloatComplex& x, const FloatComplex& y)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	462 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	463 return abs (x) <= abs (y) ? x : (xisnan (x) ? x : y);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	464 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	465
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	466 FloatComplex
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	467 xmax (const FloatComplex& x, const FloatComplex& y)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	468 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	469 return abs (x) >= abs (y) ? x : (xisnan (x) ? x : y);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	470 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7740 diff changeset	471
9812 f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	472 Complex
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	473 rc_acos (double x)
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	474 {
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	475 return fabs (x) > 1.0 ? acos (Complex (x)) : Complex (acos (x));
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	476 }
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	477
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	478 FloatComplex
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	479 rc_acos (float x)
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	480 {
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	481 return fabsf (x) > 1.0f ? acos (FloatComplex (x)) : FloatComplex (acosf (x));
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	482 }
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	483
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	484 Complex
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	485 rc_acosh (double x)
9812 f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	486 {
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	487 return x < 1.0 ? acosh (Complex (x)) : Complex (acosh (x));
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	488 }
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	489
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	490 FloatComplex
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	491 rc_acosh (float x)
9812 f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	492 {
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	493 return x < 1.0f ? acosh (FloatComplex (x)) : FloatComplex (acoshf (x));
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	494 }
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	495
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	496 Complex
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	497 rc_asin (double x)
9812 f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	498 {
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	499 return fabs (x) > 1.0 ? asin (Complex (x)) : Complex (asin (x));
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	500 }
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	501
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	502 FloatComplex
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	503 rc_asin (float x)
9812 f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	504 {
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	505 return fabsf (x) > 1.0f ? asin (FloatComplex (x)) : FloatComplex (asinf (x));
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	506 }
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	507
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	508 Complex
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	509 rc_atanh (double x)
9812 f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	510 {
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	511 return fabs (x) > 1.0 ? atanh (Complex (x)) : Complex (atanh (x));
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	512 }
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	513
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	514 FloatComplex
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	515 rc_atanh (float x)
9812 f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	516 {
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	517 return fabsf (x) > 1.0f ? atanh (FloatComplex (x)) : FloatComplex (atanhf (x));
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	518 }
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	519
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	520 Complex
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	521 rc_log (double x)
9812 f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	522 {
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	523 const double pi = 3.14159265358979323846;
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	524 return x < 0.0 ? Complex (log (-x), pi) : Complex (log (x));
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	525 }
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	526
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	527 FloatComplex
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	528 rc_log (float x)
9812 f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	529 {
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	530 const float pi = 3.14159265358979323846f;
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	531 return x < 0.0f ? FloatComplex (logf (-x), pi) : FloatComplex (logf (x));
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	532 }
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	533
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	534 Complex
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	535 rc_log2 (double x)
9812 f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	536 {
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	537 const double pil2 = 4.53236014182719380962; // = pi / log(2)
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	538 return x < 0.0 ? Complex (xlog2 (-x), pil2) : Complex (xlog2 (x));
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	539 }
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	540
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	541 FloatComplex
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	542 rc_log2 (float x)
9812 f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	543 {
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	544 const float pil2 = 4.53236014182719380962f; // = pi / log(2)
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	545 return x < 0.0f ? FloatComplex (xlog2 (-x), pil2) : FloatComplex (xlog2 (x));
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	546 }
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	547
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	548 Complex
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	549 rc_log10 (double x)
9812 f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	550 {
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	551 const double pil10 = 1.36437635384184134748; // = pi / log(10)
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	552 return x < 0.0 ? Complex (log10 (-x), pil10) : Complex (log10 (x));
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	553 }
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	554
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	555 FloatComplex
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	556 rc_log10 (float x)
9812 f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	557 {
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	558 const float pil10 = 1.36437635384184134748f; // = pi / log(10)
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	559 return x < 0.0f ? FloatComplex (log10 (-x), pil10) : FloatComplex (log10f (x));
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	560 }
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	561
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	562 Complex
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	563 rc_sqrt (double x)
9812 f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	564 {
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	565 return x < 0.0 ? Complex (0.0, sqrt (-x)) : Complex (sqrt (x));
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	566 }
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	567
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	568 FloatComplex
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	569 rc_sqrt (float x)
9812 f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	570 {
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	571 return x < 0.0f ? FloatComplex (0.0f, sqrtf (-x)) : FloatComplex (sqrtf (x));
f80c566bc751 improve unary mapper system Jaroslav Hajek <highegg@gmail.com> parents: 9319 diff changeset	572 }
11010 9478b216752e simplify more array tests Jaroslav Hajek <highegg@gmail.com> parents: 10521 diff changeset	573
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	574 bool
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	575 xnegative_sign (double x)
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	576 {
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	577 return __lo_ieee_signbit (x);
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	578 }
11010 9478b216752e simplify more array tests Jaroslav Hajek <highegg@gmail.com> parents: 10521 diff changeset	579
11212 ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	580 bool
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	581 xnegative_sign (float x)
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	582 {
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	583 return __lo_ieee_float_signbit (x);
ce27d6f4e134 use templates and inline for more lo-mappers functionos John W. Eaton <jwe@octave.org> parents: 11211 diff changeset	584 }
11211 2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	585
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	586 // Convert X to the nearest integer value. Should not pass NaN to
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	587 // this function.
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	588
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	589 // Sometimes you need a large integer, but not always.
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	590
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	591 octave_idx_type
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	592 NINTbig (double x)
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	593 {
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	594 if (x > std::numeric_limits<octave_idx_type>::max ())
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	595 return std::numeric_limits<octave_idx_type>::max ();
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	596 else if (x < std::numeric_limits<octave_idx_type>::min ())
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	597 return std::numeric_limits<octave_idx_type>::min ();
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	598 else
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	599 return static_cast<octave_idx_type> ((x > 0) ? (x + 0.5) : (x - 0.5));
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	600 }
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	601
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	602 octave_idx_type
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	603 NINTbig (float x)
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	604 {
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	605 if (x > std::numeric_limits<octave_idx_type>::max ())
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	606 return std::numeric_limits<octave_idx_type>::max ();
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	607 else if (x < std::numeric_limits<octave_idx_type>::min ())
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	608 return std::numeric_limits<octave_idx_type>::min ();
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	609 else
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	610 return static_cast<octave_idx_type> ((x > 0) ? (x + 0.5) : (x - 0.5));
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	611 }
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	612
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	613 int
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	614 NINT (double x)
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	615 {
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	616 if (x > std::numeric_limits<int>::max ())
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	617 return std::numeric_limits<int>::max ();
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	618 else if (x < std::numeric_limits<int>::min ())
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	619 return std::numeric_limits<int>::min ();
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	620 else
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	621 return static_cast<int> ((x > 0) ? (x + 0.5) : (x - 0.5));
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	622 }
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	623
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	624 int
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	625 NINT (float x)
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	626 {
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	627 if (x > std::numeric_limits<int>::max ())
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	628 return std::numeric_limits<int>::max ();
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	629 else if (x < std::numeric_limits<int>::min ())
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	630 return std::numeric_limits<int>::min ();
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	631 else
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	632 return static_cast<int> ((x > 0) ? (x + 0.5) : (x - 0.5));
2554b4a0806e use templates for some lo-mappers functions John W. Eaton <jwe@octave.org> parents: 11209 diff changeset	633 }

Mercurial > octave-nkf

annotate liboctave/lo-mappers.cc @ 11525:ab231f944252