spielwiese
|
Last change
on this file since 1ebec3 was
ef52d6,
checked in by Oleksandr Motsak <motsak@…>, 11 years ago
|
|
Update testsuite wrt current master (320c4f5e64cad778b4ada9556a9b0a6c71cab83a)
NOTE: some wrong/bad tests were removed...
TODO: add the reset of these updates after updating LIB/
|
-
Property mode set to
100644
|
|
File size:
2.4 KB
|
| Line | |
|---|
| 1 | LIB "tst.lib"; |
|---|
| 2 | tst_init(); |
|---|
| 3 | |
|---|
| 4 | // tst_status(); |
|---|
| 5 | |
|---|
| 6 | |
|---|
| 7 | // killall(); |
|---|
| 8 | |
|---|
| 9 | LIB "brnoeth.lib"; |
|---|
| 10 | |
|---|
| 11 | // special CURVE over F_64 (p=2) |
|---|
| 12 | |
|---|
| 13 | ring r=2,(x,y),lp; |
|---|
| 14 | poly f=y9+y8+xy6+x2y3+y2+x3; |
|---|
| 15 | ideal I=f,jacob(f); |
|---|
| 16 | closed_points(I); |
|---|
| 17 | I=f,x2+x,y2+y; |
|---|
| 18 | closed_points(I); |
|---|
| 19 | I=f,x4+x,y4+y; |
|---|
| 20 | closed_points(I); |
|---|
| 21 | I=f,x8+x,y8+y; |
|---|
| 22 | closed_points(I); |
|---|
| 23 | I=f,x16+x,y16+y; |
|---|
| 24 | closed_points(I); |
|---|
| 25 | I=f,x32+x,y32+y; |
|---|
| 26 | closed_points(I); |
|---|
| 27 | I=f,x64+x,y64+y; |
|---|
| 28 | closed_points(I); |
|---|
| 29 | list CURVE=Adj_div(f); |
|---|
| 30 | CURVE=NSplaces(1..6,CURVE); |
|---|
| 31 | CURVE=extcurve(6,CURVE); |
|---|
| 32 | CURVE; |
|---|
| 33 | def affR=CURVE[1][1]; |
|---|
| 34 | setring affR; |
|---|
| 35 | listvar(); |
|---|
| 36 | Aff_SLocus; |
|---|
| 37 | Inf_Points; |
|---|
| 38 | def projR=CURVE[1][2]; |
|---|
| 39 | setring projR; |
|---|
| 40 | CHI; |
|---|
| 41 | Weierstrass(4,10,CURVE); |
|---|
| 42 | tst_status(); |
|---|
| 43 | Weierstrass(3,10,CURVE); |
|---|
| 44 | tst_status(); |
|---|
| 45 | Weierstrass(2,10,CURVE); |
|---|
| 46 | tst_status(); |
|---|
| 47 | intvec DD=1,1,2,1,0,1; |
|---|
| 48 | BrillNoether(DD,CURVE); |
|---|
| 49 | DD=2,0,3; |
|---|
| 50 | BrillNoether(DD,CURVE); |
|---|
| 51 | def RP=CURVE[1][5]; |
|---|
| 52 | setring RP; |
|---|
| 53 | size(POINTS); |
|---|
| 54 | setring r; |
|---|
| 55 | killall(); |
|---|
| 56 | |
|---|
| 57 | // SUZUKI curve over F_8 (p=2) |
|---|
| 58 | |
|---|
| 59 | ring r=2,(x,y),lp; |
|---|
| 60 | poly f=x10+x3+y8+y; |
|---|
| 61 | ideal I=f,jacob(f); |
|---|
| 62 | closed_points(I); |
|---|
| 63 | I=f,x2+x,y2+y; |
|---|
| 64 | closed_points(I); |
|---|
| 65 | I=f,x4+x,y4+y; |
|---|
| 66 | closed_points(I); |
|---|
| 67 | I=f,x8+x,y8+y; |
|---|
| 68 | closed_points(I); |
|---|
| 69 | I=f,x16+x,y16+y; |
|---|
| 70 | closed_points(I); |
|---|
| 71 | I=f,x32+x,y32+y; |
|---|
| 72 | closed_points(I); |
|---|
| 73 | I=f,x64+x,y64+y; |
|---|
| 74 | closed_points(I); |
|---|
| 75 | list SUZUKI=Adj_div(f); |
|---|
| 76 | SUZUKI=NSplaces(1..3,SUZUKI); |
|---|
| 77 | SUZUKI=extcurve(3,SUZUKI); |
|---|
| 78 | SUZUKI; |
|---|
| 79 | def affR=SUZUKI[1][1]; |
|---|
| 80 | setring affR; |
|---|
| 81 | listvar(); |
|---|
| 82 | Aff_SLocus; |
|---|
| 83 | Inf_Points; |
|---|
| 84 | def projR=SUZUKI[1][2]; |
|---|
| 85 | setring projR; |
|---|
| 86 | CHI; |
|---|
| 87 | Weierstrass(1,20,SUZUKI); |
|---|
| 88 | tst_status(); |
|---|
| 89 | intvec DD=2,1,0,0,3; |
|---|
| 90 | DD[25]=3; |
|---|
| 91 | BrillNoether(DD,SUZUKI); |
|---|
| 92 | tst_status(); |
|---|
| 93 | def RP=SUZUKI[1][5]; |
|---|
| 94 | setring RP; |
|---|
| 95 | size(POINTS); |
|---|
| 96 | setring r; |
|---|
| 97 | killall(); |
|---|
| 98 | |
|---|
| 99 | // special HYPERELLIPTIC curve over F_16 (p=2) |
|---|
| 100 | |
|---|
| 101 | ring r=2,(x,y),lp; |
|---|
| 102 | intvec G,D,F; |
|---|
| 103 | poly f=x5+y2+y; |
|---|
| 104 | list EC=Adj_div(f); |
|---|
| 105 | EC=NSplaces(1..4,EC); |
|---|
| 106 | EC=extcurve(4,EC); |
|---|
| 107 | def R4=EC[1][4]; |
|---|
| 108 | setring R4; |
|---|
| 109 | G[5]=4; |
|---|
| 110 | D=1..33; |
|---|
| 111 | matrix C=AGcode_Omega(G,D,EC); |
|---|
| 112 | print(C); |
|---|
| 113 | // the code C=AGcode_Omega(G,D) is of type [33,22,>=10] over F_16; |
|---|
| 114 | // the basic decoding algorithm can correct up to 3 errors; |
|---|
| 115 | // for it, we need an auxiliary divisor F of degree 5 : |
|---|
| 116 | EC=NSplaces(1..5,EC); |
|---|
| 117 | F[14]=1; |
|---|
| 118 | list SV=prepSV(G,D,F,EC); |
|---|
| 119 | matrix y[1][33]; |
|---|
| 120 | y[1,1]=1; |
|---|
| 121 | y[1,2]=a; |
|---|
| 122 | y[1,3]=1; |
|---|
| 123 | print(decodeSV(y,SV)); |
|---|
| 124 | y[1,1]=0; |
|---|
| 125 | y[1,2]=a; |
|---|
| 126 | y[1,10]=a+1; |
|---|
| 127 | y[1,3]=0; |
|---|
| 128 | y[1,21]=1; |
|---|
| 129 | print(decodeSV(y,SV)); |
|---|
| 130 | y[1,4]=a; |
|---|
| 131 | y[1,2]=0; |
|---|
| 132 | y[1,12]=a; |
|---|
| 133 | y[1,10]=0; |
|---|
| 134 | y[1,21]=0; |
|---|
| 135 | y[1,29]=a+1; |
|---|
| 136 | print(decodeSV(y,SV)); |
|---|
| 137 | matrix H=dual_code(C); |
|---|
| 138 | print(H); |
|---|
| 139 | list SC=sys_code(C); |
|---|
| 140 | print(SC[1]); |
|---|
| 141 | print(SC[2]); |
|---|
| 142 | SC[3]; |
|---|
| 143 | list SH=sys_code(H); |
|---|
| 144 | print(SH[1]); |
|---|
| 145 | print(SH[2]); |
|---|
| 146 | SH[3]; |
|---|
| 147 | setring r; |
|---|
| 148 | killall(); |
|---|
| 149 | |
|---|
| 150 | |
|---|
| 151 | tst_status(1); |
|---|
| 152 | $ |
|---|
Note: See
TracBrowser
for help on using the repository browser.
