@@ -167,10 +167,10 @@ |
||
167 | 167 | self::init(); |
168 | 168 | if($m->dtype == \Np\matrix::FLOAT){ |
169 | 169 | return self::$ffi_lapack->LAPACKE_sgels( $matLayout, $trans, $m->row, $m->col, $b->col, $m->data, |
170 | - $m->col, $b->data, $b->col ); |
|
170 | + $m->col, $b->data, $b->col ); |
|
171 | 171 | } |
172 | 172 | return self::$ffi_lapack->LAPACKE_dgels( $matLayout, $trans, $m->row, $m->col, $b->col, $m->data, |
173 | - $m->col, $b->data, $b->col ); |
|
173 | + $m->col, $b->data, $b->col ); |
|
174 | 174 | } |
175 | 175 | |
176 | 176 | } |
@@ -163,14 +163,14 @@ |
||
163 | 163 | * @param string $trans |
164 | 164 | * @return type |
165 | 165 | */ |
166 | - public static function gels(\Np\matrix $m, \Np\matrix|\Np\vector $b, int $matLayout = self::ROW_MAJOR,string $trans = 'N') { |
|
166 | + public static function gels(\Np\matrix $m, \Np\matrix | \Np\vector $b, int $matLayout = self::ROW_MAJOR, string $trans = 'N') { |
|
167 | 167 | self::init(); |
168 | - if($m->dtype == \Np\matrix::FLOAT){ |
|
169 | - return self::$ffi_lapack->LAPACKE_sgels( $matLayout, $trans, $m->row, $m->col, $b->col, $m->data, |
|
170 | - $m->col, $b->data, $b->col ); |
|
168 | + if ($m->dtype == \Np\matrix::FLOAT) { |
|
169 | + return self::$ffi_lapack->LAPACKE_sgels($matLayout, $trans, $m->row, $m->col, $b->col, $m->data, |
|
170 | + $m->col, $b->data, $b->col); |
|
171 | 171 | } |
172 | - return self::$ffi_lapack->LAPACKE_dgels( $matLayout, $trans, $m->row, $m->col, $b->col, $m->data, |
|
173 | - $m->col, $b->data, $b->col ); |
|
172 | + return self::$ffi_lapack->LAPACKE_dgels($matLayout, $trans, $m->row, $m->col, $b->col, $m->data, |
|
173 | + $m->col, $b->data, $b->col); |
|
174 | 174 | } |
175 | 175 | |
176 | 176 | } |
@@ -24,7 +24,7 @@ discard block |
||
24 | 24 | * @param \Np\matrix|\Np\vector $d |
25 | 25 | * @return matrix|vector |
26 | 26 | */ |
27 | - public function dot(matrix|vector $d): matrix|vector { |
|
27 | + public function dot(matrix | vector $d): matrix | vector { |
|
28 | 28 | if ($this instanceof matrix) { |
29 | 29 | if ($d instanceof self) { |
30 | 30 | return $this->dotMatrix($d); |
@@ -44,7 +44,7 @@ discard block |
||
44 | 44 | * @return \Np\matrix |
45 | 45 | */ |
46 | 46 | protected function dotMatrix(matrix $matrix): matrix { |
47 | - if ($this->checkDtype($this, $matrix) && $this->checkDimensions($this,$matrix)) { |
|
47 | + if ($this->checkDtype($this, $matrix) && $this->checkDimensions($this, $matrix)) { |
|
48 | 48 | $ar = self::factory($this->row, $matrix->col, $this->dtype); |
49 | 49 | blas::gemm($this, $matrix, $ar); |
50 | 50 | return $ar; |
@@ -69,7 +69,7 @@ discard block |
||
69 | 69 | * Compute the multiplicative inverse of the matrix. |
70 | 70 | * @return matrix|null |
71 | 71 | */ |
72 | - public function inverse(): matrix|null { |
|
72 | + public function inverse(): matrix | null { |
|
73 | 73 | if (!$this->isSquare()) { |
74 | 74 | self::_err('Error::invalid Size of matrix!'); |
75 | 75 | } |
@@ -93,7 +93,7 @@ discard block |
||
93 | 93 | * Compute the (Moore-Penrose) pseudo inverse of the general matrix. |
94 | 94 | * @return matrix|null |
95 | 95 | */ |
96 | - public function pseudoInverse(): matrix|null { |
|
96 | + public function pseudoInverse(): matrix | null { |
|
97 | 97 | $k = min($this->row, $this->col); |
98 | 98 | $s = vector::factory($k, $this->dtype); |
99 | 99 | $u = self::factory($this->row, $this->row, $this->dtype); |
@@ -24,9 +24,9 @@ discard block |
||
24 | 24 | * @email [email protected] |
25 | 25 | * @copyright (c) 2020-2021, Shubham Chaudhary |
26 | 26 | */ |
27 | -class matrix extends nd{ |
|
27 | +class matrix extends nd { |
|
28 | 28 | |
29 | - use ops,linAlg; |
|
29 | + use ops, linAlg; |
|
30 | 30 | /** |
31 | 31 | * create empty 2d matrix for given data type |
32 | 32 | * @param int $row num of rows |
@@ -310,7 +310,7 @@ discard block |
||
310 | 310 | * @param int|float|matrix|vector $m |
311 | 311 | * @return matrix|vector |
312 | 312 | */ |
313 | - public function multiply(int|float|matrix|vector $m): matrix|vector { |
|
313 | + public function multiply(int | float | matrix | vector $m): matrix | vector { |
|
314 | 314 | if ($m instanceof self) { |
315 | 315 | return $this->multiplyMatrix($m); |
316 | 316 | } else if ($m instanceof vector) { |
@@ -376,7 +376,7 @@ discard block |
||
376 | 376 | * @param int|float|matrix|vector $m |
377 | 377 | * @return matrix |
378 | 378 | */ |
379 | - public function sum(int|float|matrix|vector $m): matrix { |
|
379 | + public function sum(int | float | matrix | vector $m): matrix { |
|
380 | 380 | if ($m instanceof self) { |
381 | 381 | return $this->sumMatrix($m); |
382 | 382 | } elseif ($m instanceof vector) { |
@@ -386,7 +386,7 @@ discard block |
||
386 | 386 | } |
387 | 387 | } |
388 | 388 | |
389 | - protected function sumScalar(int|float $s): matrix { |
|
389 | + protected function sumScalar(int | float $s): matrix { |
|
390 | 390 | $ar = self::factory($this->row, $this->col, $this->dtype); |
391 | 391 | for ($i = 0; $i < $this->ndim; ++$i) { |
392 | 392 | $ar->data[$i] = $this->data[$i] + $s; |
@@ -421,7 +421,7 @@ discard block |
||
421 | 421 | * @param int|float|matrix|vector $d matrix|$scalar to subtract this matrix |
422 | 422 | * @return \Np\matrix |
423 | 423 | */ |
424 | - public function subtract(int|float|matrix|vector $d): matrix { |
|
424 | + public function subtract(int | float | matrix | vector $d): matrix { |
|
425 | 425 | if ($d instanceof self) { |
426 | 426 | return $this->subtractMatrix($d); |
427 | 427 | } elseif ($d instanceof vector) { |
@@ -431,7 +431,7 @@ discard block |
||
431 | 431 | } |
432 | 432 | } |
433 | 433 | |
434 | - protected function subtractScalar(int|float $s): matrix { |
|
434 | + protected function subtractScalar(int | float $s): matrix { |
|
435 | 435 | $ar = self::factory($this->row, $this->col, $this->dtype); |
436 | 436 | for ($i = 0; $i < $this->ndim; ++$i) { |
437 | 437 | $ar->data[$i] = $this->data[$i] - $s; |
@@ -460,7 +460,7 @@ discard block |
||
460 | 460 | * @return matrix |
461 | 461 | */ |
462 | 462 | protected function subtractVector(vector $v): matrix { |
463 | - if ($this->checkDimensions($v, $this) && $this->checkDtype($this,$v)) { |
|
463 | + if ($this->checkDimensions($v, $this) && $this->checkDtype($this, $v)) { |
|
464 | 464 | $ar = self::factory($this->row, $this->col, $this->dtype); |
465 | 465 | for ($i = 0; $i < $this->row; ++$i) { |
466 | 466 | for ($j = 0; $j < $this->col; ++$j) { |
@@ -493,7 +493,7 @@ discard block |
||
493 | 493 | * @param int|float|matrix $d |
494 | 494 | * @return matrix |
495 | 495 | */ |
496 | - public function divide(int|float|matrix|vector $d): matrix { |
|
496 | + public function divide(int | float | matrix | vector $d): matrix { |
|
497 | 497 | if ($d instanceof self) { |
498 | 498 | return $this->divideMatrix($d); |
499 | 499 | } elseif ($d instanceof vector) { |
@@ -525,7 +525,7 @@ discard block |
||
525 | 525 | } |
526 | 526 | } |
527 | 527 | |
528 | - protected function divideScalar(int|float $s): matrix { |
|
528 | + protected function divideScalar(int | float $s): matrix { |
|
529 | 529 | $ar = self::factory($this->row, $this->col, $this->dtype); |
530 | 530 | for ($i = 0; $i < $this->ndim; ++$i) { |
531 | 531 | $ar->data[$i] = $this->data[$i] / $s; |
@@ -540,7 +540,7 @@ discard block |
||
540 | 540 | * @param int|float|matrix $m |
541 | 541 | * @return matrix |
542 | 542 | */ |
543 | - public function pow(int|float|matrix|vector $d): matrix { |
|
543 | + public function pow(int | float | matrix | vector $d): matrix { |
|
544 | 544 | if ($d instanceof self) { |
545 | 545 | return $this->powMatrix($d); |
546 | 546 | } else if ($d instanceof vector) { |
@@ -572,7 +572,7 @@ discard block |
||
572 | 572 | } |
573 | 573 | } |
574 | 574 | |
575 | - protected function powScalar(int|float $s): matrix { |
|
575 | + protected function powScalar(int | float $s): matrix { |
|
576 | 576 | $ar = $this->copy(); |
577 | 577 | for ($i = 0; $i < $this->ndim; ++$i) { |
578 | 578 | $ar->data[$i] **= $s; |
@@ -585,7 +585,7 @@ discard block |
||
585 | 585 | * @param int|float|matrix|vector $d |
586 | 586 | * @return matrix |
587 | 587 | */ |
588 | - public function mod(int|float|matrix|vector $d): matrix { |
|
588 | + public function mod(int | float | matrix | vector $d): matrix { |
|
589 | 589 | if ($d instanceof self) { |
590 | 590 | $this->modMatrix($d); |
591 | 591 | } else if ($d instanceof vector) { |
@@ -617,7 +617,7 @@ discard block |
||
617 | 617 | } |
618 | 618 | } |
619 | 619 | |
620 | - protected function modScalar(int|float $s): matrix { |
|
620 | + protected function modScalar(int | float $s): matrix { |
|
621 | 621 | $ar = $this->copy(); |
622 | 622 | for ($i = 0; $i < $this->ndim; ++$i) { |
623 | 623 | $ar->data[$i] %= $s; |
@@ -729,7 +729,7 @@ discard block |
||
729 | 729 | * @param int|float $scalar |
730 | 730 | * @return matrix |
731 | 731 | */ |
732 | - public function scale(int|float $scalar): matrix { |
|
732 | + public function scale(int | float $scalar): matrix { |
|
733 | 733 | if ($scalar == 0) { |
734 | 734 | return self::zeros($this->row, $this->col, $this->dtype); |
735 | 735 | } |
@@ -747,7 +747,7 @@ discard block |
||
747 | 747 | * @param int $row |
748 | 748 | * @param int|float $c |
749 | 749 | */ |
750 | - public function scaleRow(int $row, int|float $c) { |
|
750 | + public function scaleRow(int $row, int | float $c) { |
|
751 | 751 | for ($i = 0; $i < $this->col; ++$i) { |
752 | 752 | $this->data[$row * $this->col + $i] *= $c; |
753 | 753 | } |
@@ -758,7 +758,7 @@ discard block |
||
758 | 758 | * @param int $col |
759 | 759 | * @param int|float $c |
760 | 760 | */ |
761 | - public function scaleCol(int $col, int|float $c) { |
|
761 | + public function scaleCol(int $col, int | float $c) { |
|
762 | 762 | for ($i = 0; $i < $this->row; ++$i) { |
763 | 763 | $this->data[$i * $this->col + $col] *= $c; |
764 | 764 | } |
@@ -769,14 +769,14 @@ discard block |
||
769 | 769 | * @param int|float $c |
770 | 770 | * @param bool $lDig |
771 | 771 | */ |
772 | - public function scaleDigonalCol(int|float $c, bool $lDig = true) { |
|
773 | - if($lDig){ |
|
774 | - for ($i = 0; $i < $this->row ; ++$i) { |
|
772 | + public function scaleDigonalCol(int | float $c, bool $lDig = true) { |
|
773 | + if ($lDig) { |
|
774 | + for ($i = 0; $i < $this->row; ++$i) { |
|
775 | 775 | $this->data[$i * $this->col + $i] *= $c; |
776 | 776 | } |
777 | 777 | } |
778 | - else{ |
|
779 | - for ($i = $this->row; $i > 0 ; --$i) { |
|
778 | + else { |
|
779 | + for ($i = $this->row; $i > 0; --$i) { |
|
780 | 780 | $this->data[$i * $this->col - $i] *= $c; |
781 | 781 | } |
782 | 782 | } |
@@ -823,7 +823,7 @@ discard block |
||
823 | 823 | * @return matrix |
824 | 824 | */ |
825 | 825 | public function joinRight(matrix $m): matrix { |
826 | - if ($this->row != $m->row && !$this->checkDtype($this,$m)) { |
|
826 | + if ($this->row != $m->row && !$this->checkDtype($this, $m)) { |
|
827 | 827 | self::_err('Error::Invalid size! or DataType!'); |
828 | 828 | } |
829 | 829 | $col = $this->col + $m->col; |
@@ -889,7 +889,7 @@ discard block |
||
889 | 889 | * |
890 | 890 | * @return matrix|null |
891 | 891 | */ |
892 | - public function ref(): matrix|null { |
|
892 | + public function ref(): matrix | null { |
|
893 | 893 | return ref::factory($this); |
894 | 894 | } |
895 | 895 | |
@@ -898,7 +898,7 @@ discard block |
||
898 | 898 | * |
899 | 899 | * @return matrix|null |
900 | 900 | */ |
901 | - public function cholesky(): matrix|null { |
|
901 | + public function cholesky(): matrix | null { |
|
902 | 902 | return cholesky::factory($this); |
903 | 903 | } |
904 | 904 | |
@@ -973,7 +973,7 @@ discard block |
||
973 | 973 | * @param bool $dignoal |
974 | 974 | * @return void |
975 | 975 | */ |
976 | - public function setData(int|float|array $data): void { |
|
976 | + public function setData(int | float | array $data): void { |
|
977 | 977 | |
978 | 978 | if (is_array($data) && is_array($data[0])) { |
979 | 979 | $f = $this->flattenArray($data); |
@@ -1004,7 +1004,7 @@ discard block |
||
1004 | 1004 | * @return object |
1005 | 1005 | */ |
1006 | 1006 | public function getShape(): object { |
1007 | - return (object) ['m' => $this->row, 'n' => $this->col]; |
|
1007 | + return (object)['m' => $this->row, 'n' => $this->col]; |
|
1008 | 1008 | } |
1009 | 1009 | |
1010 | 1010 | /** |
@@ -1186,7 +1186,7 @@ discard block |
||
1186 | 1186 | * @param vector|null $mean |
1187 | 1187 | * @return vector |
1188 | 1188 | */ |
1189 | - public function variance(vector|null $mean = null): vector { |
|
1189 | + public function variance(vector | null $mean = null): vector { |
|
1190 | 1190 | if (isset($mean)) { |
1191 | 1191 | if (!$mean instanceof vector) { |
1192 | 1192 | self::_invalidArgument('mean must be a vector!'); |
@@ -1228,7 +1228,7 @@ discard block |
||
1228 | 1228 | * @param vector|null $mean |
1229 | 1229 | * @return matrix |
1230 | 1230 | */ |
1231 | - public function covariance(vector|null $mean = null): matrix { |
|
1231 | + public function covariance(vector | null $mean = null): matrix { |
|
1232 | 1232 | if (isset($mean)) { |
1233 | 1233 | if ($mean->col !== $this->row) { |
1234 | 1234 | self::_err('Err:: given mean vector dimensionality mismatched!'); |
@@ -1256,7 +1256,7 @@ discard block |
||
1256 | 1256 | * @param int|float|matrix|vector $d |
1257 | 1257 | * @return matrix |
1258 | 1258 | */ |
1259 | - public function equal(int|float|matrix|vector $d): matrix { |
|
1259 | + public function equal(int | float | matrix | vector $d): matrix { |
|
1260 | 1260 | if ($d instanceof self) { |
1261 | 1261 | return $this->equalMatrix($d); |
1262 | 1262 | } elseif ($d instanceof vector) { |
@@ -1288,7 +1288,7 @@ discard block |
||
1288 | 1288 | } |
1289 | 1289 | } |
1290 | 1290 | |
1291 | - protected function equalScalar(int|float $s): matrix { |
|
1291 | + protected function equalScalar(int | float $s): matrix { |
|
1292 | 1292 | $ar = self::factory($this->row, $this->col, $this->dtype); |
1293 | 1293 | for ($i = 0; $i < $this->ndim; ++$i) { |
1294 | 1294 | $ar->data[$i] = $this->data[$i] == $s ? 1 : 0; |
@@ -1301,7 +1301,7 @@ discard block |
||
1301 | 1301 | * @param int|float|matrix|vector $d |
1302 | 1302 | * @return matrix |
1303 | 1303 | */ |
1304 | - public function greater(int|float|matrix|vector $d): matrix { |
|
1304 | + public function greater(int | float | matrix | vector $d): matrix { |
|
1305 | 1305 | if ($d instanceof self) { |
1306 | 1306 | return $this->greaterMatrix($d); |
1307 | 1307 | } elseif ($d instanceof vector) { |
@@ -1312,7 +1312,7 @@ discard block |
||
1312 | 1312 | } |
1313 | 1313 | |
1314 | 1314 | protected function greaterMatrix(matrix $m): matrix { |
1315 | - if ($this->checkShape($this, $m) && $this->checkDtype($this,$m)) { |
|
1315 | + if ($this->checkShape($this, $m) && $this->checkDtype($this, $m)) { |
|
1316 | 1316 | $ar = self::factory($this->row, $this->col, $this->dtype); |
1317 | 1317 | for ($i = 0; $i < $this->ndim; ++$i) { |
1318 | 1318 | $ar->data[$i] = $this->data[$i] > $m->data[$i] ? 1 : 0; |
@@ -1322,7 +1322,7 @@ discard block |
||
1322 | 1322 | } |
1323 | 1323 | |
1324 | 1324 | protected function greaterVector(vector $v): matrix { |
1325 | - if ($this->checkDimensions($v, $this) && $this->checkDtype($this,$v)) { |
|
1325 | + if ($this->checkDimensions($v, $this) && $this->checkDtype($this, $v)) { |
|
1326 | 1326 | $ar = self::factory($this->row, $this->col, $this->dtype); |
1327 | 1327 | for ($i = 0; $i < $this->row; ++$i) { |
1328 | 1328 | for ($j = 0; $j < $this->col; ++$j) { |
@@ -1333,7 +1333,7 @@ discard block |
||
1333 | 1333 | } |
1334 | 1334 | } |
1335 | 1335 | |
1336 | - protected function greaterScalar(int|float $s): matrix { |
|
1336 | + protected function greaterScalar(int | float $s): matrix { |
|
1337 | 1337 | $ar = self::factory($this->row, $this->col, $this->dtype); |
1338 | 1338 | for ($i = 0; $i < $this->ndim; ++$i) { |
1339 | 1339 | $ar->data[$i] = $this->data[$i] > $s ? 1 : 0; |
@@ -1346,10 +1346,10 @@ discard block |
||
1346 | 1346 | * @param int|float|matrix $m |
1347 | 1347 | * @return matrix |
1348 | 1348 | */ |
1349 | - public function less(int|float|matrix $m): matrix { |
|
1349 | + public function less(int | float | matrix $m): matrix { |
|
1350 | 1350 | $ar = self::factory($this->row, $this->col, $this->dtype); |
1351 | 1351 | if ($m instanceof self) { |
1352 | - if ($this->checkShape($this,$m)) { |
|
1352 | + if ($this->checkShape($this, $m)) { |
|
1353 | 1353 | for ($i = 0; $i < $this->ndim; ++$i) { |
1354 | 1354 | $ar->data[$i] = $this->data[$i] < $m->data[$i] ? 1 : 0; |
1355 | 1355 | } |
@@ -1397,7 +1397,7 @@ discard block |
||
1397 | 1397 | } |
1398 | 1398 | |
1399 | 1399 | public function __toString() { |
1400 | - return (string) $this->printMatrix(); |
|
1400 | + return (string)$this->printMatrix(); |
|
1401 | 1401 | } |
1402 | 1402 | |
1403 | 1403 | private function flattenArray(array $ar) { |