Difference between revisions of "ApCoCoA-1:Weyl.WMul"

From ApCoCoAWiki
Line 15: Line 15:
  
 
<example>
 
<example>
A1::=QQ[x,d]; --Define appropriate ring
+
A1::=QQ[x,d]; -- Define appropriate ring
 
Use A1;
 
Use A1;
 
F:=x; G:=d;
 
F:=x; G:=d;
Line 27: Line 27:
 
x^4d - 2x^3d + 4x^3 + xd^2 - 6x^2 - 2d^2 + d
 
x^4d - 2x^3d + 4x^3 + xd^2 - 6x^2 - 2d^2 + d
 
-------------------------------
 
-------------------------------
If you want to multiply Weyl polynomials that are not in normal form say for example F=d^2x^3-2dx^2+7 and G=2d^3x-5xd+3, then  
+
-- If you want to multiply Weyl polynomials that are not in normal form say for example F=d^2x^3-2dx^2+7 and G=2d^3x-5xd+3, then  
first convert them into normal form before multiplication.
+
-- first convert them into normal form before multiplication.
 
-------------------------------
 
-------------------------------
 
F:=Weyl.WNormalForm([[d^2,x^3],[-2d,x^2],[7]]);
 
F:=Weyl.WNormalForm([[d^2,x^3],[-2d,x^2],[7]]);

Revision as of 14:37, 28 April 2009

Weyl.WMul

Computes the product F*G of Weyl polynomial F and G in normal form.

Syntax

Weyl.WMul(F:POLY,G:POLY):POLY

Description

This method multiplies two Weyl polynomials F and G and returns F*G as a Weyl polynomial in normal form.

  • @param F A Weyl polynomial.

  • @param G A Weyl polynomial.

  • @return The product F*G in normal form.

Example

A1::=QQ[x,d];	-- Define appropriate ring
Use A1;
F:=x; G:=d;
Weyl.WMul(F,G);
xd
-------------------------------
Weyl.WeylMul(G,F);
xd + 1
-------------------------------
Weyl.WMul(Weyl.WMul(G,F)-2G,F^3+G);
x^4d - 2x^3d + 4x^3 + xd^2 - 6x^2 - 2d^2 + d
-------------------------------
-- If you want to multiply Weyl polynomials that are not in normal form say for example F=d^2x^3-2dx^2+7 and G=2d^3x-5xd+3, then 
-- first convert them into normal form before multiplication.
-------------------------------
F:=Weyl.WNormalForm([[d^2,x^3],[-2d,x^2],[7]]);
F;
x^3d^2 + 4x^2d + 2x + 7
-------------------------------
G:=Weyl.WNormalForm([[2d^3,x],[-5x,d],[3]]);
G;
2xd^3 - 5xd + 6d^2 + 3
-------------------------------
Weyl.WMul(F,G);
2x^4d^5 - 5x^4d^3 + 18x^3d^4 - 27x^3d^2 + 36x^2d^3 + 14xd^3 - 18x^2d + 12xd^2 - 35xd + 42d^2 + 6x + 21
-------------------------------
Weyl.WMul(G,F);
2x^4d^5 - 5x^4d^3 + 32x^3d^4 - 32x^3d^2 + 148x^2d^3 + 14xd^3 - 38x^2d + 216xd^2 - 35xd + 42d^2 - 4x + 72d + 21
-------------------------------
Weyl.WMul(Weyl.WNormalForm([[d^2,x^3],[-2d,x^2],[7]]),Weyl.WNormalForm([[2d^3,x],[-5x,d],[3]]));
2x^4d^5 - 5x^4d^3 + 18x^3d^4 - 27x^3d^2 + 36x^2d^3 + 14xd^3 - 18x^2d + 12xd^2 - 35xd + 42d^2 + 6x + 21
-------------------------------

Example

A2::=ZZ/(101)[x[1..2],y[1..2]];
Use A2;
Weyl.WMul(x[1]^11,y[1]^11);
x[1]^11y[1]^11
-------------------------------
Weyl.WMul(y[1]^11,x[1]^11);
x[1]^11y[1]^11 + 20x[1]^10y[1]^10 - 10x[1]^9y[1]^9 + 33x[1]^8y[1]^8 + 23x[1]^7y[1]^7 - 17x[1]^6y[1]^6 - x[1]^5y[1]^5 - 
18x[1]^4y[1]^4 - 36x[1]^3y[1]^3 - 36x[1]^2y[1]^2 + 26x[1]y[1] - 16

-------------------------------
Weyl.WMul(y[1]^21,x[1]^21);
x[1]^21y[1]^21 + 37x[1]^20y[1]^20 + 27x[1]^19y[1]^19 + 17x[1]^18y[1]^18 - 37x[1]^17y[1]^17 + 43x[1]^16y[1]^16 - 17x[1]^15y[1]^15 - 
27x[1]^14y[1]^14 - 5x[1]^13y[1]^13 - 49x[1]^12y[1]^12 - 39x[1]^11y[1]^11 - 25x[1]^10y[1]^10 - 40x[1]^9y[1]^9 + 46x[1]^8y[1]^8 - 
35x[1]^7y[1]^7 - 47x[1]^6y[1]^6 - 30x[1]^5y[1]^5 + 45x[1]^4y[1]^4 + 40x[1]^3y[1]^3 + 3x[1]^2y[1]^2 + 41x[1]y[1] + 26

-------------------------------
--this computation is too slow with Weyl.WeylMul(F,G)

See also

Weyl.WNormalForm