mirror of
https://github.com/zeldaret/oot-vc.git
synced 2026-10-01 15:32:12 -04:00
Yanis
84 lines
2.0 KiB
C
84 lines
2.0 KiB
C
#ifndef _MATH_H_
|
|
#define _MATH_H_
|
|
|
|
#include "intrinsics.h"
|
|
#include "macros.h"
|
|
#include "revolution/types.h"
|
|
|
|
#define DONT_INLINE_SQRT
|
|
|
|
#define M_PI 3.1415926535897932
|
|
#define M_SQRT3 1.7320499420166016
|
|
|
|
extern int __float_nan[];
|
|
extern int __float_huge[];
|
|
|
|
#define NAN (*(f32*)__float_nan)
|
|
#define INFINITY (*(f32*)__float_huge)
|
|
|
|
// f64 bit-twiddling macros
|
|
#define __HI(x) (((s32*)&x)[0])
|
|
#define __LO(x) (((s32*)&x)[1])
|
|
|
|
f64 ldexp(f64 x, int exp);
|
|
f64 pow(f64 x, f64 y);
|
|
f64 ceil(f64 x);
|
|
f64 floor(f64 x);
|
|
f64 copysign(f64 x, f64 y);
|
|
|
|
f64 sin(f64 x);
|
|
f64 cos(f64 x);
|
|
|
|
static inline f32 sinf(f32 x) { return (f32)sin((f64)x); }
|
|
|
|
static inline f32 cosf(f32 x) { return (f32)cos((f64)x); }
|
|
|
|
f32 tanf(f32 x);
|
|
f32 log10f(f32);
|
|
|
|
static inline f64 fabs(f64 x) { return __fabs(x); }
|
|
|
|
// In reality, these are "weak" functions which all have C++ names (except scalbn).
|
|
// We fake it by defining them as strong C functions instead.
|
|
f32 sin__Ff(f32 x);
|
|
f32 cos__Ff(f32 x);
|
|
f64 scalbn(f64 x, int n);
|
|
f64 fabs__Fd(f64 x);
|
|
f32 fabsf__Ff(f32 x);
|
|
|
|
#define FP_NAN 1
|
|
#define FP_INFINITE 2
|
|
#define FP_ZERO 3
|
|
#define FP_NORMAL 4
|
|
#define FP_SUBNORMAL 5
|
|
|
|
int __fpclassifyd__Fd(f64 x);
|
|
|
|
#define fpclassify(x) __fpclassifyd__Fd(x)
|
|
#define isnormal(x) (fpclassify(x) == FP_NORMAL)
|
|
#define isnan(x) (fpclassify(x) == FP_NAN)
|
|
#define isinf(x) (fpclassify(x) == FP_INFINITE)
|
|
#define isfinite(x) (fpclassify(x) > FP_INFINITE)
|
|
|
|
f64 sqrt(f64 x);
|
|
|
|
static inline f32 sqrtf(f32 x) { return (f32)sqrt((f64)x); }
|
|
|
|
static inline f32 _inv_sqrtf(f32 x) {
|
|
const f32 _half = .5f;
|
|
const f32 _three = 3.0f;
|
|
|
|
if (x > 0.0f) {
|
|
f32 guess = __frsqrte((f64)x); /* returns an approximation to */
|
|
guess = _half * guess * (_three - guess * guess * x); /* now have 8 sig bits */
|
|
guess = _half * guess * (_three - guess * guess * x); /* now have 16 sig bits */
|
|
guess = _half * guess * (_three - guess * guess * x); /* now have >24 sig bits */
|
|
return guess;
|
|
} else if (x) {
|
|
return NAN;
|
|
}
|
|
return INFINITY;
|
|
}
|
|
|
|
#endif
|