/*
 * ====================================================
 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
 *
 * Developed at SunPro, a Sun Microsystems, Inc. business.
 * Permission to use, copy, modify, and distribute this
 * software is freely granted, provided that this notice
 * is preserved.
 * ====================================================
 */

/*
 * from: @(#)fdlibm.h 5.1 93/09/24
 * $NetBSD: math_private.h,v 1.34 2024/07/17 12:00:13 riastradh Exp $
 */

/* $Id$ */
#ifndef __NBSD_MATH__
#define __NBSD_MATH__

#include <Mw/BaseTypes.h>

typedef union _ieee_double_shape_type {
	double value;
	struct
	{
		MwU32 lsw;
		MwU32 msw;
	} parts;
	struct
	{
		MwU64 w;
	} xparts;
} ieee_double_shape_type;

/* Get two 32-bit integers from a double.  */

#define EXTRACT_WORDS(ix0, ix1, d) \
	do { \
		ieee_double_shape_type ew_u; \
		ew_u.value = (d); \
		if(nbsd_endian() == 'L') { \
			(ix0) = ew_u.parts.msw; \
			(ix1) = ew_u.parts.lsw; \
		} else { \
			(ix0) = ew_u.parts.lsw; \
			(ix1) = ew_u.parts.msw; \
		} \
	} while(0)

/* Get the more significant 32-bit integer from a double.  */

#define GET_HIGH_WORD(i, d) \
	do { \
		ieee_double_shape_type gh_u; \
		gh_u.value = (d); \
		(i)	   = nbsd_endian() == 'L' ? gh_u.parts.msw : gh_u.parts.lsw; \
	} while(0)

/* Get the less significant 32-bit integer from a double.  */

#define GET_LOW_WORD(i, d) \
	do { \
		ieee_double_shape_type gl_u; \
		gl_u.value = (d); \
		(i)	   = nbsd_endian() == 'L' ? gl_u.parts.lsw : gl_u.parts.msw; \
	} while(0)

/* Set the more significant 32 bits of a double from an integer.  */

#define SET_HIGH_WORD(d, v) \
	do { \
		ieee_double_shape_type sh_u; \
		sh_u.value = (d); \
		if(nbsd_endian() == 'L') { \
			sh_u.parts.msw = (v); \
		} else { \
			sh_u.parts.lsw = (v); \
		} \
		(d) = sh_u.value; \
	} while(0)

/* Set the less significant 32 bits of a double from an integer.  */

#define SET_LOW_WORD(d, v) \
	do { \
		ieee_double_shape_type sl_u; \
		sl_u.value = (d); \
		if(nbsd_endian() == 'L') { \
			sl_u.parts.lsw = (v); \
		} else { \
			sl_u.parts.msw = (v); \
		} \
		(d) = sl_u.value; \
	} while(0)

double nbsd_pow(double a, double b);
char   nbsd_endian(void);

#endif