nx_matrix.nx

buildroot/runtime/nx_matrix.nx

20257 B522 linesdepth 3pulls 4 transitivereach 4 importersview sourcekind librarytopic matrix

docs dependencies structs consts functions

about

nx_matrix.nx -- general matrix algebra (i64, row-major). Closes the gap that nx_kalman / nx_cam / nx_pc currently bridge with inline matrix arithmetic. Provides: + alloc / set / get / copy + identity / zeros / fill + transpose + multiply (C = A * B) + 2x2 + 3x3 determinant + 2x2 + 3x3 inverse (in Q14) Entries are pure i64. For fractional operations (inverse), gains are emitted in Q14 fixed-point so subsequent multiplications can be divided by NX_MATRIX_Q. genealogy_id: classical_linear_algebra + cramer_1750_determinant lineage_id: row_major_matrix + cofactor_expansion

dependencies 2 imports · 4 importers

imports: syscalls.nxnx_isqrt.nx

imported by: nx_appliedmath_gate.nxnx_matrix_eig_gate.nxnx_matrix_test.nxnx_swcompare_bench.nx

structs

30	struct Matrix

consts

28	const NX_MATRIX_Q: i64 = 16384 // Q14
159	const NXM_EIG_OK: i64 = 0
160	const NXM_EIG_ERR_NOT_SQUARE: i64 = 0 - 1
161	const NXM_EIG_ERR_EMPTY: i64 = 0 - 2
162	const NXM_EIG_ERR_ASYMMETRIC: i64 = 0 - 3 // Jacobi is only valid for symmetric input; refuse, never guess
163	const NXM_EIG_SWEEPS_DEFAULT: i64 = 12
166	const NXM_EIG_TOL: i64 = 1
318	const NXM_SVD_OK: i64 = 0
319	const NXM_SVD_ERR_EMPTY: i64 = 0 - 2
320	const NXM_SVD_ERR_DIMS: i64 = 0 - 4 // u must be m x n and v must be n x n

functions

38	func nx_matrix_alloc(rows: i64, cols: i64) -> *Matrix called by 7: main nx_matrix_transpose nx_matrix_eig_sym nx_matrix_svd main main+1
51	func nx_matrix_set(m: *Matrix, r: i64, c: i64, v: i64) -> i64 called by 12: main nx_matrix_identity nx_matrix_transpose nx_matrix_multiply nx_matrix_eig_sym nxm_ata+6
56	func nx_matrix_get(m: *Matrix, r: i64, c: i64) -> i64 called by 16: am_meq nx_matrix_transpose nx_matrix_multiply nx_matrix_is_symmetric nx_matrix_eig_sym nx_matrix_residual+10
60	func nx_matrix_fill(m: *Matrix, v: i64) -> i64 called by 2: nx_matrix_identity main
72	func nx_matrix_identity(m: *Matrix, unit_value: i64) -> i64 called by 3: main nx_matrix_eig_sym main calls 2: nx_matrix_fill nx_matrix_set
87	func nx_matrix_transpose(a: Matrix) -> Matrix called by 3: main main main calls 3: nx_matrix_alloc nx_matrix_set nx_matrix_get
106	func nx_matrix_multiply(a: Matrix, b: Matrix, c: *Matrix) -> i64 called by 3: main main main calls 2: nx_matrix_get nx_matrix_set
168	func nxm_fmul(a: i64, b: i64) -> i64 { return (a * b) / NX_MATRIX_Q } called by 6: nx_matrix_eig_sym nx_matrix_residual nxm_ata nx_matrix_svd nx_matrix_svd_residual main
170	func nxm_fdiv(a: i64, b: i64) -> i64 called by 2: nx_matrix_eig_sym nx_matrix_svd
177	func nxm_fsqrt(x: i64) -> i64 called by 2: nx_matrix_eig_sym nx_matrix_svd calls 1: nx_isqrt
182	func nxm_abs(v: i64) -> i64 { if v < 0 { return 0 - v } return v } called by 4: nx_matrix_is_symmetric nx_matrix_eig_sym nx_matrix_residual nx_matrix_svd_residual
186	func nx_matrix_is_symmetric(m: *Matrix) -> i64 called by 2: nx_matrix_eig_sym main calls 2: nxm_abs nx_matrix_get
202	func nx_matrix_eig_sym(a: Matrix, evals: i64, v: *Matrix, sweeps: i64) -> i64 called by 2: nx_matrix_svd main calls 9: nx_matrix_is_symmetric nx_matrix_alloc nx_matrix_set nx_matrix_get nx_matrix_identity nxm_abs+3
278	func nx_matrix_residual(a: Matrix, v: Matrix, evals: *i64) -> i64 called by 1: main calls 3: nxm_fmul nx_matrix_get nxm_abs
325	func nxm_ata(a: Matrix, out: Matrix) -> i64 called by 1: nx_matrix_svd calls 3: nxm_fmul nx_matrix_get nx_matrix_set
343	func nx_matrix_svd(a: Matrix, u: Matrix, sv: i64, v: Matrix, sweeps: i64) -> i64 called by 1: main calls 8: nx_matrix_alloc nxm_ata nx_matrix_eig_sym nxm_fsqrt nx_matrix_get nx_matrix_set+2
410	func nx_matrix_svd_residual(a: Matrix, u: Matrix, sv: i64, v: Matrix) -> i64 called by 1: main calls 3: nxm_fmul nx_matrix_get nxm_abs
433	func nx_matrix_det_2x2(m: *Matrix) -> i64 called by 5: main nx_matrix_inv_2x2 main main main calls 1: nx_matrix_get
440	func nx_matrix_det_3x3(m: *Matrix) -> i64 called by 3: nx_matrix_inv_3x3 main main calls 1: nx_matrix_get
467	func nx_matrix_inv_2x2(m: Matrix, out: Matrix) -> i64 called by 1: main calls 3: nx_matrix_det_2x2 nx_matrix_get nx_matrix_set
486	func nx_matrix_inv_3x3(m: Matrix, out: Matrix) -> i64 called by 1: main calls 3: nx_matrix_det_3x3 nx_matrix_get nx_matrix_set