octave: liboctave/lo-mappers.cc annotate

annotate liboctave/lo-mappers.cc @ 11523:fd0a3ac60b0e

update copyright notices

author	John W. Eaton <jwe@octave.org>
date	Fri, 14 Jan 2011 05:47:45 -0500
parents	110e570e5f8d
children	ab231f944252

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

Mercurial > octave

annotate liboctave/lo-mappers.cc @ 11523:fd0a3ac60b0e