7.5 Entrywise Operations and Arithmetic
procedure
M : (Matrix Number) N : (Matrix Number)
procedure
M : (Matrix Number) N : (Matrix Number)
procedure
M : (Matrix Number) N : (Matrix Number)
Matrix addition, subtraction and products respectively.
For matrix addition and subtraction all matrices must have the same shape.
For matrix product the number of columns of one matrix must equal the number of rows in the following matrix.
Examples:
> (define A (matrix ([1 2] [3 4])))
> (define B (matrix ([5 6] [7 8])))
> (define C (matrix ([ 9 10 11] [12 13 14])))
> (matrix+ A B) (array #[#[6 8] #[10 12]])
> (matrix- A B) (array #[#[-4 -4] #[-4 -4]])
> (matrix* A C) (array #[#[33 36 39] #[75 82 89]])
Computes (matrix* M ...) with n arguments, but more efficiently.
M must be a square-matrix? and n must be nonnegative.
Examples:
> (matrix* (matrix-expt (matrix [[1 1] [1 0]]) 100) (col-matrix [0 1])) (array #[#[354224848179261915075] #[218922995834555169026]])
> (->col-matrix (list (fibonacci 100) (fibonacci 99))) (array #[#[354224848179261915075] #[218922995834555169026]])
Computes the matrix zM, a matrix of the same shape as M
where each entry in M is multiplied with z.
Example:
> (matrix-scale (matrix [[1 2] [3 4]]) 2) (array #[#[2 4] #[6 8]])
procedure
(matrix-map f M) → (Matrix R)
f : (A -> R) M : (Matrix A) (matrix-map f M0 M1 N ...) → (Matrix R) f : (A B Ts ... -> R) M0 : (Matrix A) M1 : (Matrix B) N : (Matrix Ts)
Like array-map, but requires at least one matrix argument and never broadcasts.
Examples:
> (matrix-map sqr (matrix [[1 2] [3 4]])) (array #[#[1 4] #[9 16]])
> (matrix-map + (matrix [[1 2] [3 4]]) (matrix [[5 6] [7 8]])) (array #[#[6 8] #[10 12]])
Returns #t when its arguments are the same size and are equal entrywise.
See matrix-relative-error and matrix-absolute-error for equality testing that is tolerant to floating-point error.