xref: /qemu/target/riscv/fpu_helper.c (revision 727385c4) 
1 /*
2  * RISC-V FPU Emulation Helpers for QEMU.
3  *
4  * Copyright (c) 2016-2017 Sagar Karandikar, sagark@eecs.berkeley.edu
5  *
6  * This program is free software; you can redistribute it and/or modify it
7  * under the terms and conditions of the GNU General Public License,
8  * version 2 or later, as published by the Free Software Foundation.
9  *
10  * This program is distributed in the hope it will be useful, but WITHOUT
11  * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
12  * FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License for
13  * more details.
14  *
15  * You should have received a copy of the GNU General Public License along with
16  * this program.  If not, see <http://www.gnu.org/licenses/>.
17  */
18 
19 #include "qemu/osdep.h"
20 #include "cpu.h"
21 #include "qemu/host-utils.h"
22 #include "exec/exec-all.h"
23 #include "exec/helper-proto.h"
24 #include "fpu/softfloat.h"
25 #include "internals.h"
26 
27 target_ulong riscv_cpu_get_fflags(CPURISCVState *env)
28 {
29     int soft = get_float_exception_flags(&env->fp_status);
30     target_ulong hard = 0;
31 
32     hard |= (soft & float_flag_inexact) ? FPEXC_NX : 0;
33     hard |= (soft & float_flag_underflow) ? FPEXC_UF : 0;
34     hard |= (soft & float_flag_overflow) ? FPEXC_OF : 0;
35     hard |= (soft & float_flag_divbyzero) ? FPEXC_DZ : 0;
36     hard |= (soft & float_flag_invalid) ? FPEXC_NV : 0;
37 
38     return hard;
39 }
40 
41 void riscv_cpu_set_fflags(CPURISCVState *env, target_ulong hard)
42 {
43     int soft = 0;
44 
45     soft |= (hard & FPEXC_NX) ? float_flag_inexact : 0;
46     soft |= (hard & FPEXC_UF) ? float_flag_underflow : 0;
47     soft |= (hard & FPEXC_OF) ? float_flag_overflow : 0;
48     soft |= (hard & FPEXC_DZ) ? float_flag_divbyzero : 0;
49     soft |= (hard & FPEXC_NV) ? float_flag_invalid : 0;
50 
51     set_float_exception_flags(soft, &env->fp_status);
52 }
53 
54 void helper_set_rounding_mode(CPURISCVState *env, uint32_t rm)
55 {
56     int softrm;
57 
58     if (rm == 7) {
59         rm = env->frm;
60     }
61     switch (rm) {
62     case 0:
63         softrm = float_round_nearest_even;
64         break;
65     case 1:
66         softrm = float_round_to_zero;
67         break;
68     case 2:
69         softrm = float_round_down;
70         break;
71     case 3:
72         softrm = float_round_up;
73         break;
74     case 4:
75         softrm = float_round_ties_away;
76         break;
77     default:
78         riscv_raise_exception(env, RISCV_EXCP_ILLEGAL_INST, GETPC());
79     }
80 
81     set_float_rounding_mode(softrm, &env->fp_status);
82 }
83 
84 static uint64_t do_fmadd_s(CPURISCVState *env, uint64_t rs1, uint64_t rs2,
85                            uint64_t rs3, int flags)
86 {
87     float32 frs1 = check_nanbox_s(rs1);
88     float32 frs2 = check_nanbox_s(rs2);
89     float32 frs3 = check_nanbox_s(rs3);
90     return nanbox_s(float32_muladd(frs1, frs2, frs3, flags, &env->fp_status));
91 }
92 
93 uint64_t helper_fmadd_s(CPURISCVState *env, uint64_t frs1, uint64_t frs2,
94                         uint64_t frs3)
95 {
96     return do_fmadd_s(env, frs1, frs2, frs3, 0);
97 }
98 
99 uint64_t helper_fmadd_d(CPURISCVState *env, uint64_t frs1, uint64_t frs2,
100                         uint64_t frs3)
101 {
102     return float64_muladd(frs1, frs2, frs3, 0, &env->fp_status);
103 }
104 
105 uint64_t helper_fmsub_s(CPURISCVState *env, uint64_t frs1, uint64_t frs2,
106                         uint64_t frs3)
107 {
108     return do_fmadd_s(env, frs1, frs2, frs3, float_muladd_negate_c);
109 }
110 
111 uint64_t helper_fmsub_d(CPURISCVState *env, uint64_t frs1, uint64_t frs2,
112                         uint64_t frs3)
113 {
114     return float64_muladd(frs1, frs2, frs3, float_muladd_negate_c,
115                           &env->fp_status);
116 }
117 
118 uint64_t helper_fnmsub_s(CPURISCVState *env, uint64_t frs1, uint64_t frs2,
119                          uint64_t frs3)
120 {
121     return do_fmadd_s(env, frs1, frs2, frs3, float_muladd_negate_product);
122 }
123 
124 uint64_t helper_fnmsub_d(CPURISCVState *env, uint64_t frs1, uint64_t frs2,
125                          uint64_t frs3)
126 {
127     return float64_muladd(frs1, frs2, frs3, float_muladd_negate_product,
128                           &env->fp_status);
129 }
130 
131 uint64_t helper_fnmadd_s(CPURISCVState *env, uint64_t frs1, uint64_t frs2,
132                          uint64_t frs3)
133 {
134     return do_fmadd_s(env, frs1, frs2, frs3,
135                       float_muladd_negate_c | float_muladd_negate_product);
136 }
137 
138 uint64_t helper_fnmadd_d(CPURISCVState *env, uint64_t frs1, uint64_t frs2,
139                          uint64_t frs3)
140 {
141     return float64_muladd(frs1, frs2, frs3, float_muladd_negate_c |
142                           float_muladd_negate_product, &env->fp_status);
143 }
144 
145 uint64_t helper_fadd_s(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
146 {
147     float32 frs1 = check_nanbox_s(rs1);
148     float32 frs2 = check_nanbox_s(rs2);
149     return nanbox_s(float32_add(frs1, frs2, &env->fp_status));
150 }
151 
152 uint64_t helper_fsub_s(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
153 {
154     float32 frs1 = check_nanbox_s(rs1);
155     float32 frs2 = check_nanbox_s(rs2);
156     return nanbox_s(float32_sub(frs1, frs2, &env->fp_status));
157 }
158 
159 uint64_t helper_fmul_s(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
160 {
161     float32 frs1 = check_nanbox_s(rs1);
162     float32 frs2 = check_nanbox_s(rs2);
163     return nanbox_s(float32_mul(frs1, frs2, &env->fp_status));
164 }
165 
166 uint64_t helper_fdiv_s(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
167 {
168     float32 frs1 = check_nanbox_s(rs1);
169     float32 frs2 = check_nanbox_s(rs2);
170     return nanbox_s(float32_div(frs1, frs2, &env->fp_status));
171 }
172 
173 uint64_t helper_fmin_s(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
174 {
175     float32 frs1 = check_nanbox_s(rs1);
176     float32 frs2 = check_nanbox_s(rs2);
177     return nanbox_s(env->priv_ver < PRIV_VERSION_1_11_0 ?
178                     float32_minnum(frs1, frs2, &env->fp_status) :
179                     float32_minimum_number(frs1, frs2, &env->fp_status));
180 }
181 
182 uint64_t helper_fmax_s(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
183 {
184     float32 frs1 = check_nanbox_s(rs1);
185     float32 frs2 = check_nanbox_s(rs2);
186     return nanbox_s(env->priv_ver < PRIV_VERSION_1_11_0 ?
187                     float32_maxnum(frs1, frs2, &env->fp_status) :
188                     float32_maximum_number(frs1, frs2, &env->fp_status));
189 }
190 
191 uint64_t helper_fsqrt_s(CPURISCVState *env, uint64_t rs1)
192 {
193     float32 frs1 = check_nanbox_s(rs1);
194     return nanbox_s(float32_sqrt(frs1, &env->fp_status));
195 }
196 
197 target_ulong helper_fle_s(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
198 {
199     float32 frs1 = check_nanbox_s(rs1);
200     float32 frs2 = check_nanbox_s(rs2);
201     return float32_le(frs1, frs2, &env->fp_status);
202 }
203 
204 target_ulong helper_flt_s(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
205 {
206     float32 frs1 = check_nanbox_s(rs1);
207     float32 frs2 = check_nanbox_s(rs2);
208     return float32_lt(frs1, frs2, &env->fp_status);
209 }
210 
211 target_ulong helper_feq_s(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
212 {
213     float32 frs1 = check_nanbox_s(rs1);
214     float32 frs2 = check_nanbox_s(rs2);
215     return float32_eq_quiet(frs1, frs2, &env->fp_status);
216 }
217 
218 target_ulong helper_fcvt_w_s(CPURISCVState *env, uint64_t rs1)
219 {
220     float32 frs1 = check_nanbox_s(rs1);
221     return float32_to_int32(frs1, &env->fp_status);
222 }
223 
224 target_ulong helper_fcvt_wu_s(CPURISCVState *env, uint64_t rs1)
225 {
226     float32 frs1 = check_nanbox_s(rs1);
227     return (int32_t)float32_to_uint32(frs1, &env->fp_status);
228 }
229 
230 target_ulong helper_fcvt_l_s(CPURISCVState *env, uint64_t rs1)
231 {
232     float32 frs1 = check_nanbox_s(rs1);
233     return float32_to_int64(frs1, &env->fp_status);
234 }
235 
236 target_ulong helper_fcvt_lu_s(CPURISCVState *env, uint64_t rs1)
237 {
238     float32 frs1 = check_nanbox_s(rs1);
239     return float32_to_uint64(frs1, &env->fp_status);
240 }
241 
242 uint64_t helper_fcvt_s_w(CPURISCVState *env, target_ulong rs1)
243 {
244     return nanbox_s(int32_to_float32((int32_t)rs1, &env->fp_status));
245 }
246 
247 uint64_t helper_fcvt_s_wu(CPURISCVState *env, target_ulong rs1)
248 {
249     return nanbox_s(uint32_to_float32((uint32_t)rs1, &env->fp_status));
250 }
251 
252 uint64_t helper_fcvt_s_l(CPURISCVState *env, target_ulong rs1)
253 {
254     return nanbox_s(int64_to_float32(rs1, &env->fp_status));
255 }
256 
257 uint64_t helper_fcvt_s_lu(CPURISCVState *env, target_ulong rs1)
258 {
259     return nanbox_s(uint64_to_float32(rs1, &env->fp_status));
260 }
261 
262 target_ulong helper_fclass_s(uint64_t rs1)
263 {
264     float32 frs1 = check_nanbox_s(rs1);
265     return fclass_s(frs1);
266 }
267 
268 uint64_t helper_fadd_d(CPURISCVState *env, uint64_t frs1, uint64_t frs2)
269 {
270     return float64_add(frs1, frs2, &env->fp_status);
271 }
272 
273 uint64_t helper_fsub_d(CPURISCVState *env, uint64_t frs1, uint64_t frs2)
274 {
275     return float64_sub(frs1, frs2, &env->fp_status);
276 }
277 
278 uint64_t helper_fmul_d(CPURISCVState *env, uint64_t frs1, uint64_t frs2)
279 {
280     return float64_mul(frs1, frs2, &env->fp_status);
281 }
282 
283 uint64_t helper_fdiv_d(CPURISCVState *env, uint64_t frs1, uint64_t frs2)
284 {
285     return float64_div(frs1, frs2, &env->fp_status);
286 }
287 
288 uint64_t helper_fmin_d(CPURISCVState *env, uint64_t frs1, uint64_t frs2)
289 {
290     return env->priv_ver < PRIV_VERSION_1_11_0 ?
291             float64_minnum(frs1, frs2, &env->fp_status) :
292             float64_minimum_number(frs1, frs2, &env->fp_status);
293 }
294 
295 uint64_t helper_fmax_d(CPURISCVState *env, uint64_t frs1, uint64_t frs2)
296 {
297     return env->priv_ver < PRIV_VERSION_1_11_0 ?
298             float64_maxnum(frs1, frs2, &env->fp_status) :
299             float64_maximum_number(frs1, frs2, &env->fp_status);
300 }
301 
302 uint64_t helper_fcvt_s_d(CPURISCVState *env, uint64_t rs1)
303 {
304     return nanbox_s(float64_to_float32(rs1, &env->fp_status));
305 }
306 
307 uint64_t helper_fcvt_d_s(CPURISCVState *env, uint64_t rs1)
308 {
309     float32 frs1 = check_nanbox_s(rs1);
310     return float32_to_float64(frs1, &env->fp_status);
311 }
312 
313 uint64_t helper_fsqrt_d(CPURISCVState *env, uint64_t frs1)
314 {
315     return float64_sqrt(frs1, &env->fp_status);
316 }
317 
318 target_ulong helper_fle_d(CPURISCVState *env, uint64_t frs1, uint64_t frs2)
319 {
320     return float64_le(frs1, frs2, &env->fp_status);
321 }
322 
323 target_ulong helper_flt_d(CPURISCVState *env, uint64_t frs1, uint64_t frs2)
324 {
325     return float64_lt(frs1, frs2, &env->fp_status);
326 }
327 
328 target_ulong helper_feq_d(CPURISCVState *env, uint64_t frs1, uint64_t frs2)
329 {
330     return float64_eq_quiet(frs1, frs2, &env->fp_status);
331 }
332 
333 target_ulong helper_fcvt_w_d(CPURISCVState *env, uint64_t frs1)
334 {
335     return float64_to_int32(frs1, &env->fp_status);
336 }
337 
338 target_ulong helper_fcvt_wu_d(CPURISCVState *env, uint64_t frs1)
339 {
340     return (int32_t)float64_to_uint32(frs1, &env->fp_status);
341 }
342 
343 target_ulong helper_fcvt_l_d(CPURISCVState *env, uint64_t frs1)
344 {
345     return float64_to_int64(frs1, &env->fp_status);
346 }
347 
348 target_ulong helper_fcvt_lu_d(CPURISCVState *env, uint64_t frs1)
349 {
350     return float64_to_uint64(frs1, &env->fp_status);
351 }
352 
353 uint64_t helper_fcvt_d_w(CPURISCVState *env, target_ulong rs1)
354 {
355     return int32_to_float64((int32_t)rs1, &env->fp_status);
356 }
357 
358 uint64_t helper_fcvt_d_wu(CPURISCVState *env, target_ulong rs1)
359 {
360     return uint32_to_float64((uint32_t)rs1, &env->fp_status);
361 }
362 
363 uint64_t helper_fcvt_d_l(CPURISCVState *env, target_ulong rs1)
364 {
365     return int64_to_float64(rs1, &env->fp_status);
366 }
367 
368 uint64_t helper_fcvt_d_lu(CPURISCVState *env, target_ulong rs1)
369 {
370     return uint64_to_float64(rs1, &env->fp_status);
371 }
372 
373 target_ulong helper_fclass_d(uint64_t frs1)
374 {
375     return fclass_d(frs1);
376 }
377