582 lines
		
	
	
		
			22 KiB
		
	
	
	
		
			C
		
	
	
	
	
	
			
		
		
	
	
			582 lines
		
	
	
		
			22 KiB
		
	
	
	
		
			C
		
	
	
	
	
	
| /*
 | |
|  * Linux/PA-RISC Project (http://www.parisc-linux.org/)
 | |
|  *
 | |
|  * Floating-point emulation code
 | |
|  *  Copyright (C) 2001 Hewlett-Packard (Paul Bame) <bame@debian.org>
 | |
|  *
 | |
|  *    This program is free software; you can redistribute it and/or modify
 | |
|  *    it under the terms of the GNU General Public License as published by
 | |
|  *    the Free Software Foundation; either version 2, or (at your option)
 | |
|  *    any later version.
 | |
|  *
 | |
|  *    This program is distributed in the hope that it will be useful,
 | |
|  *    but WITHOUT ANY WARRANTY; without even the implied warranty of
 | |
|  *    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 | |
|  *    GNU General Public License for more details.
 | |
|  *
 | |
|  *    You should have received a copy of the GNU General Public License
 | |
|  *    along with this program; if not, write to the Free Software
 | |
|  *    Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 | |
|  */
 | |
| /*
 | |
|  * BEGIN_DESC
 | |
|  * 
 | |
|  *  File: 
 | |
|  *      @(#)	pa/spmath/float.h		$Revision: 1.1 $
 | |
|  * 
 | |
|  *  Purpose:
 | |
|  *      <<please update with a synopis of the functionality provided by this file>>
 | |
|  * 
 | |
|  *  BE header:  no
 | |
|  *
 | |
|  *  Shipped:  yes
 | |
|  *	/usr/conf/pa/spmath/float.h
 | |
|  *
 | |
|  * END_DESC  
 | |
| */
 | |
| 
 | |
| #ifdef __NO_PA_HDRS
 | |
|     PA header file -- do not include this header file for non-PA builds.
 | |
| #endif
 | |
| 
 | |
| #include "fpbits.h"
 | |
| #include "hppa.h"
 | |
| /*
 | |
|  * Want to pick up the FPU capability flags, not the PDC structures.
 | |
|  * 'LOCORE' isn't really true in this case, but we don't want the C structures
 | |
|  * so it suits our purposes
 | |
|  */
 | |
| #define LOCORE
 | |
| #include "fpu.h"
 | |
| 
 | |
| /*
 | |
|  * Declare the basic structures for the 3 different
 | |
|  * floating-point precisions.
 | |
|  *        
 | |
|  * Single number  
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  * |s|       exp     |               mantissa                      |
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  */
 | |
| #define	Sall(object) (object)
 | |
| #define	Ssign(object) Bitfield_extract( 0,  1,object)
 | |
| #define	Ssignedsign(object) Bitfield_signed_extract( 0,  1,object)
 | |
| #define	Sexponent(object) Bitfield_extract( 1,  8,object)
 | |
| #define	Smantissa(object) Bitfield_mask( 9, 23,object)
 | |
| #define	Ssignaling(object) Bitfield_extract( 9,  1,object)
 | |
| #define	Ssignalingnan(object) Bitfield_extract( 1,  9,object)
 | |
| #define	Shigh2mantissa(object) Bitfield_extract( 9,  2,object)
 | |
| #define	Sexponentmantissa(object) Bitfield_mask( 1, 31,object)
 | |
| #define	Ssignexponent(object) Bitfield_extract( 0,  9,object)
 | |
| #define	Shidden(object) Bitfield_extract( 8,  1,object)
 | |
| #define	Shiddenoverflow(object) Bitfield_extract( 7,  1,object)
 | |
| #define	Shiddenhigh7mantissa(object) Bitfield_extract( 8,  8,object)
 | |
| #define	Shiddenhigh3mantissa(object) Bitfield_extract( 8,  4,object)
 | |
| #define	Slow(object) Bitfield_mask( 31,  1,object)
 | |
| #define	Slow4(object) Bitfield_mask( 28,  4,object)
 | |
