2 NetWinder Floating Point Emulator
3 (c) Rebel.COM, 1998,1999
5 Direct questions, comments to Scott Bambrough <scottb@netwinder.org>
7 This program is free software; you can redistribute it and/or modify
8 it under the terms of the GNU General Public License as published by
9 the Free Software Foundation; either version 2 of the License, or
10 (at your option) any later version.
12 This program is distributed in the hope that it will be useful,
13 but WITHOUT ANY WARRANTY; without even the implied warranty of
14 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
15 GNU General Public License for more details.
17 You should have received a copy of the GNU General Public License
18 along with this program; if not, see <http://www.gnu.org/licenses/>.
21 #include "qemu/osdep.h"
26 //#include "fpmodule.h"
27 //#include "fpmodule.inl"
29 //#include <asm/system.h>
32 FPA11
* qemufpa
= NULL
;
33 CPUARMState
* user_registers
;
35 /* Reset the FPA11 chip. Called to initialize and reset the emulator. */
39 FPA11
*fpa11
= GET_FPA11();
41 /* initialize the register type array */
44 fpa11
->fType
[i
] = typeNone
;
47 /* FPSR: set system id to FP_EMULATOR, set AC, clear all other bits */
48 fpa11
->fpsr
= FP_EMULATOR
| BIT_AC
;
50 /* FPCR: set SB, AB and DA bits, clear all others */
52 fpa11
->fpcr
= MASK_RESET
;
56 * Real FPA11 hardware does not handle NaNs, but always takes an
57 * exception for them to be software-emulated (ARM7500FE datasheet
58 * section 10.4). There is no documented architectural requirement
59 * for NaN propagation rules and it will depend on how the OS
60 * level software emulation opted to do it. We here use prop_s_ab
61 * which matches the later VFP hardware choice and how QEMU's
62 * fpa11 emulation has worked in the past. The real Linux kernel
63 * does something slightly different: arch/arm/nwfpe/softfloat-specialize
64 * propagateFloat64NaN() has the curious behaviour that it prefers
65 * the QNaN over the SNaN, but if both are QNaN it picks A and
66 * if both are SNaN it picks B. In theory we could add this as
67 * a NaN propagation rule, but in practice FPA11 emulation is so
68 * close to totally dead that it's not worth trying to match it at
71 set_float_2nan_prop_rule(float_2nan_prop_s_ab
, &fpa11
->fp_status
);
74 void SetRoundingMode(const unsigned int opcode
)
77 FPA11
*fpa11
= GET_FPA11();
80 fpa11
->fpcr
&= ~MASK_ROUNDING_MODE
;
82 switch (opcode
& MASK_ROUNDING_MODE
)
85 case ROUND_TO_NEAREST
:
86 rounding_mode
= float_round_nearest_even
;
88 fpa11
->fpcr
|= ROUND_TO_NEAREST
;
92 case ROUND_TO_PLUS_INFINITY
:
93 rounding_mode
= float_round_up
;
95 fpa11
->fpcr
|= ROUND_TO_PLUS_INFINITY
;
99 case ROUND_TO_MINUS_INFINITY
:
100 rounding_mode
= float_round_down
;
102 fpa11
->fpcr
|= ROUND_TO_MINUS_INFINITY
;
107 rounding_mode
= float_round_to_zero
;
109 fpa11
->fpcr
|= ROUND_TO_ZERO
;
113 set_float_rounding_mode(rounding_mode
, &fpa11
->fp_status
);
116 void SetRoundingPrecision(const unsigned int opcode
)
118 FloatX80RoundPrec rounding_precision
;
119 FPA11
*fpa11
= GET_FPA11();
121 fpa11
->fpcr
&= ~MASK_ROUNDING_PRECISION
;
123 switch (opcode
& MASK_ROUNDING_PRECISION
) {
125 rounding_precision
= floatx80_precision_s
;
127 fpa11
->fpcr
|= ROUND_SINGLE
;
132 rounding_precision
= floatx80_precision_d
;
134 fpa11
->fpcr
|= ROUND_DOUBLE
;
139 rounding_precision
= floatx80_precision_x
;
141 fpa11
->fpcr
|= ROUND_EXTENDED
;
146 rounding_precision
= floatx80_precision_x
;
149 set_floatx80_rounding_precision(rounding_precision
, &fpa11
->fp_status
);
152 /* Emulate the instruction in the opcode. */
153 /* ??? This is not thread safe. */
154 unsigned int EmulateAll(unsigned int opcode
, FPA11
* qfpa
, CPUARMState
* qregs
)
156 unsigned int nRc
= 0;
157 // unsigned long flags;
160 // save_flags(flags); sti();
162 /* Check that this is really an FPA11 instruction: the coprocessor
163 * field in bits [11:8] must be 1 or 2.
165 cp
= (opcode
>> 8) & 0xf;
166 if (cp
!= 1 && cp
!= 2) {
171 user_registers
=qregs
;
174 fprintf(stderr
,"emulating FP insn 0x%08x, PC=0x%08x\n",
175 opcode
, qregs
[ARM_REG_PC
]);
179 if (fpa11
->initflag
== 0) /* good place for __builtin_expect */
182 SetRoundingMode(ROUND_TO_NEAREST
);
183 SetRoundingPrecision(ROUND_EXTENDED
);
187 set_float_exception_flags(0, &fpa11
->fp_status
);
189 if (TEST_OPCODE(opcode
,MASK_CPRT
))
191 //fprintf(stderr,"emulating CPRT\n");
192 /* Emulate conversion opcodes. */
193 /* Emulate register transfer opcodes. */
194 /* Emulate comparison opcodes. */
195 nRc
= EmulateCPRT(opcode
);
197 else if (TEST_OPCODE(opcode
,MASK_CPDO
))
199 //fprintf(stderr,"emulating CPDO\n");
200 /* Emulate monadic arithmetic opcodes. */
201 /* Emulate dyadic arithmetic opcodes. */
202 nRc
= EmulateCPDO(opcode
);
204 else if (TEST_OPCODE(opcode
,MASK_CPDT
))
206 //fprintf(stderr,"emulating CPDT\n");
207 /* Emulate load/store opcodes. */
208 /* Emulate load/store multiple opcodes. */
209 nRc
= EmulateCPDT(opcode
);
213 /* Invalid instruction detected. Return FALSE. */
217 // restore_flags(flags);
218 if(nRc
== 1 && get_float_exception_flags(&fpa11
->fp_status
))
220 //printf("fef 0x%x\n",float_exception_flags);
221 nRc
= -get_float_exception_flags(&fpa11
->fp_status
);
224 //printf("returning %d\n",nRc);
229 unsigned int EmulateAll1(unsigned int opcode
)
231 switch ((opcode
>> 24) & 0xf)
235 if ((opcode
>> 20) & 0x1)
237 switch ((opcode
>> 8) & 0xf)
239 case 0x1: return PerformLDF(opcode
); break;
240 case 0x2: return PerformLFM(opcode
); break;
246 switch ((opcode
>> 8) & 0xf)
248 case 0x1: return PerformSTF(opcode
); break;
249 case 0x2: return PerformSFM(opcode
); break;
257 return EmulateCPDO(opcode
);
259 return EmulateCPRT(opcode
);