|
|
|
last edited 7 years ago by rrogers |
| 1 2 3 4 5 6 7 8 9 10 11 | ||
|
Editor: rrogers
Time: 2018/08/15 22:10:35 GMT+0 |
||
| Note: | ||
removed:
---
---
-
-- Functions to facilitate Transforming Polynomial Generating Functions into
--
--
-- 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) |