1: 2: 3: 4: 5: 6: 7: 8: 9: 10: 11: 12: 13: 14: 15: 16: 17: 18: 19: 20: 21: 22: 23: 24: 25: 26: 27: 28: 29: 30: 31: 32: 33: 34: 35: 36: 37: 38: 39: 40: 41: 42: 43: 44: 45: 46: 47: 48: 49: 50: 51: 52: 53: 54: 55: 56: 57: 58: 59: 60: 61: 62: 63: 64: 65: 66: 67: 68: 69: 70: 71: 72: 73: 74: 75: 76: 77: 78: 79: 80: 81: 82: 83: 84: 85: 86: 87: 88: 89: 90: 91: 92: 93: 94: 95: 96: 97: 98: 99: 100: 101: 102: 103: 104: 105: 106: 107: 108: 109: 110: 111: 112: 113: 114: 115: 116: 117: 118: 119: 120: 121: 122: 123: 124: 125: 126: 127: 128: 129: 130: 131: 132: 133: 134: 135: 136: 137: 138: 139: 140: 141: 142: 143: 144: 145: 146: 147: 148: 149: 150: 151: 152: 153: 154: 155: 156: 157: 158: 159: 160: 161: 162: 163: 164: 165: 166: 167: 168: 169: 170: 171: 172: 173: 174: 175: 176: 177: 178: 179: 180: 181: 182: 183: 184: 185: 186: 187: 188: 189: 190: 191: 192: 193: 194: 195: 196: 197: 198: 199: 200: 201: 202: 203: 204: 205: 206: 207: 208: 209: 210: 211: 212: 213: 214: 215: 216: 217: 218: 219: 220: 221: 222: 223: 224: 225: 226: 227: 228: 229: 230: 231: 232: 233: 234: 235: 236: 237: 238: 239: 240:
<?php
namespace Chippyash\Math\Matrix\Computation\Mul;
use Chippyash\Math\Matrix\Computation\AbstractComputation;
use Chippyash\Matrix\Transformation\Colslice;
use Chippyash\Math\Matrix\NumericMatrix;
use Chippyash\Math\Matrix\Special\Zeros as ZMatrix;
use Chippyash\Math\Matrix\Exceptions\ComputationException;
use Chippyash\Math\Matrix\Exceptions\UndefinedComputationException;
use Chippyash\Math\Matrix\Traits\CreateCorrectMatrixType;
use Chippyash\Math\Matrix\Traits\CreateCorrectScalarType;
use Chippyash\Math\Matrix\Traits\AssertMatrixIsNumeric;
use Chippyash\Matrix\Traits\AssertParameterIsMatrix;
use Chippyash\Type\TypeFactory;
use Chippyash\Math\Type\Calculator;
class Matrix extends AbstractComputation
{
use CreateCorrectMatrixType;
use CreateCorrectScalarType;
use AssertMatrixIsNumeric;
use AssertParameterIsMatrix;
public function compute(NumericMatrix $mA, $extra = null)
{
$this->assertParameterIsMatrix($extra, 'Parameter is not a matrix')
->assertMatrixIsNumeric($extra, 'Parameter is not a numeric matrix');
if ($mA->is('empty') || $extra->is('empty')) {
return $this->createCorrectMatrixType($mA, []);
}
$this->checkCompatibility($mA, $extra);
$product = $this->doComputation($mA, $extra);
return $this->massageProduct($mA, $extra, $product);
}
protected function doComputation(NumericMatrix $mA, NumericMatrix $mB)
{
$size = max([$mA->columns(), $mA->rows(), $mB->columns(), $mB->rows()]);
$product = (new ZMatrix())->create([$size, $size])->toArray();
$dA = $mA->toArray();
$dB = $mB->toArray();
$zero = $this->createCorrectScalarType($mA, 0);
$calc = new Calculator();
for ($i = 0; $i < $size; $i++) {
for ($k = 0; $k < $size; $k++) {
for ($j = 0; $j < $size; $j++) {
$a = isset($dA[$i][$k]) ? $dA[$i][$k] : $zero;
$b = isset($dB[$k][$j]) ? $dB[$k][$j] : $zero;
$product[$i][$j] = $calc->add($product[$i][$j], $calc->mul($a, $b));
}
}
}
return $product;
}
protected function checkCompatibility(NumericMatrix $mA, NumericMatrix $mB)
{
if ($mA->is('rowvector')) {
$this->checkRowVectorCompatibility($mA, $mB);
}
if ($mA->is('columnvector')) {
$this->checkColumnVectorCompatibility($mA, $mB);
}
if ($mA->is('square')) {
$this->checkSquareMatrixCompatibility($mA, $mB);
}
if ($mA->is('rectangle')) {
$this->checkRectangleMatrixCompatibility($mA, $mB);
}
}
protected function checkRowVectorCompatibility(NumericMatrix $mRV, NumericMatrix $mB)
{
if ($mB->is('columnvector') && ($mRV->columns() != $mB->rows())) {
throw new ComputationException('Two matrices cannot be multiplied: mA->columns != mB->rows');
}
}
protected function checkColumnVectorCompatibility(NumericMatrix $mCV, NumericMatrix $mB)
{
if ($mB->is('rowvector') && ($mCV->rows() != $mB->columns())) {
throw new ComputationException('Two matrices cannot be multiplied: mA->rows != mB->columns');
}
if ($mB->is('square')) {
throw new UndefinedComputationException('columnVector x square');
}
}
protected function checkSquareMatrixCompatibility(NumericMatrix $mSq, NumericMatrix $mB)
{
if ($mB->is('columnvector') && ($mSq->columns() != $mB->rows())) {
throw new ComputationException('Two matrices cannot be multiplied: mA->columns != mB->rows');
}
if ($mB->is('square') && ($mSq->vertices() != $mB->vertices())) {
throw new ComputationException('Two matrices cannot be multiplied: mA->columns != mB->rows');
}
}
protected function checkRectangleMatrixCompatibility(NumericMatrix $mA, NumericMatrix $mB)
{
if ($mB->is('rectangle') && ($mA->columns() != $mB->rows())) {
throw new ComputationException('Two matrices cannot be multiplied: mA->columns != mB->rows');
}
}
protected function massageProduct(NumericMatrix $mA, NumericMatrix $mB, array $product)
{
if ($mA->is('rowvector') && $mB->is('columnvector')) {
return $this->createCorrectMatrixType($mA, [$product[0][0]]);
}
if ($mA->is('square')) {
return $this->massageSquare($mA, $mB, $product);
}
if ($mA->is('rectangle')) {
return $this->massageRectangle($mA, $mB, $product);
}
return $this->createCorrectMatrixType($mA, $product);
}
protected function massageSquare(NumericMatrix $mA, NumericMatrix $mB, array $product)
{
if ($mB->is('columnvector')) {
$fC = new Colslice();
return $fC($this->createCorrectMatrixType($mA, $product), array(1, 1));
}
return $this->createCorrectMatrixType($mA, $product);
}
protected function massageRectangle(NumericMatrix $mA, NumericMatrix $mB, array $product)
{
if ($mB->is('rectangle') && ($mA->rows() < $mB->rows())) {
return $this->createCorrectMatrixType(
$mA,
[[$product[0][0], $product[0][1]],
[$product[1][0], $product[1][1]]
]
);
}
return $this->createCorrectMatrixType($mA, $product);
}
}