octave-nkf: liboctave/oct-inttypes.cc annotate

annotate liboctave/oct-inttypes.cc @ 9003:0631d397fbe0

replace lo_ieee_isnan by xisnan, add missing includes

author	Jaroslav Hajek <highegg@gmail.com>
date	Sat, 21 Mar 2009 20:46:21 +0100
parents	eb63fbe60fab
children	3b1908b58662

rev	line source
5072 97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	1 /*
97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	2
8920 eb63fbe60fab update copyright notices John W. Eaton <jwe@octave.org> parents: 8864 diff changeset	3 Copyright (C) 2008, 2009 Jaroslav Hajek
7017 a1dbe9d80eee [project @ 2007-10-12 21:27:11 by jwe] jwe parents: 7016 diff changeset	4 Copyright (C) 2004, 2005, 2006, 2007 John W. Eaton
5072 97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	5
97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	6 This file is part of Octave.
97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	7
97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	8 Octave is free software; you can redistribute it and/or modify it
97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	9 under the terms of the GNU General Public License as published by the
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6108 diff changeset	10 Free Software Foundation; either version 3 of the License, or (at your
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6108 diff changeset	11 option) any later version.
5072 97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	12
97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	13 Octave is distributed in the hope that it will be useful, but WITHOUT
97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	14 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	15 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	16 for more details.
97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	17
97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	18 You should have received a copy of the GNU General Public License
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6108 diff changeset	19 along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6108 diff changeset	20 <http://www.gnu.org/licenses/>.
5072 97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	21
97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	22 */
97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	23
97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	24 #ifdef HAVE_CONFIG_H
97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	25 #include <config.h>
97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	26 #endif
97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	27
8169 66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	28 #include "lo-error.h"
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	29
5072 97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	30 #include "oct-inttypes.h"
97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	31
8185 69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	32 template<class T>
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	33 const octave_int<T> octave_int<T>::zero (static_cast<T> (0));
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	34
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	35 template<class T>
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	36 const octave_int<T> octave_int<T>::one (static_cast<T> (1));
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	37
8169 66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	38 // define type names.
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	39 #define DECLARE_OCTAVE_INT_TYPENAME(TYPE, TYPENAME) \
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	40 template <> \
8319 c374691576f6 Fix for MSVC compilation Michael Goffioul <michael.goffioul@gmail.com> parents: 8293 diff changeset	41 OCTAVE_API const char * \
8169 66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	42 octave_int<TYPE>::type_name () { return TYPENAME; }
5072 97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	43
8169 66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	44 DECLARE_OCTAVE_INT_TYPENAME (int8_t, "int8")
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	45 DECLARE_OCTAVE_INT_TYPENAME (int16_t, "int16")
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	46 DECLARE_OCTAVE_INT_TYPENAME (int32_t, "int32")
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	47 DECLARE_OCTAVE_INT_TYPENAME (int64_t, "int64")
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	48 DECLARE_OCTAVE_INT_TYPENAME (uint8_t, "uint8")
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	49 DECLARE_OCTAVE_INT_TYPENAME (uint16_t, "uint16")
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	50 DECLARE_OCTAVE_INT_TYPENAME (uint32_t, "uint32")
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	51 DECLARE_OCTAVE_INT_TYPENAME (uint64_t, "uint64")
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	52
8185 69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	53 #ifndef OCTAVE_INT_USE_LONG_DOUBLE
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	54
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	55 // Define comparison operators
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	56
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	57 template <class xop>
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	58 OCTAVE_API bool
8856 ab4db66e286f workaround gcc 4.3 explicit instantiation bug in octave_int_cmp_op Jaroslav Hajek <highegg@gmail.com> parents: 8797 diff changeset	59 octave_int_cmp_op::emulate_mop (uint64_t x, double y)
8185 69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	60 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	61 static const double xxup = std::numeric_limits<uint64_t>::max ();
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	62 // This converts to the nearest double. Unless there's an equality, the
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	63 // result is clear.
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	64 double xx = x;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	65 if (xx != y)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	66 return xop::op (xx, y);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	67 else
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	68 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	69 // If equality occured we compare as integers.
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	70 if (xx == xxup)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	71 return xop::gtval;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	72 else
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	73 return xop::op (x, static_cast<uint64_t> (xx));
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	74 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	75 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	76
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	77 template <class xop>
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	78 OCTAVE_API bool
8856 ab4db66e286f workaround gcc 4.3 explicit instantiation bug in octave_int_cmp_op Jaroslav Hajek <highegg@gmail.com> parents: 8797 diff changeset	79 octave_int_cmp_op::emulate_mop (int64_t x, double y)
8169 66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	80 {
8185 69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	81 static const double xxup = std::numeric_limits<int64_t>::max ();
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	82 static const double xxlo = std::numeric_limits<int64_t>::min ();
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	83 // This converts to the nearest double. Unless there's an equality, the
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	84 // result is clear.
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	85 double xx = x;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	86 if (xx != y)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	87 return xop::op (xx, y);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	88 else
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	89 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	90 // If equality occured we compare as integers.
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	91 if (xx == xxup)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	92 return xop::gtval;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	93 else if (xx == xxlo)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	94 return xop::ltval;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	95 else
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	96 return xop::op (x, static_cast<int64_t> (xx));
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	97 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	98
8169 66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	99 }
5072 97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	100
8185 69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	101 // We define double-int operations by reverting the operator
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	102
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	103 // A trait class reverting the operator
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	104 template <class xop>
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	105 class rev_op
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	106 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	107 public:
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	108 typedef xop op;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	109 };
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	110
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	111 #define DEFINE_REVERTED_OPERATOR(OP1,OP2) \
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	112 template <> \
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	113 class rev_op<octave_int_cmp_op::OP1> \
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	114 { \
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	115 public: \
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	116 typedef octave_int_cmp_op::OP2 op; \
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	117 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	118
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	119 DEFINE_REVERTED_OPERATOR(lt,gt);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	120 DEFINE_REVERTED_OPERATOR(gt,lt);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	121 DEFINE_REVERTED_OPERATOR(le,ge);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	122 DEFINE_REVERTED_OPERATOR(ge,le);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	123
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	124 template <class xop>
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	125 OCTAVE_API bool
8856 ab4db66e286f workaround gcc 4.3 explicit instantiation bug in octave_int_cmp_op Jaroslav Hajek <highegg@gmail.com> parents: 8797 diff changeset	126 octave_int_cmp_op::emulate_mop (double x, uint64_t y)
8185 69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	127 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	128 typedef typename rev_op<xop>::op rop;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	129 return mop<rop> (y, x);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	130 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	131
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	132 template <class xop>
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	133 OCTAVE_API bool
8856 ab4db66e286f workaround gcc 4.3 explicit instantiation bug in octave_int_cmp_op Jaroslav Hajek <highegg@gmail.com> parents: 8797 diff changeset	134 octave_int_cmp_op::emulate_mop (double x, int64_t y)
8185 69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	135 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	136 typedef typename rev_op<xop>::op rop;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	137 return mop<rop> (y, x);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	138 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	139
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	140
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	141 // Define handlers for int64 multiplication
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	142
8169 66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	143 template <>
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	144 uint64_t
8185 69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	145 octave_int_arith_base<uint64_t, false>::mul (uint64_t x, uint64_t y)
8169 66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	146 {
8185 69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	147 // Get upper words
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	148 uint64_t ux = x >> 32, uy = y >> 32;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	149 uint64_t res;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	150 if (ux)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	151 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	152 if (uy)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	153 goto overflow;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	154 else
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	155 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	156 uint64_t ly = static_cast<uint32_t> (y), uxly = ux*ly;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	157 if (uxly >> 32)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	158 goto overflow;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	159 uxly <<= 32; // never overflows
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	160 uint64_t lx = static_cast<uint32_t> (x), lxly = lx*ly;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	161 res = add (uxly, lxly);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	162 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	163 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	164 else if (uy)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	165 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	166 uint64_t lx = static_cast<uint32_t> (x), uylx = uy*lx;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	167 if (uylx >> 32)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	168 goto overflow;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	169 uylx <<= 32; // never overflows
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	170 uint64_t ly = static_cast<uint32_t> (y), lylx = ly*lx;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	171 res = add (uylx, lylx);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	172 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	173 else
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	174 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	175 uint64_t lx = static_cast<uint32_t> (x);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	176 uint64_t ly = static_cast<uint32_t> (y);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	177 res = lx*ly;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	178 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	179
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	180 return res;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	181
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	182 overflow:
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	183 ftrunc = true;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	184 return max_val ();
8169 66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	185 }
8185 69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	186
8169 66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	187 template <>
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	188 int64_t
8185 69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	189 octave_int_arith_base<int64_t, true>::mul (int64_t x, int64_t y)
8169 66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	190 {
8185 69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	191 // The signed case is far worse. The problem is that
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	192 // even if neither integer fits into signed 32-bit range, the result may
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	193 // still be OK. Uh oh.
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	194
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	195 // Essentially, what we do is compute sign, multiply absolute values
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	196 // (as above) and impose the sign.
8333 9238637cb81c style fixes John W. Eaton <jwe@octave.org> parents: 8319 diff changeset	197 // FIXME -- can we do something faster if we HAVE_FAST_INT_OPS?
8185 69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	198
8293 ad5bb02d267a workaround missing std::abs(int64_t) in MSVC Jaroslav Hajek <highegg@gmail.com> parents: 8202 diff changeset	199 uint64_t usx = octave_int_abs (x), usy = octave_int_abs (y);
8185 69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	200 bool positive = (x < 0) == (y < 0);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	201
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	202 // Get upper words
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	203 uint64_t ux = usx >> 32, uy = usy >> 32;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	204 uint64_t res;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	205 if (ux)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	206 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	207 if (uy)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	208 goto overflow;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	209 else
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	210 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	211 uint64_t ly = static_cast<uint32_t> (usy), uxly = ux*ly;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	212 if (uxly >> 32)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	213 goto overflow;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	214 uxly <<= 32; // never overflows
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	215 uint64_t lx = static_cast<uint32_t> (usx), lxly = lx*ly;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	216 res = uxly + lxly;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	217 if (res < uxly)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	218 goto overflow;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	219 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	220 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	221 else if (uy)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	222 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	223 uint64_t lx = static_cast<uint32_t> (usx), uylx = uy*lx;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	224 if (uylx >> 32)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	225 goto overflow;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	226 uylx <<= 32; // never overflows
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	227 uint64_t ly = static_cast<uint32_t> (usy), lylx = ly*lx;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	228 res = uylx + lylx;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	229 if (res < uylx)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	230 goto overflow;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	231 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	232 else
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	233 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	234 uint64_t lx = static_cast<uint32_t> (usx);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	235 uint64_t ly = static_cast<uint32_t> (usy);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	236 res = lx*ly;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	237 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	238
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	239 if (positive)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	240 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	241 if (res > static_cast<uint64_t> (max_val ()))
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	242 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	243 ftrunc = true;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	244 return max_val ();
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	245 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	246 else
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	247 return static_cast<int64_t> (res);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	248 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	249 else
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	250 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	251 if (res > static_cast<uint64_t> (-min_val ()))
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	252 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	253 ftrunc = true;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	254 return min_val ();
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	255 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	256 else
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	257 return -static_cast<int64_t> (res);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	258 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	259
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	260
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	261 overflow:
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	262 ftrunc = true;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	263 return positive ? max_val () : min_val ();
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	264
8169 66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	265 }
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	266
8185 69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	267 #define INT_DOUBLE_BINOP_DECL(OP,SUFFIX) \
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	268 template <> \
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	269 OCTAVE_API octave_ ## SUFFIX \
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	270 operator OP (const octave_ ## SUFFIX & x, const double& y)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	271
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	272 #define DOUBLE_INT_BINOP_DECL(OP,SUFFIX) \
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	273 template <> \
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	274 OCTAVE_API octave_ ## SUFFIX \
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	275 operator OP (const double& x, const octave_ ## SUFFIX & y)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	276
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	277 INT_DOUBLE_BINOP_DECL (+, uint64)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	278 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	279 return (y < 0) ? x - octave_uint64(-y) : x + octave_uint64(y);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	280 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	281
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	282 DOUBLE_INT_BINOP_DECL (+, uint64)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	283 { return y + x; }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	284
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	285 INT_DOUBLE_BINOP_DECL (+, int64)
8169 66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	286 {
8185 69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	287 if (fabs (y) < static_cast<double> (octave_int64::max ()))
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	288 return x + octave_int64 (y);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	289 else
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	290 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	291 // If the number is within the int64 range (the most common case,
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	292 // probably), the above will work as expected. If not, it's more
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	293 // complicated - as long as y is within _twice_ the signed range, the
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	294 // result may still be an integer. An instance of such an operation is
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	295 // 3262 + (1+intmin('int64')) that should yield int64(2*62) + 1. So
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	296 // what we do is to try to convert y/2 and add it twice. Note that if y/2
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	297 // overflows, the result must overflow as well, and that y/2 cannot be a
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	298 // fractional number.
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	299 octave_int64 y2 (y / 2);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	300 return (x + y2) + y2;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	301 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	302 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	303
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	304 DOUBLE_INT_BINOP_DECL (+, int64)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	305 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	306 return y + x;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	307 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	308
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	309 INT_DOUBLE_BINOP_DECL (-, uint64)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	310 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	311 return x + (-y);
8169 66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	312 }
5072 97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	313
8185 69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	314 DOUBLE_INT_BINOP_DECL (-, uint64)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	315 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	316 if (x <= static_cast<double> (octave_uint64::max ()))
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	317 return octave_uint64(x) - y;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	318 else
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	319 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	320 // Again a trick to get the corner cases right. Things like
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	321 // 3263 - intmax('uint64') should produce the correct result, i.e.
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	322 // int64(2**63) + 1.
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	323 const double p2_64 = std::pow (2.0, 64);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	324 if (y.bool_value ())
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	325 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	326 const uint64_t p2_64my = (~y.value ()) + 1; // Equals 2**64 - y
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	327 return octave_uint64 (x - p2_64) + octave_uint64 (p2_64my);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	328 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	329 else
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	330 return octave_uint64 (p2_64);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	331 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	332 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	333
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	334 INT_DOUBLE_BINOP_DECL (-, int64)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	335 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	336 return x + (-y);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	337 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	338
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	339 DOUBLE_INT_BINOP_DECL (-, int64)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	340 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	341 static const bool twosc = (std::numeric_limits<int64_t>::min ()
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	342 < -std::numeric_limits<int64_t>::max ());
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	343 // In case of symmetric integers (not two's complement), this will probably
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	344 // be eliminated at compile time.
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	345 if (twosc && y.value () == std::numeric_limits<int64_t>::min ())
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	346 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	347 return octave_int64 (x + std::pow(2.0, 63));
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	348 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	349 else
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	350 return x + (-y);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	351 }
5072 97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	352
8185 69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	353 // NOTE:
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	354 // Emulated mixed multiplications are tricky due to possible precision loss.
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	355 // Here, after sorting out common cases for speed, we follow the strategy
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	356 // of converting the double number into the form sign * 64-bit integer* 2**exponent,
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	357 // multiply the 64-bit integers to get a 128-bit number, split that number into 32-bit words
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	358 // and form 4 double-valued summands (none of which loases precision), then convert these
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	359 // into integers and sum them. Though it is not immediately obvious, this should work
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	360 // even w.r.t. rounding (none of the summands lose precision).
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	361
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	362 // Multiplies two unsigned 64-bit ints to get a 128-bit number represented
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	363 // as four 32-bit words.
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	364 static void
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	365 umul128 (uint64_t x, uint64_t y, uint32_t w[4])
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	366 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	367 uint64_t lx = static_cast<uint32_t> (x), ux = x >> 32;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	368 uint64_t ly = static_cast<uint32_t> (y), uy = y >> 32;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	369 uint64_t a = lx * ly;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	370 w[0] = a; a >>= 32;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	371 uint64_t uxly = uxly, uylx = uylx;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	372 a += static_cast<uint32_t> (uxly); uxly >>= 32;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	373 a += static_cast<uint32_t> (uylx); uylx >>= 32;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	374 w[1] = a; a >>= 32;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	375 uint64_t uxuy = ux * uy;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	376 a += uxly; a += uylx; a += uxuy;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	377 w[2] = a; a >>= 32;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	378 w[3] = a;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	379 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	380
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	381 // Splits a double into bool sign, unsigned 64-bit mantissa and int exponent
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	382 static void
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	383 dblesplit (double x, bool& sign, uint64_t& mtis, int& exp)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	384 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	385 sign = x < 0; x = fabs (x);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	386 x = frexp (x, &exp);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	387 exp -= 52;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	388 mtis = static_cast<uint64_t> (ldexp (x, 52));
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	389 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	390
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	391 // Gets a double number from a 32-bit unsigned integer mantissa, exponent and sign.
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	392 static double
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	393 dbleget (bool sign, uint32_t mtis, int exp)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	394 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	395 double x = ldexp (static_cast<double> (mtis), exp);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	396 return sign ? -x : x;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	397 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	398
8169 66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	399
8185 69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	400 INT_DOUBLE_BINOP_DECL (*, uint64)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	401 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	402 if (y >= 0 && y < octave_uint64::max () && y == xround (y))
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	403 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	404 return x * octave_uint64 (static_cast<uint64_t> (y));
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	405 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	406 else if (y == 0.5)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	407 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	408 return x / octave_uint64 (static_cast<uint64_t> (2));
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	409 }
9003 0631d397fbe0 replace lo_ieee_isnan by xisnan, add missing includes Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	410 else if (y < 0 \|\| xisnan (x) \|\| xisinf (x))
8185 69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	411 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	412 return octave_uint64 (x.value () * y);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	413 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	414 else
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	415 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	416 bool sign;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	417 uint64_t my;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	418 int e;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	419 dblesplit (y, sign, my, e);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	420 uint32_t w[4];
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	421 umul128 (x.value (), my, w);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	422 octave_uint64 res = octave_uint64::zero;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	423 for (short i = 0; i < 4; i++)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	424 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	425 res += octave_uint64 (dbleget (sign, w[i], e));
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	426 e += 32;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	427 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	428 return res;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	429 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	430 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	431
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	432 DOUBLE_INT_BINOP_DECL (*, uint64)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	433 { return y * x; }
8169 66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	434
8185 69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	435 INT_DOUBLE_BINOP_DECL (*, int64)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	436 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	437 if (fabs (y) < octave_int64::max () && y == xround (y))
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	438 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	439 return x * octave_int64 (static_cast<int64_t> (y));
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	440 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	441 else if (fabs (y) == 0.5)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	442 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	443 return x / octave_int64 (static_cast<uint64_t> (4*y));
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	444 }
9003 0631d397fbe0 replace lo_ieee_isnan by xisnan, add missing includes Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	445 else if (xisnan (x) \|\| xisinf (x))
8185 69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	446 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	447 return octave_int64 (x.value () * y);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	448 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	449 else
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	450 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	451 bool sign;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	452 uint64_t my;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	453 int e;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	454 dblesplit (y, sign, my, e);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	455 uint32_t w[4];
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	456 sign = (sign != (x.value () < 0));
8293 ad5bb02d267a workaround missing std::abs(int64_t) in MSVC Jaroslav Hajek <highegg@gmail.com> parents: 8202 diff changeset	457 umul128 (octave_int_abs (x.value ()), my, w);
8185 69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	458 octave_int64 res = octave_int64::zero;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	459 for (short i = 0; i < 4; i++)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	460 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	461 res += octave_int64 (dbleget (sign, w[i], e));
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	462 e += 32;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	463 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	464 return res;
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	465 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	466 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	467
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	468 DOUBLE_INT_BINOP_DECL (*, int64)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	469 { return y * x; }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	470
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	471 DOUBLE_INT_BINOP_DECL (/, uint64)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	472 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	473 return octave_uint64 (x / static_cast<double> (y));
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	474 }
5072 97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	475
8185 69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	476 DOUBLE_INT_BINOP_DECL (/, int64)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	477 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	478 return octave_int64 (x / static_cast<double> (y));
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	479 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	480
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	481 INT_DOUBLE_BINOP_DECL (/, uint64)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	482 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	483 if (y >= 0 && y < octave_uint64::max () && y == xround (y))
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	484 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	485 return x / octave_uint64 (y);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	486 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	487 else
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	488 return x * (1.0/y);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	489 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	490
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	491 INT_DOUBLE_BINOP_DECL (/, int64)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	492 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	493 if (fabs (y) < octave_int64::max () && y == xround (y))
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	494 {
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	495 return x / octave_int64 (y);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	496 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	497 else
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	498 return x * (1.0/y);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	499 }
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	500
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	501 #define INSTANTIATE_INT64_DOUBLE_CMP_OP0(OP,T1,T2) \
8856 ab4db66e286f workaround gcc 4.3 explicit instantiation bug in octave_int_cmp_op Jaroslav Hajek <highegg@gmail.com> parents: 8797 diff changeset	502 template bool \
ab4db66e286f workaround gcc 4.3 explicit instantiation bug in octave_int_cmp_op Jaroslav Hajek <highegg@gmail.com> parents: 8797 diff changeset	503 octave_int_cmp_op::emulate_mop<octave_int_cmp_op::OP> (T1 x, T2 y)
8185 69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	504
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	505 #define INSTANTIATE_INT64_DOUBLE_CMP_OP(OP) \
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	506 INSTANTIATE_INT64_DOUBLE_CMP_OP0(OP, double, int64_t); \
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	507 INSTANTIATE_INT64_DOUBLE_CMP_OP0(OP, double, uint64_t); \
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	508 INSTANTIATE_INT64_DOUBLE_CMP_OP0(OP, int64_t, double); \
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	509 INSTANTIATE_INT64_DOUBLE_CMP_OP0(OP, uint64_t, double)
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	510
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	511 INSTANTIATE_INT64_DOUBLE_CMP_OP(lt);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	512 INSTANTIATE_INT64_DOUBLE_CMP_OP(le);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	513 INSTANTIATE_INT64_DOUBLE_CMP_OP(gt);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	514 INSTANTIATE_INT64_DOUBLE_CMP_OP(ge);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	515 INSTANTIATE_INT64_DOUBLE_CMP_OP(eq);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	516 INSTANTIATE_INT64_DOUBLE_CMP_OP(ne);
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	517
69c5cce38c29 implement 64-bit arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8169 diff changeset	518 #endif
8169 66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	519
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	520 //template <class T>
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	521 //bool
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	522 //xisnan (const octave_int<T>&)
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	523 //{
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	524 // return false;
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	525 //}
5072 97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	526
8169 66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	527 template <class T>
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	528 octave_int<T>
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	529 pow (const octave_int<T>& a, const octave_int<T>& b)
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	530 {
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	531 octave_int<T> retval;
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	532
8864 4d328b8979c8 simplify integer power Jaroslav Hajek <highegg@gmail.com> parents: 8856 diff changeset	533 octave_int<T> zero = static_cast<T> (0);
4d328b8979c8 simplify integer power Jaroslav Hajek <highegg@gmail.com> parents: 8856 diff changeset	534 octave_int<T> one = static_cast<T> (1);
8169 66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	535
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	536 if (b == zero \|\| a == one)
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	537 retval = one;
8864 4d328b8979c8 simplify integer power Jaroslav Hajek <highegg@gmail.com> parents: 8856 diff changeset	538 else if (b < zero)
8169 66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	539 {
8864 4d328b8979c8 simplify integer power Jaroslav Hajek <highegg@gmail.com> parents: 8856 diff changeset	540 if (a == -one)
8169 66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	541 retval = (b.value () % 2) ? a : one;
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	542 else
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	543 retval = zero;
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	544 }
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	545 else
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	546 {
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	547 octave_int<T> a_val = a;
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	548 T b_val = b; // no need to do saturation on b
5072 97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	549
8169 66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	550 retval = a;
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	551
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	552 b_val -= 1;
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	553
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	554 while (b_val != 0)
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	555 {
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	556 if (b_val & 1)
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	557 retval = retval * a_val;
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	558
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	559 b_val = b_val >> 1;
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	560
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	561 if (b_val)
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	562 a_val = a_val * a_val;
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	563 }
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	564 }
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	565
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	566 return retval;
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	567 }
5072 97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	568
8169 66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	569 template <class T>
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	570 octave_int<T>
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	571 pow (const double& a, const octave_int<T>& b)
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	572 { return octave_int<T> (pow (a, b.double_value ())); }
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	573
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	574 template <class T>
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	575 octave_int<T>
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	576 pow (const octave_int<T>& a, const double& b)
8195 ec3a55bd876b fix integer exponentiation Jaroslav Hajek <highegg@gmail.com> parents: 8192 diff changeset	577 {
ec3a55bd876b fix integer exponentiation Jaroslav Hajek <highegg@gmail.com> parents: 8192 diff changeset	578 return ((b >= 0 && b < std::numeric_limits<T>::digits && b == xround (b))
ec3a55bd876b fix integer exponentiation Jaroslav Hajek <highegg@gmail.com> parents: 8192 diff changeset	579 ? pow (a, octave_int<T> (static_cast<T> (b)))
ec3a55bd876b fix integer exponentiation Jaroslav Hajek <highegg@gmail.com> parents: 8192 diff changeset	580 : octave_int<T> (pow (a.double_value (), b)));
ec3a55bd876b fix integer exponentiation Jaroslav Hajek <highegg@gmail.com> parents: 8192 diff changeset	581 }
5072 97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	582
8169 66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	583 template <class T>
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	584 octave_int<T>
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	585 powf (const float& a, const octave_int<T>& b)
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	586 { return octave_int<T> (pow (a, b.float_value ())); }
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	587
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	588 template <class T>
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	589 octave_int<T>
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	590 powf (const octave_int<T>& a, const float& b)
8195 ec3a55bd876b fix integer exponentiation Jaroslav Hajek <highegg@gmail.com> parents: 8192 diff changeset	591 {
ec3a55bd876b fix integer exponentiation Jaroslav Hajek <highegg@gmail.com> parents: 8192 diff changeset	592 return ((b >= 0 && b < std::numeric_limits<T>::digits && b == xround (b))
ec3a55bd876b fix integer exponentiation Jaroslav Hajek <highegg@gmail.com> parents: 8192 diff changeset	593 ? pow (a, octave_int<T> (static_cast<T> (b)))
8319 c374691576f6 Fix for MSVC compilation Michael Goffioul <michael.goffioul@gmail.com> parents: 8293 diff changeset	594 : octave_int<T> (pow (a.double_value (), static_cast<double> (b))));
8195 ec3a55bd876b fix integer exponentiation Jaroslav Hajek <highegg@gmail.com> parents: 8192 diff changeset	595 }
5072 97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	596
97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	597 #define INSTANTIATE_INTTYPE(T) \
6108 143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 6105 diff changeset	598 template class OCTAVE_API octave_int<T>; \
143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 6105 diff changeset	599 template OCTAVE_API octave_int<T> pow (const octave_int<T>&, const octave_int<T>&); \
8169 66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	600 template OCTAVE_API octave_int<T> pow (const double&, const octave_int<T>&); \
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	601 template OCTAVE_API octave_int<T> pow (const octave_int<T>&, const double&); \
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	602 template OCTAVE_API octave_int<T> powf (const float&, const octave_int<T>&); \
66bc6f9b4f72 rewrite integer arithmetics and conversions Jaroslav Hajek <highegg@gmail.com> parents: 7521 diff changeset	603 template OCTAVE_API octave_int<T> powf (const octave_int<T>&, const float&); \
6108 143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 6105 diff changeset	604 template OCTAVE_API octave_int<T> \
5072 97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	605 bitshift (const octave_int<T>&, int, const octave_int<T>&); \
97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	606
5828 22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5307 diff changeset	607 INSTANTIATE_INTTYPE (int8_t);
22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5307 diff changeset	608 INSTANTIATE_INTTYPE (int16_t);
22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5307 diff changeset	609 INSTANTIATE_INTTYPE (int32_t);
22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5307 diff changeset	610 INSTANTIATE_INTTYPE (int64_t);
5072 97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	611
5828 22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5307 diff changeset	612 INSTANTIATE_INTTYPE (uint8_t);
22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5307 diff changeset	613 INSTANTIATE_INTTYPE (uint16_t);
22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5307 diff changeset	614 INSTANTIATE_INTTYPE (uint32_t);
22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5307 diff changeset	615 INSTANTIATE_INTTYPE (uint64_t);
5072 97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	616
97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	617
8192 9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	618 // Tests follow.
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	619
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	620 /*
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	621
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	622 %!assert(intmax("int64")/intmin("int64"),int64(-1))
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	623 %!assert(intmin("int64")/int64(-1),intmax("int64"))
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	624 %!assert(int64(2**63),intmax("int64"))
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	625 %!test
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	626 %! wstate = warning("query", "Octave:int-convert-overflow");
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	627 %! warning("on", "Octave:int-convert-overflow");
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	628 %! fail("int64(2*63)","warning",".")
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	629 %! warning(wstate.state, "Octave:int-convert-overflow");
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	630 %!assert(uint64(2**64),intmax("uint64"))
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	631 %!test
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	632 %! a = 1.9*2^61; b = uint64(a); b++; assert(b > a)
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	633 %!test
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	634 %! a = -1.9*2^61; b = int64(a); b++; assert(b > a)
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	635 %!test
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	636 %! a = int64(-2*60) + 2; assert(1.25a == (5*a)/4)
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	637 %!test
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	638 %! a = uint64(2*61) + 2; assert(1.25a == (5*a)/4)
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	639 %!assert(int32(2**31+0.5),intmax('int32'))
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	640 %!test
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	641 %! wstate = warning("query", "Octave:int-convert-overflow");
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	642 %! warning("on", "Octave:int-convert-overflow");
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	643 %! fail("int32(2*31+0.5)","warning",".")
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	644 %! warning(wstate.state, "Octave:int-convert-overflow");
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	645 %!assert(int32(-2**31-0.5),intmin('int32'))
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	646 %!test
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	647 %! wstate = warning("query", "Octave:int-convert-overflow");
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	648 %! warning("on", "Octave:int-convert-overflow");
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	649 %! fail("int32(-2*31-0.5)","warning",".")
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	650 %! warning(wstate.state, "Octave:int-convert-overflow");
8195 ec3a55bd876b fix integer exponentiation Jaroslav Hajek <highegg@gmail.com> parents: 8192 diff changeset	651 %!assert((int64(262)+1)1, int64(2**62)+1)
ec3a55bd876b fix integer exponentiation Jaroslav Hajek <highegg@gmail.com> parents: 8192 diff changeset	652 %!assert((int64(230)+1)2, int64(260+231) + 1)
8192 9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	653 */
9a0a66f650b1 tests for int64 arithmetics Jaroslav Hajek <highegg@gmail.com> parents: 8185 diff changeset	654
5072 97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	655 /*
97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	656 ;;; Local Variables: ***
97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	657 ;;; mode: C++ ***
97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	658 ;;; End: ***
97b62f0c1bee [project @ 2004-11-09 05:51:30 by jwe] jwe parents: diff changeset	659 */

Mercurial > octave-nkf

annotate liboctave/oct-inttypes.cc @ 9003:0631d397fbe0