#### D.4.35.6 syModStd

Procedure from library `symodstd.lib` (see symodstd_lib).

Usage:
  syModStd(I,sigma); I ideal, sigma intvec

Assume:
  size(sigma) = nvars(basering) and sigma(I) = I. If size(#) > 0, then # contains either 1, 2 or 4 integers such that
  - #[1] is the number of available processors for the computation,
  - #[2] is an optional parameter for the exactness of the computation, if #[2] = 1, the procedure computes a standard basis for sure,
  - #[3] is the number of primes until the first lifting,
  - #[4] is the constant number of primes between two liftings until the computation stops.

Return:
  ideal, a standard basis of I if no warning appears;

Note:
  The procedure computes a standard basis of the ideal I (over the rational numbers) by using modular methods and the fact that I is invariant under the variable permutation sigma.
  By default the procedure computes a standard basis of I for sure, but if the optional parameter #[2] = 0, it computes a standard basis of I with high probability.
  The procedure distinguishes between different variants for the standard basis computation in positive characteristic depending on the ordering of the basering, the parameter #[2] and if the ideal I is homogeneous.
  - variant = 1, if I is homogeneous,
  - variant = 2, if I is not homogeneous, 1-block-ordering,
  - variant = 3, if I is not homogeneous, complicated ordering (lp or > 1 block),
  - variant = 4, if I is not homogeneous, ordering lp, dim(I) = 0.

Example:

    LIB "symodstd.lib";
    ring R1 = 0, (x,y,z), dp;
    ideal I;
    I[1] = -2xyz4+xz5+xz;
    I[2] = -2xyz4+yz5+yz;
    intvec sigma = 2,1,3;
    ideal sI = syModStd(I,sigma);
    sI;
    ==> sI[1]=x2yz-xy2z
    ==> sI[2]=xz5-yz5+xz-yz
    ==> sI[3]=xyz4-1/2yz5-1/2yz
    ==> sI[4]=y2z5-1/2yz6-xyz+y2z-1/2yz2
    ring R2 = 0, x(1..4), dp;
    ideal I = cyclic(4);
    I;
    ==> I[1]=x(1)+x(2)+x(3)+x(4)
    ==> I[2]=x(1)*x(2)+x(2)*x(3)+x(1)*x(4)+x(3)*x(4)
    ==> I[3]=x(1)*x(2)*x(3)+x(1)*x(2)*x(4)+x(1)*x(3)*x(4)+x(2)*x(3)*x(4)
    ==> I[4]=x(1)*x(2)*x(3)*x(4)-1
    intvec pi = 2,3,4,1;
    ideal sJ1 = syModStd(I,pi,1);
    ideal sJ2 = syModStd(I,pi,1,0);
    size(reduce(sJ1,sJ2));
    ==> 0
    size(reduce(sJ2,sJ1));
    ==> 0

User manual for Singular version 4.4.1, 2025, generated by texi2html.