Jingga/phpOMS
1
0
Fork 0
mirror of https://github.com/Karaka-Management/phpOMS.git synced 2026-01-30 17:58:39 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

fixes #158

Browse Source
This commit is contained in:
Dennis Eichhorn 2018-08-09 18:43:03 +02:00
parent 175be313d2
commit 24659e324f
2 changed files with 83 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
tests/Math/Matrix/CholeskyDecompositionTest.php
View File

@ -36,11 +36,23 @@ class CholeskyDecompositionTest extends \PHPUnit\Framework\TestCase
[-1, 1, 3], [-1, 1, 3],
], $cholesky->getL()->toArray(), '', 0.2); ], $cholesky->getL()->toArray(), '', 0.2);
self::assertTrue($cholesky->isSpd());
}
public function testSolve()
{
$A = new Matrix();
$A->setMatrix([
[25, 15, -5],
[15, 17, 0],
[-5, 0, 11],
]);
$cholesky = new CholeskyDecomposition($A);
$vec = new Vector(); $vec = new Vector();
$vec->setMatrix([[40], [49], [28]]); $vec->setMatrix([[40], [49], [28]]);
self::assertEquals([[1], [2], [3]], $cholesky->solve($vec)->toArray(), '', 0.2); self::assertEquals([[1], [2], [3]], $cholesky->solve($vec)->toArray(), '', 0.2);
self::assertTrue($cholesky->isSpd());
} }
/** /**

69
tests/Math/Matrix/QRDecompositionTest.php Normal file
View File

@ -0,0 +1,69 @@
<?php
/**
* Orange Management
*
* PHP Version 7.2
*
* @package TBD
* @copyright Dennis Eichhorn
* @license OMS License 1.0
* @version 1.0.0
* @link http://website.orange-management.de
*/
namespace phpOMS\tests\Math\Matrix;
use phpOMS\Math\Matrix\Matrix;
use phpOMS\Math\Matrix\Vector;
use phpOMS\Math\Matrix\QRDecomposition;
class QRDecompositionTest extends \PHPUnit\Framework\TestCase
{
/**public function testDecomposition()
{
$A = new Matrix();
$A->setMatrix([
[12, -51, 4],
[6, 167, -68],
[-4, 24, -41],
]);
$QR = new QRDecomposition($A);
self::assertTrue($QR->isFullRank());
self::assertEquals([
[12, -69, -58 / 5],
[6, 158, 6 / 5],
[-4, 30, -33],
], $QR->getH()->toArray(), '', 0.2);
self::assertEquals([
[6 / 7, -69 / 175, -58 / 175],
[3 / 7, 158 / 175, 6 / 175],
[-2 / 7, 6 / 35, -33 / 35],
], $QR->getQ()->toArray(), '', 0.2);
self::assertEquals([
[14, 21, -14],
[0, 175, -70],
[0, 0, 35],
], $QR->getR()->toArray(), '', 0.2);
}*/
public function testSolve()
{
$A = new Matrix();
$A->setMatrix([
[25, 15, -5],
[15, 17, 0],
[-5, 0, 11],
]);
$QR = new QRDecomposition($A);
$vec = new Vector();
$vec->setMatrix([[40], [49], [28]]);
self::assertEquals([[1], [2], [3]], $QR->solve($vec)->toArray(), '', 0.2);
}
}