Changeset 6d1466 in git


Ignore:
Timestamp:
Feb 24, 1999, 2:40:06 PM (25 years ago)
Author:
Wilfred Pohl <pohl@…>
Branches:
(u'spielwiese', '8e0ad00ce244dfd0756200662572aef8402f13d5')
Children:
a3836d9cf2acddc950056d641aefb9462400d7cd
Parents:
34cb4bea84226d98334d714b43a4556d8f1c5487
Message:
new bareiss


git-svn-id: file:///usr/local/Singular/svn/trunk@2868 2c84dea3-7e68-4137-9b89-c4e89433aadc
Location:
Tst
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • Tst/Old/barei.res.gz.uu

    r34cb4b r6d1466  
    11begin 644 barei.res.gz
    2 M'XL("$&WJ34``V)A<F5I+G)E<P"U6DMOX\@1ONM7,(,<))/TL%]\C&,?DEP6
    3 M"'+9W`Q[H;$TLT8L6Q`U"Y)!_GNJGU7-AR$'&,&6FE75U?7XJKI%ZM=__?V7
    4 M?R9)PNZ2SY]7O[HK?@U7R7G?GI/VZ?1\/"??WD[)U^UI_]RVR=/;X;!]W05I
    5 M<4WGRNODN/V^?]F_?C__GMPFK(#736"KN^3T_/H].0%+\*(0V;K+^FS89+LC
    6 M2I4@=;,"&Y+DZ??M:?MTWI^>V_/S4_+%SK*\UQ^'K_M3\O8M^6-[:C7/TLWK
    7 MZ\O;T[^U:T!_.^WV9MG=D4C@ZTORNCWL6SWLDCX9DAD]G.KY6S"UNDL.V_/I
    8 MN4L.[%X\P-]ME_%,Z/\NDQG3PPY=JT&>W:Q`F&4,9.V(/]QR.X+Y0H^XYG([
    9 M`J[H[!#84H^$9C,[XFZ*T%S>A:6:N^0%HI8<#NS6)6]]8)M@"RON-.]F=<\>
    10 MOJQ@PF_6IKSL1,J;+N>EI\(2A1^+,#8VYJ)+2W_-S6R>,NXI*.U,=F-N3'47
    11 MQNE[;JTP$0LV`C*/;R]]<MSMS[<'':('8R5ZP>^2YV_KE_UV]_2V_[;6@IN_
    12 M%)O`!WS^)UP`/!.K:YVSS94>H:;R.OEON*BNC5P._S9J16!!#C_E%[P^W:PN
    13 M$0MZ&UJ%O+@.R.+WZN&^A$@"JEB60R1YEC<PY)F":)59KH<`M@+&FL>U$,@&
    14 MWS@+8.`(!HY@X!#&(V`;G(7EV(-V.$N2Z?\2F2M8CQ>9D'4FN,P<N82/G(&%
    15 M`FQB-7R6F5/",YUL5)*+8(SX.1'F,HJPP@@+7;L2ZK$'ZSH^0&P'WF<]'S@#
    16 M"EQU$'48#:(7&+0R1%5@5`6)JFX/!S$ML8[W<BC33@V<U)@F`T6E8,&TW%@F
    17 MX](#Z4&0TH-KL)>67M##R5Q3B%U/*K$@A1@J5-KI6)8:89BB&O$B+%ZH4QEU
    18 MQ6$A61F#,]=8G;69NS!2*X@_MN6"C`U2M#U^_>9"B%R*D0`2`6V1H$2P:]N!
    19 M6@9[%N]$SU,N.M6GK)!=7Z:J'$B#%]PUK);K+>ZSF9#SSV(H>PE]$5R0*"R\
    20 ML`#AG'>RAX#)0:&$]!(2)*K/=5?U8LBQL(7R`@H$)`15@!X(+4J47J+4^_%G
    21 MX$%_!\F!R%1>I@*9+NW3`7FUY]6:QU,(`)W:>':CV4/:]6F/LV6!&Z1T&V2K
    22 M"ZG5]=0J/:JRML[:!N=`LSI`E&""J187:!MUFP%5=E`Y1L#@UXQ`-8FAIG&W
    23 M34I;("0\EB;<7%L2VFU[98K).FH)NI;Z%+SK`P(E#\4OL?@E%K\4NOCEI/AU
    24 M#NNT[LJ^,>_UD$N35@6?@!+W"8@Q],I<U:ED)?C&"G-9#>ZZL0,!C`%B4T%-
    25 M*?C4T^TGAXIL>%?W*N=%K6?82SGHSPJ*U-%K.Q"@,I>EEC#6E?9"&&*EIVFJ
    26 MTA^J+P?8[DNM&A14SGSX4$:M3L7`TZX<=$:,99"8*N5<PJ4$EAV)09B1`)K(
    27 M&>.&QMV(.YK0<A<<2$*>R)E$ZQ:Y3K\T[V)P.-)8$.:=`W)]RN?/+@7IF!1)
    28 MM'?Z;FEVWP`$>6FONKQ;8;^2BFYJLL1-3>EZ`V-U=4&--5E;FHI#A%8!P@HA
    29 MK`B$29=7>"I8>A'>BFL0Y/`&7:]T":U-.GE&Y=[31Z/L+?II?5\5_KL)@X/6
    30 MNENSZ^MFL\G6NR,&1#%L9Z5K9X=M][S;;U_6Y&RML#F4&-F2"`B,;&DC:]O#
    31 M*!RSX=%$V&[7Y>:J6U>;M%M+/6HVF2?6FJ@<T8AKP8DN+3A#;,+I6?G]Y\AN
    32 MK:'QZ5OY[>?(;\U)F;I8CD[F;/.G6[SB&URE(F=T59LS.DN/N,F4!3F:EY"#
    33 M_4N[#]><S"Z%G9W3V4K/#N?W$DKDPG*\N!J#[LJ#B`<0,?`TV[5H3HT8JMRY
    34 M,YD#4=D$C%1`!66;K%MSR%BW%N9=;K*5SG.FTZX)E7FO#;G19%:8=S,3#/&J
    35 MJR+@LT)\5KAVQ1"?U3OX=-UX`J*%LC9D#=W:8M.C%$S+/!5,3AV,(S+3Y,J3
    36 MK2;CY605Z_8<F<V326#XQ;WZP^"H!":^AJPO-8]*AMC7/N\)2?S*)#XQB4],
    37 MYE<D\Z9^-64I]2JDOL;4UV3Y$E-?O]OTWTVRW8X_#(P&(5`1"#0(@7H,@?\_
    38 MV=5/3#:I\D8G6Y)DDR_>%19Y8Y/--X_V6Y%QP?EJQI9N'';D5;A@KMK#W`:G
    39 M^N!EMF"N7,CT^-%_`<.P8KT9^9K(NZK-K'Q8RB<*LN&\JK&_-`BR!MVN27]I
    40 M)OT%0>)&2*"XLL>[;$%,H)AP8O:<1LX?U#]TSSF'C@6KN3[2-[@-6DFKX%'8
    41 MJ#V*/)!]L!U'IA/.H[1+T#E!([\B&)B3<'IYD*@V5+>(^79U$4T5@>6FCE2+
    42 M\>+$YD=._?&K>5Y.5XM=B2<N2LQ:Y/C".UMOYD*!_";X8,(?S0G2,F3X)]UZ
    43 MJB6Y4<ST_3P5-?]L;B.H52@2QN9V@H3N!*O)(<!N!6X[7)&]()MJFFXJD:85
    44 M544U%;:.O,5X1XS1N\[4*7-/C$WO.VL[KMQ1UIUL3?IR;>:5.RB,&,*FL/)0
    45 M2K5C7C90<ZI:!0`8<:?:2'JJ'ANZ"*O584D6-(T8W"M!6YQYRN,X#^:%.L[I
    46 M4J47]'YPLH21I2'R=9'&`2)DMSX&#<@Y#5'$28-ZC(:=P&>B1SAH.SG`,6)I
    47 M/(?-A);,<3Y[-S#+DO:V:&U%.29]++:8G"HXB4DS6IHF(3:7$]P1#H5-G`Q,
    48 ME.DS"!QL[I'N,HBF07%%:2,$EZ[AY3&$"9GFLG+D=)++:J0'0T#FT(SY/CT*
    49 M98UZ1@&K1VNC(\UHSCB7C)25#V5L$^*>LM1H[0*UR1A<@8XQ;<;:!*G!$0OA
    50 M33:PG(([IM.V9#:K=-)*/#FT"$_+XVY`R'(:!*."4RN(/.U*_FL6-DX"(Q[G
    51 MS+<63EH+6F\/":.Z&W/&Z.(Q[B2&@<>X4Y0S+9QH#B*>3!AAHL$)%!*$3!%!
    52 M#E@4#20<:`\%*B==0HWIDBQ`62G%SX@5Y\U7*'&2=&XVM]2$->[##*V(&O%8
    53 M6U2F<1CJ698<%5WL;`S?&!KE/(]-X1UG&_&*A1;A"/MJ&9/SX`EJR$<]G:A@
    54 M8RLU_?U;U?1A'LO4Z)GZS,ECE(]FDX[/#'D<?NP8LUM=-2).^WX^+3"ZKY`<
    55 MIS/0QZ_J4:^:;4=7D]YRM=`+HCQ%[896UABSZ3R&Z!."A0/)Z"R2OG=$J$?$
    56 MN6W]:FY3QPTAG]LW\[F-,9W=^_*Y#IO.M=#(FG%`XU8USF.07:B].BZ[^;K@
    57 MT6,;^F!;C1Y41\UPM`G&-8FFCE)L'^C%NW*DR1\S@^]$%=IL'P-2-6RD9I2!
    58 M<4\D!M''\>BQ=-5?A9J7T3SIL$I.[&HD((@`HS=2O8!^KC\R1Q(3E`UZL_&7
    59 M)G2PEK\V,0#5_EJ2J<KP0"W]U0`XF.'WZTM_///AYS>U^[W3']O3>OORLC$_
    60 M(#OQF7N-]\5#DESIV_3A=U[Z"^F\G%:9)>WSL/^2<)2?$;?R[DN^2N"?J&\6
    61 MQ!?4SXA'ZB6HET1]_3'U,^(3]8*HKSZF?D8\4B^\]:>9*%K)*#O'Q2SJ9T\H
    62 MMJ@M$CL<RH^Y,R,^<<=$ZS0C../.X:`^9L",^,0`(GZ0'U,_(S[OGWFUB]B,
    63 MHMPN8BP66\1*++:8@UAL,5*QV*+'L9BX3.PR=+:7HG-FU?>RMVCDN-B<^AEK
    64 MWU._Z!QI="66X&Y_OLC)I?:YU&XO<-)`]!^__'5!LCV;![/M>?NZVYYVUR_/
    65 M7_W6T11W9N*G<WO6]$\W,/@-),\_6GT_]L\KEMS=)4B\?OMQ7OT/E/LN&CTM
    66 "````
     2M'XL("(L!U#8``V)A<F5I+G)E<P"U6DMO&[L5WNM7L$$7EF;&'K[FD<!>M-U<
     3MH.CF=F<X@6(IN<:5;$&C!#,J^M][2`[)PX=<>Q$C-H8\'P_/^QR-\ON___';
     4MOP@A]([<W"Q^GU?L&E;DM!U.9'@\/AU.Y-O+D7Q='[=/PT`>7_;[]?/&H?DU
     5M/BNNR6']?;O;/G\__4%N":WAYY,CRSMR?'K^3HY`JLNKL9S*\[+<'#RB`<2G
     6M!=Q/R.,?Z^/Z\;0]/@VGIT?RD=1F__G'_NOV2%Z^D9_KXP#[W.SKGZ^[E\<_
     7ME4JP_W+<;/5UFP-"^)^/Y'F]WP[J<203.9,,'X;Y_-V)V=Z1_?IT?!K)GM[S
     8M!_AW.Y:LY.IW+$5)U>/HU>H`3S\M`$Q+"ECSQ!YNF7F"\UP],45EY@FH?#2/
     9M0!;JB2LR-4]L/L(5E8WNJOZ.[,!B9+^GM[/3KO9TZ62A]9VB?5K<TX>/"SCP
     10MY5[=.:Z^;Y^O^+(0^H$M"ZX?Z-)@X+I*.%`_6E0_;U7=O&/@(%/5C-SB66^/
     11M5JR9GQ;`4PM`E=6<>!",/[>/)PBYGYOMZ7:O[*,,[!5@=^3PLIO(0='A-U2/
     12MWY&G;U>[[7KS^++]=J5`R[_<NK5BNEPZ-`3L?]Q"7H,T"E"H4[.8GG-[3?Z[
     13MJ-VRNR;;W;!UZQZQ8K5B]6%[/+X</S@.C-D+JMP%3*@+[`)RY4/UAA]@_Q:8
     14MX]O@9&>@DPUD=B\?[IN'VQK<0<L*8HV555\J[T@(YZ:LU"/$=@W/BL84"+!>
     15M@<[%'O.QQY"&$)P'2"5P&5Q''X`"+.J2$/6[@#\5;_1Z0>&9=G!GP^QI7O\:
     16MD_"@_G'F3<)5;@O(UPF$&MD9C'%F4SFQ,Z.P`ZL1S`1/9SYQIR7GS@S<FP$Y
     17MF@N5@CQ*0>`?9]S()GX6-HDL`.<8(.29R1AR*;NX]`[@Q@%*+65Q8IQ@?Y0S
     18MB*+E]Y5<I;U=7>MN:-[HI+=ZR;L)BB[V$^2?+@,#A6["1CZQ@O%13@6MQ3@U
     19MA6S.J`3S?JX:`P,XO]$'*G;#S\TD"LI`!>'`HK9@#N"*C6(")<59>@2U"`&(
     20M]J8;6S!)Y7-!V"HU2``(93'@PR;F$=PB&M4I;X`V\@*09X01%M,"9BRFXNQI
     21MTM(Z16,%&``?;2RY5^1S,4[%A$ZC%B;F%C:H4!Y41`]2/;7ET)5#[\^H-@96
     22M@@.ZC<V&-E8W'I#->&8&`/%;FR<5J-Z&:H_-C4R8/H?,8_;X?%8W/*VV6>F<
     23M,(J:#=5XIP*TFUP$"M\!A4\_X=-/Z@XHHO0SRMBVIG5R"Z6:76@-[<(UQ,D^
     24MV/PKDN9)H1]._7RRZAE$C.NUK.Z`:%MH(1H53A;:@2C2D6!ACU6,"3C&O3R4
     25M:EG=1J7IMO4J,PN\X$Y:KZ)>L(GCA:LJD.E.8_=TMD^^+*ETZ-QD,#93CQ>=
     26M8R=TVDB_A&0,EY"?&-QBFKM!V59,M,;$]AQ2^V#-`>TVE"N8NQ:6#1)*+YF3
     27M0CFNFYR(RG'":Z>)XHR6+53($-P%:^X%4V[OO(.4I9J`Q#&R1?<HJ$0K.34.
     28MVBB:=Z8(32B0*R065:6K<VPU-MX>0&%X=6Z=I2`<!;!DX08_<[S!`>&B2P6L
     29M<!'L-EB(X(B':VQJ2G&-329CH[!C(PSUJQ4AQ^UF^^WI6<WQ"@%[[FPZ4*IJ
     30MD9X[1.?>-6I*/&K*UT=-&8Z:,AHU)1XUF^RHV;P^:C9XU&S>.FJ^HVG[MMT$
     31M`V>#!DZIV@Z41M5DH-7TY=#HQN.E5/WF9?-CMX5Z?;OWW;?Q)5[Z$K_W^K6U
     32MGW.DF7-0CRGQ0'/Y9V%;C4+K0EKH"EJ8TEG-(X3G>RX7MAFI(TQE5*7J=:$K
     33M$$1YU>G@9O@6*S)$\)]/NYWZ_&NWV*^:I5IN/XG36_5!_(I>7_?+97FU.2`C
     34M"C\B-/.(L%^/3YOM>G>%/G*UTGFC\=YH$*#QWFCFJ1.B,>L&LP?T9EDJT&J\
     35MDE!EH)>M])ZC]X[>.3KL?=%#@;NY=67A0&_-[?BS9-O9U&<F\;'4?93>%"?W
     36M@?G4[FJ4CQU5F7>@Q8'%6=?Q(*T[$:9U)S&;QK"I,FPZE+Q=_^;D?7/N6MY]
     37M;:.$N2BAH'>Y&9PP/?5!TLX?ET@N2GKF@J#5`4"5`]FR)#H6X*]80OZ`MTOM
     38M:-AH]=].;VMWTUK_U2=!$,?:?^!J?0"VZ&[A`["U`0C7D7PA<#'6S7*ME"R5
     39M"[<204",<H:`7!:C=L-([-]>9-_OI\;[H`,'7$K4OG5FZ*P+"/+!;!3M`Z*=
     40ML$!.(-H)B\M>\)_^.^^%#EV//OUWJ`R\ZH6Y%)!<+<"8WF.">J`G4O!#Z7S%
     41MO:^8]Y4`C*W%M'[KFX;W.XO6*&%ZY2R!G.7?E=#:YTMOG,66GW4K(EK\E=4'
     42MGLV^]LR\O7`+.B>..]O[H[U%&`NNC&\U9H9KUSM\A_$=PG?.VNK97=7.1XFW
     43MK$_5W@=)C]36[T;Z^-V(N]I]$''"N\]B[CJ](]2.,TCR+M-S_.RF4\_S,_=<
     44ME3ZKF)50Z-:B67ROM<UJ=ELL=40)WNA88QJN3A9AY(-5*#6Z+3JA-?+2FY.A
     45ME)K67U2CC1GG%>TC"\U:":.ON3>Q187DYBFM72+)9K(UA#&6,&$1JA>:*U#0
     46MIL-G<<$"M2)%+G8`=PW#(KT"S%HL5`P)G3$LMAO->*]-39_`P@`0J7N\V#PP
     47M-0L8\)SN63;_#Y@5.G*O-5*WS%DR@?4^?D1T49@8"!"^'2V%JSZ_YJT[K5&'
     48MIE2]:)=!BRXS[9K6OEU3FNO7!/?K13(U$=L559%>H(Y=IIS2UA]P6F!6F%.-
     49M/\/0VH\`%'_[A)7J]==/\?=/NBT[YTE;KTV%4<_2EG-<+TTAUPU=H./,3`H!
     50M/QT1AE$USVN8LYRS<GYM9B>Z&G47B7,>O>SR`CB6P@T7,MZUHH37SX?G3)%F
     51MCK0IY;N9X]C'RC5H5Z)1%>]:RW"LFA:*.@6";>$XVVW-A#L5@FTK&&[3V)#1
     52MOG0">IL62+C4UB+P-JH/B>-SM,#P6/0*FS]1-6P@H<60QG7LXB:A>=MU&5H5
     53M.0<7S<`]F!`X*`R>5="]P_N%I])LP*:T>=L5]C!RD5EI5J$^9*C=$<P>T:$N
     54M.32';]0XPE"5::CYJ&]S1)]GF=DQ,@?+YEQ`BW*L"\]53NLN/!0XK,LS[,(9
     55M0;H@9&G"".=IF1+Q)Z@,,4GD2(,HG1/]D`8X4L.T=A23V#)*;&_8I';)RW3G
     56M2%PW99BL,3'(""]\'YL*^\T3V0I[NTV/![F;IS/+.9%ZUF?V!DYBBD0**5Z:
     57M)CH6A$N7\L36C^I@%=H^I<Z2HGRH8PO$"5/'L2]P9M2QVR/_I`"O<UK=*5(M
     58MRX%F?)AR"-(PTR2$3YN\>&F[F&.(ABKF.H,W<$JEH0%2$;SX[25`$-X7^U+4
     59M@W"`H\]'F7MS/8KBPIMI7]X@C:TV<2WP5@FH.*#:D%HE`=7FF7M[IAQPO$0=
     60M,7)6ES"/?/%:TT0!E7"(`RK3)3.#5Q3,V7X9SZC9CAH/3!$Y,P$E73>)MP@1
     61MSS8Y'\?O)R()O'=P<*43"Z9&<TE27(/)!%,S$PCFFYE!PO+KB[U1.FGT0;@G
     62MS7R.UW1(P-4I'5:BTG,!$"=#?MQ!0WDBO(_7=+Z)RD2.@T_E]'@4K7UR'`=K
     63M2L6Q&EDOFM=?&\!PLH6Z>;]FR`)=GD$$@V$>$4:-JU6)C5!W3<EQ<TP1<7>D
     64MB:!!>[QPQ:7)TANS>PTAXH*3LU54U'+1VKP*H6FF9D/.)UM277)CL)>PR5(+
     65M9X&$;1&UWY0OC=5"Y.!=E\1ON^B;O^9XY_?+5/TO7O4BZ.?Z>+7>[<S_ICCB
     66MK[WMSWW]0,A*?<OH_O>U>CV4QRF6)1F>SMN/A'E\!F[P\RLW2>`7L>\OP"^P
     67MS\`#]@+8"\2^>Q_[##QASQ'[]GWL,_"`/;?2'S-6-,C`.X>+7E1?I'O816[F
     68M^WBD4/,^A3+P1"%MKW_^]K<+R.&DO]8>3NOGS?JXN=X]?761R^[TP0^GX:3V
     69B/WR"AR^`//T8U,O9ORXHN;LC?O/ZY<=I\3\)>&E-2C$`````
    6770`
    6871end

  • Tst/Old/barei.tst

    r34cb4b r6d1466  
    33//
    44pagelength = 10000;
    5 ring r = 32003,(x,y,z),dp;
     5ring r = 0,(x,y,z),dp;
    66r;
    77matrix m1[3][3]=x,2,3,2,3x,4,1,3,2x;
    … …  
    99list mm1=bareiss(m1);
    1010mm1;
    11 poly pdet=mm1[1][1,1];
    12 if(leadcoef(pdet)<0)
     11vector vdet=mm1[1][3];
     12poly pdet=det(m1);
     13if(leadcoef(pdet)!=leadcoef(vdet))
    1314{
    14   pdet=(-1)*pdet;
     15  vdet+pdet*gen(3);
    1516}
    16 pdet-det(m1);
     17else
     18{
     19  "error";
     20  vdet-pdet*gen(3);
     21}
    1722"------------------------------------";
    1823//
    … …  
    4146list mm4=bareiss(m4);
    4247mm4;
     48vector vdet=mm4[1][3];
     49poly pdet=det(m4);
     50if(leadcoef(pdet)!=leadcoef(vdet))
     51{
     52  vdet+pdet*gen(3);
     53}
     54else
     55{
     56  "error";
     57  vdet-pdet*gen(3);
     58}
    4359"---------------------------------------";
    4460//
    4561matrix m5[3][2]=s8,s7,s9,s6,s5,0;
    46 list mm5=bareiss(m5);
     62module mm=m5;
     63list mm5=bareiss(mm);
    4764print(mm5[1]);
     65kill r;
    4866"-------------------------------------";
    4967ring r1=0,(x(1..9)),(dp);
    … …  
    5169list mm6=bareiss(m6);
    5270print(mm6[1]);
    53 poly p1=mm6[1][1,1];
     71vector p1=mm6[1][3];
    5472poly p2=det(m6);
    5573if(leadcoef(p1)!=leadcoef(p2))
    5674{
    57   p1+p2;
     75  p1+p2*gen(3);
    5876}
    5977else
    6078{
    61   p1-p2;
     79  p1-p2*gen(3);
    6280}
    6381"--------------------------------------";
    … …  
    7694print(m9);
    7795list mm9=bareiss(m9);
    78 print(mm9[1]);
    79 mm9[1][1,1];
     96mm9;
    8097"------------------------------------";
    8198matrix m11[5][5]=maxideal(1),maxideal(1);

  • Tst/Short/barei_s.res.gz.uu

    r34cb4b r6d1466  
    11begin 644 barei_s.res.gz
    2 M'XL("!4[IC4``V)A<F5I7W,N<F5S`+5:2V_CR!&^ZU=T!CE()FFS7WS,Q#XD
    3 MN2P0Y+*Y&?9"8VEFA5BV(6H6%(/\]U0_JYJD##G`"+;4K*JNKL=7U2U2O_[K
    4 M[[_\DS'&[]C-S>)7?R6NX8H=M]V1=4^'W=N1?7L]L*_KPW;7=>SI=;]?OVRB
    5 MM+RF<]4U^\<O?V6?CMWQ^GGW]=.7R-%W#(B_[5YVQ^7JR\)\LCM'ZX[KXX_N
    6 M^O7',4I7=RR.ZVOVMOZ^?=Z^?#_^SFX9+^&%BIL[=MB]?&<'8$E1EC)?]ODI
    7 M'U;YY@VE6I#ZL@"_&'OZ?7U8/QVWAUUWW#VQSVZ6X[W\V'_='MCK-_;'^M`9
    8 MGJ/;U]?GUZ=_FW`!_?6PV=IE-V]$`E^?V<MZO^W,L&<G-K`9/8+J^5LPE9=W
    9 M;+\^'G8]V_-[^0!_MWTN<FG^^USEW`S[Z!J'[.WYEP4(\YR#K!N)AUOA1C!?
    10 MFI$P7.%&P)6]&P);F9$T;.Y&PD^1ABOZN)2X8\\0-;;?\UL/B.6>K]`6>6=X
    11 M7Q;W_.'S`ES\S=E45+W,1-L7H@I46*(,8QG'UL9"]ED5KH6=+3(N`@6EO<E^
    12 M+*RI_L(Z?2^<%39BT49UQ]Y>GT_L;;,]WNY-B!ZLE>@%8'7W;?F\76^>7K??
    13 MED9P]9=R%?G5-?M/O`!X,J=K6?#5E1FAIO::_3=<B/+:RA7P[Z)61A;D\%-Q
    14 MP0L*ZA*QJ%?0RA90JP%9XEX_W%<024`5SPN(I,B+%H8BUQ"M*B_,$,!6PMCP
    15 MA!$"V>B;4!$,`L$@$`P"PO@&V`9G83G^8!S.&9O^GR,+#>N),I>JR:50N2=7
    16 M\%%PL%""3;R!SRKW2D1NDHU*"AF-J7Y2A.LDP@U&6)K:55"/)["N%P/$=A"G
    17 M_"0&P8$"5SU$'4:#/$D,6ANC*C&J$J,J37O8RVF)]>*DABKK]2!(C1DR4'0&
    18 M%DS+C><J+3V0'B0I/;@&>VGI13V"S+6%V)]()9:D$&.%*C<=R](@+*9(<L2+
    19 M='BA3N74%8\%MK`&Y[ZQ>FMS?V&E%A!_;,LE&5ND&'O"^N)"B%R*D0@2"6V1
    20 MH$3"+FD[4,=ASQ*]/(E,R%Z?,EZJ_E1ENAI(@Y?:-ZQ.F"WNQDXHQ(T<JI."
    21 MO@@N*!2N@K`$X4+TZ@0!4X-&B3I(*)"H;YJ^/LFAP,*631#0(*`@J!+T0&A1
    22 MH@T2E=F/;X`'_1TD!Y1199"I0:;/3MF`/!YXC>&)#`)`IXK`;@U[R/I3=B*S
    23 M)6Z0RF^0G2FDSM13I\VHSKLF[UJ<`\UJ#U&"";9:?*!=U%T&=-5#Y5@!BU\[
    24 M`M4DAH8F_#:I7(&0\#B:]'-=21BWW94M)N>H(YA:.F7@W2DB4.E8_`J+7V'Q
    25 MJ\H4OYH4O\EADS5]=6KM>S,4RJ95PR>@Q'\"8BR]ME=-IG@%OO'27M:#OV[=
    26 M0`)C@-C44%,:/LUT]RF@(EO1-R==B+(Q,]RE&LQG#47JZ8T;2%!9J,I(6.LJ
    27 M=R$ML3;3#%6;#WVJ!MCN*Z,:%-3>?/C05JU)Q2"ROAI,1JQED)@Z$T+!I0*6
    28 M&\E!VI$$FBPX%Y8F_$AXFC1R%QQ(8I[(F<3HEH5)O[+O<O`X,EB0]ET`<D/*
    29 MY\\N)>F8%$FT=X9N:7??"(3ZTEYU>;?"?J4:NJFI%C<U;>H-C#75!376YEUE
    30 M*RXB5)<1PAHAK!'"FG1YC:>"<R_"6P@#@@+>H.M5/J&-3:?(J=Q[^FB4@T4_
    31 MK>]K&;Z;<#AH+?LEO[YN5ZM\N7DC`5'8SBK?SO;K?K?9KI^7Y&RML3E4&-F*
    32 M"%08V<I%UK6'43AFPV.(L-TNJ]55OZQ76;]49M2N\D!L#%%[HA4W@A-=1G"&
    33 MV,;3LP[[SQN_=8:FIV\=MI\W<6M/RM3%=G0RYZL_W>*56,55JI*<T2MNS^@\
    34 M>\--II+D:%Y!#K;/W39>:SJ[<K,+.KLQL^/YO8(2N;`<+Z[&H+LN`XA$!!$'
    35 M3_--%\VI.6*H]N=.-@>B6D2,U$`%9:N\7PK(6+^4]EVM\H7)<V[2;@BU?6\L
    36 MN35D7MIW.Q,,B:IEQ&>-^*S)V@KQ6;^#3]^-)R`Z4]:6;*#;.&P&E()I>:""
    37 MR9F'<4+FAEP'LM-DO9RLXMR>(_-Y,@F,OKA7?QP<%2:^@:R?:QYU'6/?A+PS
    38 MDOB%33RSB6<V\PN2>5N_AG(N]4U,?8.I;\CR+::^>;?IOYMDMQU_&!@M0J`F
    39 M$&@1`LT8`O]WLIORYR6[(57>FF0KDFSRQ;O!(F]=LL7JT7TKLBYX7^W8T:W#
    40 MGKR(%]Q7>YS;XM00O-P5S)4/F1D_AB]@&%:L-RO?$'E?M;F3CTN%1+%X'[#!
    41 M_M(BR%KB-NDO[:2_($C\"`D45^YXEY\1DR@FO9@[IY'S!_4/W?/.H6/1:FV.
    42 M]"UN@T[2*7B4+FJ/LHCD$&S/4=F$\ZC<$G1.U"BN"`;F)+Q>$27J%=4M4[Y;
    43 M7293963YJ2/5<KPXL?E14'_":H%7T-525]*)9R5F+?)\&9QM5G.A0'X;?;#A
    44 M3^9$:14S_)-N/34UN5',S?T\G33_?&XC:)I8))S/[02,[@2+R2'`;05^.UR0
    45 MO2"?:IIN*HFF!55%-96NCH+%>$>,T[O.Q*G6WA/CT_O.QHXK?Y3U)UN;OL*8
    46 M>>4/"B.&="FL`Y0RXUB0C=2"JM81`%;<J[:2@6K&EB[C:DU<DD=-(X8(2M`6
    47 M;YX.."ZB>;&."[I4%02#'X(L865IB$)=9&F`"-FOCT$#<D%#E'"RJ!ZCX2:(
    48 MF>@1#MI.#G"<6)K.X3.A)7.\S\$-S+*BO2U96U..31]/+2:G"D%BTHZ6IDE(
    49 MS14$=X1#89,F`Q-E^PP"!YM[HKN*HEE47%/:",&5;WA%"F%"IKFL/3F;Y+(>
    50 MZ<$0D#DT8Z%/CT+9H)Y1P)K1VNA(.YHSSB4G915"F=J$N*<L/5J[1&TJ!5>D
    51 M8TS;L39):G#$0GB3#:R@X$[IM"W9S2J;M))`CBTBT(JT&Q"RF@;!JA#4"B)/
    52 MNU+XFH6-D\!(I#D+K460UH+6NT/"J.[&G#&Z1(H[A6$0*>XTY4P+)YF#B"<3
    53 M1IAH<0*%!"%31)`#%D4#"0?:0X$J2)?08[HB"U!61O$S8J5Y"Q5*G"2=F\\M
    54 M-6&-^S!'*Y)&/-:6E&D:AF:6I49%ESJ;PC>%1C7/XU-XI]E&O&*A)3C"OEJE
    55 MY")Z@AJ*44\G*OC82D-__U8U?9C'<SUZICYS\ACEHUUEXS-#D88?.\;L5E>/
    56 MB-.^7TP+C.XK),?9#/3QJWK2JV;;T=6DMUR=Z05)GI)V0RMKC-EL'D/T"<&9
    57 M`\GH+)*]=T1H1L2Y;?UJ;E/'#:&8VS>+N8TQF]W[BKD.F\VUT,2:<4#35C7.
    58 M8Y0]4WM-6G;S=2&2QS;TP;8>/:A.FN%H$TQK$DT=I=@]T$MWY413.&9&WXDJ
    59 MM-D]!J1J^$C-*`/CGD@,HH_CT6/EJ[^.-:^2><ICE9S8]4A`$@%.;Z0&`?-<
    60 M?V2.(B9H%_1V%2YMZ&"M<&UC`*K#M2)3M>6!6OJK`7`PC]^OVTM_///AYS>M
    61 M_[W3'^O#<OW\O+(_(#N(F7N-]^4#8U?F-GW\G9?Y0CHO9U3FK-L-V\],H/R,
    62 MN)/W7_(U@W^BOCTC?D;]C'BB7H%Z1=0W'U,_(SY1+XGZ^F/J9\03]3)8?YB)
    63 MHI-,LO-V-HOFV1.*G=66B.WWU<?<F1&?N&.C=9@1G'%GO]<?,V!&?&(`$=^K
    64 MCZF?$9_WS[ZZL]A,HMR=Q5@J=A8KJ=C9'*1B9R.5BIWU.!63EXE=AL[N4G3.
    65 MK/I>]LX:.2XVKW[&VO?4GW6.-+H*2W"S/5[DY+GV>:[=7N"DA:CY0?.\9'>T
    66 I#V:[X_IELSYLS"^><__+Y[B%2/H;9W,G]L\+/O/#Y_\!CT=ZV8LM````
     2M'XL("'(!U#8``V)A<F5I7W,N<F5S`,5:36\<.9*]UZ_(;<Q!<E;*9/"[#>DP
     3MNY<!%GN9O1ER0[:J>X31%U3J1B87^]_W160FR2J5>]2'QLIVB4$&@\&(%\$@
     4MRW__[__XVW]U7:>ONH\?-W]?*+H`U;WN]J_=_MO+W?-K]_/32_?UYF5WM]]W
     5MWYX>'FX>;PNWN6CGVHON/__VU^Z'U_WKQ?W=UQ\^E1%WU:'SI[O'N]>S\T\;
     6M_MU=S7W[UYO77_<73[^^%FY_U95VN.B>;W[9W>\>?WG]1W?9:86?*CA>=2]W
     7MC[]T+Q@RI)39GHW;:9O/M[?/E2N!Z],&^^JZ;_^X>;GY]KI[N=N_WGWK?IQG
     8MS6./OSY\W;UT3S]WO]V\['EL[I>?K_=/W_[)YD+_T\OM3I:]?6XXZL^/W>/-
     9MPV[/S;&;NMR=D$.MG']?5=7JJGNX>7VY&[L'_=E<X^_EN*6MX7_CUFXU-\>R
     10M-0WO/>A/&S#KK0;OW*+K2YI;F&^X13Q*<PNC9IR;&+;<,CRLYQ8M4PR/TEB6
     11MHJON'E;K'A[TY0*(LP=]7G4Q5SSV:?-97_^XP19_^LQKCA]^V3V>F?/>2H/.
     12M>R,-?3[S8+G!%J8TKEQIZ1KBTC.S0Z?!CV;EI[1.'<@OK0UDB@*:K5;4LU?=
     13M;[MOKX#S;[>[U\L'M@\;N&X`,'U^NI^Z9Q['O\/M`99W/Y_=[VYNOSWM?CYC
     14MIO-_NRPT"ST_+]P`[O\4(EY`&V;H>=:B9I%,ZJ+[WXTJI+[H=O?[7:&I$46&
     15M1?VP>WEY>JD!1FY=8#BY0.`%5@(Q\\^[^WL.B;4+`?+#\(X?K/@>ME6N46UN
     16M,=A7"595`O7L&V*UZFJH1@!]=M>?_?6E@A_U=@!(:3ND+;O5(0[\=N`F@D*A
     17MS6/$3."MTAB33[>_WN^ZD2X?J`[8@F8J:!ZIT<,QFND(S9\A7R%6"G0_*RS)
     18MBL3K@D_TH5=[NBY(Y,"E(AI(@C:O+S>/^^>G_>[LH5TW7'U?I_@]G>A8)PV=
     19M#+%>HI\_TL[&;;08)6S&-'I*Q,"B;NO+HND-7JSZ<_!B#\XB2]_#RWV#%VLJ
     20M7@QG3(LL.&&+(V48/=.TG2B39K=M\P@,H97-9*H$VT@PK5N:$+*NX,74[&=:
     21M#L_.,<?.836@S'3DH)%8*^+?V;"N1PX:'8^;81V?7"8W"-_J*_'ZUI;U!3BM
     22M<HUN`,XS;(F,!B/I:XS,>LF1M/R2GPT3/':Z?S(9L%G4Z<%75GB+$O=>E+P7
     23M)@4G[@`GKN#D\KM9Q9DEM>\UP$2CF:@G`R/W6MEQ\KWSN3E6G5W9B0N+CS)A
     24MH(\F^\GVFK!Q6YG78V-OP#S0:">8QF97.?S*8<$1/L8QP)!#35,NK`P.#);M
     25M##F`2.6(*X?G*N@CQD;3@S,W/&GE">`9^ZG/9<RK=2SR&/4P0#/5ZW4X\7#N
     26MQZF?FME-4K9+6;+G0-IS/.T=M\)V'[?[5.>8@CD+1RQVGHT^.\#Y,=.V15EC
     27MO.U&-P,'/XU]#F9O9,NG9JW[Y3'L#/N;LLQ8=:VG@:T!9"M^O)P&]CBX47-]
     28M;UNBRV%.MAJ%RY2&1/"_[4E%D+&WGN$P1$AP\*B;[$!D,60R]5J+4`@<I&O@
     29MS5OY-'E9D"UAY)/@N&6KQ\<18R[BTT,!_HP9:S$,>4V@>OD-A$M_$"J*SG;2
     30M2LB0%SKE93-:Y8%W0]GA-T^??U/VO,V('?$V[;SKB&7FW;N\],>\V"%E-@0X
     31M1#L_$R;/UL$T[G59K.1S[YGJ&2UA45\LQV(909F&T4,5M$E^Y2`VM6"DI66R
     32M6>QLLQ$[VV)QGF86VX.OYES.N*6B]/ZXHK1'%:4/;RK*%E3QCU24/C5E8%"_
     33M6U$&.J@H@SFL*(-M1;E3%64(OUM1AM14E/'=B?[]J;XF^WB0["-VMF8BQYD(
     34ML<5Y!]DG;?=><E'1,IH2UJZ&M:O;B+:>BVX^%[^77+[_LVFSBP1I+T':ST$Z
     35M+(='E9NWFS4)\11BX`X<][T$%T`X1($@M:NL*KLW!VWT?]9!&X%?6*=[O*R'
     36M68Q+W^/EXY=Z@L2T7L6U',-G^N+B\?'\X"*>FMNM__QXC;^7#S?CW>WNYOZL
     37MN6HE7?SFJ]]\PT`E\IXUXLYSW#W6N$OK<?],<]2U<^U1U.DVYIZI1EQR39@D
     38M?RI,DMSLGG7_3!(BCV4=K50-$:VPHR8"M6KN=%J96<AP0HAK`UDK*/'.2'MW
     39MH!79874@%0=J8@?NJSJQ.C`LQ79WPH-:I1)8@8M-C&U'W.&WW8@TPI_V'%%P
     40MYKC;2T>0SRC=B;NUDD^9"456T?Q$LH`C5'"$NC:_B:Q!'99BEY?K3H?SW+N1
     41MI46O#ZS+P'I^F#LK"]38+BS0:^7AWI_DN:2H0._.B'_83[JY\$0X@(_X4R[0
     42M-;?%U05=XX/%*.*#3IRP:9S0B1,VW_="O0_%ZH78+.^K%^+J!?,OO,`V]HN"
     43M'UBS@17^,'<V/*GRQ,+#G=WLAVWQE:F^HNHK"YZN:!K^1&<U`9/86;9Q%C76
     44MJO&29F?1^1<Y4#I1_\.Z'[3G?O',TKTIA%X"I\Q-=6I:.68+?IA]*SP+N[B^
     45M\,>6/S;\L5B;VV6IL$PM9Y6F&JJI@B35;5,3JFD%2=W7X?VT1&'=R$_RFKD]
     46M9H@-`[]G+M&Y+0RKJLQ@9P:&A/RR)8J);SVB%]SV:3-O=C;#%S.O\&5^B+2,
     47MKV5TM>S"8!>&_@T#AD2)$Q+*,O2A<?_O,"Z+T3%C.&\7-"?99LW,@2!SS+$(
     48M.EK/?$>Q9G=?Z(0!5DV.6(96D\.]GQ3SKQA/*KVPF2,CQ?-3EGS#ELJF6\\>
     49M2"AS5X:")_.G/)MI:AZPM.8W4W=P)&Q/'0_D2NAI?>I\Z-KS8?/FE.[6+,Q)
     50M8=.<$-NWDMX>-0>2-JVH5I)J*U]-OGX!T7X#T6XJ'-:%6I[XW75EB$>%X<'L
     51M],[*4!O5E&_\KOVF,M3&'%:&]=:AC6LJ0WX,;BM#$UK1\;`R;(38@R\-M(42
     52M\[5`;U^H<KV[$/F#5P/-[Z[+J]OA%V^'+V]:7E?G%_B'VNGF.\3MI;>UT\^=
     53M7R]C[5LN(+=?+V^'K[5[N8/<U1[X[N>GE[.[R]M/=U?JT]TP%&^YUEM.O/7P
     54M^>Z:'[K$JG=56V<:U_`KX"+SZUNAKA7J5Z&7_-%/L]S^:R,YMI(/M&7IC63?
     55MJNOU2<G#;2/:MTI[VXH^DMSJ[(]TSHOD5G"KLR\ZOS5&:%4.^J3@OE4YM"J'
     56M]?7M95_#NN$%6%[V',BW"&1YS?K_^S-Z?J!R?HH]X9(>)^T'AUN[FLCV08T&
     57MMWSB'K*35<RC_63YX2U.SO<3IGL[>IIR;RT$R5N:-M:.UD\ZY)[(0DJ<R.5A
     58M;JO)F-P+#]%D=1YXIH;\E'L6-CG,T\G1Z!6_&@PH.P.-COAY#E3P9K1VTA%$
     59MLB:,QD_D,_4VA)'B9/A1#!F$W[759"&AIP1"T^1$@%?\J,;377!A=/PP9WJ7
     60M$C:@)JVS&4+P:;18,8%(WKK1V(D"LX6(5?QD'$9P@XQIU'&R_/R6(@REU>0T
     61MO[;!5E@0\_NHDF<;VVQ[$W42.U.V,%345HRM,(18YG5HH@C*11O8ZL:#T3EG
     62MQ?)69CDR;'\1X2?/D[&V9?.[[(:@G!/[\_.E@?;B`)U=[YW2[`%*&(DI&':!
     63M"2"2<TJ<``%P#BS.CF`!2.%DV?H^^]X')[:W:">K9]L3"%C"S]97V0_>P\8P
     64MOHG\<)JB8>-;+Y/DG=2Q`.S*L\5##D.,K#,,[D"8&+P8W.0`!+G9X!J$@V9L
     65M;Y-`:).([6U9@`Z)7Z,<SV=$Y]@[CX5AZ\CM%(DMK7D`MSD3Q=)6*.NUF)IR
     66M'$A1$%,K#/F0&.F6):#\GT3L`O2<8!ICV-`)[9B,V!F@QP"*;2.&=J`,[F1B
     67M:`."D@IB:)U33SKQ0[!-,F=!.66M!A^C%TLKIJ)W8FI@GJGD$IL:H.>7;`6H
     68MLJTMCP7R3HPM4JR#4RS#'HPKP$W6FK\G4D8,KD%BU\&RQ8%W4$9$>@8\LP(T
     69M1HSNA+21=Z`8\TPB10O264Z!-W1!H)D$,[/1H0PAG6`>&UTQI0%2`3C0#I*4
     70M6!UH9R(A]41&N^;DHQCJ:*VXAAH(/1N2%H,;)J,E+P:'&@@^I8.8'$`':9UV
     71M;',@G6?B$$]L=9:S(AKK6NQ[,;-E*@;`A\U,3'F3G)A9,>4B8@)F!L!!H5"#
     72MNRTCG.<M4,9BKD_!1S&L8PJIQXAAH3,",%HM=M4\9@QB!F8%M#&&5(#,0HQM
     73MC$ENSMH7&#,Q>Q<&!:9Y"`F3S6EY"&`U)`8E'@NXFK)!%1,Z<EH4"6NFSCH4
     74M`&/YT*<4(IL3:`8%^"6QIF/*L_G8FF:>!K4]PUF'DINQ:*R@Q:J(FA0T6Y,A
     75M#$K!>5"#(2SQ..=GQG!D?'..L")ES<18+%68:J:B`1^LR:#%F+>"4@9MZF.B
     76MV9J.*2S`MA09)>GBF%`%E(`H*)N2V%`Q$945$P*@I&!1;^0[',^4-5[PB>::
     77M8%TF76`(4((B`KC8;EHH!@GL!DPR)_(:?\,40)14"ME4<`<4@B*;C%B*OU\'
     78ML!.G)V(88M!$SG",0B'FK`F)IB`-N",^IV3MR,!CTG*8LG4TDV3D-3]Q6P[_
     79MC,._Y,C(U`(N((WL8#V?%9ZAQHR1H06<H1V5$I2AN1[\_#U_28>)OV9?T01L
     80M@4J&Y/0G^1\!P"T.)<_8`E5.]TR^9#TN`GP!$%<!B(0`),`J7`8P^),<\5P(
     81M^')Z0V(HF0T``K5"!@`"%3R.)A@%`*+0&\EJ@`_:Y62&P%@S%V,F5IA`K\CH
     82M4@GV8-!$235*(!/+V0N)J28H!DJJT-!,(2H3QQ0#)=63U6>C2AH".(RJ<``Z
     83M0`*4LF_%Q'INHFC1)=4`#J"*_P$'H_L8+(<'P`!"BKYLJ!Z+$51-*(``#XK3
     84M`0!N4PQJDDFEQ,O&E&,/)9`I20-.Q]#J9C@=8TO=!@&VG&;PLK$U,<#-QA;'
     85MPLT8+)58-JX>4ER^N1+^7+^YQ9=<OKE25T%UOYXZ<*3Q)<+A2*;$=7"C\:52
     86M@K:A'";PG`EK&,-Q(-:R!^K%<D2PI^(:J>RG6*H8*)1*WF?7I"46V2]I+KNS
     87M5:4F05VHUEP./U@U!YPPE8HZ6UV*C,3$DJ%A>:MK78SBD4JU`&N#6O(N;`VB
     88ME+8F6U,/?9B72<FG,*XUM3;-UJXG-U>W=DZ27-K:4EZZ;%TY@&%!4&MYB/WX
     89M<I+":-:7`@^:AW(BPDXVS#>1;$N)%KDM)YMTKW>*;$N1E:2]7`0HNU(=8?,@
     90MUN(=)6\M:K!;ITO!;;.;RQ%LSU&MDK,S:UWKL[-K*8JRVO%%*>.BQ,6C-.2"
     91LDUU8;B$0%)>K`A9.@]C1JS[[^FX;?/M?B_G1ZR\;?>+_&_\?H4K;YP(M````
    6792`
    6893end

  • Tst/Short/barei_s.tst

    r34cb4b r6d1466  
    1212list mm1=bareiss(m1);
    1313mm1;
    14 poly pdet=mm1[1][1,1];
    15 if(leadcoef(pdet)<0)
     14vector vdet=mm1[1][3];
     15poly pdet=det(m1);
     16if(leadcoef(pdet)!=leadcoef(vdet))
    1617{
    17   pdet=(-1)*pdet;
     18  vdet+pdet*gen(3);
    1819}
    19 pdet-det(m1);
     20else
     21{
     22  "error";
     23  vdet-pdet*gen(3);
     24}
     25kill r;
    2026"------------------------------------";
    2127//
     28ring r = 0,(x,y,z),(c,dp);
    2229matrix m2[5][6]=0,2,1,-3,12,-9,1,2,5,3,6,-1,2,4,10,6,12,-2,-3,2,1;
    23 list mm2=bareiss(m2);
    24 print(mm2[1]);
     30module x2=m2;
     31list mm2=bareiss(x2);
     32mm2;
     33x2=transpose(m2);
     34mm2=bareiss(x2);
     35mm2;
     36kill r;
    2537"------------------------------------";
    2638//
     39ring r = 0,(x,y,z),(c,lp);
    2740matrix m3[3][4]=xy,0,x2z,1,z2y,y2z21,0,1,zx3,1,0,z3y3;
    28 list mm3=bareiss(m3);
     41matrix m33=transpose(m3);
     42list mm3=bareiss(m33);
    2943mm3;
     44mm3=bareiss(m3);
    3045print(mm3[1]);
     46kill r;
    3147"-------------------------------------";
    32  //
     48//
     49ring r=0,(x,y,z),(c,dp);
    3350poly s1 = 2x3y2+23x5y+104xy6+56z2x;
    3451poly s2 = 3/2x3y2-2/3z6y4+12x2z4;
    … …  
    4158poly s9 = xz+xy+yz;
    4259matrix m4[3][3]=s1,0,s3,1,s5,0,s7,s8,s9;
    43 m4;
     60print(m4);
    4461list mm4=bareiss(m4);
    4562mm4;
     63vector vdet=mm4[1][3];
     64poly pdet=det(m4);
     65if(leadcoef(pdet)!=leadcoef(vdet))
     66{
     67  vdet+pdet*gen(3);
     68}
     69else
     70{
     71  "error";
     72  vdet-pdet*gen(3);
     73}
    4674"---------------------------------------";
    4775//
    … …  
    4977list mm5=bareiss(m5);
    5078print(mm5[1]);
     79kill r;
    5180"-------------------------------------";
    52 ring r1=0,(x(1..9)),(dp);
    53 matrix m6[3][3]=maxideal(1);
     81int n=5;
     82int nn=n^2;
     83ring r1=0,(x(1..nn)),dp;
     84matrix m6[n][n]=maxideal(1);
    5485list mm6=bareiss(m6);
    55 print(mm6[1]);
    56 poly p1=mm6[1][1,1];
     86vector p1=mm6[1][n];
    5787poly p2=det(m6);
    5888if(leadcoef(p1)!=leadcoef(p2))
    5989{
    60   p1+p2;
     90  "error";
     91  p1+p2*gen(n);
    6192}
    6293else
    6394{
    64   p1-p2;
     95  p1-p2*gen(n);
    6596}
    6697"--------------------------------------";
    … …  
    80111list mm9=bareiss(m9);
    81112print(mm9[1]);
    82 mm9[1][1,1];
     113mm9[1][4];
    83114"------------------------------------";
    84115matrix m11[5][5]=maxideal(1),maxideal(1);
    85116print(m11);
    86117list mm11=bareiss(m11);
    87 mm11;
     118vector p1=mm11[1][5];
     119poly p2=det(m11);
     120if(leadcoef(p1)!=leadcoef(p2))
     121{
     122  "error";
     123  p1+p2*gen(5);
     124}
     125else
     126{
     127  p1-p2*gen(5);
     128}
     129kill r1,r2;
    88130"-------------------------------------";
    89 listvar(all);
     131ring r=32003,(x,y,z),(c,dp);
     132module m;
     133int d=64;
     134int b=8;
     135int db=d-b;
     136int i;
     137for(i=d;i>0;i--)
     138{
     139  m[i] = x*gen(i);
     140}
     141for(i=db;i>0;i--)
     142{
     143  m[i]=m[i]+y*gen(i+b);
     144}
     145for(i=d;i>db;i--)
     146{
     147  m[i]=m[i]+y*gen(i-db);
     148}
     149for(i=d;i>b;i--)
     150{
     151  m[i]=m[i]+z*gen(i-b);
     152}
     153for(i=b;i>0;i--)
     154{
     155  m[i]=m[i]+z*gen(i+db);
     156}
     157list rs=bareiss(m);
     158rs[1][d];
    90159tst_status(1);$
Note: See TracChangeset for help on using the changeset viewer.