git.gag.com Git - debian/gnuradio/blob - gnuradio-core/src/lib/general/gr_math.cc

   1 /* -*- c++ -*- */
   2 /*
   3  * Copyright 2003 Free Software Foundation, Inc.
   4  *
   5  * This file is part of GNU Radio
   6  *
   7  * GNU Radio is free software; you can redistribute it and/or modify
   8  * it under the terms of the GNU General Public License as published by
   9  * the Free Software Foundation; either version 2, or (at your option)
  10  * any later version.
  11  *
  12  * GNU Radio is distributed in the hope that it will be useful,
  13  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  14  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  15  * GNU General Public License for more details.
  16  *
  17  * You should have received a copy of the GNU General Public License
  18  * along with GNU Radio; see the file COPYING.  If not, write to
  19  * the Free Software Foundation, Inc., 51 Franklin Street,
  20  * Boston, MA 02110-1301, USA.
  21  */
  22
  23 #ifdef HAVE_CONFIG_H
  24 #include <config.h>
  25 #endif
  26
  27 #include <gr_math.h>
  28 #include <math.h>
  29
  30 /*
  31  * Greatest Common Divisor, using Euclid's algorithm.
  32  * [There are faster algorithms.  See Knuth 4.5.2 if you care]
  33  */
  34
  35 long
  36 gr_gcd (long m, long n)
  37 {
  38   if (m < 0)
  39     m = -m;
  40
  41   if (n < 0)
  42     n = -n;
  43
  44   while (n != 0){
  45     long        t = m % n;
  46     m = n;
  47     n = t;
  48   }
  49
  50   return m;
  51 }
  52
  53
  54 /*
  55  * These really need some configure hacking to figure out the right answer.
  56  * As a stop gap, try for a macro, and if not that, then try std::
  57  */
  58
  59 // returns a non-zero value if value is "not-a-number" (NaN), and 0 otherwise
  60
  61 #if defined(isnan) || !defined(CXX_HAS_STD_ISNAN)
  62
  63 int
  64 gr_isnan (double value)
  65 {
  66   return isnan (value);
  67 }
  68
  69 #else
  70
  71 int
  72 gr_isnan (double value)
  73 {
  74   return std::isnan (value);
  75 }
  76
  77 #endif
  78
  79 // returns a non-zero value if the value of x has its sign bit set.
  80 //
  81 // This  is  not  the  same  as `x < 0.0', because IEEE 754 floating point
  82 // allows zero to be signed.  The comparison `-0.0 < 0.0'  is  false,  but
  83 // `gr_signbit (-0.0)' will return a nonzero value.
  84
  85 #ifdef signbit
  86
  87 int
  88 gr_signbit (double x)
  89 {
  90   return signbit (x);
  91 }
  92
  93 #else
  94
  95 int
  96 gr_signbit (double x)
  97 {
  98   return std::signbit (x);
  99 }
 100
 101
 102 #endif