|
|
last edited 6 years ago by rrogers |
1 2 3 4 5 6 7 8 9 10 11 | ||
Editor: rrogers
Time: 2018/08/15 22:11:00 GMT+0 |
||
Note: |
added:
--
--
--
-- -- -- Lower strict Jordan -- When used on the left of variables x it divides by x -- with 1/x =0 -- When used on the right of the coefficient matrix it shifts -- the columns to the left by one. -- Really a horizontal shift of the coefficient array and -- Should be applied A*SJ_lower*[x...]^T -- SJ_lower(dim) == J : Matrix(Polynomial(Fraction(Integer))) := new(dim,dim, 0) for i in 1..(dim-1) repeat J(i+1, i):=1 J
-- -- -- Upper strict Jordan -- When used on the left of variables x it multiplies by x -- With last entry 0 -- When used on the right of coefficient matrix it shifts -- the columns to the right by one. -- SJ_upper(dim) == J : Matrix(FPFI) := new(dim,dim, 0) for i in 1..(dim-1) repeat J(i, i+1):=1 J
-- -- SJ_upper(4)
Compiling function SJ_upper with type PositiveInteger -> Matrix( Fraction(Polynomial(Fraction(Integer))))
(1) |