| #define	Slow31(object) Bitfield_mask( 1, 31,object)
 | |
| #define	Shigh31(object) Bitfield_extract( 0, 31,object)
 | |
| #define	Ssignedhigh31(object) Bitfield_signed_extract( 0, 31,object)
 | |
| #define	Shigh4(object) Bitfield_extract( 0,  4,object)
 | |
| #define	Sbit24(object) Bitfield_extract( 24,  1,object)
 | |
| #define	Sbit28(object) Bitfield_extract( 28,  1,object)
 | |
| #define	Sbit29(object) Bitfield_extract( 29,  1,object)
 | |
| #define	Sbit30(object) Bitfield_extract( 30,  1,object)
 | |
| #define	Sbit31(object) Bitfield_mask( 31,  1,object)
 | |
| 
 | |
| #define Deposit_ssign(object,value) Bitfield_deposit(value,0,1,object)
 | |
| #define Deposit_sexponent(object,value) Bitfield_deposit(value,1,8,object)
 | |
| #define Deposit_smantissa(object,value) Bitfield_deposit(value,9,23,object)
 | |
| #define Deposit_shigh2mantissa(object,value) Bitfield_deposit(value,9,2,object)
 | |
| #define Deposit_sexponentmantissa(object,value) \
 | |
|     Bitfield_deposit(value,1,31,object)
 | |
| #define Deposit_ssignexponent(object,value) Bitfield_deposit(value,0,9,object)
 | |
| #define Deposit_slow(object,value) Bitfield_deposit(value,31,1,object)
 | |
| #define Deposit_shigh4(object,value) Bitfield_deposit(value,0,4,object)
 | |
| 
 | |
| #define	Is_ssign(object) Bitfield_mask( 0,  1,object)
 | |
| #define	Is_ssignaling(object) Bitfield_mask( 9,  1,object)
 | |
| #define	Is_shidden(object) Bitfield_mask( 8,  1,object)
 | |
| #define	Is_shiddenoverflow(object) Bitfield_mask( 7,  1,object)
 | |
| #define	Is_slow(object) Bitfield_mask( 31,  1,object)
 | |
| #define	Is_sbit24(object) Bitfield_mask( 24,  1,object)
 | |
| #define	Is_sbit28(object) Bitfield_mask( 28,  1,object)
 | |
| #define	Is_sbit29(object) Bitfield_mask( 29,  1,object)
 | |
| #define	Is_sbit30(object) Bitfield_mask( 30,  1,object)
 | |
| #define	Is_sbit31(object) Bitfield_mask( 31,  1,object)
 | |
| 
 | |
| /* 
 | |
|  * Double number.
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  * |s|       exponent      |          mantissa part 1              |
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  *
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  * |                    mantissa part 2                            |
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  */
 | |
| #define Dallp1(object) (object)
 | |
| #define Dsign(object) Bitfield_extract( 0,  1,object)
 | |
| #define Dsignedsign(object) Bitfield_signed_extract( 0,  1,object)
 | |
| #define Dexponent(object) Bitfield_extract( 1,  11,object)
 | |
| #define Dmantissap1(object) Bitfield_mask( 12, 20,object)
 | |
| #define Dsignaling(object) Bitfield_extract( 12,  1,object)
 | |
| #define Dsignalingnan(object) Bitfield_extract( 1,  12,object)
 | |
| #define Dhigh2mantissa(object) Bitfield_extract( 12,  2,object)
 | |
| #define Dexponentmantissap1(object) Bitfield_mask( 1, 31,object)
 | |
| #define Dsignexponent(object) Bitfield_extract( 0, 12,object)
 | |
| #define Dhidden(object) Bitfield_extract( 11,  1,object)
 | |
| #define Dhiddenoverflow(object) Bitfield_extract( 10,  1,object)
 | |
| #define Dhiddenhigh7mantissa(object) Bitfield_extract( 11,  8,object)
 | |
| #define Dhiddenhigh3mantissa(object) Bitfield_extract( 11,  4,object)
 | |
| #define Dlowp1(object) Bitfield_mask( 31,  1,object)
 | |
