Ore_algebra[rand_skew_poly] - random skew polynomial generator
|
Calling Sequence
|
|
rand_skew_poly(vars, eqns, Alg)
|
|
Parameters
|
|
vars
|
-
|
indeterminate or list or set of indeterminates
|
eqns
|
-
|
(optional) equations specifying properties where option=value
|
Alg
|
-
|
Ore algebra table
|
|
|
|
|
Description
|
|
•
|
A rand_skew_poly command generates a random polynomial in the variables vars from the algebra Alg. It is useful for generating test problems for debugging, testing and demonstration purposes. Several options can be specified, determining the form of the polynomial. This allows for general expressions with certain properties.
|
•
|
The syntax of rand_skew_poly (described below) mimics that of randpoly.
|
•
|
The first argument vars specifies the variables in which the polynomial is to be generated. If vars is a single variable, a univariate polynomial in that variable is generated. If vars is a list or set of variables, then a multivariate polynomial is generated.
|
•
|
The remaining arguments eqns are equations of the form option=value.
|
•
|
The possible options (and their default values) are:
|
Option
|
Use
|
Default Value
|
|
|
|
coeffs
|
generate the coefficients
|
rand(-99..99)
|
expons
|
generate the exponents
|
rand(6)
|
terms
|
number of terms generated
|
|
degree
|
total degree for a dense polynomial
|
|
dense
|
the polynomial is to be dense
|
|
|
|
|
Important: The terms option is intended for specifying sparse polynomials, where the number of terms is often significantly smaller than the degree of the polynomial. The terms option will be overridden by the degree option. A call rand_skew_poly(x, terms=7) will return a polynomial with 6 terms, since the default degree is five. If you want a dense polynomial, use the dense option.
|
|
|
Examples
|
|
>
|
|
>
|
|
>
|
|
| (1) |
>
|
|
| (2) |
>
|
rand_skew_poly([Dx],coeffs=proc() randpoly(x) end proc,A);
|
| (3) |
>
|
|
| (4) |
>
|
|
>
|
rand_skew_poly(Sn,expons=rand(-5..5),coeffs=proc() randpoly(n) end proc,A);
|
| (5) |
|
|
Download Help Document
Was this information helpful?