Blame sysdeps/ieee754/dbl-64/ulog.h

Packit 6c4009
/*
Packit 6c4009
 * IBM Accurate Mathematical Library
Packit 6c4009
 * Written by International Business Machines Corp.
Packit 6c4009
 * Copyright (C) 2001-2018 Free Software Foundation, Inc.
Packit 6c4009
 *
Packit 6c4009
 * This program is free software; you can redistribute it and/or modify
Packit 6c4009
 * it under the terms of the GNU Lesser General Public License as published by
Packit 6c4009
 * the Free Software Foundation; either version 2.1 of the License, or
Packit 6c4009
 * (at your option) any later version.
Packit 6c4009
 *
Packit 6c4009
 * This program is distributed in the hope that it will be useful,
Packit 6c4009
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
Packit 6c4009
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
Packit 6c4009
 * GNU Lesser General Public License for more details.
Packit 6c4009
 *
Packit 6c4009
 * You should have received a copy of the GNU Lesser General Public License
Packit 6c4009
 * along with this program; if not, see <http://www.gnu.org/licenses/>.
Packit 6c4009
 */
Packit 6c4009
Packit 6c4009
/******************************************************************/
Packit 6c4009
/*                                                                */
Packit 6c4009
/* MODULE_NAME:ulog.h                                             */
Packit 6c4009
/*                                                                */
Packit 6c4009
/* common data and variables prototype and definition             */
Packit 6c4009
/******************************************************************/
Packit 6c4009
Packit 6c4009
#ifndef ULOG_H
Packit 6c4009
#define ULOG_H
Packit 6c4009
Packit 6c4009
#ifdef BIG_ENDI
Packit 6c4009
  static const number
Packit 6c4009
  /* polynomial I */
Packit 6c4009
/**/ a2             = {{0xbfe00000, 0x0001aa8f} }, /* -0.500... */
Packit 6c4009
/**/ a3             = {{0x3fd55555, 0x55588d2e} }, /*  0.333... */
Packit 6c4009
  /* polynomial II */
Packit 6c4009
/**/ b0             = {{0x3fd55555, 0x55555555} }, /*  0.333... */
Packit 6c4009
/**/ b1             = {{0xbfcfffff, 0xffffffbb} }, /* -0.249... */
Packit 6c4009
/**/ b2             = {{0x3fc99999, 0x9999992f} }, /*  0.199... */
Packit 6c4009
/**/ b3             = {{0xbfc55555, 0x556503fd} }, /* -0.166... */
Packit 6c4009
/**/ b4             = {{0x3fc24924, 0x925b3d62} }, /*  0.142... */
Packit 6c4009
/**/ b5             = {{0xbfbffffe, 0x160472fc} }, /* -0.124... */
Packit 6c4009
/**/ b6             = {{0x3fbc71c5, 0x25db58ac} }, /*  0.111... */
Packit 6c4009
/**/ b7             = {{0xbfb9a4ac, 0x11a2a61c} }, /* -0.100... */
Packit 6c4009
/**/ b8             = {{0x3fb75077, 0x0df2b591} }, /*  0.091... */
Packit 6c4009
  /* constants    */
Packit 6c4009
/**/ sqrt_2         = {{0x3ff6a09e, 0x667f3bcc} }, /* sqrt(2)   */
Packit 6c4009
/**/ h1             = {{0x3fd2e000, 0x00000000} }, /* 151/2**9  */
Packit 6c4009
/**/ h2             = {{0x3f669000, 0x00000000} }, /* 361/2**17 */
Packit 6c4009
/**/ delu           = {{0x3f700000, 0x00000000} }, /* 1/2**8    */
Packit 6c4009
/**/ delv           = {{0x3ef00000, 0x00000000} }, /* 1/2**16   */
Packit 6c4009
/**/ ln2a           = {{0x3fe62e42, 0xfefa3800} }, /* ln(2) 43 bits */
Packit 6c4009
/**/ ln2b           = {{0x3d2ef357, 0x93c76730} }, /* ln(2)-ln2a    */
Packit 6c4009
/**/ two54          = {{0x43500000, 0x00000000} }, /* 2**54         */
Packit 6c4009
/**/ u03            = {{0x3f9eb851, 0xeb851eb8} }; /* 0.03          */
Packit 6c4009
Packit 6c4009
#else
Packit 6c4009
#ifdef LITTLE_ENDI
Packit 6c4009
  static const number
Packit 6c4009
  /* polynomial I */
Packit 6c4009
/**/ a2             = {{0x0001aa8f, 0xbfe00000} }, /* -0.500... */
Packit 6c4009
/**/ a3             = {{0x55588d2e, 0x3fd55555} }, /*  0.333... */
Packit 6c4009
  /* polynomial II */
Packit 6c4009
/**/ b0             = {{0x55555555, 0x3fd55555} }, /*  0.333... */
Packit 6c4009
/**/ b1             = {{0xffffffbb, 0xbfcfffff} }, /* -0.249... */
Packit 6c4009
/**/ b2             = {{0x9999992f, 0x3fc99999} }, /*  0.199... */
Packit 6c4009
/**/ b3             = {{0x556503fd, 0xbfc55555} }, /* -0.166... */
Packit 6c4009
/**/ b4             = {{0x925b3d62, 0x3fc24924} }, /*  0.142... */
Packit 6c4009
/**/ b5             = {{0x160472fc, 0xbfbffffe} }, /* -0.124... */
Packit 6c4009
/**/ b6             = {{0x25db58ac, 0x3fbc71c5} }, /*  0.111... */
Packit 6c4009
/**/ b7             = {{0x11a2a61c, 0xbfb9a4ac} }, /* -0.100... */
Packit 6c4009
/**/ b8             = {{0x0df2b591, 0x3fb75077} }, /*  0.091... */
Packit 6c4009
  /* constants    */
Packit 6c4009
/**/ sqrt_2         = {{0x667f3bcc, 0x3ff6a09e} }, /* sqrt(2)   */
Packit 6c4009
/**/ h1             = {{0x00000000, 0x3fd2e000} }, /* 151/2**9  */
Packit 6c4009
/**/ h2             = {{0x00000000, 0x3f669000} }, /* 361/2**17 */
Packit 6c4009
/**/ delu           = {{0x00000000, 0x3f700000} }, /* 1/2**8    */
Packit 6c4009
/**/ delv           = {{0x00000000, 0x3ef00000} }, /* 1/2**16   */
Packit 6c4009
/**/ ln2a           = {{0xfefa3800, 0x3fe62e42} }, /* ln(2) 43 bits */
Packit 6c4009
/**/ ln2b           = {{0x93c76730, 0x3d2ef357} }, /* ln(2)-ln2a    */
Packit 6c4009
/**/ two54          = {{0x00000000, 0x43500000} }, /* 2**54         */
Packit 6c4009
/**/ u03            = {{0xeb851eb8, 0x3f9eb851} }; /* 0.03          */
Packit 6c4009
Packit 6c4009
#endif
Packit 6c4009
#endif
Packit 6c4009
Packit 6c4009
#define  SQRT_2    sqrt_2.d
Packit 6c4009
#define  DEL_U     delu.d
Packit 6c4009
#define  DEL_V     delv.d
Packit 6c4009
#define  LN2A      ln2a.d
Packit 6c4009
#define  LN2B      ln2b.d
Packit 6c4009
#define  U03       u03.d
Packit 6c4009
Packit 6c4009
#endif