code wiki / _hdl_build / _tn_probe.nx
_tn_probe.nx source
↩ module page · 31 lines · 1555 B
1// _tn_probe.nx -- localizer: reduction + tan sign at the pi/2 pole crossing.
2import "nx_syscalls.nx"
3import "nx_bigfloat120.nx"
4import "nx_bigfloat120_div.nx"
5import "nx_bigfloat120_trig.nx"
6import "nx_bigfloat120_tan.nx"
7
8func tp_puts(s: *u8) -> i64 { var n: i64 = 0; while s[n] != (0 as u8) { n = n + 1 } sys_write(1, s, n); return 0 }
9func tp_putn(v: i64) -> i64 { let bb: *u8 = sys_mmap(28); var m: i64 = v; if m < 0 { m = 0 - m; sys_write(1, "-" as *u8, 1) }; let t: *u8 = sys_mmap(28); var k: i64 = 0; if m == 0 { t[0] = 48; k = 1 }; while m > 0 { t[k] = 48 + (m % 10); m = m / 10; k = k + 1 }; var i: i64 = 0; while i < k { bb[i] = t[k-1-i]; i = i + 1 }; sys_write(1, bb, k); return 0 }
10func tp_puthex(v: i64) -> i64 { let bb: *u8 = sys_mmap(20); var i: i64 = 0; while i < 16 { let nib: i64 = (v >> ((15 - i) * 4)) & 15; if nib < 10 { bb[i] = 48 + nib } else { bb[i] = 55 + nib } i = i + 1 } sys_write(1, bb, 16); return 0 }
11
12func main() -> i64 {
13 let x: i64 = 0x3FF921FB54442D19
14 let xb: *i64 = bf_new()
15 bf_set_f64(xb, x)
16 let r: *i64 = bf_new()
17 let kq: *i64 = sys_mmap(16) as *i64
18 _bf_trig_reduce(r, kq, xb)
19 tp_puts("k=" as *u8); tp_putn(kq[0])
20 tp_puts(" rneg=" as *u8); tp_putn(kq[1])
21 tp_puts(" r.e=" as *u8); tp_putn(r[0])
22 tp_puts("\n" as *u8)
23 let mag: *i64 = bf_new()
24 let s: i64 = bf_tan_xb(mag, xb)
25 tp_puts("sign=" as *u8); tp_putn(s)
26 tp_puts(" mag.e=" as *u8); tp_putn(mag[0])
27 tp_puts(" out=" as *u8); tp_puthex(bf_tan_f64(x))
28 tp_puts("\n" as *u8)
29 sys_exit(0)
30 return 0
31}