| #define Dlow31p1(object) Bitfield_mask( 1, 31,object)
 | |
| #define Dhighp1(object) Bitfield_extract( 0,  1,object)
 | |
| #define Dhigh4p1(object) Bitfield_extract( 0,  4,object)
 | |
| #define Dhigh31p1(object) Bitfield_extract( 0, 31,object)
 | |
| #define Dsignedhigh31p1(object) Bitfield_signed_extract( 0, 31,object)
 | |
| #define Dbit3p1(object) Bitfield_extract( 3,  1,object)
 | |
| 
 | |
| #define Deposit_dsign(object,value) Bitfield_deposit(value,0,1,object)
 | |
| #define Deposit_dexponent(object,value) Bitfield_deposit(value,1,11,object)
 | |
| #define Deposit_dmantissap1(object,value) Bitfield_deposit(value,12,20,object)
 | |
| #define Deposit_dhigh2mantissa(object,value) Bitfield_deposit(value,12,2,object)
 | |
| #define Deposit_dexponentmantissap1(object,value) \
 | |
|     Bitfield_deposit(value,1,31,object)
 | |
| #define Deposit_dsignexponent(object,value) Bitfield_deposit(value,0,12,object)
 | |
| #define Deposit_dlowp1(object,value) Bitfield_deposit(value,31,1,object)
 | |
| #define Deposit_dhigh4p1(object,value) Bitfield_deposit(value,0,4,object)
 | |
| 
 | |
| #define Is_dsign(object) Bitfield_mask( 0,  1,object)
 | |
| #define Is_dsignaling(object) Bitfield_mask( 12,  1,object)
 | |
| #define Is_dhidden(object) Bitfield_mask( 11,  1,object)
 | |
| #define Is_dhiddenoverflow(object) Bitfield_mask( 10,  1,object)
 | |
| #define Is_dlowp1(object) Bitfield_mask( 31,  1,object)
 | |
| #define Is_dhighp1(object) Bitfield_mask( 0,  1,object)
 | |
| #define Is_dbit3p1(object) Bitfield_mask( 3,  1,object)
 | |
| 
 | |
| #define Dallp2(object) (object)
 | |
| #define Dmantissap2(object) (object)
 | |
| #define Dlowp2(object) Bitfield_mask( 31,  1,object)
 | |
| #define Dlow4p2(object) Bitfield_mask( 28,  4,object)
 | |
| #define Dlow31p2(object) Bitfield_mask( 1, 31,object)
 | |
| #define Dhighp2(object) Bitfield_extract( 0,  1,object)
 | |
| #define Dhigh31p2(object) Bitfield_extract( 0, 31,object)
 | |
| #define Dbit2p2(object) Bitfield_extract( 2,  1,object)
 | |
| #define Dbit3p2(object) Bitfield_extract( 3,  1,object)
 | |
| #define Dbit21p2(object) Bitfield_extract( 21,  1,object)
 | |
| #define Dbit28p2(object) Bitfield_extract( 28,  1,object)
 | |
| #define Dbit29p2(object) Bitfield_extract( 29,  1,object)
 | |
| #define Dbit30p2(object) Bitfield_extract( 30,  1,object)
 | |
| #define Dbit31p2(object) Bitfield_mask( 31,  1,object)
 | |
| 
 | |
| #define Deposit_dlowp2(object,value) Bitfield_deposit(value,31,1,object)
 | |
| 
 | |
| #define Is_dlowp2(object) Bitfield_mask( 31,  1,object)
 | |
| #define Is_dhighp2(object) Bitfield_mask( 0,  1,object)
 | |
| #define Is_dbit2p2(object) Bitfield_mask( 2,  1,object)
 | |
| #define Is_dbit3p2(object) Bitfield_mask( 3,  1,object)
 | |
| #define Is_dbit21p2(object) Bitfield_mask( 21,  1,object)
 | |
| #define Is_dbit28p2(object) Bitfield_mask( 28,  1,object)
 | |
| #define Is_dbit29p2(object) Bitfield_mask( 29,  1,object)
 | |
