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: 241: 242: 243: 244: 245: 246: 247: 248: 249: 250: 251: 252: 253: 254: 255: 256: 257: 258: 259: 260: 261: 262: 263: 264: 265: 266: 267: 268: 269: 270: 271: 272: 273: 274: 275: 276: 277: 278: 279: 280: 281: 282: 283: 284:
<?php
namespace Chippyash\Math\Matrix;
use Chippyash\Matrix\Matrix;
use Chippyash\Math\Matrix\Traits\ConvertNumberToNumeric;
use Chippyash\Math\Matrix\Exceptions\MathMatrixException;
use Chippyash\Math\Matrix\Interfaces\ComputatationInterface;
use Chippyash\Math\Matrix\Interfaces\DerivativeInterface;
use Chippyash\Math\Matrix\Interfaces\DecompositionInterface;
use Chippyash\Matrix\Interfaces\TransformationInterface;
use Chippyash\Type\Number\Rational\RationalTypeFactory;
use Chippyash\Type\TypeFactory;
use Chippyash\Type\Interfaces\NumericTypeInterface;
use Chippyash\Math\Type\Comparator;
class NumericMatrix extends Matrix
{
use ConvertNumberToNumeric;
const NS_NUMERIC_ATTRIBUTE = 'Chippyash\Math\Matrix\Attribute\Is';
const NS_COMPUTATION = 'Chippyash\Math\Matrix\Computation\\';
const NS_NTRANSFORMATION = 'Chippyash\Math\Matrix\Transformation\\';
const NS_DERIVATIVE = 'Chippyash\Math\Matrix\Derivative\\';
const NS_DECOMPOSITION = 'Chippyash\Math\Matrix\Decomposition\\';
public function __construct($source, $normalizeDefault = 0)
{
if ($source instanceof self) {
$this->store($source->toArray());
return;
}
if (is_array($source)) {
if (is_int($normalizeDefault)) {
$default = TypeFactory::createInt($normalizeDefault);
} elseif (is_float($normalizeDefault)) {
$default = RationalTypeFactory::fromFloat($normalizeDefault);
} elseif (!$normalizeDefault instanceof NumericTypeInterface) {
throw new MathMatrixException('NumericMatrix expects numeric default value');
} else {
$default = $normalizeDefault;
}
parent::__construct($source, false, true, $default);
} else {
throw new MathMatrixException('NumericMatrix expects NumericMatrix or array as source data');
}
}
public function test($attribute)
{
if (is_string($attribute)) {
$attribute = ucfirst(strtolower($attribute));
$class = self::NS_NUMERIC_ATTRIBUTE . $attribute;
if (class_exists($class)) {
$obj = new $class();
} else {
return parent::test($attribute);
}
} else {
$obj = $attribute;
}
return parent::test($obj);
}
public function compute(ComputatationInterface $computation, $extra = null)
{
return $computation->compute($this, $extra);
}
public function derive(DerivativeInterface $derivative, $extra = null)
{
return $derivative->derive($this, $extra);
}
public function decompose(DecompositionInterface $decomposition, $extra = null)
{
return $decomposition->decompose($this, $extra);
}
public function transform(TransformationInterface $transformation, $extra = null)
{
return new self(parent::transform($transformation, $extra)->toArray());
}
public function __invoke()
{
$numArgs = func_num_args();
if ($numArgs == 1) {
$operationName = func_get_arg(0);
$extra = null;
} elseif($numArgs == 2) {
$operationName = func_get_arg(0);
$extra = func_get_arg(1);
} else {
throw new \InvalidArgumentException(self::ERR_INVALID_INVOKE_ARG);
}
$cName = self::NS_COMPUTATION . $operationName;
if (class_exists($cName, true)) {
return $this->compute(new $cName(), $extra);
}
$tName = self::NS_NTRANSFORMATION . $operationName;
if (class_exists($tName, true)) {
return $this->transform(new $tName(), $extra);
}
$dName = self::NS_DERIVATIVE . $operationName;
if (class_exists($dName, true)) {
return $this->derive(new $dName(), $extra);
}
$dcName = self::NS_DECOMPOSITION . $operationName;
if (class_exists($dcName, true)) {
return $this->decompose(new $dcName(), $extra);
}
$tName = self::NS_TRANSFORMATION . $operationName;
if (class_exists($tName, true)) {
return $this->transform(new $tName(), $extra);
}
throw new \InvalidArgumentException(self::ERR_INVALID_OP_NAME);
}
protected function checkEntryEquality(Matrix $mB, $strict)
{
if ($strict) {
if (get_class($this) !== get_class($mB)) {
return false;
}
}
$dA = $this->toArray();
$dB = $mB->toArray();
$m = $this->rows();
$n = $this->columns();
$comp = new Comparator();
for ($i=0; $i<$m; $i++) {
for ($j=0; $j<$n; $j++) {
if ($strict) {
if ($dA[$i][$j] !== $dB[$i][$j]) {
return false;
}
} else {
if ($comp->neq($dA[$i][$j], $dB[$i][$j])) {
return false;
}
}
}
}
return true;
}
protected function store(array $data) {
foreach ($data as &$row) {
foreach ($row as &$item) {
$item = $this->convertNumberToNumeric($item);
}
}
$this->data = $data;
}
}