mcordingley/linearalgebra
Matrix math for PHP.
Requires
 php: ^7.0.0
Requires (Dev)
 phpunit/phpunit: ^5.5.3
 codeclimate/phptestreporter: devmaster
README
Standalone Linear Algebra Library for PHP
New Major Version
LinearAlgebra has just reached version 2.0. This release includes several breaking changes. Check the changelog for specifics about what has changed, but the biggest changes are that the root namespace has been capitalized and that the library now depends on PHP 7.0 or higher. A number of lesserused functions have been removed, new functions added, and the internals of the package are now much cleaner.
If you're still using PHP 5.x, the 1.0 series still exists and supports PHP 5, though it will receive no further updates.
Installation
composer require mcordingley/LinearAlgebra
Alternately, include this in your composer.json and then update:
"mcordingley/linearalgebra": "^2.0.0"
If Composer isn't an option for you, clone this repository and run buildphar.php
to generate a phar
archive that you can include into your project. PHP will autoload classes from inside the archive as needed.
Usage
Start with a use
statement for the class:
use MCordingley\LinearAlgebra\Matrix;
Then, instantiate a new instance of the matrix class like so:
$matrix = new Matrix([
[0, 1, 2],
[3, 4, 5],
[6, 7, 8]
]);
You can also generate an identity matrix with the identity
factory function:
$threeByThreeIdentityMatrix = Matrix::identity(3);
With the matrix instance, you can retrieve individual elements with get
using
the zerobased indices of the row and column that you want:
$element = $matrix>get($row, $column);
It's also possible to find out how large the matrix is with getRowCount()
and getColumnCount()
:
$rows = $matrix>getRowCount();
$columns = $matrix>getColumnCount();
You can also add, subtract, and multiply the matrix with scalar values and other matrices. All operations return a new Matrix and do not modify the underlying matrix:
$addedScalar = $matrix>addScalar(3);
$addedMatrix = $matrix>addMatrix($anotherMatrix);
$subtractedScalar = $matrix>subtractScalar(2);
$subtractedMatrix = $matrix>subtractMatrix($anotherMatrix);
$multipliedByScalar = $matrix>multiplyScalar(4);
$multipliedByMatrix = $matrix>multiplyMatrix($anotherMatrix);
Matrices can be compared with equals
to see if they're equal:
if ($matrix1>equals($matrix2)) {
// Equality for all!
}
In addition to these basic operations, the Matrix class offers other common matrix operations:
$matrix>inverse()
$matrix>adjugate()
$matrix>determinant()
$matrix>trace()
$matrix>transpose()
You can get the upper and lower triangular matrices by calling upper(bool)
and lower(bool)
. The lone argument tells
whether the main diagonal of the triangular matrix should be set to ones (true
) or the value of the parent matrix
(false
).
It's also possible to run a map over the matrix:
$squaredElements = $matrix>map(function($element, $row, $column, $matrix) {
return $element * $element
});
Submatrices may be extracted with sliceColumns($offset, $length)
and sliceRows($offset, $length)
. The semantics of
the arguments are the same as PHP's array_slice
.
Similarly, spliceColumns($offset, $length, $replacement)
and spliceRows($offset, $length, $replacement)
can be used
to create new matrices with specific rows or columns removed or replaced. Unlike the native PHP array_splice
, these
operations do not modify the matrix in place and return the removed elements, but instead return a new matrix with the
splice applied.
If you need to combine together matrices, you can do so by calling the concatenation methods:
$m1 = new Matrix([
[1,2,3],
[4,5,6],
]);
$m2 = new Matrix([
[7],
[8],
]);
$m3 = new Matrix([[3,2,1]]);
$m4 = $m1>concatenateRight($m2);
// [
// [1,2,3,7],
// [4,5,6,8],
// ]
$m5 = $m1>concatenateBottom($m3);
// [
// [1,2,3],
// [4,5,6],
// [3,2,1],
// ]
LU and LUP decomposition methods are available as separate classes and both expose lower()
and upper()
for the L
and U portions of the decompositions, respectively. The LUP decomposition additionally exposes permutationMatrix
and
permutationArray
to fetch the P component of the decomposition as well as parity
to return the total number of
pivots performed.
Changelog

2.0.0
 Drop support for PHP 5.x
 Introduce strict scalar type hints
 Drop deprecated functions and properties.
 Tighten up interface with the
final
andprivate
keywords. diagonal()
now returns a full matrix, not a vector. Rename
adjoint()
toadjugate()
for clarity.  Add
entrywise()
to compute the Hadamard product.  Add
upper()
andlower()
 Add
sliceColumns()
andsliceRows()
 Add
spliceColumns()
andspliceRows()
 Add
LU
andLUP
decompositions as classes.

1.3.2
 Deprecate
__toString()
magic method.  Deprecate
isSymmetric()
.
 Deprecate

1.3.1
 Deprecate use of the
ArrayAccess
interface.  More internal code style fixes.
 Deprecate use of the

1.3.0
 Fix typo in names of
concatenateRight()
andconcatenateBottom()
 Remove generated Phar file. Users who need it should use the
buildphar.php
script to generate one.  Refactor LUDecomposition to have a less awkward constructor.
 Split
add()
intoaddMatrix()
andaddScalar()
. Deprecateadd()
.  Split
subtract()
intosubtractMatrix()
andsubtractScalar()
. Deprecatesubtract()
.  Split
multiply()
intomultiplyMatrix()
andmultiplyScalar()
. Deprecatemultiply()
.  Add
getRowCount()
andgetColumnCount()
accessors.  Deprecate
rows
andcolumns
properties.
 Fix typo in names of

1.2.0
 Added
concatenateBottom($other)
 Added
concatencateRight($other)
 Added

1.1.0
 Added
diagonal()
.
 Added

1.0.0
 Switch to PSR4 from PSR0.
 Take
isSymmetric()
public.  Rearrange source in
Matrix.php
to be more readable and PSRcompliant.

0.9.1
 Fix several bugs with the Cholesky decomposition and inverse.

0.9.0
 Bump version up to represent that this is close to it's final form.
 Merged PR for faster
inverse
calculations  KISS
Vector
class goodbye.  Renamed
eq
toequals
.  Removed
set
function, so instantiated objects are immutable.

0.3.0
 Added the
identity
factory function  Using Cholesky decomposition for faster matrix inversions for applicable matrices
 Added
eq
function to test equality of matrices  Implemented the ArrayAccess interface
 Added the

0.2.0
 Created the Vector type
\MCordingley
namespace is now\mcordingley
 Matrix functions that return a new Matrix now return a new instance of the called class

0.1.0
 Created the Matrix type
 Scalar Addition
 Scalar Subtraction
 Scalar Multiplication
 Matrix Addition
 Matrix Subtraction
 Matrix Multiplication
 Inverse
 Adjoint
 Determinant
 Trace
 Transpose
 Submatrix
 Map