^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 1) /* SPDX-License-Identifier: GPL-2.0-or-later */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 2) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 3) * Linux/PA-RISC Project (http://www.parisc-linux.org/)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 4) *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 5) * Floating-point emulation code
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 6) * Copyright (C) 2001 Hewlett-Packard (Paul Bame) <bame@debian.org>
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 7) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 8) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 9) * BEGIN_DESC
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 10) *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 11) * File:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 12) * @(#) pa/spmath/float.h $Revision: 1.1 $
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 13) *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 14) * Purpose:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 15) * <<please update with a synopis of the functionality provided by this file>>
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 16) *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 17) * BE header: no
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 18) *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 19) * Shipped: yes
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 20) * /usr/conf/pa/spmath/float.h
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 21) *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 22) * END_DESC
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 23) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 24)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 25) #ifdef __NO_PA_HDRS
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 26) PA header file -- do not include this header file for non-PA builds.
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 27) #endif
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 28)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 29) #include "fpbits.h"
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 30) #include "hppa.h"
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 31) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 32) * Want to pick up the FPU capability flags, not the PDC structures.
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 33) * 'LOCORE' isn't really true in this case, but we don't want the C structures
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 34) * so it suits our purposes
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 35) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 36) #define LOCORE
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 37) #include "fpu.h"
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 38)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 39) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 40) * Declare the basic structures for the 3 different
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 41) * floating-point precisions.
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 42) *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 43) * Single number
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 44) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 45) * |s| exp | mantissa |
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 46) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 47) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 48) #define Sall(object) (object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 49) #define Ssign(object) Bitfield_extract( 0, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 50) #define Ssignedsign(object) Bitfield_signed_extract( 0, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 51) #define Sexponent(object) Bitfield_extract( 1, 8,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 52) #define Smantissa(object) Bitfield_mask( 9, 23,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 53) #define Ssignaling(object) Bitfield_extract( 9, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 54) #define Ssignalingnan(object) Bitfield_extract( 1, 9,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 55) #define Shigh2mantissa(object) Bitfield_extract( 9, 2,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 56) #define Sexponentmantissa(object) Bitfield_mask( 1, 31,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 57) #define Ssignexponent(object) Bitfield_extract( 0, 9,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 58) #define Shidden(object) Bitfield_extract( 8, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 59) #define Shiddenoverflow(object) Bitfield_extract( 7, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 60) #define Shiddenhigh7mantissa(object) Bitfield_extract( 8, 8,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 61) #define Shiddenhigh3mantissa(object) Bitfield_extract( 8, 4,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 62) #define Slow(object) Bitfield_mask( 31, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 63) #define Slow4(object) Bitfield_mask( 28, 4,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 64) #define Slow31(object) Bitfield_mask( 1, 31,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 65) #define Shigh31(object) Bitfield_extract( 0, 31,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 66) #define Ssignedhigh31(object) Bitfield_signed_extract( 0, 31,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 67) #define Shigh4(object) Bitfield_extract( 0, 4,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 68) #define Sbit24(object) Bitfield_extract( 24, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 69) #define Sbit28(object) Bitfield_extract( 28, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 70) #define Sbit29(object) Bitfield_extract( 29, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 71) #define Sbit30(object) Bitfield_extract( 30, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 72) #define Sbit31(object) Bitfield_mask( 31, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 73)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 74) #define Deposit_ssign(object,value) Bitfield_deposit(value,0,1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 75) #define Deposit_sexponent(object,value) Bitfield_deposit(value,1,8,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 76) #define Deposit_smantissa(object,value) Bitfield_deposit(value,9,23,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 77) #define Deposit_shigh2mantissa(object,value) Bitfield_deposit(value,9,2,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 78) #define Deposit_sexponentmantissa(object,value) \
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 79) Bitfield_deposit(value,1,31,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 80) #define Deposit_ssignexponent(object,value) Bitfield_deposit(value,0,9,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 81) #define Deposit_slow(object,value) Bitfield_deposit(value,31,1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 82) #define Deposit_shigh4(object,value) Bitfield_deposit(value,0,4,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 83)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 84) #define Is_ssign(object) Bitfield_mask( 0, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 85) #define Is_ssignaling(object) Bitfield_mask( 9, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 86) #define Is_shidden(object) Bitfield_mask( 8, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 87) #define Is_shiddenoverflow(object) Bitfield_mask( 7, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 88) #define Is_slow(object) Bitfield_mask( 31, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 89) #define Is_sbit24(object) Bitfield_mask( 24, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 90) #define Is_sbit28(object) Bitfield_mask( 28, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 91) #define Is_sbit29(object) Bitfield_mask( 29, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 92) #define Is_sbit30(object) Bitfield_mask( 30, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 93) #define Is_sbit31(object) Bitfield_mask( 31, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 94)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 95) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 96) * Double number.
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 97) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 98) * |s| exponent | mantissa part 1 |
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 99) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 100) *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 101) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 102) * | mantissa part 2 |
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 103) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 104) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 105) #define Dallp1(object) (object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 106) #define Dsign(object) Bitfield_extract( 0, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 107) #define Dsignedsign(object) Bitfield_signed_extract( 0, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 108) #define Dexponent(object) Bitfield_extract( 1, 11,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 109) #define Dmantissap1(object) Bitfield_mask( 12, 20,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 110) #define Dsignaling(object) Bitfield_extract( 12, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 111) #define Dsignalingnan(object) Bitfield_extract( 1, 12,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 112) #define Dhigh2mantissa(object) Bitfield_extract( 12, 2,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 113) #define Dexponentmantissap1(object) Bitfield_mask( 1, 31,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 114) #define Dsignexponent(object) Bitfield_extract( 0, 12,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 115) #define Dhidden(object) Bitfield_extract( 11, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 116) #define Dhiddenoverflow(object) Bitfield_extract( 10, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 117) #define Dhiddenhigh7mantissa(object) Bitfield_extract( 11, 8,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 118) #define Dhiddenhigh3mantissa(object) Bitfield_extract( 11, 4,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 119) #define Dlowp1(object) Bitfield_mask( 31, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 120) #define Dlow31p1(object) Bitfield_mask( 1, 31,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 121) #define Dhighp1(object) Bitfield_extract( 0, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 122) #define Dhigh4p1(object) Bitfield_extract( 0, 4,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 123) #define Dhigh31p1(object) Bitfield_extract( 0, 31,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 124) #define Dsignedhigh31p1(object) Bitfield_signed_extract( 0, 31,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 125) #define Dbit3p1(object) Bitfield_extract( 3, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 126)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 127) #define Deposit_dsign(object,value) Bitfield_deposit(value,0,1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 128) #define Deposit_dexponent(object,value) Bitfield_deposit(value,1,11,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 129) #define Deposit_dmantissap1(object,value) Bitfield_deposit(value,12,20,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 130) #define Deposit_dhigh2mantissa(object,value) Bitfield_deposit(value,12,2,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 131) #define Deposit_dexponentmantissap1(object,value) \
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 132) Bitfield_deposit(value,1,31,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 133) #define Deposit_dsignexponent(object,value) Bitfield_deposit(value,0,12,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 134) #define Deposit_dlowp1(object,value) Bitfield_deposit(value,31,1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 135) #define Deposit_dhigh4p1(object,value) Bitfield_deposit(value,0,4,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 136)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 137) #define Is_dsign(object) Bitfield_mask( 0, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 138) #define Is_dsignaling(object) Bitfield_mask( 12, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 139) #define Is_dhidden(object) Bitfield_mask( 11, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 140) #define Is_dhiddenoverflow(object) Bitfield_mask( 10, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 141) #define Is_dlowp1(object) Bitfield_mask( 31, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 142) #define Is_dhighp1(object) Bitfield_mask( 0, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 143) #define Is_dbit3p1(object) Bitfield_mask( 3, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 144)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 145) #define Dallp2(object) (object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 146) #define Dmantissap2(object) (object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 147) #define Dlowp2(object) Bitfield_mask( 31, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 148) #define Dlow4p2(object) Bitfield_mask( 28, 4,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 149) #define Dlow31p2(object) Bitfield_mask( 1, 31,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 150) #define Dhighp2(object) Bitfield_extract( 0, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 151) #define Dhigh31p2(object) Bitfield_extract( 0, 31,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 152) #define Dbit2p2(object) Bitfield_extract( 2, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 153) #define Dbit3p2(object) Bitfield_extract( 3, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 154) #define Dbit21p2(object) Bitfield_extract( 21, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 155) #define Dbit28p2(object) Bitfield_extract( 28, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 156) #define Dbit29p2(object) Bitfield_extract( 29, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 157) #define Dbit30p2(object) Bitfield_extract( 30, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 158) #define Dbit31p2(object) Bitfield_mask( 31, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 159)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 160) #define Deposit_dlowp2(object,value) Bitfield_deposit(value,31,1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 161)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 162) #define Is_dlowp2(object) Bitfield_mask( 31, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 163) #define Is_dhighp2(object) Bitfield_mask( 0, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 164) #define Is_dbit2p2(object) Bitfield_mask( 2, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 165) #define Is_dbit3p2(object) Bitfield_mask( 3, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 166) #define Is_dbit21p2(object) Bitfield_mask( 21, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 167) #define Is_dbit28p2(object) Bitfield_mask( 28, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 168) #define Is_dbit29p2(object) Bitfield_mask( 29, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 169) #define Is_dbit30p2(object) Bitfield_mask( 30, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 170) #define Is_dbit31p2(object) Bitfield_mask( 31, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 171)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 172) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 173) * Quad number.
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 174) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 175) * |s| exponent | mantissa part 1 |
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 176) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 177) *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 178) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 179) * | mantissa part 2 |
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 180) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 181) *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 182) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 183) * | mantissa part 3 |
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 184) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 185) *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 186) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 187) * | mantissa part 4 |
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 188) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 189) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 190) typedef struct
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 191) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 192) union
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 193) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 194) struct { unsigned qallp1; } u_qallp1;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 195) /* Not needed for now...
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 196) Bitfield_extract( 0, 1,u_qsign,qsign)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 197) Bitfield_signed_extract( 0, 1,u_qsignedsign,qsignedsign)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 198) Bitfield_extract( 1, 15,u_qexponent,qexponent)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 199) Bitfield_extract(16, 16,u_qmantissap1,qmantissap1)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 200) Bitfield_extract(16, 1,u_qsignaling,qsignaling)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 201) Bitfield_extract(1, 16,u_qsignalingnan,qsignalingnan)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 202) Bitfield_extract(16, 2,u_qhigh2mantissa,qhigh2mantissa)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 203) Bitfield_extract( 1, 31,u_qexponentmantissap1,qexponentmantissap1)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 204) Bitfield_extract( 0, 16,u_qsignexponent,qsignexponent)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 205) Bitfield_extract(15, 1,u_qhidden,qhidden)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 206) Bitfield_extract(14, 1,u_qhiddenoverflow,qhiddenoverflow)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 207) Bitfield_extract(15, 8,u_qhiddenhigh7mantissa,qhiddenhigh7mantissa)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 208) Bitfield_extract(15, 4,u_qhiddenhigh3mantissa,qhiddenhigh3mantissa)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 209) Bitfield_extract(31, 1,u_qlowp1,qlowp1)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 210) Bitfield_extract( 1, 31,u_qlow31p1,qlow31p1)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 211) Bitfield_extract( 0, 1,u_qhighp1,qhighp1)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 212) Bitfield_extract( 0, 4,u_qhigh4p1,qhigh4p1)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 213) Bitfield_extract( 0, 31,u_qhigh31p1,qhigh31p1)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 214) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 215) } quad_u1;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 216) union
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 217) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 218) struct { unsigned qallp2; } u_qallp2;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 219) /* Not needed for now...
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 220) Bitfield_extract(31, 1,u_qlowp2,qlowp2)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 221) Bitfield_extract( 1, 31,u_qlow31p2,qlow31p2)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 222) Bitfield_extract( 0, 1,u_qhighp2,qhighp2)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 223) Bitfield_extract( 0, 31,u_qhigh31p2,qhigh31p2)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 224) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 225) } quad_u2;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 226) union
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 227) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 228) struct { unsigned qallp3; } u_qallp3;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 229) /* Not needed for now...
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 230) Bitfield_extract(31, 1,u_qlowp3,qlowp3)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 231) Bitfield_extract( 1, 31,u_qlow31p3,qlow31p3)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 232) Bitfield_extract( 0, 1,u_qhighp3,qhighp3)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 233) Bitfield_extract( 0, 31,u_qhigh31p3,qhigh31p3)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 234) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 235) } quad_u3;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 236) union
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 237) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 238) struct { unsigned qallp4; } u_qallp4;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 239) /* Not need for now...
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 240) Bitfield_extract(31, 1,u_qlowp4,qlowp4)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 241) Bitfield_extract( 1, 31,u_qlow31p4,qlow31p4)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 242) Bitfield_extract( 0, 1,u_qhighp4,qhighp4)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 243) Bitfield_extract( 0, 31,u_qhigh31p4,qhigh31p4)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 244) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 245) } quad_u4;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 246) } quad_floating_point;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 247)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 248) /* Extension - An additional structure to hold the guard, round and
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 249) * sticky bits during computations.
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 250) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 251) #define Extall(object) (object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 252) #define Extsign(object) Bitfield_extract( 0, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 253) #define Exthigh31(object) Bitfield_extract( 0, 31,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 254) #define Extlow31(object) Bitfield_extract( 1, 31,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 255) #define Extlow(object) Bitfield_extract( 31, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 256)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 257) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 258) * Single extended - The upper word is just like single precision,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 259) * but one additional word of mantissa is needed.
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 260) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 261) #define Sextallp1(object) (object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 262) #define Sextallp2(object) (object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 263) #define Sextlowp1(object) Bitfield_extract( 31, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 264) #define Sexthighp2(object) Bitfield_extract( 0, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 265) #define Sextlow31p2(object) Bitfield_extract( 1, 31,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 266) #define Sexthiddenoverflow(object) Bitfield_extract( 4, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 267) #define Is_sexthiddenoverflow(object) Bitfield_mask( 4, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 268)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 269) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 270) * Double extended - The upper two words are just like double precision,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 271) * but two additional words of mantissa are needed.
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 272) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 273) #define Dextallp1(object) (object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 274) #define Dextallp2(object) (object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 275) #define Dextallp3(object) (object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 276) #define Dextallp4(object) (object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 277) #define Dextlowp2(object) Bitfield_extract( 31, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 278) #define Dexthighp3(object) Bitfield_extract( 0, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 279) #define Dextlow31p3(object) Bitfield_extract( 1, 31,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 280) #define Dexthiddenoverflow(object) Bitfield_extract( 10, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 281) #define Is_dexthiddenoverflow(object) Bitfield_mask( 10, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 282) #define Deposit_dextlowp4(object,value) Bitfield_deposit(value,31,1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 283)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 284) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 285) * Declare the basic structures for the 3 different
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 286) * fixed-point precisions.
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 287) *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 288) * Single number
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 289) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 290) * |s| integer |
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 291) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 292) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 293) typedef int sgl_integer;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 294)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 295) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 296) * Double number.
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 297) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 298) * |s| high integer |
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 299) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 300) *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 301) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 302) * | low integer |
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 303) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 304) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 305) struct dint {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 306) int wd0;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 307) unsigned int wd1;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 308) };
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 309)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 310) struct dblwd {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 311) unsigned int wd0;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 312) unsigned int wd1;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 313) };
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 314)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 315) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 316) * Quad number.
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 317) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 318) * |s| integer part1 |
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 319) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 320) *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 321) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 322) * | integer part 2 |
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 323) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 324) *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 325) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 326) * | integer part 3 |
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 327) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 328) *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 329) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 330) * | integer part 4 |
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 331) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 332) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 333)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 334) struct quadwd {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 335) int wd0;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 336) unsigned int wd1;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 337) unsigned int wd2;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 338) unsigned int wd3;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 339) };
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 340)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 341) typedef struct quadwd quad_integer;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 342)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 343)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 344) /* useful typedefs */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 345) typedef unsigned int sgl_floating_point;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 346) typedef struct dblwd dbl_floating_point;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 347) typedef struct dint dbl_integer;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 348) typedef struct dblwd dbl_unsigned;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 349)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 350) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 351) * Define the different precisions' parameters.
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 352) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 353) #define SGL_BITLENGTH 32
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 354) #define SGL_EMAX 127
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 355) #define SGL_EMIN (-126)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 356) #define SGL_BIAS 127
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 357) #define SGL_WRAP 192
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 358) #define SGL_INFINITY_EXPONENT (SGL_EMAX+SGL_BIAS+1)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 359) #define SGL_THRESHOLD 32
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 360) #define SGL_EXP_LENGTH 8
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 361) #define SGL_P 24
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 362)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 363) #define DBL_BITLENGTH 64
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 364) #define DBL_EMAX 1023
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 365) #define DBL_EMIN (-1022)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 366) #define DBL_BIAS 1023
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 367) #define DBL_WRAP 1536
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 368) #define DBL_INFINITY_EXPONENT (DBL_EMAX+DBL_BIAS+1)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 369) #define DBL_THRESHOLD 64
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 370) #define DBL_EXP_LENGTH 11
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 371) #define DBL_P 53
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 372)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 373) #define QUAD_BITLENGTH 128
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 374) #define QUAD_EMAX 16383
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 375) #define QUAD_EMIN (-16382)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 376) #define QUAD_BIAS 16383
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 377) #define QUAD_WRAP 24576
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 378) #define QUAD_INFINITY_EXPONENT (QUAD_EMAX+QUAD_BIAS+1)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 379) #define QUAD_P 113
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 380)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 381) /* Boolean Values etc. */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 382) #define FALSE 0
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 383) #define TRUE (!FALSE)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 384) #define NOT !
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 385) #define XOR ^
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 386)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 387) /* other constants */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 388) #undef NULL
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 389) #define NULL 0
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 390) #define NIL 0
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 391) #define SGL 0
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 392) #define DBL 1
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 393) #define BADFMT 2
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 394) #define QUAD 3
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 395)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 396)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 397) /* Types */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 398) typedef int boolean;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 399) typedef int FORMAT;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 400) typedef int VOID;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 401)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 402)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 403) /* Declare status register equivalent to FPUs architecture.
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 404) *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 405) * 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
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 406) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 407) * |V|Z|O|U|I|C| rsv | model | version |RM |rsv|T|r|V|Z|O|U|I|
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 408) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 409) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 410) #define Cbit(object) Bitfield_extract( 5, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 411) #define Tbit(object) Bitfield_extract( 25, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 412) #define Roundingmode(object) Bitfield_extract( 21, 2,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 413) #define Invalidtrap(object) Bitfield_extract( 27, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 414) #define Divisionbyzerotrap(object) Bitfield_extract( 28, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 415) #define Overflowtrap(object) Bitfield_extract( 29, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 416) #define Underflowtrap(object) Bitfield_extract( 30, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 417) #define Inexacttrap(object) Bitfield_extract( 31, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 418) #define Invalidflag(object) Bitfield_extract( 0, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 419) #define Divisionbyzeroflag(object) Bitfield_extract( 1, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 420) #define Overflowflag(object) Bitfield_extract( 2, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 421) #define Underflowflag(object) Bitfield_extract( 3, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 422) #define Inexactflag(object) Bitfield_extract( 4, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 423) #define Allflags(object) Bitfield_extract( 0, 5,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 424)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 425) /* Definitions relevant to the status register */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 426)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 427) /* Rounding Modes */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 428) #define ROUNDNEAREST 0
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 429) #define ROUNDZERO 1
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 430) #define ROUNDPLUS 2
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 431) #define ROUNDMINUS 3
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 432)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 433) /* Exceptions */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 434) #define NOEXCEPTION 0x0
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 435) #define INVALIDEXCEPTION 0x20
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 436) #define DIVISIONBYZEROEXCEPTION 0x10
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 437) #define OVERFLOWEXCEPTION 0x08
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 438) #define UNDERFLOWEXCEPTION 0x04
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 439) #define INEXACTEXCEPTION 0x02
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 440) #define UNIMPLEMENTEDEXCEPTION 0x01
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 441)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 442) /* New exceptions for the 2E Opcode */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 443) #define OPC_2E_INVALIDEXCEPTION 0x30
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 444) #define OPC_2E_OVERFLOWEXCEPTION 0x18
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 445) #define OPC_2E_UNDERFLOWEXCEPTION 0x0c
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 446) #define OPC_2E_INEXACTEXCEPTION 0x12
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 447)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 448) /* Declare exception registers equivalent to FPUs architecture
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 449) *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 450) * 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
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 451) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 452) * |excepttype | r1 | r2/ext | operation |parm |n| t/cond |
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 453) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 454) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 455) #define Allexception(object) (object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 456) #define Exceptiontype(object) Bitfield_extract( 0, 6,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 457) #define Instructionfield(object) Bitfield_mask( 6,26,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 458) #define Parmfield(object) Bitfield_extract( 23, 3,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 459) #define Rabit(object) Bitfield_extract( 24, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 460) #define Ibit(object) Bitfield_extract( 25, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 461)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 462) #define Set_exceptiontype(object,value) Bitfield_deposit(value, 0, 6,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 463) #define Set_parmfield(object,value) Bitfield_deposit(value, 23, 3,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 464) #define Set_exceptiontype_and_instr_field(exception,instruction,object) \
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 465) object = exception << 26 | instruction
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 466)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 467) /* Declare the condition field
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 468) *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 469) * 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
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 470) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 471) * | |G|L|E|U|X|
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 472) * +-------+-------+-------+-------+-------+-------+-------+-------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 473) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 474) #define Greaterthanbit(object) Bitfield_extract( 27, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 475) #define Lessthanbit(object) Bitfield_extract( 28, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 476) #define Equalbit(object) Bitfield_extract( 29, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 477) #define Unorderedbit(object) Bitfield_extract( 30, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 478) #define Exceptionbit(object) Bitfield_extract( 31, 1,object)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 479)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 480) /* An alias name for the status register */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 481) #define Fpustatus_register (*status)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 482)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 483) /**************************************************
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 484) * Status register referencing and manipulation. *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 485) **************************************************/
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 486)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 487) /* Rounding mode */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 488) #define Rounding_mode() Roundingmode(Fpustatus_register)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 489) #define Is_rounding_mode(rmode) \
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 490) (Roundingmode(Fpustatus_register) == rmode)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 491) #define Set_rounding_mode(value) \
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 492) Bitfield_deposit(value,21,2,Fpustatus_register)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 493)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 494) /* Boolean testing of the trap enable bits */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 495) #define Is_invalidtrap_enabled() Invalidtrap(Fpustatus_register)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 496) #define Is_divisionbyzerotrap_enabled() Divisionbyzerotrap(Fpustatus_register)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 497) #define Is_overflowtrap_enabled() Overflowtrap(Fpustatus_register)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 498) #define Is_underflowtrap_enabled() Underflowtrap(Fpustatus_register)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 499) #define Is_inexacttrap_enabled() Inexacttrap(Fpustatus_register)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 500)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 501) /* Set the indicated flags in the status register */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 502) #define Set_invalidflag() Bitfield_deposit(1,0,1,Fpustatus_register)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 503) #define Set_divisionbyzeroflag() Bitfield_deposit(1,1,1,Fpustatus_register)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 504) #define Set_overflowflag() Bitfield_deposit(1,2,1,Fpustatus_register)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 505) #define Set_underflowflag() Bitfield_deposit(1,3,1,Fpustatus_register)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 506) #define Set_inexactflag() Bitfield_deposit(1,4,1,Fpustatus_register)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 507)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 508) #define Clear_all_flags() Bitfield_deposit(0,0,5,Fpustatus_register)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 509)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 510) /* Manipulate the trap and condition code bits (tbit and cbit) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 511) #define Set_tbit() Bitfield_deposit(1,25,1,Fpustatus_register)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 512) #define Clear_tbit() Bitfield_deposit(0,25,1,Fpustatus_register)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 513) #define Is_tbit_set() Tbit(Fpustatus_register)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 514) #define Is_cbit_set() Cbit(Fpustatus_register)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 515)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 516) #define Set_status_cbit(value) \
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 517) Bitfield_deposit(value,5,1,Fpustatus_register)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 518)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 519) /*******************************
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 520) * Condition field referencing *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 521) *******************************/
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 522) #define Unordered(cond) Unorderedbit(cond)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 523) #define Equal(cond) Equalbit(cond)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 524) #define Lessthan(cond) Lessthanbit(cond)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 525) #define Greaterthan(cond) Greaterthanbit(cond)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 526) #define Exception(cond) Exceptionbit(cond)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 527)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 528)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 529) /* Defines for the extension */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 530) #define Ext_isone_sign(extent) (Extsign(extent))
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 531) #define Ext_isnotzero(extent) \
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 532) (Extall(extent))
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 533) #define Ext_isnotzero_lower(extent) \
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 534) (Extlow31(extent))
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 535) #define Ext_leftshiftby1(extent) \
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 536) Extall(extent) <<= 1
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 537) #define Ext_negate(extent) \
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 538) (int )Extall(extent) = 0 - (int )Extall(extent)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 539) #define Ext_setone_low(extent) Bitfield_deposit(1,31,1,extent)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 540) #define Ext_setzero(extent) Extall(extent) = 0
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 541)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 542) typedef int operation;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 543)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 544) /* error messages */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 545)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 546) #define NONE 0
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 547) #define UNDEFFPINST 1
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 548)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 549) /* Function definitions: opcode, opclass */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 550) #define FTEST (1<<2) | 0
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 551) #define FCPY (2<<2) | 0
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 552) #define FABS (3<<2) | 0
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 553) #define FSQRT (4<<2) | 0
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 554) #define FRND (5<<2) | 0
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 555)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 556) #define FCNVFF (0<<2) | 1
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 557) #define FCNVXF (1<<2) | 1
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 558) #define FCNVFX (2<<2) | 1
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 559) #define FCNVFXT (3<<2) | 1
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 560)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 561) #define FCMP (0<<2) | 2
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 562)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 563) #define FADD (0<<2) | 3
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 564) #define FSUB (1<<2) | 3
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 565) #define FMPY (2<<2) | 3
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 566) #define FDIV (3<<2) | 3
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 567) #define FREM (4<<2) | 3
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 568)