Changeset a1bbf6 in git


Ignore:
Timestamp:
Mar 31, 2015, 10:47:36 AM (8 years ago)
Author:
Hans Schoenemann <hannes@…>
Branches:
(u'spielwiese', '8e0ad00ce244dfd0756200662572aef8402f13d5')
Children:
2851932a1bfd822c8a77ba701ab4061fda4db650
Parents:
f6cd60bc278fb9a300668ad9a5c51e492c8b6409320644e83870d617264b4e83eea678323a62d15a
Message:
Merge branch 'tigistdereje-algemodstd' into spielwiese
Files:
4 added
4 deleted
4 edited
3 moved

Legend:

Unmodified
Added
Removed
  • Singular/LIB/nfmodstd.lib

    rf6cd60 ra1bbf6  
    11////////////////////////////////////////////////////////////////////////////////
    2 version="version algemodstd.lib 4.0.1.0 Sep_2014 ";  // $Id$
     2version="version nfmodstd.lib 4.0.1.0 Sep_2014 ";  // $Id$
    33category="Commutative Algebra";
    44info="
    55
    6 LIBRARY:   algemodstd.lib  Groebner bases of ideals in polynomial rings
     6LIBRARY:   nfmodstd.lib  Groebner bases of ideals in polynomial rings
    77                           over algebraic number fields
    88AUTHORS:   D.K. Boku       boku@mathematik.uni-kl.de
     
    3333PROCEDURES:
    3434  chinrempoly(l,m);       chinese remaindering for polynomials
    35   algemodStd(I);          standard basis of I over algebraic number field using modular methods
     35  nfmodStd(I);          standard basis of I over algebraic number field using modular methods
    3636";
    3737
     
    443443    intvec opt = option(get);
    444444    option(redSB);
    445     I = modular("Algemodstd::LiftPolyCRT", list(I), PrimeTestTask, Modstd::deleteUnluckyPrimes_std,
     445    I = modular("Nfmodstd::LiftPolyCRT", list(I), PrimeTestTask, Modstd::deleteUnluckyPrimes_std,
    446446              PtestStd, Modstd::finalTest_std,536870909);
    447447    attrib(I, "isSB", 1);
     
    452452////////////////////////////////////////////////////////////////////////////////
    453453/* main procedure */
    454 proc algemodStd(ideal I, list #)
    455 "USAGE:  algemodStd(I, #); I ideal, # optional parameters
     454proc nfmodStd(ideal I, list #)
     455"USAGE:  nfmodStd(I, #); I ideal, # optional parameters
    456456RETURN:  standard basis of I over algebraic number field
    457457NOTE: The procedure passes to @ref{modStd} if the ground field has no
     
    459459      are directly passed to @ref{modStd}.
    460460SEE ALSO: modStd
    461 EXAMPLE: example algemodStd; shows an example
     461EXAMPLE: example nfmodStd; shows an example
    462462"
    463463{
     
    512512    minpoly =a^2+1;
    513513    ideal k=(a/2+1)*x^2+2/3y, 3*x-a*y+ a/7+2;
    514     ideal I=algemodStd(k);
     514    ideal I=nfmodStd(k);
    515515    I;
    516516    ring r2 =(0,a),(x,y,z),dp;
    517517    minpoly =a^3 +2;
    518518    ideal k=(a^2+a/2)*x^2+(a^2 -2/3*a)*yz, (3*a^2+1)*zx-(a+4/7)*y+ a+2/5;
    519     ideal IJ=algemodStd(k);
     519    ideal IJ=nfmodStd(k);
    520520    IJ;
    521521    ring r3=0,(x,y),dp;// ring without parameter
    522522    ideal I = x2 + y, xy - 7y + 2x;
    523     I=algemodStd(I);
     523    I=nfmodStd(I);
    524524    I;
    525525}
  • Tst/Long/nfmodstd_l.tst

    rf6cd60 ra1bbf6  
    11LIB "tst.lib";
    2 LIB "algemodstd.lib";
     2LIB "nfmodstd.lib";
    33tst_init();
    4 proc tst_test_algemodstd(ideal I)
     4proc tst_test_nfmodstd(ideal I)
    55{
    6    ideal Jtst = algemodStd(I);
     6   ideal Jtst = nfmodStd(I);
    77   Jtst;
    88}
     
    1616          -x*w^7+y*w^7];
    1717
    18  tst_test_algemodstd(HJa);
     18 tst_test_nfmodstd(HJa);
    1919kill r;
    2020ring r=(0,a),(x,y,z,w,u,v,n),dp;
     
    2424          8*x^3+(12*a+a^2)*y^3 +x*z^2+3,
    2525          (7*a^2-a)*x^2*y^4+18*y^3*z^2+y^3*z^3];
    26 tst_test_algemodstd(I1);
     26tst_test_nfmodstd(I1);
    2727
    2828ideal HUa = [a*x+(a-1)*y+(a+2)*w+u,x*y+(a-1)*y*z+z*w + (a+2)*x*w+a*w*u,
    2929            a*x*y*z+(a+5)*y*z*w+a*x*y*u+(a+2)*x*w*u+a*z*w*u,(a-11)*x*y*z*w+(a+5)*x*y*z*u+a*x*y*w*u +
    3030            a*x*z*w*u+a*y*z*w*u,(a + 3)*x*y*z*w*u -(a+23)*v^5];
    31 tst_test_algemodstd(HUa);
     31tst_test_nfmodstd(HUa);
    3232kill r;
    3333ring r=(0,a),(x,y,z,w,u,v,n),dp;
     
    3636          + (a+1)*z*w*u+x*y*v+x*u*v+w*u*v,(a-1)*x*y*z*w + y*z*w*u +x*y*z*v +x*y*u*v+ x*w*u*v+ z*w*u*v,
    3737          x*y*z*w*u+(a+1)*x*y*z*w*v+x*y*z*u*v+x*y*w*u*v+x*z*w*u*v+y*z*w*u*v,x*y*z*w*u*v-a+2];
    38 tst_test_algemodstd(I6a);
     38tst_test_nfmodstd(I6a);
    3939kill r;
    4040ring r=(0,a),(x,y,z,w,u,v,n),dp;
     
    4343            2*w*y^3-10*y^2 - 10*w*y,(w^2-a^2+2*a)*x + 11*w*y*z-2*z, (-z^2 + 4*w^2+a^2+2)*x*z
    4444            + (4*w*z^2 + 2*w^3-10*w)*y + 4*z^2-10*w^2 + a^2+a];
    45 tst_test_algemodstd(I);
     45tst_test_nfmodstd(I);
    4646ideal I1= [(2*w*y-(2*a+3)*z^2)*x-z*y^2-z^2*x,
    4747            2*z*x^2*w + (4*w*y + 5*z^2)*x^2 + (4*z*y^2 + 4*z^2*w)*x + 2*w*y^3-10*x^2*y^2 - 10*x*w*y*z,
    4848            (w^2-(a^2+a)*y*z)*x + 11*w*y*z-2*z^3,
    4949            (-z^2 + 4*w^2 + (a^2+2)*y^2)*x*z + (4*w*z^2 + 2*w^3-10*w^3)*y + 4*z^2*x*y-10*w^2*x*z + (a^2+2)*y^3*w];
    50             tst_test_algemodstd(I1);
     50            tst_test_nfmodstd(I1);
    5151kill r;
    5252ring r=(0,a),(x,y,z,w,u),dp;
     
    5555          (a + 2)*x*w + a*w*u,a*x*y*z + (a + 5)*y*z*w + a*x*y*u + (a+2)*x*w*u + a*z*w*u,
    5656          (a-11)*x*y*z*w + (a + 5)*x*y*z*u + a*x*y*w*u + a*x*z*w*u + a*y*z*w*u,(a + 3)*x*y*z*w*u + a+23];
    57 tst_test_algemodstd(Ua);
     57tst_test_nfmodstd(Ua);
    5858kill r;
    5959
  • Tst/Long/ok_l2.lst

    rf6cd60 ra1bbf6  
    1 algemodstd_l
     1nfmodstd_l
    22bug_tr642
    33bug_tr677
  • Tst/Short.lst

    rf6cd60 ra1bbf6  
    22Short/alexpoly.tst
    33Short/algebralib.tst
    4 Short/algemodstd_s
     4Short/nfmodstd_s
    55Short/allres_s.tst
    66Short/arcAtPoint.tst
  • Tst/Short/nfmodstd_s.tst

    rf6cd60 ra1bbf6  
    11LIB "tst.lib";
    2 LIB "algemodstd.lib";
     2LIB "nfmodstd.lib";
    33tst_init();
    44
    5 proc tst_test_algemodstd(ideal I)
     5proc tst_test_nfmodstd(ideal I)
    66{
    7    ideal Jtst = algemodStd(I);
     7   ideal Jtst = nfmodStd(I);
    88   Jtst;
    99}
     
    1515          8*x^3+(12+a)*y^3+x*z^2+3,
    1616          (7-a)*x^2*y^4+18*y^3*z^2+y^3*z^3];
    17 tst_test_algemodstd(HIa);
     17tst_test_nfmodstd(HIa);
    1818ideal HJa=[(2*a+3)*x*y^4*z^2 + (a+2)*x^2*y^3*z*w-x^2*y^3*z*w + 2*x*y*z^2*w^3 + (7*a-1)*y^3*w^4,
    1919          2*x^2*y^4*z + x^2*y*z^2*w^2-a*x*y^2*z^2*w^2 + (a+11)*x^2*y*z*w^3 -12*x*w^6 + 12*y*w^6,
     
    2222          -x*w^7+y*w^7];
    2323
    24  tst_test_algemodstd(HJa);
     24 tst_test_nfmodstd(HJa);
    2525ideal HBa =[x+(a/2147483647)*x*y+(a/2147483646)*y*z,
    2626          x^2*z + (a/2147483647)*y^3 + (a+1)*y*z^2,
    2727          (2*a*x)/3-(7*a*y)/23 + (9*z)/7];
    2828
    29   tst_test_algemodstd(HBa);
     29  tst_test_nfmodstd(HBa);
    3030kill r;
    3131ring r=(0,a),(x,y,z,w,u),dp;
     
    3434        (a+3)*x*y^2*z^2 + (a+1)*x^5 + (a-11)*y^2*z^2,
    3535        (a+2)*x*y*z + (a+7)*x^3 + 12*y^3 + 1,3*x^3 +(a- 4)*y^3 + y*z^2];
    36 tst_test_algemodstd(ka);
     36tst_test_nfmodstd(ka);
    3737
    3838ideal Hka=[(a+5)*x^3*y^2*z + (a-3)*y^3*x^2*z + (a+7)*x*y^2*z^2,
     
    4040          4*x*y*z + (a+7)*x^3 + 12*y^3 + a*z^3,
    4141          3*x^3 +(a- 4)*y^3 + y*z^2];
    42 tst_test_algemodstd(Hka);
     42tst_test_nfmodstd(Hka);
    4343
    4444ideal Ua = [a*x + (a-1)*y + z + (a+2)*w + u,x*y + (a-1)*y*z + z*w +
    4545          (a + 2)*x*w + a*w*u,a*x*y*z + (a + 5)*y*z*w + a*x*y*u + (a+2)*x*w*u + a*z*w*u,
    4646          (a-11)*x*y*z*w + (a + 5)*x*y*z*u + a*x*y*w*u + a*x*z*w*u + a*y*z*w*u,(a + 3)*x*y*z*w*u + a+23];
    47 tst_test_algemodstd(Ua);
     47tst_test_nfmodstd(Ua);
    4848kill r;
    4949
  • libpolys/polys/nc/gb_hack.h

    r320644 ra1bbf6  
    5151D(ideal sca_mora(const ideal, const ideal, const intvec *, const intvec *, kStrategy, const ring _currRing))
    5252
    53 D(poly kNF(ideal, ideal, poly, int, int, const ring _currRing))
     53D(poly k_NF(ideal, ideal, poly, int, int, const ring _currRing))
    5454
    5555#endif // # ifdef PLURAL_INTERNAL_DECLARATIONS_GB_HACK
Note: See TracChangeset for help on using the changeset viewer.