lib/libc/alpha/gen/frexp.c

   1 /*      $NetBSD: frexp.c,v 1.1 1995/02/10 17:50:22 cgd Exp $    */
   2 /* $FreeBSD: src/lib/libc/alpha/gen/frexp.c,v 1.1.1.1.6.1 2000/08/21 21:09:28 jhb Exp $ */
   3
   4 /*
   5  * Copyright (c) 1994, 1995 Carnegie-Mellon University.
   6  * All rights reserved.
   7  *
   8  * Author: Chris G. Demetriou
   9  *
  10  * Permission to use, copy, modify and distribute this software and
  11  * its documentation is hereby granted, provided that both the copyright
  12  * notice and this permission notice appear in all copies of the
  13  * software, derivative works or modified versions, and any portions
  14  * thereof, and that both notices appear in supporting documentation.
  15  *
  16  * CARNEGIE MELLON ALLOWS FREE USE OF THIS SOFTWARE IN ITS "AS IS"
  17  * CONDITION.  CARNEGIE MELLON DISCLAIMS ANY LIABILITY OF ANY KIND
  18  * FOR ANY DAMAGES WHATSOEVER RESULTING FROM THE USE OF THIS SOFTWARE.
  19  *
  20  * Carnegie Mellon requests users of this software to return to
  21  *
  22  *  Software Distribution Coordinator  or  Software.Distribution@CS.CMU.EDU
  23  *  School of Computer Science
  24  *  Carnegie Mellon University
  25  *  Pittsburgh PA 15213-3890
  26  *
  27  * any improvements or extensions that they make and grant Carnegie the
  28  * rights to redistribute these changes.
  29  */
  30
  31 #include <sys/types.h>
  32 #include <machine/ieee.h>
  33 #include <math.h>
  34
  35 double
  36 frexp(value, eptr)
  37         double value;
  38         int *eptr;
  39 {
  40         union doub {
  41                 double v;
  42                 struct ieee_double s;
  43         } u;
  44
  45         if (value) {
  46                 u.v = value;
  47                 *eptr = u.s.dbl_exp - (DBL_EXP_BIAS - 1);
  48                 u.s.dbl_exp = DBL_EXP_BIAS - 1;
  49                 return(u.v);
  50         } else {
  51                 *eptr = 0;
  52                 return((double)0);
  53         }
  54 }