Changeset ec66b3f in git
- Timestamp:
- Sep 18, 2018, 3:39:45 PM (6 years ago)
- Branches:
- (u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', 'b4f17ed1d25f93d46dbe29e4b499baecc2fd51bb')
- Children:
- 1e440acd9490b8589048022b1691c2b4e7e90c3d
- Parents:
- 1adfacf6826d2b0c6660854e9fa7689b6adf7538
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/ffmodstd.lib
r1adfacf rec66b3f 1 1 /////////////////////////////////////////////////////////////////////////////// 2 //version="version ffmodstd.lib 4.1.1.0 Dec_2017 "; // $Id$ 3 version=" $ Id$ "; 2 version="version ffmodstd.lib 4.1.1.0 Sep_2018 "; // $Id$ 4 3 category="Commutative Algebra"; 5 4 info=" … … 7 6 over rational function fields 8 7 AUTHORS: D.K. Boku boku@mathematik.uni-kl.de 9 @* W. Decker decker@mathematik.uni-kl.de 10 @*C. Fieker fieker@mathematik.uni-kl.de8 W. Decker decker@mathematik.uni-kl.dei 9 C. Fieker fieker@mathematik.uni-kl.de 11 10 12 11 OVERVIEW: … … 53 52 fareypoly(g,f); univariate rational function reconstruction 54 53 polyInterpolation(l,m); univariate polynomial interpolation 55 modrationalInterpolation(l,m); modular univariate rational interpolation 54 modrationalInterpolation(l,m); modular univariate rational interpolation 56 55 BerlekampMassey(L,i); Berlekamp/Massey algorithm 57 56 modberlekampMassey(L,i) modular Berlekamp/Massey algorithm … … 355 354 the length of the sequence L is long enough. If this polynomial does not split into 356 355 linear factors, an optional parameter M = BerlekampMassey(L',i) can be provided to 357 add more elements to the sequence. 356 add more elements to the sequence. 358 357 REFERENCES: 359 358 @* [1] E. Kaltofen and W.-s. Lee: Early termination in sparse interpolation … … 405 404 if(D1 == 0) 406 405 { 407 if(j>1)//(2*l0) < j && 406 if(j>1)//(2*l0) < j && 408 407 { 409 408 // early termination is triggered here … … 413 412 B1 = var(i)*B0; 414 413 l1 = l0; 415 416 414 } 417 415 else … … 1135 1133 vL[j] = T[j][1]; 1136 1134 if(evbd < T[j][3]) 1137 { 1135 { 1138 1136 evbd = T[j][3]; 1139 1137 } … … 1384 1382 /////////////////////////////////////////////////////////////////////////////// 1385 1383 // +++++++++++++++++ std for one parameter begins here +++++++++++++++++++ 1386 1384 1387 1385 static proc test_fmodI(poly f,ideal I) 1388 1386 { … … 1491 1489 } 1492 1490 } 1493 number Nd = 0; 1491 number Nd = 0; 1494 1492 number Nn = 1; 1495 1493 list Pl, Il, L, M; … … 1575 1573 1576 1574 /////////////////////////////////////////////////////////////////////////////// 1577 1575 1578 1576 static proc firststd(ideal I, ideal cI, int in_value, list #) 1579 1577 { … … 2055 2053 } 2056 2054 v = (v0-q*v1)*N; 2057 v0 = v1; 2058 v1 = v; 2055 v0 = v1; 2056 v1 = v; 2059 2057 f = g1; 2060 2058 g1= r*N; … … 2120 2118 { 2121 2119 // check that whether the relation -L[j+t] = c_0*L[j] + ... + c_{t-1}*L[j+t-1] 2122 // for all j >=1 where t= size(L)/2 2120 // for all j >=1 where t= size(L)/2 2123 2121 int i, j; 2124 2122 poly f = fareyresult[1]; … … 2136 2134 { 2137 2135 return(0); 2138 } 2136 } 2139 2137 } 2140 2138 return(1); … … 2271 2269 * cannot be satisfied for all i 2272 2270 */ 2273 return(1); 2271 return(1); 2274 2272 } 2275 2273
Note: See TracChangeset
for help on using the changeset viewer.