| #define Is_dbit30p2(object) Bitfield_mask( 30,  1,object)
 | |
| #define Is_dbit31p2(object) Bitfield_mask( 31,  1,object)
 | |
| 
 | |
| /* 
 | |
|  * Quad number.
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  * |s|          exponent           |      mantissa part 1          |
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  *
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  * |                    mantissa part 2                            |
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  *
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  * |                    mantissa part 3                            |
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  *
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  * |                    mantissa part 4                            |
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  */
 | |
| typedef struct
 | |
|     {
 | |
|     union
 | |
| 	{
 | |
| 	struct { unsigned qallp1; } u_qallp1;
 | |
| /* Not needed for now...
 | |
| 	Bitfield_extract( 0,  1,u_qsign,qsign)
 | |
| 	Bitfield_signed_extract( 0,  1,u_qsignedsign,qsignedsign)
 | |
| 	Bitfield_extract( 1, 15,u_qexponent,qexponent)
 | |
| 	Bitfield_extract(16, 16,u_qmantissap1,qmantissap1)
 | |
| 	Bitfield_extract(16,  1,u_qsignaling,qsignaling)
 | |
| 	Bitfield_extract(1,  16,u_qsignalingnan,qsignalingnan)
 | |
| 	Bitfield_extract(16,  2,u_qhigh2mantissa,qhigh2mantissa)
 | |
| 	Bitfield_extract( 1, 31,u_qexponentmantissap1,qexponentmantissap1)
 | |
| 	Bitfield_extract( 0, 16,u_qsignexponent,qsignexponent)
 | |
| 	Bitfield_extract(15,  1,u_qhidden,qhidden)
 | |
| 	Bitfield_extract(14,  1,u_qhiddenoverflow,qhiddenoverflow)
 | |
| 	Bitfield_extract(15,  8,u_qhiddenhigh7mantissa,qhiddenhigh7mantissa)
 | |
| 	Bitfield_extract(15,  4,u_qhiddenhigh3mantissa,qhiddenhigh3mantissa)
 | |
| 	Bitfield_extract(31,  1,u_qlowp1,qlowp1)
 | |
| 	Bitfield_extract( 1, 31,u_qlow31p1,qlow31p1)
 | |
| 	Bitfield_extract( 0,  1,u_qhighp1,qhighp1)
 | |
| 	Bitfield_extract( 0,  4,u_qhigh4p1,qhigh4p1)
 | |
| 	Bitfield_extract( 0, 31,u_qhigh31p1,qhigh31p1)
 | |
|   */
 | |
| 	} quad_u1;
 | |
|     union
 | |
| 	{
 | |
| 	struct { unsigned qallp2; } u_qallp2;
 | |
|   /* Not needed for now...
 | |
| 	Bitfield_extract(31,  1,u_qlowp2,qlowp2)
 | |
| 	Bitfield_extract( 1, 31,u_qlow31p2,qlow31p2)
 | |
| 	Bitfield_extract( 0,  1,u_qhighp2,qhighp2)
 | |
| 	Bitfield_extract( 0, 31,u_qhigh31p2,qhigh31p2)
 | |
|    */
 | |
| 	} quad_u2;
 | |
|     union
 | |
| 	{
 | |
| 	struct { unsigned qallp3; } u_qallp3;
 | |
|   /* Not needed for now...
 | |
| 	Bitfield_extract(31,  1,u_qlowp3,qlowp3)
 | |
| 	Bitfield_extract( 1, 31,u_qlow31p3,qlow31p3)
 | |
| 	Bitfield_extract( 0,  1,u_qhighp3,qhighp3)
 | |
| 	Bitfield_extract( 0, 31,u_qhigh31p3,qhigh31p3)
 | |
|    */ 
 | |
| 	} quad_u3;
 | |
|     union
 | |
| 	{
 | |
| 	struct { unsigned qallp4; } u_qallp4;
 | |
|     /* Not need for now...
 | |
| 	Bitfield_extract(31,  1,u_qlowp4,qlowp4)
 | |
| 	Bitfield_extract( 1, 31,u_qlow31p4,qlow31p4)
 | |
| 	Bitfield_extract( 0,  1,u_qhighp4,qhighp4)
 | |
| 	Bitfield_extract( 0, 31,u_qhigh31p4,qhigh31p4)
 | |
|      */
 | |
| 	} quad_u4;
 | |
|     } quad_floating_point;
 | |
