^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) union float64_components {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 15) float64 f64;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 16) unsigned int i[2];
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 17) };
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 18)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 19) float64 float64_exp(float64 Fm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 20) float64 float64_ln(float64 Fm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 21) float64 float64_sin(float64 rFm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 22) float64 float64_cos(float64 rFm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 23) float64 float64_arcsin(float64 rFm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 24) float64 float64_arctan(float64 rFm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 25) float64 float64_log(float64 rFm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 26) float64 float64_tan(float64 rFm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 27) float64 float64_arccos(float64 rFm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 28) float64 float64_pow(float64 rFn, float64 rFm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 29) float64 float64_pol(float64 rFn, float64 rFm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 30)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 31) static float64 float64_rsf(struct roundingData *roundData, float64 rFn, float64 rFm)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 32) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 33) return float64_sub(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 float64 float64_rdv(struct roundingData *roundData, float64 rFn, float64 rFm)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 37) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 38) return float64_div(roundData, rFm, rFn);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 39) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 40)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 41) static float64 (*const dyadic_double[16])(struct roundingData*, float64 rFn, float64 rFm) = {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 42) [ADF_CODE >> 20] = float64_add,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 43) [MUF_CODE >> 20] = float64_mul,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 44) [SUF_CODE >> 20] = float64_sub,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 45) [RSF_CODE >> 20] = float64_rsf,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 46) [DVF_CODE >> 20] = float64_div,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 47) [RDF_CODE >> 20] = float64_rdv,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 48) [RMF_CODE >> 20] = float64_rem,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 49)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 50) /* strictly, these opcodes should not be implemented */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 51) [FML_CODE >> 20] = float64_mul,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 52) [FDV_CODE >> 20] = float64_div,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 53) [FRD_CODE >> 20] = float64_rdv,
^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 float64 float64_mvf(struct roundingData *roundData,float64 rFm)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 57) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 58) return rFm;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 59) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 60)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 61) static float64 float64_mnf(struct roundingData *roundData,float64 rFm)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 62) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 63) union float64_components u;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 64)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 65) u.f64 = rFm;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 66) #ifdef __ARMEB__
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 67) u.i[0] ^= 0x80000000;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 68) #else
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 69) u.i[1] ^= 0x80000000;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 70) #endif
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 71)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 72) return u.f64;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 73) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 74)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 75) static float64 float64_abs(struct roundingData *roundData,float64 rFm)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 76) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 77) union float64_components u;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 78)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 79) u.f64 = rFm;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 80) #ifdef __ARMEB__
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 81) u.i[0] &= 0x7fffffff;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 82) #else
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 83) u.i[1] &= 0x7fffffff;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 84) #endif
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 85)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 86) return u.f64;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 87) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 88)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 89) static float64 (*const monadic_double[16])(struct roundingData *, float64 rFm) = {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 90) [MVF_CODE >> 20] = float64_mvf,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 91) [MNF_CODE >> 20] = float64_mnf,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 92) [ABS_CODE >> 20] = float64_abs,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 93) [RND_CODE >> 20] = float64_round_to_int,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 94) [URD_CODE >> 20] = float64_round_to_int,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 95) [SQT_CODE >> 20] = float64_sqrt,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 96) [NRM_CODE >> 20] = float64_mvf,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 97) };
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 98)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 99) unsigned int DoubleCPDO(struct roundingData *roundData, const unsigned int opcode, FPREG * rFd)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 100) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 101) FPA11 *fpa11 = GET_FPA11();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 102) float64 rFm;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 103) unsigned int Fm, opc_mask_shift;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 104)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 105) Fm = getFm(opcode);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 106) if (CONSTANT_FM(opcode)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 107) rFm = getDoubleConstant(Fm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 108) } else {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 109) switch (fpa11->fType[Fm]) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 110) case typeSingle:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 111) rFm = float32_to_float64(fpa11->fpreg[Fm].fSingle);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 112) break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 113)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 114) case typeDouble:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 115) rFm = fpa11->fpreg[Fm].fDouble;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 116) break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 117)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 118) default:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 119) return 0;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 120) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 121) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 122)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 123) opc_mask_shift = (opcode & MASK_ARITHMETIC_OPCODE) >> 20;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 124) if (!MONADIC_INSTRUCTION(opcode)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 125) unsigned int Fn = getFn(opcode);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 126) float64 rFn;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 127)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 128) switch (fpa11->fType[Fn]) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 129) case typeSingle:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 130) rFn = float32_to_float64(fpa11->fpreg[Fn].fSingle);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 131) break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 132)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 133) case typeDouble:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 134) rFn = fpa11->fpreg[Fn].fDouble;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 135) break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 136)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 137) default:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 138) return 0;
^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) if (dyadic_double[opc_mask_shift]) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 142) rFd->fDouble = dyadic_double[opc_mask_shift](roundData, rFn, rFm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 143) } else {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 144) return 0;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 145) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 146) } else {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 147) if (monadic_double[opc_mask_shift]) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 148) rFd->fDouble = monadic_double[opc_mask_shift](roundData, rFm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 149) } else {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 150) return 0;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 151) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 152) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 153)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 154) return 1;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 155) }
Orange Pi5 kernel
Deprecated Linux kernel 5.10.110 for OrangePi 5/5B/5+ boards
3 Commits
0 Branches
0 Tags