Changeset ec66b3f in git


Ignore:
Timestamp:
Sep 18, 2018, 3:39:45 PM (6 years ago)
Author:
Andreas Steenpass <steenpass@…>
Branches:
(u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', 'b4f17ed1d25f93d46dbe29e4b499baecc2fd51bb')
Children:
1e440acd9490b8589048022b1691c2b4e7e90c3d
Parents:
1adfacf6826d2b0c6660854e9fa7689b6adf7538
Message:
chg: fix issues in ffmodstd found by formal checker
File:
1 edited

Legend:

Unmodified
Added
Removed
  • Singular/LIB/ffmodstd.lib

    r1adfacf rec66b3f  
    11///////////////////////////////////////////////////////////////////////////////
    2 //version="version ffmodstd.lib 4.1.1.0 Dec_2017 "; // $Id$
    3 version=" $ Id$ ";
     2version="version ffmodstd.lib 4.1.1.0 Sep_2018 "; // $Id$
    43category="Commutative Algebra";
    54info="
     
    76                           over rational function fields
    87AUTHORS:   D.K. Boku       boku@mathematik.uni-kl.de
    9 @*         W. Decker       decker@mathematik.uni-kl.de
    10 @*         C. Fieker       fieker@mathematik.uni-kl.de
     8           W. Decker       decker@mathematik.uni-kl.dei
     9           C. Fieker       fieker@mathematik.uni-kl.de
    1110
    1211OVERVIEW:
     
    5352  fareypoly(g,f);         univariate rational function reconstruction
    5453  polyInterpolation(l,m);  univariate polynomial interpolation
    55   modrationalInterpolation(l,m); modular univariate rational interpolation 
     54  modrationalInterpolation(l,m); modular univariate rational interpolation
    5655  BerlekampMassey(L,i);    Berlekamp/Massey algorithm
    5756  modberlekampMassey(L,i)  modular Berlekamp/Massey algorithm
     
    355354      the length of the sequence L is long enough. If this polynomial does not split into
    356355      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.
    358357REFERENCES:
    359358@* [1] E. Kaltofen and W.-s. Lee: Early termination in sparse interpolation
     
    405404       if(D1 == 0)
    406405       {
    407            if(j>1)//(2*l0) < j && 
     406           if(j>1)//(2*l0) < j &&
    408407           {
    409408               // early termination is triggered here
     
    413412           B1 = var(i)*B0;
    414413           l1 = l0;
    415            
    416414       }
    417415       else
     
    11351133           vL[j] = T[j][1];
    11361134           if(evbd < T[j][3])
    1137            { 
     1135           {
    11381136             evbd = T[j][3];
    11391137           }
     
    13841382///////////////////////////////////////////////////////////////////////////////
    13851383// +++++++++++++++++ std for one parameter begins here +++++++++++++++++++
    1386  
     1384
    13871385static proc test_fmodI(poly f,ideal I)
    13881386{
     
    14911489        }
    14921490    }
    1493     number Nd = 0; 
     1491    number Nd = 0;
    14941492    number Nn = 1;
    14951493    list Pl, Il, L, M;
     
    15751573
    15761574///////////////////////////////////////////////////////////////////////////////
    1577  
     1575
    15781576static proc firststd(ideal I, ideal cI, int in_value, list #)
    15791577{
     
    20552053       }
    20562054       v = (v0-q*v1)*N;
    2057        v0 = v1; 
    2058        v1 = v; 
     2055       v0 = v1;
     2056       v1 = v;
    20592057       f = g1;
    20602058       g1= r*N;
     
    21202118{
    21212119    // 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
    21232121    int i, j;
    21242122    poly f = fareyresult[1];
     
    21362134        {
    21372135            return(0);
    2138         }     
     2136        }
    21392137    }
    21402138    return(1);
     
    22712269    * cannot be satisfied for all i
    22722270    */
    2273     return(1); 
     2271    return(1);
    22742272}
    22752273
Note: See TracChangeset for help on using the changeset viewer.