frexp.c revision 1.5 
11.5Sabs/*	$NetBSD: frexp.c,v 1.5 2012/06/25 22:32:43 abs Exp $	*/
21.2Schristos
31.1Sragge/*-
41.1Sragge * Copyright (c) 1991, 1993
51.1Sragge *	The Regents of the University of California.  All rights reserved.
61.1Sragge *
71.1Sragge * Redistribution and use in source and binary forms, with or without
81.1Sragge * modification, are permitted provided that the following conditions
91.1Sragge * are met:
101.1Sragge * 1. Redistributions of source code must retain the above copyright
111.1Sragge *    notice, this list of conditions and the following disclaimer.
121.1Sragge * 2. Redistributions in binary form must reproduce the above copyright
131.1Sragge *    notice, this list of conditions and the following disclaimer in the
141.1Sragge *    documentation and/or other materials provided with the distribution.
151.4Sagc * 3. Neither the name of the University nor the names of its contributors
161.1Sragge *    may be used to endorse or promote products derived from this software
171.1Sragge *    without specific prior written permission.
181.1Sragge *
191.1Sragge * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND
201.1Sragge * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
211.1Sragge * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
221.1Sragge * ARE DISCLAIMED.  IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE
231.1Sragge * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
241.1Sragge * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
251.1Sragge * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
261.1Sragge * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
271.1Sragge * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
281.1Sragge * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
291.1Sragge * SUCH DAMAGE.
301.1Sragge */
311.1Sragge
321.2Schristos#include <sys/cdefs.h>
331.1Sragge#if defined(LIBC_SCCS) && !defined(lint)
341.2Schristos#if 0
351.2Schristosstatic char sccsid[] = "@(#)frexp.c	8.1 (Berkeley) 6/4/93";
361.2Schristos#else
371.5Sabs__RCSID("$NetBSD: frexp.c,v 1.5 2012/06/25 22:32:43 abs Exp $");
381.2Schristos#endif
391.1Sragge#endif /* LIBC_SCCS and not lint */
401.1Sragge
411.1Sragge#include <sys/types.h>
421.1Sragge#include <math.h>
431.1Sragge
441.1Sraggedouble
451.5Sabsfrexp(double value, int *eptr)
461.1Sragge{
471.1Sragge	union {
481.1Sragge                double v;
491.1Sragge                struct {
501.1Sragge			u_int u_mant1 :  7;
511.1Sragge			u_int   u_exp :  8;
521.1Sragge			u_int  u_sign :  1;
531.1Sragge			u_int u_mant2 : 16;
541.1Sragge			u_int u_mant3 : 32;
551.1Sragge                } s;
561.1Sragge        } u;
571.1Sragge
581.1Sragge	if (value) {
591.1Sragge		u.v = value;
601.1Sragge		*eptr = u.s.u_exp - 128;
611.1Sragge		u.s.u_exp = 128;
621.3Smycroft		return (u.v);
631.1Sragge	} else {
641.1Sragge		*eptr = 0;
651.3Smycroft		return (0.0);
661.1Sragge	}
671.1Sragge}
68