Home Online Manual
Top
Back: weierstr_lib
Forward: weierstrPrep
FastBack:
FastForward:
Up: weierstr_lib
Top: Singular Manual
Contents: Table of Contents
Index: Index
About: About this document

D.12.5.1 weierstrDiv

Procedure from library weierstr.lib (see weierstr_lib).

Usage:
weierstrDiv(g,f,d); g,f=poly, d=integer

Assume:
f must be general of finite order, say b, in the last ring variable, say T; if not use the procedure lastvarGeneral first

Purpose:
perform the Weierstrass division of g by f up to order d

Return:
- a list, say l, of two polynomials and an integer, such that
g = l[1]*f + l[2], deg_T(l[2]) < b, up to (including) total degree d
- l[3] is the number of iterations used
- if f is not T-general, return (0,g)

Note:
the procedure works for any monomial ordering

Theory:
the proof of Grauert-Remmert (Analytische Stellenalgebren) is used for the algorithm

Example:
 
LIB "weierstr.lib";
ring R = 0,(x,y),ds;
poly f = y - xy2 + x2;
poly g = y;
list l = weierstrDiv(g,f,10); l;"";
==> [1]:
==>    1+xy-x3+x2y2-2x4y+2x6+x3y3-3x5y2+5x7y+x4y4-5x9-4x6y3+9x8y2+x5y5
==> [2]:
==>    -x2+x5-2x8
==> [3]:
==>    5
==> 
l[1]*f + l[2];               //g = l[1]*f+l[2] up to degree 10
==> y-5x11+14x10y2+5x7y5-9x9y4-x6y7