^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) NetWinder Floating Point Emulator
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 4) (c) Rebel.COM, 1998,1999
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 5)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 6) Direct questions, comments to Scott Bambrough <scottb@netwinder.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)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 10) #include "fpa11.h"
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 11) #include "softfloat.h"
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 12) #include "fpopcode.h"
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 13)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 14) floatx80 floatx80_exp(floatx80 Fm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 15) floatx80 floatx80_ln(floatx80 Fm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 16) floatx80 floatx80_sin(floatx80 rFm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 17) floatx80 floatx80_cos(floatx80 rFm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 18) floatx80 floatx80_arcsin(floatx80 rFm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 19) floatx80 floatx80_arctan(floatx80 rFm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 20) floatx80 floatx80_log(floatx80 rFm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 21) floatx80 floatx80_tan(floatx80 rFm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 22) floatx80 floatx80_arccos(floatx80 rFm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 23) floatx80 floatx80_pow(floatx80 rFn, floatx80 rFm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 24) floatx80 floatx80_pol(floatx80 rFn, floatx80 rFm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 25)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 26) static floatx80 floatx80_rsf(struct roundingData *roundData, floatx80 rFn, floatx80 rFm)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 27) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 28) return floatx80_sub(roundData, rFm, rFn);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 29) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 30)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 31) static floatx80 floatx80_rdv(struct roundingData *roundData, floatx80 rFn, floatx80 rFm)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 32) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 33) return floatx80_div(roundData, rFm, rFn);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 34) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 35)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 36) static floatx80 (*const dyadic_extended[16])(struct roundingData*, floatx80 rFn, floatx80 rFm) = {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 37) [ADF_CODE >> 20] = floatx80_add,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 38) [MUF_CODE >> 20] = floatx80_mul,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 39) [SUF_CODE >> 20] = floatx80_sub,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 40) [RSF_CODE >> 20] = floatx80_rsf,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 41) [DVF_CODE >> 20] = floatx80_div,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 42) [RDF_CODE >> 20] = floatx80_rdv,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 43) [RMF_CODE >> 20] = floatx80_rem,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 44)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 45) /* strictly, these opcodes should not be implemented */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 46) [FML_CODE >> 20] = floatx80_mul,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 47) [FDV_CODE >> 20] = floatx80_div,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 48) [FRD_CODE >> 20] = floatx80_rdv,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 49) };
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 50)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 51) static floatx80 floatx80_mvf(struct roundingData *roundData, floatx80 rFm)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 52) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 53) return rFm;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 54) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 55)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 56) static floatx80 floatx80_mnf(struct roundingData *roundData, floatx80 rFm)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 57) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 58) rFm.high ^= 0x8000;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 59) return rFm;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 60) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 61)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 62) static floatx80 floatx80_abs(struct roundingData *roundData, floatx80 rFm)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 63) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 64) rFm.high &= 0x7fff;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 65) return rFm;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 66) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 67)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 68) static floatx80 (*const monadic_extended[16])(struct roundingData*, floatx80 rFm) = {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 69) [MVF_CODE >> 20] = floatx80_mvf,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 70) [MNF_CODE >> 20] = floatx80_mnf,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 71) [ABS_CODE >> 20] = floatx80_abs,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 72) [RND_CODE >> 20] = floatx80_round_to_int,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 73) [URD_CODE >> 20] = floatx80_round_to_int,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 74) [SQT_CODE >> 20] = floatx80_sqrt,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 75) [NRM_CODE >> 20] = floatx80_mvf,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 76) };
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 77)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 78) unsigned int ExtendedCPDO(struct roundingData *roundData, const unsigned int opcode, FPREG * rFd)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 79) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 80) FPA11 *fpa11 = GET_FPA11();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 81) floatx80 rFm;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 82) unsigned int Fm, opc_mask_shift;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 83)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 84) Fm = getFm(opcode);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 85) if (CONSTANT_FM(opcode)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 86) rFm = getExtendedConstant(Fm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 87) } else {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 88) switch (fpa11->fType[Fm]) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 89) case typeSingle:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 90) rFm = float32_to_floatx80(fpa11->fpreg[Fm].fSingle);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 91) break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 92)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 93) case typeDouble:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 94) rFm = float64_to_floatx80(fpa11->fpreg[Fm].fDouble);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 95) break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 96)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 97) case typeExtended:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 98) rFm = fpa11->fpreg[Fm].fExtended;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 99) break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 100)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 101) default:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 102) return 0;
^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)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 106) opc_mask_shift = (opcode & MASK_ARITHMETIC_OPCODE) >> 20;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 107) if (!MONADIC_INSTRUCTION(opcode)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 108) unsigned int Fn = getFn(opcode);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 109) floatx80 rFn;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 110)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 111) switch (fpa11->fType[Fn]) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 112) case typeSingle:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 113) rFn = float32_to_floatx80(fpa11->fpreg[Fn].fSingle);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 114) break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 115)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 116) case typeDouble:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 117) rFn = float64_to_floatx80(fpa11->fpreg[Fn].fDouble);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 118) break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 119)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 120) case typeExtended:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 121) rFn = fpa11->fpreg[Fn].fExtended;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 122) break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 123)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 124) default:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 125) return 0;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 126) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 127)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 128) if (dyadic_extended[opc_mask_shift]) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 129) rFd->fExtended = dyadic_extended[opc_mask_shift](roundData, rFn, rFm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 130) } else {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 131) return 0;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 132) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 133) } else {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 134) if (monadic_extended[opc_mask_shift]) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 135) rFd->fExtended = monadic_extended[opc_mask_shift](roundData, rFm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 136) } else {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 137) return 0;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 138) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 139) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 140)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 141) return 1;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 142) }
Orange Pi5 kernel
Deprecated Linux kernel 5.10.110 for OrangePi 5/5B/5+ boards
3 Commits
0 Branches
0 Tags