| 
 | |
| /* Extension - An additional structure to hold the guard, round and
 | |
|  *             sticky bits during computations.
 | |
|  */
 | |
| #define Extall(object) (object)
 | |
| #define Extsign(object) Bitfield_extract( 0,  1,object)
 | |
| #define Exthigh31(object) Bitfield_extract( 0, 31,object)
 | |
| #define Extlow31(object) Bitfield_extract( 1, 31,object)
 | |
| #define Extlow(object) Bitfield_extract( 31,  1,object)
 | |
| 
 | |
| /*
 | |
|  * Single extended - The upper word is just like single precision,
 | |
|  *                 but one additional word of mantissa is needed.
 | |
|  */
 | |
| #define Sextallp1(object) (object)
 | |
| #define Sextallp2(object) (object)
 | |
| #define Sextlowp1(object) Bitfield_extract( 31,  1,object)
 | |
| #define Sexthighp2(object) Bitfield_extract( 0,  1,object)
 | |
| #define Sextlow31p2(object) Bitfield_extract( 1, 31,object)
 | |
| #define Sexthiddenoverflow(object) Bitfield_extract( 4,  1,object)
 | |
| #define Is_sexthiddenoverflow(object) Bitfield_mask( 4,  1,object)
 | |
| 
 | |
| /*
 | |
|  * Double extended - The upper two words are just like double precision,
 | |
|  *		     but two additional words of mantissa are needed.
 | |
|  */
 | |
| #define Dextallp1(object) (object)
 | |
| #define Dextallp2(object) (object)
 | |
| #define Dextallp3(object) (object)
 | |
| #define Dextallp4(object) (object)
 | |
| #define Dextlowp2(object) Bitfield_extract( 31,  1,object)
 | |
| #define Dexthighp3(object) Bitfield_extract( 0,  1,object)
 | |
| #define Dextlow31p3(object) Bitfield_extract( 1, 31,object)
 | |
| #define Dexthiddenoverflow(object) Bitfield_extract( 10,  1,object)
 | |
| #define Is_dexthiddenoverflow(object) Bitfield_mask( 10,  1,object)
 | |
| #define Deposit_dextlowp4(object,value) Bitfield_deposit(value,31,1,object)
 | |
| 
 | |
| /*
 | |
|  * Declare the basic structures for the 3 different
 | |
|  * fixed-point precisions.
 | |
|  *        
 | |
|  * Single number  
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  * |s|                    integer                                  |
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  */
 | |
| typedef int sgl_integer;
 | |
| 
 | |
| /* 
 | |
|  * Double number.
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  * |s|                     high integer                            |
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  *
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  * |                       low integer                             |
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  */
 | |
| struct dint {
 | |
|         int  wd0;
 | |
|         unsigned int wd1;
 | |
| };
 | |
| 
 | |
| struct dblwd {
 | |
|         unsigned int wd0;
 | |
|         unsigned int wd1;
 | |
| };
 | |
| 
 | |
| /* 
 | |
|  * Quad number.
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  * |s|                  integer part1                              |
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  *
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  * |                    integer part 2                             |
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  *
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  * |                    integer part 3                             |
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  *
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  * |                    integer part 4                             |
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  */
 | |
| 
 | |
| struct quadwd {
 | |
|         int  wd0;
 | |
|         unsigned int wd1;
 | |
|         unsigned int wd2;
 | |
|         unsigned int wd3;
 | |
| };
 | |
| 
 | |
| typedef struct quadwd quad_integer;
 | |
| 
 | |
| 
 | |
| /* useful typedefs */
 | |
| typedef unsigned int sgl_floating_point;
 | |
| typedef struct dblwd dbl_floating_point;
 | |
| typedef struct dint dbl_integer;
 | |
