modules/base/math/fixed_point/f64_double.c

   1 /*
   2  *  Copyright Droids Corporation, Microb Technology, Eirbot (2005)
   3  *
   4  *  This program is free software; you can redistribute it and/or modify
   5  *  it under the terms of the GNU General Public License as published by
   6  *  the Free Software Foundation; either version 2 of the License, or
   7  *  (at your option) any later version.
   8  *
   9  *  This program is distributed in the hope that it will be useful,
  10  *  but WITHOUT ANY WARRANTY; without even the implied warranty of
  11  *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  12  *  GNU General Public License for more details.
  13  *
  14  *  You should have received a copy of the GNU General Public License
  15  *  along with this program; if not, write to the Free Software
  16  *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
  17  *
  18  *  Revision : $Id: f64_double.c,v 1.5.4.2 2008-05-10 15:06:26 zer0 Exp $
  19  *
  20  */
  21
  22 #include <stdio.h>
  23 #include <f64.h>
  24
  25 /* value of 2^^32 in float */
  26 #define POW2_32F (4294967296.0)
  27
  28
  29 /**************** double <-> f64 */
  30
  31 f64 f64_from_double(double d)
  32 {
  33         f64 f;
  34         f.f64_decimal = (uint32_t) ((d - (int32_t)d) * POW2_32F);
  35         f.f64_integer = (d < 0 && f.f64_decimal != 0 ? (int32_t)d-1 : (int32_t)d) ;
  36
  37         return f;
  38 }
  39
  40 double f64_to_double(f64 f)
  41 {
  42         double d;
  43
  44         d = f.f64_integer;
  45         d += ((double)f.f64_decimal / POW2_32F);
  46
  47         return d;
  48 }