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