| typedef struct dblwd dbl_unsigned;
 | |
| 
 | |
| /* 
 | |
|  * Define the different precisions' parameters.
 | |
|  */
 | |
| #define SGL_BITLENGTH 32
 | |
| #define SGL_EMAX 127
 | |
| #define SGL_EMIN (-126)
 | |
| #define SGL_BIAS 127
 | |
| #define SGL_WRAP 192
 | |
| #define SGL_INFINITY_EXPONENT (SGL_EMAX+SGL_BIAS+1)
 | |
| #define SGL_THRESHOLD 32
 | |
| #define SGL_EXP_LENGTH 8
 | |
| #define SGL_P 24
 | |
| 
 | |
| #define DBL_BITLENGTH 64
 | |
| #define DBL_EMAX 1023
 | |
| #define DBL_EMIN (-1022)
 | |
| #define DBL_BIAS 1023
 | |
| #define DBL_WRAP 1536
 | |
| #define DBL_INFINITY_EXPONENT (DBL_EMAX+DBL_BIAS+1)
 | |
| #define DBL_THRESHOLD 64
 | |
| #define DBL_EXP_LENGTH 11
 | |
| #define DBL_P 53
 | |
| 
 | |
| #define QUAD_BITLENGTH 128
 | |
| #define QUAD_EMAX 16383
 | |
| #define QUAD_EMIN (-16382)
 | |
| #define QUAD_BIAS 16383
 | |
| #define QUAD_WRAP 24576
 | |
| #define QUAD_INFINITY_EXPONENT (QUAD_EMAX+QUAD_BIAS+1)
 | |
| #define QUAD_P 113
 | |
| 
 | |
| /* Boolean Values etc. */
 | |
| #define FALSE 0
 | |
| #define TRUE (!FALSE)
 | |
| #define NOT !
 | |
| #define XOR ^
 | |
| 
 | |
| /* other constants */
 | |
| #undef NULL
 | |
| #define NULL 0
 | |
| #define NIL 0
 | |
| #define SGL 0
 | |
| #define DBL 1
 | |
| #define BADFMT 2
 | |
| #define QUAD 3
 | |
| 
 | |
| 
 | |
| /* Types */
 | |
| typedef int boolean;
 | |
| typedef int FORMAT;
 | |
| typedef int VOID;
 | |
| 
 | |
| 
 | |
| /* Declare status register equivalent to FPUs architecture.
 | |
|  *
 | |
|  *  0 1 2 3 4 5 6 7 8 910 1 2 3 4 5 6 7 8 920 1 2 3 4 5 6 7 8 930 1
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  * |V|Z|O|U|I|C|  rsv  |  model    | version |RM |rsv|T|r|V|Z|O|U|I|
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  */
 | |
| #define Cbit(object) Bitfield_extract( 5, 1,object)
 | |
| #define Tbit(object) Bitfield_extract( 25, 1,object)
 | |
| #define Roundingmode(object) Bitfield_extract( 21, 2,object)
 | |
| #define Invalidtrap(object) Bitfield_extract( 27, 1,object)
 | |
| #define Divisionbyzerotrap(object) Bitfield_extract( 28, 1,object)
 | |
| #define Overflowtrap(object) Bitfield_extract( 29, 1,object)
 | |
| #define Underflowtrap(object) Bitfield_extract( 30, 1,object)
 | |
| #define Inexacttrap(object) Bitfield_extract( 31, 1,object)
 | |
| #define Invalidflag(object) Bitfield_extract( 0, 1,object)
 | |
| #define Divisionbyzeroflag(object) Bitfield_extract( 1, 1,object)
 | |
| #define Overflowflag(object) Bitfield_extract( 2, 1,object)
 | |
| #define Underflowflag(object) Bitfield_extract( 3, 1,object)
 | |
| #define Inexactflag(object) Bitfield_extract( 4, 1,object)
 | |
| #define Allflags(object) Bitfield_extract( 0, 5,object)
 | |
| 
 | |
| /* Definitions relevant to the status register */
 | |
| 
 | |
| /* Rounding Modes */
 | |
