Changeset 95a42a2 in git for Singular/LIB/freegb.lib


Ignore:
Timestamp:
Mar 16, 2018, 6:11:32 PM (6 years ago)
Author:
Karim Abou Zeid <karim23697@…>
Branches:
(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
Children:
0be0394f8808f23d95c9df90daf22918ecc53e49
Parents:
3ea833075e746ca13c85f94a5894c0c7fd4e6a56
Message:
Fix lpNF bug where degbound must be larger than neccessary
File:
1 edited

Legend:

Unmodified
Added
Removed
  • Singular/LIB/freegb.lib

    r3ea833 r95a42a2  
    34813481
    34823482static proc getlpCoeffs(poly q, poly p)
    3483 "
    3484 "
    34853483{list R; poly m; intvec cq,t,lv,rv,bla;
    34863484 int n = attrib(basering,"lV"); int d = attrib(basering,"uptodeg");
     
    35053503"
    35063504{poly l,r,qt; int i;
    3507  g = shiftPoly(g,s);
    3508  list K = getlpCoeffs(lead(g),lead(p));
     3505 list K = getlpCoeffs(lead(shiftPoly(g,s)), lead(p));
    35093506 l = K[1]; r = K[2];
    35103507 kill K;
    35113508 for (i = 1; i <= size(g); i++)
    3512  {qt = qt + lpMult(lpMult(l,g[i]),r);
     3509 {
     3510   qt = qt + lpMult(lpMult(l,g[i]),r);
    35133511 }
    35143512 return(p - leadcoef(p)*normalize(qt));
Note: See TracChangeset for help on using the changeset viewer.