sortix-mirror/libm/arch/i387/machine/npx.h
Jonas 'Sortie' Termansen 5980be9b3c Add Sortix Math Library.
This work is based in part on code from NetBSD libm, libc and kernel.

The library is partly public domain and partly BSD-style licensed.
2013-12-17 14:30:39 +01:00

198 lines
6.6 KiB
C

/* $NetBSD: npx.h,v 1.25 2010/07/31 21:47:54 joerg Exp $ */
/*-
* Copyright (c) 1990 The Regents of the University of California.
* All rights reserved.
*
* This code is derived from software contributed to Berkeley by
* William Jolitz.
*
* Redistribution and use in source and binary forms, with or without
* modification, are permitted provided that the following conditions
* are met:
* 1. Redistributions of source code must retain the above copyright
* notice, this list of conditions and the following disclaimer.
* 2. Redistributions in binary form must reproduce the above copyright
* notice, this list of conditions and the following disclaimer in the
* documentation and/or other materials provided with the distribution.
* 3. Neither the name of the University nor the names of its contributors
* may be used to endorse or promote products derived from this software
* without specific prior written permission.
*
* THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND
* ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
* IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
* ARE DISCLAIMED. IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE
* FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
* DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
* OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
* HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
* OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
* SUCH DAMAGE.
*
* @(#)npx.h 5.3 (Berkeley) 1/18/91
*/
/*
* 287/387 NPX Coprocessor Data Structures and Constants
* W. Jolitz 1/90
*/
#ifndef _I386_NPX_H_
#define _I386_NPX_H_
/* Environment information of floating point unit */
struct env87 {
long en_cw; /* control word (16bits) */
long en_sw; /* status word (16bits) */
long en_tw; /* tag word (16bits) */
long en_fip; /* floating point instruction pointer */
__uint16_t en_fcs; /* floating code segment selector */
__uint16_t en_opcode; /* opcode last executed (11 bits ) */
long en_foo; /* floating operand offset */
long en_fos; /* floating operand segment selector */
};
/* Contents of each floating point accumulator */
struct fpacc87 {
#ifdef dontdef /* too unportable */
__uint32_t fp_mantlo; /* mantissa low (31:0) */
__uint32_t fp_manthi; /* mantissa high (63:32) */
int fp_exp:15; /* exponent */
int fp_sgn:1; /* mantissa sign */
#else
__uint8_t fp_bytes[10];
#endif
};
/* Floating point context */
struct save87 {
struct env87 sv_env; /* floating point control/status */
struct fpacc87 sv_ac[8]; /* accumulator contents, 0-7 */
#ifndef dontdef
__uint32_t sv_ex_sw; /* status word for last exception (was pad) */
__uint32_t sv_ex_tw; /* tag word for last exception (was pad) */
__uint8_t sv_pad[8 * 2 - 2 * 4]; /* bogus historical padding */
#endif
};
/* Environment of FPU/MMX/SSE/SSE2. */
struct envxmm {
/*0*/ __uint16_t en_cw; /* FPU Control Word */
__uint16_t en_sw; /* FPU Status Word */
__uint8_t en_tw; /* FPU Tag Word (abridged) */
__uint8_t en_rsvd0;
__uint16_t en_opcode; /* FPU Opcode */
__uint32_t en_fip; /* FPU Instruction Pointer */
__uint16_t en_fcs; /* FPU IP selector */
__uint16_t en_rsvd1;
/*16*/ __uint32_t en_foo; /* FPU Data pointer */
__uint16_t en_fos; /* FPU Data pointer selector */
__uint16_t en_rsvd2;
__uint32_t en_mxcsr; /* MXCSR Register State */
__uint32_t en_rsvd3;
};
/* FPU regsters in the extended save format. */
struct fpaccxmm {
__uint8_t fp_bytes[10];
__uint8_t fp_rsvd[6];
};
/* SSE/SSE2 registers. */
struct xmmreg {
__uint8_t sse_bytes[16];
};
/* FPU/MMX/SSE/SSE2 context */
struct savexmm {
struct envxmm sv_env; /* control/status context */
struct fpaccxmm sv_ac[8]; /* ST/MM regs */
struct xmmreg sv_xmmregs[8]; /* XMM regs */
__uint8_t sv_rsvd[16 * 14];
/* 512-bytes --- end of hardware portion of save area */
__uint32_t sv_ex_sw; /* saved SW from last exception */
__uint32_t sv_ex_tw; /* saved TW from last exception */
} __attribute__((__aligned__(16)));
union savefpu {
struct save87 sv_87;
struct savexmm sv_xmm;
};
/*
* The i387 defaults to Intel extended precision mode and round to nearest,
* with all exceptions masked.
*/
#define __INITIAL_NPXCW__ 0x037f
/* NetBSD uses IEEE double precision. */
#define __NetBSD_NPXCW__ 0x127f
/* FreeBSD leaves some exceptions unmasked as well. */
#define __FreeBSD_NPXCW__ 0x1272
/* iBCS2 goes a bit further and leaves the underflow exception unmasked. */
#define __iBCS2_NPXCW__ 0x0262
/* Linux just uses the default control word. */
#define __Linux_NPXCW__ 0x037f
/* SVR4 uses the same control word as iBCS2. */
#define __SVR4_NPXCW__ 0x0262
/*
* The default MXCSR value at reset is 0x1f80, IA-32 Instruction
* Set Reference, pg. 3-369.
*/
#define __INITIAL_MXCSR__ 0x1f80
/*
* 80387 control word bits
*/
#define EN_SW_INVOP 0x0001 /* Invalid operation */
#define EN_SW_DENORM 0x0002 /* Denormalized operand */
#define EN_SW_ZERODIV 0x0004 /* Divide by zero */
#define EN_SW_OVERFLOW 0x0008 /* Overflow */
#define EN_SW_UNDERFLOW 0x0010 /* Underflow */
#define EN_SW_PRECLOSS 0x0020 /* Loss of precision */
#define EN_SW_DATACHAIN 0x0080 /* Data chain exception */
#define EN_SW_CTL_PREC 0x0300 /* Precision control */
#define EN_SW_CTL_ROUND 0x0c00 /* Rounding control */
#define EN_SW_CTL_INF 0x1000 /* Infinity control */
/*
* The standard control word from finit is 0x37F, giving:
* round to nearest
* 64-bit precision
* all exceptions masked.
*
* Now we want:
* affine mode (if we decide to support 287's)
* round to nearest
* 53-bit precision
* all exceptions masked.
*
* 64-bit precision often gives bad results with high level languages
* because it makes the results of calculations depend on whether
* intermediate values are stored in memory or in FPU registers.
*
* The iBCS control word has underflow, overflow, zero divide, and invalid
* operation exceptions unmasked. But that causes an unexpected exception
* in the test program 'paranoia' and makes denormals useless (DBL_MIN / 2
* underflows). It doesn't make a lot of sense to trap underflow without
* trapping denormals.
*/
#ifdef _KERNEL
void probeintr(void);
void probetrap(void);
int npx586bug1(int, int);
void npxinit(struct cpu_info *);
void process_xmm_to_s87(const struct savexmm *, struct save87 *);
void process_s87_to_xmm(const struct save87 *, struct savexmm *);
struct lwp;
int npxtrap(struct lwp *);
#endif
#endif /* !_I386_NPX_H_ */