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: 

<?php
/*
 * Math-Matrix library
 *
 * @author Ashley Kitson <akitson@zf4.biz>
 * @copyright Ashley Kitson, UK, 2014
 * @licence GPL V3 or later : http://www.gnu.org/licenses/gpl.html
 * @link http://en.wikipedia.org/wiki/Matrix_(mathematics)
 */
namespace Chippyash\Math\Matrix;

use Chippyash\Type\Number\Complex\ComplexType;
use Chippyash\Type\Number\Complex\ComplexTypeFactory;
use Chippyash\Type\Number\Rational\RationalType;
use Chippyash\Type\Number\Rational\RationalTypeFactory;
use Chippyash\Math\Matrix\Exceptions\MathMatrixException;
use Chippyash\Type\Number\IntType;

/**
 * Static factory to create the various standard numerical matrices
 *
 */
abstract class MatrixFactory
{

    /**
     *
     * @param string $type
     * @param array $data
     * @return \Chippyash\Math\Matrix\NumericMatrix
     */
    public static function create($type, array $data)
    {
        switch (strtolower($type)) {
            case 'complex':
                return self::createComplex($data);
            case 'rational':
                return self::createRational($data);
            case 'numeric':
            default:
                return self::createNumeric($data);
        }
    }

    /**
     * Create and return a complex number matrix
     * $data elements are either:
     *  - a ComplexType
     *  - string representations of complex number
     *  - a 2 item array representing r & i e.g. [2,-4] = '2-4i'
     *
     * @param array $data
     *
     * @return \Chippyash\Math\Matrix\ComplexMatrix
     *
     * @throws \Chippyash\Math\Matrix\Exceptions\MathMatrixException
     */
    public static function createComplex(array $data)
    {
        foreach ($data as &$row) {
            foreach ($row as &$item) {
                if (!$item instanceof ComplexType) {
                    if (is_array($item) && count($item) == 2) {
                        $item = ComplexTypeFactory::create($item[0], $item[1]);
                    } elseif (is_string($item)) {
                        try {
                            $item = ComplexTypeFactory::fromString($item);
                        } catch (\InvalidArgumentException $e) {
                            throw new MathMatrixException('Invalid item type for Complex Matrix');
                        }
                    } else {
                        throw new MathMatrixException('Invalid item type for Complex Matrix');
                    }
                }
            }
        }
        return new ComplexMatrix($data);
    }

    /**
     * Create and return a rational number matrix
     * $data elements are either:
     *  - a RationalType
     *  - string representations of rational number
     *  - a PHP float
     *  - a 2 item array representing numerator & denominator e.g. [2,-4] = '-2/4'
     *
     * @param array $data
     *
     * @return \Chippyash\Math\Matrix\RationalMatrix
     *
     * @throws \Chippyash\Math\Matrix\Exceptions\MathMatrixException
     */
    public static function createRational(array $data)
    {
        foreach ($data as &$row) {
            foreach ($row as &$item) {
                if (!$item instanceof RationalType) {
                    if (is_array($item) && count($item) == 2) {
                        $item = RationalTypeFactory::create($item[0], $item[1]);
                    } elseif (is_string($item)) {
                        try {
                            $item = RationalTypeFactory::fromString($item);
                        } catch (\InvalidArgumentException $e) {
                            throw new MathMatrixException('Invalid item type for Rational Matrix');
                        }
                    } elseif(is_float($item)) {
                        $item = RationalTypeFactory::fromFloat($item);
                    } else {
                        throw new MathMatrixException('Invalid item type for Rational Matrix');
                    }
                }
            }
        }

        return new RationalMatrix($data);
    }

    /**
     * Create and return a numeric value matrix
     *
     * @param array $data
     * @return \Chippyash\Math\Matrix\NumericMatrix
     */
    public static function createNumeric(array $data)
    {
        return new NumericMatrix($data);
    }

    /**
     * Construct a complete matrix whose entries are a result of a function
     *
     * The function must accept two parameters
     *  e.g. $function($row, $col) {return $row - $col;}
     *
     * $row and $col are 1 based
     *
     * @param callable $fn
     * @param IntType $rows Number of required rows
     * @param IntType $cols Number or required columns
     * @param string $type ['numeric'|'rational'|'complex']
     *
     * @return \Chippyash\Math\Matrix\ComplexMatrix|\Chippyash\Math\Matrix\NumericMatrix|\Chippyash\Math\Matrix\RationalMatrix
     *
     * @throws \InvalidArgumentException
     */
    public static function createFromFunction(callable $fn, IntType $rows, IntType $cols, $type = 'numeric')
    {
        if ($rows() < 1) {
            throw new \InvalidArgumentException('$rows must be >= 1');
        }
        if ($cols() < 1) {
            throw new \InvalidArgumentException('$cols must be >= 1');
        }

        $source = array();
        $rc = $rows();
        $cc = $cols();
        for ($r = 0; $r < $rc; $r++) {
            for ($c = 0; $c < $cc; $c++) {
                $source[$r][$c] = $fn($r + 1, $c + 1);
            }
        }

        return self::create($type, $source);
    }
    
    /**
     * Create a matrix representation of a complex number
     * For z = a+bi
     * Returns [[a, -b]
     *          [b, a]]
     * 
     * @param \Chippyash\Type\Number\Complex\ComplexType $c
     * @return \Chippyash\Math\Matrix\RationalMatrix
     */
    public static function createFromComplex(ComplexType $c)
    {
        $a = clone $c->r();
        $b = clone $c->i();
        $b2 = clone $c->i();
        $bi = new RationalType($b2->numerator()->negate(), $b->denominator());
        return self::createRational([[$a, $bi],[$b, $a]]);
    }
}