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:
<?php
namespace Chippyash\Math\Matrix\Derivative;
use Chippyash\Math\Matrix\Derivative\AbstractDerivative;
use Chippyash\Math\Matrix\NumericMatrix;
use Chippyash\Math\Matrix\Exceptions\UndefinedComputationException;
use Chippyash\Matrix\Traits\AssertMatrixIsSquare;
use Chippyash\Math\Matrix\Derivative\Strategy\Determinant\Laplace;
use Chippyash\Math\Matrix\Derivative\Strategy\Determinant\Lu;
use Chippyash\Math\Matrix\Interfaces\TuningInterface;
use Chippyash\Type\String\StringType;
class Determinant extends AbstractDerivative implements TuningInterface
{
use AssertMatrixIsSquare;
const METHOD_AUTO = 0;
const METHOD_LAPLACE = 1;
const METHOD_LU = 2;
static protected $luLimit = 20;
protected $method;
public function __construct($method = self::METHOD_AUTO)
{
$this->method = $method;
}
public function derive(NumericMatrix $mA, $extra = null)
{
$this->assertMatrixIsSquare($mA, 'No determinant for non-square matrix');
return $this->getDeterminant($mA);
}
public function tune(StringType $name, $value)
{
if ($name() != 'luLimit') {
throw new \InvalidArgumentException("{$name} is unknown for tuning");
}
$ret = self::$luLimit;
self::$luLimit = $value;
return $ret;
}
protected function getDeterminant(NumericMatrix $mA)
{
switch ($this->method) {
case self::METHOD_AUTO;
if ($mA->rows() <= self::$luLimit) {
$strategy = new Lu();
} else {
throw new UndefinedComputationException('No available strategy found to determine the determinant');
}
break;
case self::METHOD_LAPLACE:
$strategy = new Laplace();
break;
case self::METHOD_LU;
$strategy = new Lu();
break;
default:
throw new UndefinedComputationException('Unknown determinant computation method');
}
return $strategy->determinant($mA);
}
}