/* SPDX-License-Identifier: GPL-2.0-or-later */ /* * Linux/PA-RISC Project (http://www.parisc-linux.org/) * * Floating-point emulation code * Copyright (C) 2001 Hewlett-Packard (Paul Bame) <bame@debian.org> */ /* * 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