Changeset af598e in git
- Timestamp:
- May 4, 2011, 7:50:23 PM (12 years ago)
- Branches:
- (u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', '0604212ebb110535022efecad887940825b97c3f')
- Children:
- 107b6a63f83a561264dff51b21ba9a08c2b381f8
- Parents:
- 02c42db0b19db524ec869930bffc30346b16abed
- git-author:
- Oleksandr Motsak <motsak@mathematik.uni-kl.de>2011-05-04 19:50:23+02:00
- git-committer:
- Mohamed Barakat <mohamed.barakat@rwth-aachen.de>2011-11-09 12:31:35+01:00
- Location:
- libpolys/polys
- Files:
-
- 3 edited
- 2 moved
Legend:
- Unmodified
- Added
- Removed
-
libpolys/polys/Makefile.am
r02c42d raf598e 37 37 pDebug.cc pInline0.cc polys0.cc prCopy.cc \ 38 38 kbuckets.cc sbuckets.cc weight.cc weight0.c simpleideals.cc matpol.cc \ 39 ${USE_P_PROCS_STATIC_CC} ${USE_P_PROCS_DYNAMIC_CC} templates/mod_raw.cc 39 ${USE_P_PROCS_STATIC_CC} ${USE_P_PROCS_DYNAMIC_CC} templates/mod_raw.cc \ 40 clapsing.cc 40 41 41 42 BUILT_SOURCES = templates/p_Procs.inc … … 51 52 templates/p_Procs_Dynamic.h templates/p_Procs_Impl.h templates/p_Procs_Set.h templates/p_Procs_Static.h \ 52 53 monomials/p_polys.h monomials/polys-impl.h monomials/maps.h polys.h prCopy.h prCopyMacros.h \ 53 kbuckets.h sbuckets.h simpleideals.h weight.h \54 matpol.h54 kbuckets.h sbuckets.h simpleideals.h weight.h matpol.h \ 55 clapsing.h 55 56 56 57 P_PROCS_CPPFLAGS_COMMON = -DHAVE_CONFIG_H -DDYNAMIC_VERSION -
libpolys/polys/clapsing.cc
r02c42d raf598e 9 9 10 10 //#define FACTORIZE2_DEBUG 11 #include <kernel/mod2.h> 11 #include "config.h" 12 #include <misc/auxiliary.h> 13 14 15 TODO(Martin, Please adapt the following code for the use in SW) 16 #ifdef HAVE_FACTORY 17 18 #define SI_DONT_HAVE_GLOBAL_VARS 12 19 #include <omalloc/omalloc.h> 13 #ifdef HAVE_FACTORY 14 #define SI_DONT_HAVE_GLOBAL_VARS 15 #include <kernel/structs.h> 16 #include <kernel/clapsing.h> 17 #include <kernel/numbers.h> 18 #include <kernel/ring.h> 19 #include <kernel/ideals.h> 20 #include <kernel/ffields.h> 20 #include <coeffs/numbers.h> 21 #include <coeffs/coeffs.h> 22 23 // #include <kernel/ffields.h> 24 // #include <kernel/clapconv.h> 25 // #include <libfac/factor.h> 26 21 27 #include <factory/factory.h> 22 #include <kernel/clapconv.h> 23 #include <libfac/factor.h> 24 #include <kernel/ring.h> 28 29 #include "clapsing.h" 30 #include "monomials/ring.h" 31 #include "simpleideals.h" 32 25 33 26 34 void out_cf(const char *s1,const CanonicalForm &f,const char *s2); … … 1652 1660 return mm; 1653 1661 } 1654 #endif 1655 1662 /* 1656 1663 napoly singclap_alglcm ( napoly f, napoly g ) 1657 1664 { … … 1735 1742 Off(SW_RATIONAL); 1736 1743 } 1744 */ 1737 1745 1738 1746 #if 0 -
libpolys/polys/clapsing.h
r02c42d raf598e 12 12 #define INCL_FACTORYSING_H 13 13 14 #include <kernel/structs.h> 15 #include <kernel/polys.h> 16 #include <kernel/ideals.h> 17 #include <kernel/intvec.h> 18 #include <kernel/matpol.h> 19 #include <kernel/longtrans.h> 14 // #include <kernel/structs.h> 15 #include <polys/monomials/p_polys.h> 16 #include <polys/monomials/ring.h> 17 #include <polys/simpleideals.h> 18 #include <misc/intvec.h> 19 #include <polys/matpol.h> 20 //#include <kernel/longtrans.h> 20 21 21 22 poly singclap_gcd ( poly f, poly g ); 22 23 poly singclap_gcd_r ( poly f, poly g, const ring r ); 24 /* 23 25 napoly singclap_alglcm ( napoly f, napoly g ); 24 26 void singclap_algdividecontent ( napoly f, napoly g, napoly &ff, napoly &gg ); 27 */ 25 28 26 29 poly singclap_resultant ( poly f, poly g , poly x); … … 53 56 poly singclap_det( const matrix m ); 54 57 int singclap_det_i( intvec * m ); 55 58 /* 56 59 BOOLEAN jjRESULTANT(leftv res, leftv u, leftv v, leftv w); 57 60 BOOLEAN jjCHARSERIES(leftv res, leftv u); … … 59 62 BOOLEAN jjIS_SQR_FREE(leftv res, leftv u); 60 63 #endif 64 */ 61 65 62 66 #endif /* INCL_FACTORYSING_H */ -
libpolys/polys/monomials/p_polys.cc
r02c42d raf598e 20 20 #include <coeffs/longrat.h> // ??? 21 21 22 #include <polys/weight.h>23 #include <polys/simpleideals.h>24 25 #include <polys/monomials/ring.h>26 #include <polys/monomials/p_polys.h>22 #include "weight.h" 23 #include "simpleideals.h" 24 25 #include "monomials/ring.h" 26 #include "monomials/p_polys.h" 27 27 28 28 // #include <???/ideals.h> … … 34 34 35 35 #ifdef HAVE_PLURAL 36 #include <polys/nc/nc.h> 37 #include <polys/nc/sca.h> 38 #endif 39 40 #include <polys/coeffrings.h> 36 #include "nc/nc.h" 37 #include "nc/sca.h" 38 #endif 39 40 #include "coeffrings.h" 41 #include "clapsing.h" 42 41 43 42 44 /*************************************************************** -
libpolys/polys/simpleideals.cc
r02c42d raf598e 11 11 #include "config.h" 12 12 #include <misc/auxiliary.h> 13 14 #include <omalloc/omalloc.h> 15 13 16 #include <misc/options.h> 14 #include <omalloc/omalloc.h>15 #include <polys/monomials/p_polys.h>16 #include <polys/weight.h>17 #include <polys/matpol.h>18 17 #include <misc/intvec.h> 19 #include <polys/simpleideals.h> 18 19 #include <coeffs/longrat.h> 20 21 #include "monomials/p_polys.h" 22 #include "weight.h" 23 #include "matpol.h" 24 #include "simpleideals.h" 20 25 #include "sbuckets.h" 21 26 … … 1664 1669 1665 1670 #ifdef HAVE_FACTORY 1671 #if 0 1666 1672 poly id_GCD(poly f, poly g, const ring r) 1667 1673 { … … 1670 1676 ideal I=idInit(2,1); I->m[0]=f; I->m[1]=g; 1671 1677 intvec *w = NULL; 1678 1672 1679 ideal S=idSyzygies(I,testHomog,&w); 1680 1673 1681 if (w!=NULL) delete w; 1674 1682 poly gg=pTakeOutComp(&(S->m[0]),2); … … 1680 1688 } 1681 1689 #endif 1690 #endif 1682 1691 1683 1692 /*2 … … 1689 1698 */ 1690 1699 #ifdef HAVE_FACTORY 1691 ideal id ChineseRemainder(ideal *xx, number *q, int rl, const ring r)1700 ideal id_ChineseRemainder(ideal *xx, number *q, int rl, const ring R) 1692 1701 { 1693 1702 int cnt=IDELEMS(xx[0])*xx[0]->nrows; … … 1708 1717 h=xx[j]->m[i]; 1709 1718 if ((h!=NULL) 1710 &&((r==NULL)||(p LmCmp(r,h)==-1)))1719 &&((r==NULL)||(p_LmCmp(r,h,R)==-1))) 1711 1720 r=h; 1712 1721 } 1713 1722 if (r==NULL) break; 1714 h=p Head(r);1723 h=p_Head(r, R); 1715 1724 for(j=rl-1;j>=0;j--) 1716 1725 { 1717 1726 hh=xx[j]->m[i]; 1718 if ((hh!=NULL) && (p LmCmp(r,hh)==0))1727 if ((hh!=NULL) && (p_LmCmp(r,hh, R)==0)) 1719 1728 { 1720 x[j]=p GetCoeff(hh);1721 hh=p LmFreeAndNext(hh);1729 x[j]=p_GetCoeff(hh, R); 1730 hh=p_LmFreeAndNext(hh, R); 1722 1731 xx[j]->m[i]=hh; 1723 1732 } 1724 1733 else 1725 x[j]=nlInit(0, r); 1726 } 1727 number n=nlChineseRemainder(x,q,rl); 1734 x[j]=nlInit(0, R->cf); // is R->cf really n_Q??? 1735 } 1736 1737 number n=nlChineseRemainder(x,q,rl); // kernel/clapconv.cc 1738 1728 1739 for(j=rl-1;j>=0;j--) 1729 1740 { 1730 1741 x[j]=NULL; // nlInit(0...) takes no memory 1731 1742 } 1732 if (nlIsZero(n )) pDelete(&h);1743 if (nlIsZero(n, R->cf)) p_Delete(&h, R); 1733 1744 else 1734 1745 { 1735 p SetCoeff(h,n);1746 p_SetCoeff(h,n, R); 1736 1747 //Print("new mon:");pWrite(h); 1737 res_p=p Add(res_p,h);1748 res_p=p_Add_q(res_p, h, R); 1738 1749 } 1739 1750 } … … 1741 1752 } 1742 1753 omFree(x); 1743 for(i=rl-1;i>=0;i--) id Delete(&(xx[i]));1754 for(i=rl-1;i>=0;i--) id_Delete(&(xx[i]), R); 1744 1755 omFree(xx); 1745 1756 return result;
Note: See TracChangeset
for help on using the changeset viewer.