target/loongarch: Implement vmskltz/vmskgez/vmsknz
[qemu/kevin.git] / target / riscv / fpu_helper.c
blob5dd14d83908f53c81a5201a6420e79323394f518
1 /*
2 * RISC-V FPU Emulation Helpers for QEMU.
4 * Copyright (c) 2016-2017 Sagar Karandikar, sagark@eecs.berkeley.edu
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.
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.
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/>.
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"
27 target_ulong riscv_cpu_get_fflags(CPURISCVState *env)
29 int soft = get_float_exception_flags(&env->fp_status);
30 target_ulong hard = 0;
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;
38 return hard;
41 void riscv_cpu_set_fflags(CPURISCVState *env, target_ulong hard)
43 int soft = 0;
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;
51 set_float_exception_flags(soft, &env->fp_status);
54 void helper_set_rounding_mode(CPURISCVState *env, uint32_t rm)
56 int softrm;
58 if (rm == RISCV_FRM_DYN) {
59 rm = env->frm;
61 switch (rm) {
62 case RISCV_FRM_RNE:
63 softrm = float_round_nearest_even;
64 break;
65 case RISCV_FRM_RTZ:
66 softrm = float_round_to_zero;
67 break;
68 case RISCV_FRM_RDN:
69 softrm = float_round_down;
70 break;
71 case RISCV_FRM_RUP:
72 softrm = float_round_up;
73 break;
74 case RISCV_FRM_RMM:
75 softrm = float_round_ties_away;
76 break;
77 default:
78 riscv_raise_exception(env, RISCV_EXCP_ILLEGAL_INST, GETPC());
81 set_float_rounding_mode(softrm, &env->fp_status);
84 void helper_set_rounding_mode_chkfrm(CPURISCVState *env, uint32_t rm)
86 int softrm;
88 /* Always validate frm, even if rm != DYN. */
89 if (unlikely(env->frm >= 5)) {
90 riscv_raise_exception(env, RISCV_EXCP_ILLEGAL_INST, GETPC());
92 if (rm == RISCV_FRM_DYN) {
93 rm = env->frm;
95 switch (rm) {
96 case RISCV_FRM_RNE:
97 softrm = float_round_nearest_even;
98 break;
99 case RISCV_FRM_RTZ:
100 softrm = float_round_to_zero;
101 break;
102 case RISCV_FRM_RDN:
103 softrm = float_round_down;
104 break;
105 case RISCV_FRM_RUP:
106 softrm = float_round_up;
107 break;
108 case RISCV_FRM_RMM:
109 softrm = float_round_ties_away;
110 break;
111 case RISCV_FRM_ROD:
112 softrm = float_round_to_odd;
113 break;
114 default:
115 g_assert_not_reached();
118 set_float_rounding_mode(softrm, &env->fp_status);
121 static uint64_t do_fmadd_h(CPURISCVState *env, uint64_t rs1, uint64_t rs2,
122 uint64_t rs3, int flags)
124 float16 frs1 = check_nanbox_h(env, rs1);
125 float16 frs2 = check_nanbox_h(env, rs2);
126 float16 frs3 = check_nanbox_h(env, rs3);
127 return nanbox_h(env, float16_muladd(frs1, frs2, frs3, flags,
128 &env->fp_status));
131 static uint64_t do_fmadd_s(CPURISCVState *env, uint64_t rs1, uint64_t rs2,
132 uint64_t rs3, int flags)
134 float32 frs1 = check_nanbox_s(env, rs1);
135 float32 frs2 = check_nanbox_s(env, rs2);
136 float32 frs3 = check_nanbox_s(env, rs3);
137 return nanbox_s(env, float32_muladd(frs1, frs2, frs3, flags,
138 &env->fp_status));
141 uint64_t helper_fmadd_s(CPURISCVState *env, uint64_t frs1, uint64_t frs2,
142 uint64_t frs3)
144 return do_fmadd_s(env, frs1, frs2, frs3, 0);
147 uint64_t helper_fmadd_d(CPURISCVState *env, uint64_t frs1, uint64_t frs2,
148 uint64_t frs3)
150 return float64_muladd(frs1, frs2, frs3, 0, &env->fp_status);
153 uint64_t helper_fmadd_h(CPURISCVState *env, uint64_t frs1, uint64_t frs2,
154 uint64_t frs3)
156 return do_fmadd_h(env, frs1, frs2, frs3, 0);
159 uint64_t helper_fmsub_s(CPURISCVState *env, uint64_t frs1, uint64_t frs2,
160 uint64_t frs3)
162 return do_fmadd_s(env, frs1, frs2, frs3, float_muladd_negate_c);
165 uint64_t helper_fmsub_d(CPURISCVState *env, uint64_t frs1, uint64_t frs2,
166 uint64_t frs3)
168 return float64_muladd(frs1, frs2, frs3, float_muladd_negate_c,
169 &env->fp_status);
172 uint64_t helper_fmsub_h(CPURISCVState *env, uint64_t frs1, uint64_t frs2,
173 uint64_t frs3)
175 return do_fmadd_h(env, frs1, frs2, frs3, float_muladd_negate_c);
178 uint64_t helper_fnmsub_s(CPURISCVState *env, uint64_t frs1, uint64_t frs2,
179 uint64_t frs3)
181 return do_fmadd_s(env, frs1, frs2, frs3, float_muladd_negate_product);
184 uint64_t helper_fnmsub_d(CPURISCVState *env, uint64_t frs1, uint64_t frs2,
185 uint64_t frs3)
187 return float64_muladd(frs1, frs2, frs3, float_muladd_negate_product,
188 &env->fp_status);
191 uint64_t helper_fnmsub_h(CPURISCVState *env, uint64_t frs1, uint64_t frs2,
192 uint64_t frs3)
194 return do_fmadd_h(env, frs1, frs2, frs3, float_muladd_negate_product);
197 uint64_t helper_fnmadd_s(CPURISCVState *env, uint64_t frs1, uint64_t frs2,
198 uint64_t frs3)
200 return do_fmadd_s(env, frs1, frs2, frs3,
201 float_muladd_negate_c | float_muladd_negate_product);
204 uint64_t helper_fnmadd_d(CPURISCVState *env, uint64_t frs1, uint64_t frs2,
205 uint64_t frs3)
207 return float64_muladd(frs1, frs2, frs3, float_muladd_negate_c |
208 float_muladd_negate_product, &env->fp_status);
211 uint64_t helper_fnmadd_h(CPURISCVState *env, uint64_t frs1, uint64_t frs2,
212 uint64_t frs3)
214 return do_fmadd_h(env, frs1, frs2, frs3,
215 float_muladd_negate_c | float_muladd_negate_product);
218 uint64_t helper_fadd_s(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
220 float32 frs1 = check_nanbox_s(env, rs1);
221 float32 frs2 = check_nanbox_s(env, rs2);
222 return nanbox_s(env, float32_add(frs1, frs2, &env->fp_status));
225 uint64_t helper_fsub_s(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
227 float32 frs1 = check_nanbox_s(env, rs1);
228 float32 frs2 = check_nanbox_s(env, rs2);
229 return nanbox_s(env, float32_sub(frs1, frs2, &env->fp_status));
232 uint64_t helper_fmul_s(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
234 float32 frs1 = check_nanbox_s(env, rs1);
235 float32 frs2 = check_nanbox_s(env, rs2);
236 return nanbox_s(env, float32_mul(frs1, frs2, &env->fp_status));
239 uint64_t helper_fdiv_s(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
241 float32 frs1 = check_nanbox_s(env, rs1);
242 float32 frs2 = check_nanbox_s(env, rs2);
243 return nanbox_s(env, float32_div(frs1, frs2, &env->fp_status));
246 uint64_t helper_fmin_s(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
248 float32 frs1 = check_nanbox_s(env, rs1);
249 float32 frs2 = check_nanbox_s(env, rs2);
250 return nanbox_s(env, env->priv_ver < PRIV_VERSION_1_11_0 ?
251 float32_minnum(frs1, frs2, &env->fp_status) :
252 float32_minimum_number(frs1, frs2, &env->fp_status));
255 uint64_t helper_fmax_s(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
257 float32 frs1 = check_nanbox_s(env, rs1);
258 float32 frs2 = check_nanbox_s(env, rs2);
259 return nanbox_s(env, env->priv_ver < PRIV_VERSION_1_11_0 ?
260 float32_maxnum(frs1, frs2, &env->fp_status) :
261 float32_maximum_number(frs1, frs2, &env->fp_status));
264 uint64_t helper_fsqrt_s(CPURISCVState *env, uint64_t rs1)
266 float32 frs1 = check_nanbox_s(env, rs1);
267 return nanbox_s(env, float32_sqrt(frs1, &env->fp_status));
270 target_ulong helper_fle_s(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
272 float32 frs1 = check_nanbox_s(env, rs1);
273 float32 frs2 = check_nanbox_s(env, rs2);
274 return float32_le(frs1, frs2, &env->fp_status);
277 target_ulong helper_flt_s(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
279 float32 frs1 = check_nanbox_s(env, rs1);
280 float32 frs2 = check_nanbox_s(env, rs2);
281 return float32_lt(frs1, frs2, &env->fp_status);
284 target_ulong helper_feq_s(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
286 float32 frs1 = check_nanbox_s(env, rs1);
287 float32 frs2 = check_nanbox_s(env, rs2);
288 return float32_eq_quiet(frs1, frs2, &env->fp_status);
291 target_ulong helper_fcvt_w_s(CPURISCVState *env, uint64_t rs1)
293 float32 frs1 = check_nanbox_s(env, rs1);
294 return float32_to_int32(frs1, &env->fp_status);
297 target_ulong helper_fcvt_wu_s(CPURISCVState *env, uint64_t rs1)
299 float32 frs1 = check_nanbox_s(env, rs1);
300 return (int32_t)float32_to_uint32(frs1, &env->fp_status);
303 target_ulong helper_fcvt_l_s(CPURISCVState *env, uint64_t rs1)
305 float32 frs1 = check_nanbox_s(env, rs1);
306 return float32_to_int64(frs1, &env->fp_status);
309 target_ulong helper_fcvt_lu_s(CPURISCVState *env, uint64_t rs1)
311 float32 frs1 = check_nanbox_s(env, rs1);
312 return float32_to_uint64(frs1, &env->fp_status);
315 uint64_t helper_fcvt_s_w(CPURISCVState *env, target_ulong rs1)
317 return nanbox_s(env, int32_to_float32((int32_t)rs1, &env->fp_status));
320 uint64_t helper_fcvt_s_wu(CPURISCVState *env, target_ulong rs1)
322 return nanbox_s(env, uint32_to_float32((uint32_t)rs1, &env->fp_status));
325 uint64_t helper_fcvt_s_l(CPURISCVState *env, target_ulong rs1)
327 return nanbox_s(env, int64_to_float32(rs1, &env->fp_status));
330 uint64_t helper_fcvt_s_lu(CPURISCVState *env, target_ulong rs1)
332 return nanbox_s(env, uint64_to_float32(rs1, &env->fp_status));
335 target_ulong helper_fclass_s(CPURISCVState *env, uint64_t rs1)
337 float32 frs1 = check_nanbox_s(env, rs1);
338 return fclass_s(frs1);
341 uint64_t helper_fadd_d(CPURISCVState *env, uint64_t frs1, uint64_t frs2)
343 return float64_add(frs1, frs2, &env->fp_status);
346 uint64_t helper_fsub_d(CPURISCVState *env, uint64_t frs1, uint64_t frs2)
348 return float64_sub(frs1, frs2, &env->fp_status);
351 uint64_t helper_fmul_d(CPURISCVState *env, uint64_t frs1, uint64_t frs2)
353 return float64_mul(frs1, frs2, &env->fp_status);
356 uint64_t helper_fdiv_d(CPURISCVState *env, uint64_t frs1, uint64_t frs2)
358 return float64_div(frs1, frs2, &env->fp_status);
361 uint64_t helper_fmin_d(CPURISCVState *env, uint64_t frs1, uint64_t frs2)
363 return env->priv_ver < PRIV_VERSION_1_11_0 ?
364 float64_minnum(frs1, frs2, &env->fp_status) :
365 float64_minimum_number(frs1, frs2, &env->fp_status);
368 uint64_t helper_fmax_d(CPURISCVState *env, uint64_t frs1, uint64_t frs2)
370 return env->priv_ver < PRIV_VERSION_1_11_0 ?
371 float64_maxnum(frs1, frs2, &env->fp_status) :
372 float64_maximum_number(frs1, frs2, &env->fp_status);
375 uint64_t helper_fcvt_s_d(CPURISCVState *env, uint64_t rs1)
377 return nanbox_s(env, float64_to_float32(rs1, &env->fp_status));
380 uint64_t helper_fcvt_d_s(CPURISCVState *env, uint64_t rs1)
382 float32 frs1 = check_nanbox_s(env, rs1);
383 return float32_to_float64(frs1, &env->fp_status);
386 uint64_t helper_fsqrt_d(CPURISCVState *env, uint64_t frs1)
388 return float64_sqrt(frs1, &env->fp_status);
391 target_ulong helper_fle_d(CPURISCVState *env, uint64_t frs1, uint64_t frs2)
393 return float64_le(frs1, frs2, &env->fp_status);
396 target_ulong helper_flt_d(CPURISCVState *env, uint64_t frs1, uint64_t frs2)
398 return float64_lt(frs1, frs2, &env->fp_status);
401 target_ulong helper_feq_d(CPURISCVState *env, uint64_t frs1, uint64_t frs2)
403 return float64_eq_quiet(frs1, frs2, &env->fp_status);
406 target_ulong helper_fcvt_w_d(CPURISCVState *env, uint64_t frs1)
408 return float64_to_int32(frs1, &env->fp_status);
411 target_ulong helper_fcvt_wu_d(CPURISCVState *env, uint64_t frs1)
413 return (int32_t)float64_to_uint32(frs1, &env->fp_status);
416 target_ulong helper_fcvt_l_d(CPURISCVState *env, uint64_t frs1)
418 return float64_to_int64(frs1, &env->fp_status);
421 target_ulong helper_fcvt_lu_d(CPURISCVState *env, uint64_t frs1)
423 return float64_to_uint64(frs1, &env->fp_status);
426 uint64_t helper_fcvt_d_w(CPURISCVState *env, target_ulong rs1)
428 return int32_to_float64((int32_t)rs1, &env->fp_status);
431 uint64_t helper_fcvt_d_wu(CPURISCVState *env, target_ulong rs1)
433 return uint32_to_float64((uint32_t)rs1, &env->fp_status);
436 uint64_t helper_fcvt_d_l(CPURISCVState *env, target_ulong rs1)
438 return int64_to_float64(rs1, &env->fp_status);
441 uint64_t helper_fcvt_d_lu(CPURISCVState *env, target_ulong rs1)
443 return uint64_to_float64(rs1, &env->fp_status);
446 target_ulong helper_fclass_d(uint64_t frs1)
448 return fclass_d(frs1);
451 uint64_t helper_fadd_h(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
453 float16 frs1 = check_nanbox_h(env, rs1);
454 float16 frs2 = check_nanbox_h(env, rs2);
455 return nanbox_h(env, float16_add(frs1, frs2, &env->fp_status));
458 uint64_t helper_fsub_h(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
460 float16 frs1 = check_nanbox_h(env, rs1);
461 float16 frs2 = check_nanbox_h(env, rs2);
462 return nanbox_h(env, float16_sub(frs1, frs2, &env->fp_status));
465 uint64_t helper_fmul_h(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
467 float16 frs1 = check_nanbox_h(env, rs1);
468 float16 frs2 = check_nanbox_h(env, rs2);
469 return nanbox_h(env, float16_mul(frs1, frs2, &env->fp_status));
472 uint64_t helper_fdiv_h(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
474 float16 frs1 = check_nanbox_h(env, rs1);
475 float16 frs2 = check_nanbox_h(env, rs2);
476 return nanbox_h(env, float16_div(frs1, frs2, &env->fp_status));
479 uint64_t helper_fmin_h(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
481 float16 frs1 = check_nanbox_h(env, rs1);
482 float16 frs2 = check_nanbox_h(env, rs2);
483 return nanbox_h(env, env->priv_ver < PRIV_VERSION_1_11_0 ?
484 float16_minnum(frs1, frs2, &env->fp_status) :
485 float16_minimum_number(frs1, frs2, &env->fp_status));
488 uint64_t helper_fmax_h(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
490 float16 frs1 = check_nanbox_h(env, rs1);
491 float16 frs2 = check_nanbox_h(env, rs2);
492 return nanbox_h(env, env->priv_ver < PRIV_VERSION_1_11_0 ?
493 float16_maxnum(frs1, frs2, &env->fp_status) :
494 float16_maximum_number(frs1, frs2, &env->fp_status));
497 uint64_t helper_fsqrt_h(CPURISCVState *env, uint64_t rs1)
499 float16 frs1 = check_nanbox_h(env, rs1);
500 return nanbox_h(env, float16_sqrt(frs1, &env->fp_status));
503 target_ulong helper_fle_h(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
505 float16 frs1 = check_nanbox_h(env, rs1);
506 float16 frs2 = check_nanbox_h(env, rs2);
507 return float16_le(frs1, frs2, &env->fp_status);
510 target_ulong helper_flt_h(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
512 float16 frs1 = check_nanbox_h(env, rs1);
513 float16 frs2 = check_nanbox_h(env, rs2);
514 return float16_lt(frs1, frs2, &env->fp_status);
517 target_ulong helper_feq_h(CPURISCVState *env, uint64_t rs1, uint64_t rs2)
519 float16 frs1 = check_nanbox_h(env, rs1);
520 float16 frs2 = check_nanbox_h(env, rs2);
521 return float16_eq_quiet(frs1, frs2, &env->fp_status);
524 target_ulong helper_fclass_h(CPURISCVState *env, uint64_t rs1)
526 float16 frs1 = check_nanbox_h(env, rs1);
527 return fclass_h(frs1);
530 target_ulong helper_fcvt_w_h(CPURISCVState *env, uint64_t rs1)
532 float16 frs1 = check_nanbox_h(env, rs1);
533 return float16_to_int32(frs1, &env->fp_status);
536 target_ulong helper_fcvt_wu_h(CPURISCVState *env, uint64_t rs1)
538 float16 frs1 = check_nanbox_h(env, rs1);
539 return (int32_t)float16_to_uint32(frs1, &env->fp_status);
542 target_ulong helper_fcvt_l_h(CPURISCVState *env, uint64_t rs1)
544 float16 frs1 = check_nanbox_h(env, rs1);
545 return float16_to_int64(frs1, &env->fp_status);
548 target_ulong helper_fcvt_lu_h(CPURISCVState *env, uint64_t rs1)
550 float16 frs1 = check_nanbox_h(env, rs1);
551 return float16_to_uint64(frs1, &env->fp_status);
554 uint64_t helper_fcvt_h_w(CPURISCVState *env, target_ulong rs1)
556 return nanbox_h(env, int32_to_float16((int32_t)rs1, &env->fp_status));
559 uint64_t helper_fcvt_h_wu(CPURISCVState *env, target_ulong rs1)
561 return nanbox_h(env, uint32_to_float16((uint32_t)rs1, &env->fp_status));
564 uint64_t helper_fcvt_h_l(CPURISCVState *env, target_ulong rs1)
566 return nanbox_h(env, int64_to_float16(rs1, &env->fp_status));
569 uint64_t helper_fcvt_h_lu(CPURISCVState *env, target_ulong rs1)
571 return nanbox_h(env, uint64_to_float16(rs1, &env->fp_status));
574 uint64_t helper_fcvt_h_s(CPURISCVState *env, uint64_t rs1)
576 float32 frs1 = check_nanbox_s(env, rs1);
577 return nanbox_h(env, float32_to_float16(frs1, true, &env->fp_status));
580 uint64_t helper_fcvt_s_h(CPURISCVState *env, uint64_t rs1)
582 float16 frs1 = check_nanbox_h(env, rs1);
583 return nanbox_s(env, float16_to_float32(frs1, true, &env->fp_status));
586 uint64_t helper_fcvt_h_d(CPURISCVState *env, uint64_t rs1)
588 return nanbox_h(env, float64_to_float16(rs1, true, &env->fp_status));
591 uint64_t helper_fcvt_d_h(CPURISCVState *env, uint64_t rs1)
593 float16 frs1 = check_nanbox_h(env, rs1);
594 return float16_to_float64(frs1, true, &env->fp_status);