| #define ROUNDNEAREST 0
 | |
| #define ROUNDZERO    1
 | |
| #define ROUNDPLUS    2
 | |
| #define ROUNDMINUS   3
 | |
| 
 | |
| /* Exceptions */
 | |
| #define NOEXCEPTION		0x0
 | |
| #define INVALIDEXCEPTION	0x20
 | |
| #define DIVISIONBYZEROEXCEPTION	0x10
 | |
| #define OVERFLOWEXCEPTION	0x08
 | |
| #define UNDERFLOWEXCEPTION	0x04
 | |
| #define INEXACTEXCEPTION	0x02
 | |
| #define UNIMPLEMENTEDEXCEPTION	0x01
 | |
| 
 | |
| /* New exceptions for the 2E Opcode */
 | |
| #define OPC_2E_INVALIDEXCEPTION     0x30
 | |
| #define OPC_2E_OVERFLOWEXCEPTION    0x18
 | |
| #define OPC_2E_UNDERFLOWEXCEPTION   0x0c
 | |
| #define OPC_2E_INEXACTEXCEPTION     0x12
 | |
| 
 | |
| /* Declare exception registers equivalent to FPUs architecture 
 | |
|  *
 | |
|  *  0 1 2 3 4 5 6 7 8 910 1 2 3 4 5 6 7 8 920 1 2 3 4 5 6 7 8 930 1
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  * |excepttype |  r1     | r2/ext  |  operation  |parm |n| t/cond  |
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  */
 | |
| #define Allexception(object) (object)
 | |
| #define Exceptiontype(object) Bitfield_extract( 0, 6,object)
 | |
| #define Instructionfield(object) Bitfield_mask( 6,26,object)
 | |
| #define Parmfield(object) Bitfield_extract( 23, 3,object)
 | |
| #define Rabit(object) Bitfield_extract( 24, 1,object)
 | |
| #define Ibit(object) Bitfield_extract( 25, 1,object)
 | |
| 
 | |
| #define Set_exceptiontype(object,value) Bitfield_deposit(value, 0, 6,object)
 | |
| #define Set_parmfield(object,value) Bitfield_deposit(value, 23, 3,object)
 | |
| #define Set_exceptiontype_and_instr_field(exception,instruction,object) \
 | |
|     object = exception << 26 | instruction
 | |
| 
 | |
| /* Declare the condition field
 | |
|  *
 | |
|  *  0 1 2 3 4 5 6 7 8 910 1 2 3 4 5 6 7 8 920 1 2 3 4 5 6 7 8 930 1
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  * |                                                     |G|L|E|U|X|
 | |
|  * +-------+-------+-------+-------+-------+-------+-------+-------+
 | |
|  */
 | |
| #define Greaterthanbit(object) Bitfield_extract( 27, 1,object)
 | |
| #define Lessthanbit(object) Bitfield_extract( 28, 1,object)
 | |
| #define Equalbit(object) Bitfield_extract( 29, 1,object)
 | |
| #define Unorderedbit(object) Bitfield_extract( 30, 1,object)
 | |
| #define Exceptionbit(object) Bitfield_extract( 31, 1,object)
 | |
| 
 | |
| /* An alias name for the status register */
 | |
| #define Fpustatus_register (*status)
 | |
| 
 | |
| /**************************************************
 | |
|  * Status register referencing and manipulation.  *
 | |
|  **************************************************/
 | |
| 
 | |
| /* Rounding mode */
 | |
| #define Rounding_mode()  Roundingmode(Fpustatus_register)
 | |
| #define Is_rounding_mode(rmode) \
 | |
|     (Roundingmode(Fpustatus_register) == rmode)
 | |
| #define Set_rounding_mode(value) \
 | |
|     Bitfield_deposit(value,21,2,Fpustatus_register)
 | |
| 
 | |
| /* Boolean testing of the trap enable bits */
 | |
| #define Is_invalidtrap_enabled() Invalidtrap(Fpustatus_register)
 | |
| #define Is_divisionbyzerotrap_enabled() Divisionbyzerotrap(Fpustatus_register)
 | |
