Files
ss/include/nw4r/math/math_arithmetic.h
T
Elijah Thomas 9c3c480b24 g3d source (#123)
* g3d_calcvtx

GetData seems to have changed -> dwarf says r is a local and using ofs_to_ptr didnt work

* g3d_light and g3d_fog

sdata2 splits and func ordering

* g3d_scnproc

* g3d_init

* g3d_scnmdl

* g3d_scnmdlsmpl

* g3d_scnroot

* g3d_scnobj

* g3d_res* progress

* g3d_resmdl OK

* g3d_restev OK

* g3d_resmat OK

* g3d_resvtx and g3d_restex OK

* g3d_resnode OK

* g3d_resanm OK

* g3d_resanmchr Progress

* the rest of g3d_res* OK

* g3d_anmvis OK

* g3d_anmclr OK

* Some Splitting

* more OK, Inline Issue in g3d_anmtexsrt

* g3d_obj, g3d_anmobj, g3d_gpu, g3d_tmem, g3d_cpu OK

* g3d_state OK

* g3d/dcc OK

* Include fixup

* More Fixups

* g3d_camera OK

* g3d_draw OK

* g3d_calcworld OK

* g3d_calcworld actually OK

* g3d_workmem, g3d_dcc OK

* g3d_calcview OK

* g3d_anmtexsrt OK with DONT_INLINE

* g3d_transform OK (Feels Cheaty)

* g3d_resanmchr OK

* g3d_draw1mat1shp Close

* g3d_draw1mat1shp OK (Thanks Lago!). Ran symbol applying script
2025-03-16 11:26:15 -04:00

140 lines
2.4 KiB
C++

#ifndef NW4R_MATH_ARITHMETIC_H
#define NW4R_MATH_ARITHMETIC_H
#include "math.h"
#include "nw4r/types_nw4r.h" // IWYU pragma: export
#include "rvl/OS.h" // IWYU pragma: export
#define NW4R_MATH_QNAN (-(0.0f / 0.0f))
#define NW4R_MATH_FLT_MAX 3.402823466e+38f
namespace nw4r {
namespace math {
namespace detail {
f32 FExp(f32 x);
f32 FLog(f32 x);
} // namespace detail
f32 FrSqrt(f32 x);
inline f32 FAbs(register f32 x) {
// clang-format off
asm {
fabs x, x
}
// clang-format on
return x;
}
inline f32 FCeil(f32 x) {
return ceilf(x);
}
inline f32 FExp(f32 x) {
return detail::FExp(x);
}
inline f32 FFloor(f32 x) {
return floorf(x);
}
inline f32 FInv(register f32 x) {
register f32 work0, work1, work2, work3;
// clang-format off
asm {
fmr work1, x // x
fres work0, work1 // 1/x
// Refine estimate
ps_add work2, work0, work0 // 2/x
ps_mul work3, work0, work0 // 1/x^2
ps_nmsub work0, work1, work3, work2 // -(x * 1/x^2 - 2/x)
}
// clang-format on
return work0;
}
inline f32 FMod(f32 x, f32 y) {
return fmodf(x, y);
}
inline f32 FModf(f32 x, f32 *y) {
return modff(x, y);
}
inline f32 FSqrt(f32 x) {
return (x <= 0.0f) ? 0.0f : x * FrSqrt(x);
}
inline f32 FLog(f32 x) {
if (x > 0.0f) {
return detail::FLog(x);
}
return NW4R_MATH_QNAN;
}
inline f32 FSelect(register f32 value, register f32 ge_zero, register f32 lt_zero) {
register f32 ret;
// clang-format off
asm {
fsel ret, value, ge_zero, lt_zero
}
// clang-format on
return ret;
}
inline f32 U16ToF32(u16 arg) {
f32 ret;
OSu16tof32(&arg, &ret);
return ret;
}
inline u16 F32ToU16(f32 arg) {
u16 ret;
OSf32tou16(&arg, &ret);
return ret;
}
inline f32 S16ToF32(s16 arg) {
f32 ret;
OSs16tof32(&arg, &ret);
return ret;
}
inline s16 F32ToS16(f32 arg) {
s16 ret;
OSf32tos16(&arg, &ret);
return ret;
}
inline u32 F32AsU32(f32 arg) {
return *reinterpret_cast<u32 *>(&arg);
}
inline f32 U32AsF32(u32 arg) {
return *reinterpret_cast<f32 *>(&arg);
}
inline s32 FGetExpPart(f32 x) {
s32 s = F32AsU32(x);
return ((s >> 23) & 0xFF) - 0x7F;
}
inline f32 FGetMantPart(f32 x) {
u32 u = F32AsU32(x);
return U32AsF32((u & 0x807FFFFF) | 0x3F800000);
}
} // namespace math
} // namespace nw4r
#endif