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: 

<?php
/**
 * Markov
 *
 * @author Ashley Kitson
 * @copyright Ashley Kitson, 2016, UK
 * @license GPL V3+ See LICENSE.md
 */

namespace Chippyash\Math\Matrix\Attribute;


use Chippyash\Math\Matrix\Derivative\Sum;
use Chippyash\Math\Type\Comparator;
use Chippyash\Matrix\Exceptions\MatrixException;
use Chippyash\Matrix\Interfaces\AttributeInterface;
use Chippyash\Matrix\Matrix;
use Chippyash\Math\Matrix\Traits\AssertMatrixIsNumeric;
use Chippyash\Matrix\Traits\AssertMatrixIsSquare;

/**
 * IsMarkov Matrix Attribute
 * Is it A NumericMatrix?
 * Is it Square?
 * Are all row sums equal i.e. is the probability graph complete?
 * Are there at least 2 rows and columns?
 */
class IsMarkov implements AttributeInterface
{
    use AssertMatrixIsNumeric;
    use AssertMatrixIsSquare;

    /**
     * Does the matrix exhibit this attribute
     *
     * @param Matrix $mA
     *
     * @return boolean
     */
    public function is(Matrix $mA)
    {
        try {
            $this->assertMatrixIsNumeric($mA)
                 ->assertMatrixIsSquare($mA);
        } catch (MatrixException $e) {
            return false;
        }
        //we've already assured it is square, now make sure we at least have a
        // cyclic probability
        if ($mA->rows() < 2 || $mA->columns() < 2) {
            return false;
        }

        $sumDerivative = new Sum();
        $firstRow = $mA('Rowslice', [1])->derive($sumDerivative);
        $comp = new Comparator();
        //check that each row has same sum
        foreach(range(2, $mA->rows()) as $row) {
            if ($comp->neq($firstRow, $mA('RowSlice', [$row])->derive($sumDerivative)) != 0) {
                return false;
            }
        }

        return true;
    }
}