| #define Is_overflowtrap_enabled() Overflowtrap(Fpustatus_register)
 | |
| #define Is_underflowtrap_enabled() Underflowtrap(Fpustatus_register)
 | |
| #define Is_inexacttrap_enabled() Inexacttrap(Fpustatus_register)
 | |
| 
 | |
| /* Set the indicated flags in the status register */
 | |
| #define Set_invalidflag() Bitfield_deposit(1,0,1,Fpustatus_register)
 | |
| #define Set_divisionbyzeroflag() Bitfield_deposit(1,1,1,Fpustatus_register)
 | |
| #define Set_overflowflag() Bitfield_deposit(1,2,1,Fpustatus_register)
 | |
| #define Set_underflowflag() Bitfield_deposit(1,3,1,Fpustatus_register)
 | |
| #define Set_inexactflag() Bitfield_deposit(1,4,1,Fpustatus_register)
 | |
| 
 | |
| #define Clear_all_flags() Bitfield_deposit(0,0,5,Fpustatus_register)
 | |
| 
 | |
| /* Manipulate the trap and condition code bits (tbit and cbit) */
 | |
| #define Set_tbit() Bitfield_deposit(1,25,1,Fpustatus_register)
 | |
| #define Clear_tbit() Bitfield_deposit(0,25,1,Fpustatus_register)
 | |
| #define Is_tbit_set() Tbit(Fpustatus_register)
 | |
| #define Is_cbit_set() Cbit(Fpustatus_register)
 | |
| 
 | |
| #define Set_status_cbit(value)  \
 | |
|         Bitfield_deposit(value,5,1,Fpustatus_register)
 | |
| 
 | |
| /*******************************
 | |
|  * Condition field referencing *
 | |
|  *******************************/
 | |
| #define Unordered(cond) Unorderedbit(cond)
 | |
| #define Equal(cond) Equalbit(cond)
 | |
| #define Lessthan(cond) Lessthanbit(cond)
 | |
| #define Greaterthan(cond) Greaterthanbit(cond)
 | |
| #define Exception(cond) Exceptionbit(cond)
 | |
| 
 | |
| 
 | |
| /* Defines for the extension */
 | |
| #define Ext_isone_sign(extent) (Extsign(extent))
 | |
| #define Ext_isnotzero(extent) \
 | |
|     (Extall(extent))
 | |
| #define Ext_isnotzero_lower(extent) \
 | |
|     (Extlow31(extent))
 | |
| #define Ext_leftshiftby1(extent) \
 | |
|     Extall(extent) <<= 1
 | |
| #define Ext_negate(extent) \
 | |
|     (int )Extall(extent) = 0 - (int )Extall(extent)
 | |
| #define Ext_setone_low(extent) Bitfield_deposit(1,31,1,extent)
 | |
| #define Ext_setzero(extent) Extall(extent) = 0
 | |
| 
 | |
| typedef int operation;
 | |
| 
 | |
| /* error messages */
 | |
| 
 | |
| #define		NONE		0
 | |
| #define		UNDEFFPINST	1
 | |
| 
 | |
| /* Function definitions: opcode, opclass */
 | |
| #define FTEST	(1<<2) | 0
 | |
| #define FCPY	(2<<2) | 0
 | |
| #define FABS	(3<<2) | 0
 | |
| #define FSQRT   (4<<2) | 0
 | |
| #define FRND    (5<<2) | 0
 | |
| 
 | |
| #define FCNVFF	(0<<2) | 1
 | |
| #define FCNVXF	(1<<2) | 1
 | |
| #define FCNVFX	(2<<2) | 1
 | |
| #define FCNVFXT	(3<<2) | 1
 | |
| 
 | |
| #define FCMP    (0<<2) | 2
 | |
| 
 | |
| #define FADD	(0<<2) | 3
 | |
| #define FSUB	(1<<2) | 3
 | |
| #define FMPY	(2<<2) | 3
 | |
| #define FDIV	(3<<2) | 3
 | |
| #define FREM	(4<<2) | 3
 | |
| 
 | 
