mirror of
https://github.com/Ed94/Odin.git
synced 2026-08-01 04:10:07 +00:00
Add @(require_results) to core:math procedures
This commit is contained in:
gingerBill
+391
-249
File diff suppressed because it is too large
Load Diff
+27
-23
@@ -5,54 +5,58 @@ import "core:intrinsics"
|
|||||||
|
|
||||||
@(default_calling_convention="none")
|
@(default_calling_convention="none")
|
||||||
foreign _ {
|
foreign _ {
|
||||||
@(link_name="llvm.sin.f16")
|
@(link_name="llvm.sin.f16", require_results)
|
||||||
sin_f16 :: proc(θ: f16) -> f16 ---
|
sin_f16 :: proc(θ: f16) -> f16 ---
|
||||||
@(link_name="llvm.sin.f32")
|
@(link_name="llvm.sin.f32", require_results)
|
||||||
sin_f32 :: proc(θ: f32) -> f32 ---
|
sin_f32 :: proc(θ: f32) -> f32 ---
|
||||||
@(link_name="llvm.sin.f64")
|
@(link_name="llvm.sin.f64", require_results)
|
||||||
sin_f64 :: proc(θ: f64) -> f64 ---
|
sin_f64 :: proc(θ: f64) -> f64 ---
|
||||||
|
|
||||||
@(link_name="llvm.cos.f16")
|
@(link_name="llvm.cos.f16", require_results)
|
||||||
cos_f16 :: proc(θ: f16) -> f16 ---
|
cos_f16 :: proc(θ: f16) -> f16 ---
|
||||||
@(link_name="llvm.cos.f32")
|
@(link_name="llvm.cos.f32", require_results)
|
||||||
cos_f32 :: proc(θ: f32) -> f32 ---
|
cos_f32 :: proc(θ: f32) -> f32 ---
|
||||||
@(link_name="llvm.cos.f64")
|
@(link_name="llvm.cos.f64", require_results)
|
||||||
cos_f64 :: proc(θ: f64) -> f64 ---
|
cos_f64 :: proc(θ: f64) -> f64 ---
|
||||||
|
|
||||||
@(link_name="llvm.pow.f16")
|
@(link_name="llvm.pow.f16", require_results)
|
||||||
pow_f16 :: proc(x, power: f16) -> f16 ---
|
pow_f16 :: proc(x, power: f16) -> f16 ---
|
||||||
@(link_name="llvm.pow.f32")
|
@(link_name="llvm.pow.f32", require_results)
|
||||||
pow_f32 :: proc(x, power: f32) -> f32 ---
|
pow_f32 :: proc(x, power: f32) -> f32 ---
|
||||||
@(link_name="llvm.pow.f64")
|
@(link_name="llvm.pow.f64", require_results)
|
||||||
pow_f64 :: proc(x, power: f64) -> f64 ---
|
pow_f64 :: proc(x, power: f64) -> f64 ---
|
||||||
|
|
||||||
@(link_name="llvm.fmuladd.f16")
|
@(link_name="llvm.fmuladd.f16", require_results)
|
||||||
fmuladd_f16 :: proc(a, b, c: f16) -> f16 ---
|
fmuladd_f16 :: proc(a, b, c: f16) -> f16 ---
|
||||||
@(link_name="llvm.fmuladd.f32")
|
@(link_name="llvm.fmuladd.f32", require_results)
|
||||||
fmuladd_f32 :: proc(a, b, c: f32) -> f32 ---
|
fmuladd_f32 :: proc(a, b, c: f32) -> f32 ---
|
||||||
@(link_name="llvm.fmuladd.f64")
|
@(link_name="llvm.fmuladd.f64", require_results)
|
||||||
fmuladd_f64 :: proc(a, b, c: f64) -> f64 ---
|
fmuladd_f64 :: proc(a, b, c: f64) -> f64 ---
|
||||||
|
|
||||||
@(link_name="llvm.exp.f16")
|
@(link_name="llvm.exp.f16", require_results)
|
||||||
exp_f16 :: proc(x: f16) -> f16 ---
|
exp_f16 :: proc(x: f16) -> f16 ---
|
||||||
@(link_name="llvm.exp.f32")
|
@(link_name="llvm.exp.f32", require_results)
|
||||||
exp_f32 :: proc(x: f32) -> f32 ---
|
exp_f32 :: proc(x: f32) -> f32 ---
|
||||||
@(link_name="llvm.exp.f64")
|
@(link_name="llvm.exp.f64", require_results)
|
||||||
exp_f64 :: proc(x: f64) -> f64 ---
|
exp_f64 :: proc(x: f64) -> f64 ---
|
||||||
}
|
}
|
||||||
|
|
||||||
|
@(require_results)
|
||||||
sqrt_f16 :: proc "contextless" (x: f16) -> f16 {
|
sqrt_f16 :: proc "contextless" (x: f16) -> f16 {
|
||||||
return intrinsics.sqrt(x)
|
return intrinsics.sqrt(x)
|
||||||
}
|
}
|
||||||
|
@(require_results)
|
||||||
sqrt_f32 :: proc "contextless" (x: f32) -> f32 {
|
sqrt_f32 :: proc "contextless" (x: f32) -> f32 {
|
||||||
return intrinsics.sqrt(x)
|
return intrinsics.sqrt(x)
|
||||||
}
|
}
|
||||||
|
@(require_results)
|
||||||
sqrt_f64 :: proc "contextless" (x: f64) -> f64 {
|
sqrt_f64 :: proc "contextless" (x: f64) -> f64 {
|
||||||
return intrinsics.sqrt(x)
|
return intrinsics.sqrt(x)
|
||||||
}
|
}
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
@(require_results)
|
||||||
ln_f64 :: proc "contextless" (x: f64) -> f64 {
|
ln_f64 :: proc "contextless" (x: f64) -> f64 {
|
||||||
// The original C code, the long comment, and the constants
|
// The original C code, the long comment, and the constants
|
||||||
// below are from FreeBSD's /usr/src/lib/msun/src/e_log.c
|
// below are from FreeBSD's /usr/src/lib/msun/src/e_log.c
|
||||||
@@ -154,14 +158,14 @@ ln_f64 :: proc "contextless" (x: f64) -> f64 {
|
|||||||
return k*LN2_HI - ((hfsq - (s*(hfsq+R) + k*LN2_LO)) - f)
|
return k*LN2_HI - ((hfsq - (s*(hfsq+R) + k*LN2_LO)) - f)
|
||||||
}
|
}
|
||||||
|
|
||||||
ln_f16 :: proc "contextless" (x: f16) -> f16 { return #force_inline f16(ln_f64(f64(x))) }
|
@(require_results) ln_f16 :: proc "contextless" (x: f16) -> f16 { return #force_inline f16(ln_f64(f64(x))) }
|
||||||
ln_f32 :: proc "contextless" (x: f32) -> f32 { return #force_inline f32(ln_f64(f64(x))) }
|
@(require_results) ln_f32 :: proc "contextless" (x: f32) -> f32 { return #force_inline f32(ln_f64(f64(x))) }
|
||||||
ln_f16le :: proc "contextless" (x: f16le) -> f16le { return #force_inline f16le(ln_f64(f64(x))) }
|
@(require_results) ln_f16le :: proc "contextless" (x: f16le) -> f16le { return #force_inline f16le(ln_f64(f64(x))) }
|
||||||
ln_f16be :: proc "contextless" (x: f16be) -> f16be { return #force_inline f16be(ln_f64(f64(x))) }
|
@(require_results) ln_f16be :: proc "contextless" (x: f16be) -> f16be { return #force_inline f16be(ln_f64(f64(x))) }
|
||||||
ln_f32le :: proc "contextless" (x: f32le) -> f32le { return #force_inline f32le(ln_f64(f64(x))) }
|
@(require_results) ln_f32le :: proc "contextless" (x: f32le) -> f32le { return #force_inline f32le(ln_f64(f64(x))) }
|
||||||
ln_f32be :: proc "contextless" (x: f32be) -> f32be { return #force_inline f32be(ln_f64(f64(x))) }
|
@(require_results) ln_f32be :: proc "contextless" (x: f32be) -> f32be { return #force_inline f32be(ln_f64(f64(x))) }
|
||||||
ln_f64le :: proc "contextless" (x: f64le) -> f64le { return #force_inline f64le(ln_f64(f64(x))) }
|
@(require_results) ln_f64le :: proc "contextless" (x: f64le) -> f64le { return #force_inline f64le(ln_f64(f64(x))) }
|
||||||
ln_f64be :: proc "contextless" (x: f64be) -> f64be { return #force_inline f64be(ln_f64(f64(x))) }
|
@(require_results) ln_f64be :: proc "contextless" (x: f64be) -> f64be { return #force_inline f64be(ln_f64(f64(x))) }
|
||||||
ln :: proc{
|
ln :: proc{
|
||||||
ln_f16, ln_f16le, ln_f16be,
|
ln_f16, ln_f16le, ln_f16be,
|
||||||
ln_f32, ln_f32le, ln_f32be,
|
ln_f32, ln_f32le, ln_f32be,
|
||||||
|
|||||||
@@ -7,46 +7,47 @@ foreign import "odin_env"
|
|||||||
|
|
||||||
@(default_calling_convention="c")
|
@(default_calling_convention="c")
|
||||||
foreign odin_env {
|
foreign odin_env {
|
||||||
@(link_name="sin")
|
@(link_name="sin", require_results)
|
||||||
sin_f64 :: proc(θ: f64) -> f64 ---
|
sin_f64 :: proc(θ: f64) -> f64 ---
|
||||||
@(link_name="cos")
|
@(link_name="cos", require_results)
|
||||||
cos_f64 :: proc(θ: f64) -> f64 ---
|
cos_f64 :: proc(θ: f64) -> f64 ---
|
||||||
@(link_name="pow")
|
@(link_name="pow", require_results)
|
||||||
pow_f64 :: proc(x, power: f64) -> f64 ---
|
pow_f64 :: proc(x, power: f64) -> f64 ---
|
||||||
@(link_name="fmuladd")
|
@(link_name="fmuladd", require_results)
|
||||||
fmuladd_f64 :: proc(a, b, c: f64) -> f64 ---
|
fmuladd_f64 :: proc(a, b, c: f64) -> f64 ---
|
||||||
@(link_name="ln")
|
@(link_name="ln", require_results)
|
||||||
ln_f64 :: proc(x: f64) -> f64 ---
|
ln_f64 :: proc(x: f64) -> f64 ---
|
||||||
@(link_name="exp")
|
@(link_name="exp", require_results)
|
||||||
exp_f64 :: proc(x: f64) -> f64 ---
|
exp_f64 :: proc(x: f64) -> f64 ---
|
||||||
}
|
}
|
||||||
|
|
||||||
|
@(require_results)
|
||||||
sqrt_f64 :: proc "contextless" (x: f64) -> f64 {
|
sqrt_f64 :: proc "contextless" (x: f64) -> f64 {
|
||||||
return intrinsics.sqrt(x)
|
return intrinsics.sqrt(x)
|
||||||
}
|
}
|
||||||
|
|
||||||
sqrt_f16 :: proc "c" (x: f16) -> f16 { return f16(sqrt_f64(f64(x))) }
|
@(require_results) sqrt_f16 :: proc "c" (x: f16) -> f16 { return f16(sqrt_f64(f64(x))) }
|
||||||
sin_f16 :: proc "c" (θ: f16) -> f16 { return f16(sin_f64(f64(θ))) }
|
@(require_results) sin_f16 :: proc "c" (θ: f16) -> f16 { return f16(sin_f64(f64(θ))) }
|
||||||
cos_f16 :: proc "c" (θ: f16) -> f16 { return f16(cos_f64(f64(θ))) }
|
@(require_results) cos_f16 :: proc "c" (θ: f16) -> f16 { return f16(cos_f64(f64(θ))) }
|
||||||
pow_f16 :: proc "c" (x, power: f16) -> f16 { return f16(pow_f64(f64(x), f64(power))) }
|
@(require_results) pow_f16 :: proc "c" (x, power: f16) -> f16 { return f16(pow_f64(f64(x), f64(power))) }
|
||||||
fmuladd_f16 :: proc "c" (a, b, c: f16) -> f16 { return f16(fmuladd_f64(f64(a), f64(a), f64(c))) }
|
@(require_results) fmuladd_f16 :: proc "c" (a, b, c: f16) -> f16 { return f16(fmuladd_f64(f64(a), f64(a), f64(c))) }
|
||||||
ln_f16 :: proc "c" (x: f16) -> f16 { return f16(ln_f64(f64(x))) }
|
@(require_results) ln_f16 :: proc "c" (x: f16) -> f16 { return f16(ln_f64(f64(x))) }
|
||||||
exp_f16 :: proc "c" (x: f16) -> f16 { return f16(exp_f64(f64(x))) }
|
@(require_results) exp_f16 :: proc "c" (x: f16) -> f16 { return f16(exp_f64(f64(x))) }
|
||||||
|
|
||||||
sqrt_f32 :: proc "c" (x: f32) -> f32 { return f32(sqrt_f64(f64(x))) }
|
@(require_results) sqrt_f32 :: proc "c" (x: f32) -> f32 { return f32(sqrt_f64(f64(x))) }
|
||||||
sin_f32 :: proc "c" (θ: f32) -> f32 { return f32(sin_f64(f64(θ))) }
|
@(require_results) sin_f32 :: proc "c" (θ: f32) -> f32 { return f32(sin_f64(f64(θ))) }
|
||||||
cos_f32 :: proc "c" (θ: f32) -> f32 { return f32(cos_f64(f64(θ))) }
|
@(require_results) cos_f32 :: proc "c" (θ: f32) -> f32 { return f32(cos_f64(f64(θ))) }
|
||||||
pow_f32 :: proc "c" (x, power: f32) -> f32 { return f32(pow_f64(f64(x), f64(power))) }
|
@(require_results) pow_f32 :: proc "c" (x, power: f32) -> f32 { return f32(pow_f64(f64(x), f64(power))) }
|
||||||
fmuladd_f32 :: proc "c" (a, b, c: f32) -> f32 { return f32(fmuladd_f64(f64(a), f64(a), f64(c))) }
|
@(require_results) fmuladd_f32 :: proc "c" (a, b, c: f32) -> f32 { return f32(fmuladd_f64(f64(a), f64(a), f64(c))) }
|
||||||
ln_f32 :: proc "c" (x: f32) -> f32 { return f32(ln_f64(f64(x))) }
|
@(require_results) ln_f32 :: proc "c" (x: f32) -> f32 { return f32(ln_f64(f64(x))) }
|
||||||
exp_f32 :: proc "c" (x: f32) -> f32 { return f32(exp_f64(f64(x))) }
|
@(require_results) exp_f32 :: proc "c" (x: f32) -> f32 { return f32(exp_f64(f64(x))) }
|
||||||
|
|
||||||
ln_f16le :: proc "contextless" (x: f16le) -> f16le { return #force_inline f16le(ln_f64(f64(x))) }
|
@(require_results) ln_f16le :: proc "contextless" (x: f16le) -> f16le { return #force_inline f16le(ln_f64(f64(x))) }
|
||||||
ln_f16be :: proc "contextless" (x: f16be) -> f16be { return #force_inline f16be(ln_f64(f64(x))) }
|
@(require_results) ln_f16be :: proc "contextless" (x: f16be) -> f16be { return #force_inline f16be(ln_f64(f64(x))) }
|
||||||
ln_f32le :: proc "contextless" (x: f32le) -> f32le { return #force_inline f32le(ln_f64(f64(x))) }
|
@(require_results) ln_f32le :: proc "contextless" (x: f32le) -> f32le { return #force_inline f32le(ln_f64(f64(x))) }
|
||||||
ln_f32be :: proc "contextless" (x: f32be) -> f32be { return #force_inline f32be(ln_f64(f64(x))) }
|
@(require_results) ln_f32be :: proc "contextless" (x: f32be) -> f32be { return #force_inline f32be(ln_f64(f64(x))) }
|
||||||
ln_f64le :: proc "contextless" (x: f64le) -> f64le { return #force_inline f64le(ln_f64(f64(x))) }
|
@(require_results) ln_f64le :: proc "contextless" (x: f64le) -> f64le { return #force_inline f64le(ln_f64(f64(x))) }
|
||||||
ln_f64be :: proc "contextless" (x: f64be) -> f64be { return #force_inline f64be(ln_f64(f64(x))) }
|
@(require_results) ln_f64be :: proc "contextless" (x: f64be) -> f64be { return #force_inline f64be(ln_f64(f64(x))) }
|
||||||
ln :: proc{
|
ln :: proc{
|
||||||
ln_f16, ln_f16le, ln_f16be,
|
ln_f16, ln_f16le, ln_f16be,
|
||||||
ln_f32, ln_f32le, ln_f32be,
|
ln_f32, ln_f32le, ln_f32be,
|
||||||
|
|||||||
+14
-12
@@ -117,13 +117,14 @@ erf :: proc{
|
|||||||
erf_f64,
|
erf_f64,
|
||||||
}
|
}
|
||||||
|
|
||||||
erf_f16 :: proc "contextless" (x: f16) -> f16 { return f16(erf_f64(f64(x))) }
|
@(require_results) erf_f16 :: proc "contextless" (x: f16) -> f16 { return f16(erf_f64(f64(x))) }
|
||||||
erf_f16le :: proc "contextless" (x: f16le) -> f16le { return f16le(erf_f64(f64(x))) }
|
@(require_results) erf_f16le :: proc "contextless" (x: f16le) -> f16le { return f16le(erf_f64(f64(x))) }
|
||||||
erf_f16be :: proc "contextless" (x: f16be) -> f16be { return f16be(erf_f64(f64(x))) }
|
@(require_results) erf_f16be :: proc "contextless" (x: f16be) -> f16be { return f16be(erf_f64(f64(x))) }
|
||||||
erf_f32 :: proc "contextless" (x: f32) -> f32 { return f32(erf_f64(f64(x))) }
|
@(require_results) erf_f32 :: proc "contextless" (x: f32) -> f32 { return f32(erf_f64(f64(x))) }
|
||||||
erf_f32le :: proc "contextless" (x: f32le) -> f32le { return f32le(erf_f64(f64(x))) }
|
@(require_results) erf_f32le :: proc "contextless" (x: f32le) -> f32le { return f32le(erf_f64(f64(x))) }
|
||||||
erf_f32be :: proc "contextless" (x: f32be) -> f32be { return f32be(erf_f64(f64(x))) }
|
@(require_results) erf_f32be :: proc "contextless" (x: f32be) -> f32be { return f32be(erf_f64(f64(x))) }
|
||||||
|
|
||||||
|
@(require_results)
|
||||||
erf_f64 :: proc "contextless" (x: f64) -> f64 {
|
erf_f64 :: proc "contextless" (x: f64) -> f64 {
|
||||||
erx :: 0h3FEB0AC160000000
|
erx :: 0h3FEB0AC160000000
|
||||||
// Coefficients for approximation to erf in [0, 0.84375]
|
// Coefficients for approximation to erf in [0, 0.84375]
|
||||||
@@ -268,13 +269,14 @@ erfc :: proc{
|
|||||||
erfc_f64,
|
erfc_f64,
|
||||||
}
|
}
|
||||||
|
|
||||||
erfc_f16 :: proc "contextless" (x: f16) -> f16 { return f16(erfc_f64(f64(x))) }
|
@(require_results) erfc_f16 :: proc "contextless" (x: f16) -> f16 { return f16(erfc_f64(f64(x))) }
|
||||||
erfc_f16le :: proc "contextless" (x: f16le) -> f16le { return f16le(erfc_f64(f64(x))) }
|
@(require_results) erfc_f16le :: proc "contextless" (x: f16le) -> f16le { return f16le(erfc_f64(f64(x))) }
|
||||||
erfc_f16be :: proc "contextless" (x: f16be) -> f16be { return f16be(erfc_f64(f64(x))) }
|
@(require_results) erfc_f16be :: proc "contextless" (x: f16be) -> f16be { return f16be(erfc_f64(f64(x))) }
|
||||||
erfc_f32 :: proc "contextless" (x: f32) -> f32 { return f32(erfc_f64(f64(x))) }
|
@(require_results) erfc_f32 :: proc "contextless" (x: f32) -> f32 { return f32(erfc_f64(f64(x))) }
|
||||||
erfc_f32le :: proc "contextless" (x: f32le) -> f32le { return f32le(erfc_f64(f64(x))) }
|
@(require_results) erfc_f32le :: proc "contextless" (x: f32le) -> f32le { return f32le(erfc_f64(f64(x))) }
|
||||||
erfc_f32be :: proc "contextless" (x: f32be) -> f32be { return f32be(erfc_f64(f64(x))) }
|
@(require_results) erfc_f32be :: proc "contextless" (x: f32be) -> f32be { return f32be(erfc_f64(f64(x))) }
|
||||||
|
|
||||||
|
@(require_results)
|
||||||
erfc_f64 :: proc "contextless" (x: f64) -> f64 {
|
erfc_f64 :: proc "contextless" (x: f64) -> f64 {
|
||||||
erx :: 0h3FEB0AC160000000
|
erx :: 0h3FEB0AC160000000
|
||||||
// Coefficients for approximation to erf in [0, 0.84375]
|
// Coefficients for approximation to erf in [0, 0.84375]
|
||||||
|
|||||||
@@ -65,7 +65,7 @@ package math
|
|||||||
// The polynomial is valid for 33 <= x <= 172; larger values are only used
|
// The polynomial is valid for 33 <= x <= 172; larger values are only used
|
||||||
// in reciprocal and produce denormalized floats. The lower precision there
|
// in reciprocal and produce denormalized floats. The lower precision there
|
||||||
// masks any imprecision in the polynomial.
|
// masks any imprecision in the polynomial.
|
||||||
@(private="file")
|
@(private="file", require_results)
|
||||||
stirling :: proc "contextless" (x: f64) -> (f64, f64) {
|
stirling :: proc "contextless" (x: f64) -> (f64, f64) {
|
||||||
@(static) gamS := [?]f64{
|
@(static) gamS := [?]f64{
|
||||||
+7.87311395793093628397e-04,
|
+7.87311395793093628397e-04,
|
||||||
@@ -93,6 +93,7 @@ stirling :: proc "contextless" (x: f64) -> (f64, f64) {
|
|||||||
return y1, SQRT_TWO_PI * w * y2
|
return y1, SQRT_TWO_PI * w * y2
|
||||||
}
|
}
|
||||||
|
|
||||||
|
@(require_results)
|
||||||
gamma_f64 :: proc "contextless" (x: f64) -> f64 {
|
gamma_f64 :: proc "contextless" (x: f64) -> f64 {
|
||||||
is_neg_int :: proc "contextless" (x: f64) -> bool {
|
is_neg_int :: proc "contextless" (x: f64) -> bool {
|
||||||
if x < 0 {
|
if x < 0 {
|
||||||
@@ -210,14 +211,14 @@ gamma_f64 :: proc "contextless" (x: f64) -> f64 {
|
|||||||
}
|
}
|
||||||
|
|
||||||
|
|
||||||
gamma_f16 :: proc "contextless" (x: f16) -> f16 { return f16(gamma_f64(f64(x))) }
|
@(require_results) gamma_f16 :: proc "contextless" (x: f16) -> f16 { return f16(gamma_f64(f64(x))) }
|
||||||
gamma_f16le :: proc "contextless" (x: f16le) -> f16le { return f16le(gamma_f64(f64(x))) }
|
@(require_results) gamma_f16le :: proc "contextless" (x: f16le) -> f16le { return f16le(gamma_f64(f64(x))) }
|
||||||
gamma_f16be :: proc "contextless" (x: f16be) -> f16be { return f16be(gamma_f64(f64(x))) }
|
@(require_results) gamma_f16be :: proc "contextless" (x: f16be) -> f16be { return f16be(gamma_f64(f64(x))) }
|
||||||
gamma_f32 :: proc "contextless" (x: f32) -> f32 { return f32(gamma_f64(f64(x))) }
|
@(require_results) gamma_f32 :: proc "contextless" (x: f32) -> f32 { return f32(gamma_f64(f64(x))) }
|
||||||
gamma_f32le :: proc "contextless" (x: f32le) -> f32le { return f32le(gamma_f64(f64(x))) }
|
@(require_results) gamma_f32le :: proc "contextless" (x: f32le) -> f32le { return f32le(gamma_f64(f64(x))) }
|
||||||
gamma_f32be :: proc "contextless" (x: f32be) -> f32be { return f32be(gamma_f64(f64(x))) }
|
@(require_results) gamma_f32be :: proc "contextless" (x: f32be) -> f32be { return f32be(gamma_f64(f64(x))) }
|
||||||
gamma_f64le :: proc "contextless" (x: f64le) -> f64le { return f64le(gamma_f64(f64(x))) }
|
@(require_results) gamma_f64le :: proc "contextless" (x: f64le) -> f64le { return f64le(gamma_f64(f64(x))) }
|
||||||
gamma_f64be :: proc "contextless" (x: f64be) -> f64be { return f64be(gamma_f64(f64(x))) }
|
@(require_results) gamma_f64be :: proc "contextless" (x: f64be) -> f64be { return f64be(gamma_f64(f64(x))) }
|
||||||
|
|
||||||
gamma :: proc{
|
gamma :: proc{
|
||||||
gamma_f16, gamma_f16le, gamma_f16be,
|
gamma_f16, gamma_f16le, gamma_f16be,
|
||||||
|
|||||||
@@ -80,7 +80,9 @@ package math
|
|||||||
//
|
//
|
||||||
|
|
||||||
|
|
||||||
|
@(require_results)
|
||||||
lgamma_f64 :: proc "contextless" (x: f64) -> (lgamma: f64, sign: int) {
|
lgamma_f64 :: proc "contextless" (x: f64) -> (lgamma: f64, sign: int) {
|
||||||
|
@(require_results)
|
||||||
sin_pi :: proc "contextless" (x: f64) -> f64 {
|
sin_pi :: proc "contextless" (x: f64) -> f64 {
|
||||||
if x < 0.25 {
|
if x < 0.25 {
|
||||||
return -sin(PI * x)
|
return -sin(PI * x)
|
||||||
@@ -345,14 +347,14 @@ lgamma_f64 :: proc "contextless" (x: f64) -> (lgamma: f64, sign: int) {
|
|||||||
}
|
}
|
||||||
|
|
||||||
|
|
||||||
lgamma_f16 :: proc "contextless" (x: f16) -> (lgamma: f16, sign: int) { r, s := lgamma_f64(f64(x)); return f16(r), s }
|
@(require_results) lgamma_f16 :: proc "contextless" (x: f16) -> (lgamma: f16, sign: int) { r, s := lgamma_f64(f64(x)); return f16(r), s }
|
||||||
lgamma_f32 :: proc "contextless" (x: f32) -> (lgamma: f32, sign: int) { r, s := lgamma_f64(f64(x)); return f32(r), s }
|
@(require_results) lgamma_f32 :: proc "contextless" (x: f32) -> (lgamma: f32, sign: int) { r, s := lgamma_f64(f64(x)); return f32(r), s }
|
||||||
lgamma_f16le :: proc "contextless" (x: f16le) -> (lgamma: f16le, sign: int) { r, s := lgamma_f64(f64(x)); return f16le(r), s }
|
@(require_results) lgamma_f16le :: proc "contextless" (x: f16le) -> (lgamma: f16le, sign: int) { r, s := lgamma_f64(f64(x)); return f16le(r), s }
|
||||||
lgamma_f16be :: proc "contextless" (x: f16be) -> (lgamma: f16be, sign: int) { r, s := lgamma_f64(f64(x)); return f16be(r), s }
|
@(require_results) lgamma_f16be :: proc "contextless" (x: f16be) -> (lgamma: f16be, sign: int) { r, s := lgamma_f64(f64(x)); return f16be(r), s }
|
||||||
lgamma_f32le :: proc "contextless" (x: f32le) -> (lgamma: f32le, sign: int) { r, s := lgamma_f64(f64(x)); return f32le(r), s }
|
@(require_results) lgamma_f32le :: proc "contextless" (x: f32le) -> (lgamma: f32le, sign: int) { r, s := lgamma_f64(f64(x)); return f32le(r), s }
|
||||||
lgamma_f32be :: proc "contextless" (x: f32be) -> (lgamma: f32be, sign: int) { r, s := lgamma_f64(f64(x)); return f32be(r), s }
|
@(require_results) lgamma_f32be :: proc "contextless" (x: f32be) -> (lgamma: f32be, sign: int) { r, s := lgamma_f64(f64(x)); return f32be(r), s }
|
||||||
lgamma_f64le :: proc "contextless" (x: f64le) -> (lgamma: f64le, sign: int) { r, s := lgamma_f64(f64(x)); return f64le(r), s }
|
@(require_results) lgamma_f64le :: proc "contextless" (x: f64le) -> (lgamma: f64le, sign: int) { r, s := lgamma_f64(f64(x)); return f64le(r), s }
|
||||||
lgamma_f64be :: proc "contextless" (x: f64be) -> (lgamma: f64be, sign: int) { r, s := lgamma_f64(f64(x)); return f64be(r), s }
|
@(require_results) lgamma_f64be :: proc "contextless" (x: f64be) -> (lgamma: f64be, sign: int) { r, s := lgamma_f64(f64(x)); return f64be(r), s }
|
||||||
|
|
||||||
lgamma :: proc{
|
lgamma :: proc{
|
||||||
lgamma_f16, lgamma_f16le, lgamma_f16be,
|
lgamma_f16, lgamma_f16le, lgamma_f16be,
|
||||||
|
|||||||
@@ -90,15 +90,16 @@ log1p :: proc {
|
|||||||
log1p_f64le,
|
log1p_f64le,
|
||||||
log1p_f64be,
|
log1p_f64be,
|
||||||
}
|
}
|
||||||
log1p_f16 :: proc "contextless" (x: f16) -> f16 { return f16(log1p_f64(f64(x))) }
|
@(require_results) log1p_f16 :: proc "contextless" (x: f16) -> f16 { return f16(log1p_f64(f64(x))) }
|
||||||
log1p_f32 :: proc "contextless" (x: f32) -> f32 { return f32(log1p_f64(f64(x))) }
|
@(require_results) log1p_f32 :: proc "contextless" (x: f32) -> f32 { return f32(log1p_f64(f64(x))) }
|
||||||
log1p_f16le :: proc "contextless" (x: f16le) -> f16le { return f16le(log1p_f64(f64(x))) }
|
@(require_results) log1p_f16le :: proc "contextless" (x: f16le) -> f16le { return f16le(log1p_f64(f64(x))) }
|
||||||
log1p_f16be :: proc "contextless" (x: f16be) -> f16be { return f16be(log1p_f64(f64(x))) }
|
@(require_results) log1p_f16be :: proc "contextless" (x: f16be) -> f16be { return f16be(log1p_f64(f64(x))) }
|
||||||
log1p_f32le :: proc "contextless" (x: f32le) -> f32le { return f32le(log1p_f64(f64(x))) }
|
@(require_results) log1p_f32le :: proc "contextless" (x: f32le) -> f32le { return f32le(log1p_f64(f64(x))) }
|
||||||
log1p_f32be :: proc "contextless" (x: f32be) -> f32be { return f32be(log1p_f64(f64(x))) }
|
@(require_results) log1p_f32be :: proc "contextless" (x: f32be) -> f32be { return f32be(log1p_f64(f64(x))) }
|
||||||
log1p_f64le :: proc "contextless" (x: f64le) -> f64le { return f64le(log1p_f64(f64(x))) }
|
@(require_results) log1p_f64le :: proc "contextless" (x: f64le) -> f64le { return f64le(log1p_f64(f64(x))) }
|
||||||
log1p_f64be :: proc "contextless" (x: f64be) -> f64be { return f64be(log1p_f64(f64(x))) }
|
@(require_results) log1p_f64be :: proc "contextless" (x: f64be) -> f64be { return f64be(log1p_f64(f64(x))) }
|
||||||
|
|
||||||
|
@(require_results)
|
||||||
log1p_f64 :: proc "contextless" (x: f64) -> f64 {
|
log1p_f64 :: proc "contextless" (x: f64) -> f64 {
|
||||||
SQRT2_M1 :: 0h3fda827999fcef34 // sqrt(2)-1
|
SQRT2_M1 :: 0h3fda827999fcef34 // sqrt(2)-1
|
||||||
SQRT2_HALF_M1 :: 0hbfd2bec333018866 // sqrt(2)/2-1
|
SQRT2_HALF_M1 :: 0hbfd2bec333018866 // sqrt(2)/2-1
|
||||||
|
|||||||
Reference in New Issue
Block a user