Changeset ed6b03 in git


Ignore:
Timestamp:
Jun 8, 2022, 2:05:49 PM (2 years ago)
Author:
Hans Schoenemann <hannes@…>
Branches:
(u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', 'b4f17ed1d25f93d46dbe29e4b499baecc2fd51bb')
Children:
762afe1b4c00a799ba78e7d2cf486e090c6a5333
Parents:
1ead6f7a96c0b74c88cce4a314fd1d652c7fb0af
Message:
compiler warnings (cohomo.cc)
File:
1 edited

Legend:

Unmodified
Added
Removed
  • Singular/dyn_modules/cohomo/cohomo.cc

    r1ead6f7 red6b03  
    656656int pvert(poly p)
    657657{
    658   int i, j, vert=0;
     658  int i, vert=0;
    659659  for(i=currRing->N;i>0;i--)
    660660  {
     
    781781ideal sfreemon(ideal h,int deg)
    782782{
    783   int i,j,t;
     783  int j;
    784784  ideal temp;
    785785  temp=idInit(1,1);
     
    808808ideal id_sfmon(ideal h)
    809809{
    810   ideal asfmons,sfmons,mons,p;
     810  ideal asfmons,sfmons,mons;
    811811  int j, vert=idvert(h);
    812812  mons=id_MaxIdeal(1, currRing);
     
    856856bool IsInX(poly p,ideal X)
    857857{
    858   int i,j;
     858  int i;
    859859  for(i=0;i<IDELEMS(X);i++)
    860860  {
     
    878878{
    879879  ideal h,qrh;
    880   int i;
    881880  h=idAdd(h1,h2);
    882881  qrh=scKBase(deg,h);
     
    910909ideal idsrRing(ideal h)
    911910{
    912   int max,i,j,n;
     911  int i,n;
    913912  ideal pp,qq,rsr,ppp,hc=idCopy(h);
    914913  for(i=1;i<=rVar(currRing);i++)
     
    990989ideal IsSimplex(ideal h)
    991990{
    992   int i,j,ifbreak=0,max=id_maxdeg(h);
     991  int i,max=id_maxdeg(h);
    993992  poly e=pOne();
    994993  ideal id_re, id_so=idCopy(h);
     
    11071106  poly e=pOne();
    11081107  ideal h2=id_complement(h), aset=idInit(1,1);
    1109   int i,j,deg1=pTotaldegree(S);
     1108  int i,deg1=pTotaldegree(S);
    11101109  int tdeg=deg1+ddeg;
    11111110  if(tdeg!=0)
     
    14201419std::vector<int> ofindbases1(int num, int vnum, std::vector<int> bset,std::vector<std::vector<int> > gset)
    14211420{
    1422   int i,j,m;
    14231421  std::vector<std::vector<int> > goodset;
    14241422  std::vector<int> fvars=freevars(num,   bset,  gset), oset, base;
     
    14441442std::vector<std::vector<int> > ofindbases(int num,  std::vector<int> bset,std::vector<std::vector<int> > gset)
    14451443{
    1446   int i,j,m;
     1444  int i,m;
    14471445  std::vector<std::vector<int> > bases;
    14481446  std::vector<int> fvars=freevars(num,   bset,  gset), base1;
     
    14781476std::vector<std::vector<int> > eli2(int num, std::vector<int> bset,std::vector<std::vector<int> > gset)
    14791477{
    1480   int i,j;
    14811478  std::vector<int> badset;
    14821479  std::vector<std::vector<int> > goodset, solve;
     
    15721569ideal p_a(ideal h)
    15731570{
    1574   poly e=pOne(), p;
     1571  poly p;
    15751572  int i,j,deg=0,deg0;
    15761573  ideal aset=idCopy(h),ia,h1=idsrRing(h);
     
    16921689  std::vector<int> vert=vertset(lk), bv;
    16931690  res=b_subsets(vert);
    1694   int i, j, nu=res.size(), adg=pTotaldegree(a);
     1691  int i, adg=pTotaldegree(a);
    16951692  poly e=pOne();
    16961693  ideal idd=idInit(1,1);
     
    18581855poly pMake3(std::vector<int> vbase)
    18591856{
    1860   int n=vbase.size(),co=1;
     1857  int co=1;
    18611858  poly p,q=0;
    18621859  for(int i=0;i<3;i++)
     
    18941891void equmab(int num)
    18951892{
    1896   int i,j;
     1893  int i;
    18971894  //Print("There are %d new variables for equations solving.\n",num);
    18981895  ring r=currRing;
     
    21232120{
    21242121  //ring r=currRing;
    2125   int i;
    21262122  //assume (LIB "presolve.lib");
    21272123  sleftv a;a.Init();
     
    21532149std::vector<int> numfree(ideal h)
    21542150{
    2155   int i,j,num=0;
     2151  int i,j;
    21562152  std::vector<int> fvar;
    21572153  for(j=1;j<=currRing->N;j++)
     
    23262322std::vector<std::vector<int> > subspacet(std::vector<std::vector<int> > mv, std::vector<int> bv,std::vector<std::vector<int> > ntvs)
    23272323{
    2328   int i,j;
    23292324  std::vector<int> alset=findalpha(mv,bv), subase;
    23302325  std::vector<std::vector<int> > subases;
    2331   for(i=0;i<alset.size();i++)
     2326  for(unsigned i=0;i<alset.size();i++)
    23322327  {
    23332328    subase=subspacet1(alset[i],ntvs);
     
    25242519//returns true if pv union qv union av minus bv is in hvs
    25252520//hvs is simplicial complex
    2526 bool nabtconditionv(std::vector<std::vector<int> > hvs,  std::vector<int> pv, std::vector<int> qv, std::vector<int> av,  std::vector<int> bv)
     2521static bool nabtconditionv(std::vector<std::vector<int> > hvs,  std::vector<int> pv, std::vector<int> qv)
    25272522{
    25282523  std::vector<int> v1;
     
    25462541    {
    25472542      var.clear();
    2548       if(nabtconditionv(hvs, Nv[i], Nv[j], av, bv))
     2543      if(nabtconditionv(hvs, Nv[i], Nv[j]))
    25492544      {
    25502545        var.push_back(i);
     
    27862781      for(j=i+1;j<n;j++)
    27872782      {
    2788         if(nabtconditionv(hvs,nv[i],nv[j],av,bv))
     2783        if(nabtconditionv(hvs,nv[i],nv[j]))
    27892784        {
    27902785          good=listsinsertlist(good,i+1,j+1);
     
    28372832{
    28382833  ideal bi=findb(h),ai;
    2839   int mm=0,index=0;
     2834  int mm=0;
    28402835  id_print(bi);
    28412836  poly a,b;
     
    29862981ideal genst(ideal h, poly a, poly b)
    29872982{
    2988   int i,j;
    29892983  std::vector<std::vector<int> > hvs=supports(h),mv,mts;
    29902984  std::vector<int> av=support1(a), bv=support1(b);
     
    29932987  std::vector<std::vector<poly> > pvs=idMakei(mv,mts);
    29942988  ideal gens=idInit(1,1);
    2995   for(i=0;i<pvs.size();i++)
     2989  for(unsigned i=0;i<pvs.size();i++)
    29962990  {
    29972991    idInsertPoly(gens,pvs[i][0]);
     
    31223116        pWrite(b);
    31233117        intvec *solve=gradedpiece2n(h, a, b);
     3118        delete solve;
    31243119        gp++;
    31253120      }
     
    32263221      for(j=i+1;j<n;j++)
    32273222      {
    3228         if(nabtconditionv(hvs,nv[i],nv[j],av,bv))
     3223        if(nabtconditionv(hvs,nv[i],nv[j]))
    32293224        {
    32303225          good=listsinsertlist(good,i+1,j+1);
     
    32683263{
    32693264  t_start=clock();
    3270   int i,j,co;
    32713265  poly e=pOne();
    32723266  std::vector<int> av=support1(a),bv=support1(b),index, em;
     
    34223416intvec * gradedpiece2nl(ideal h,poly a,poly b)
    34233417{
    3424   int i,j,t;
    34253418  poly e=pOne();
    34263419  std::vector<int> av=support1(a), bv=support1(b), em;
     
    35043497std::vector<std::vector<int> > triface(poly p, int vert)
    35053498{
    3506   int i;
    35073499  std::vector<int> vec, fv=support1(p);
    35083500  std::vector<std::vector<int> > fvs0, fvs;
     
    35113503  fvs0=b_subsets(fv);
    35123504  fvs0=vsMinusv(fvs0,fv);
    3513   for(i=0;i<fvs0.size();i++)
     3505  for(unsigned i=0;i<fvs0.size();i++)
    35143506  {
    35153507    vec=fvs0[i];
     
    35193511  return (fvs);
    35203512}
    3521 
    3522 
    3523 
    3524 
    3525 
    3526 
    35273513
    35283514// the size of p's support must be 3
     
    35743560std::vector<int> commonedge(poly p, poly q)
    35753561{
    3576   int i,j;
    35773562  std::vector<int> ev, fv1= support1(p), fv2= support2(q);
    3578   for(i=0;i<fv1.size();i++)
     3563  for(unsigned i=0;i<fv1.size();i++)
    35793564  {
    35803565    if(IsinL(fv1[i], fv2))
     
    35883573{
    35893574  intvec *m;
    3590   int i,j;
     3575  int i;
    35913576  std::vector<int> dg=commonedge(p, q);
    35923577  int lg=dg.size();
     
    36063591std::vector<std::vector<int> > tetraface(poly p, poly q, int vert)
    36073592{
    3608   int i;
    36093593  std::vector<int> ev=commonedge(p, q), vec, fv1=support1(p), fv2=support1(q);
    36103594  std::vector<std::vector<int> > fvs1, fvs2, fvs;
     
    36173601  fvs2=vsUnion(fvs1, fvs2);
    36183602  fvs2=vsMinusv(fvs2, ev);
    3619   for(i=0;i<fvs2.size();i++)
     3603  for(unsigned i=0;i<fvs2.size();i++)
    36203604  {
    36213605    vec=fvs2[i];
     
    36303614ideal triangulations2(ideal h, poly p, poly q, int vert)
    36313615{
    3632   int i,j;
    36333616  std::vector<int> ev, fv1=support1(p), fv2=support1(q);
    36343617  std::vector<std::vector<int> > vecs=supports(h), vs1;
     
    36493632std::vector<std::vector<int> > penface(poly p, poly q, poly g, int vert)
    36503633{
    3651   int i, en=0;
     3634  int en=0;
    36523635  std::vector<int> ev1=commonedge(p, q), ev2=commonedge(p, g), ev3=commonedge(q, g), ind, vec, fv1=support1(p), fv2=support1(q), fv3=support1(g);
    36533636  std::vector<std::vector<int> > fvs1, fvs2, fvs3, fvs, evec;
     
    36553638  evec.push_back(ev2);
    36563639  evec.push_back(ev3);
    3657   for(i=0;i<evec.size();i++)
     3640  for(unsigned i=0;i<evec.size();i++)
    36583641  {
    36593642    if(evec[i].size()==2)
     
    36743657    fvs3=vsUnion(fvs3, fvs2);
    36753658    fvs3=vsUnion(fvs3, fvs1);
    3676     for(i=0;i<evec.size();i++)
     3659    for(unsigned i=0;i<evec.size();i++)
    36773660    {
    36783661      if(evec[i].size()==2)
     
    36813664      }
    36823665    }
    3683     for(i=0;i<fvs3.size();i++)
     3666    for(unsigned i=0;i<fvs3.size();i++)
    36843667    {
    36853668      vec=fvs3[i];
     
    36953678ideal triangulations3(ideal h, poly p, poly q, poly g, int vert)
    36963679{
    3697   int i,j;
    36983680  std::vector<int> ev1=commonedge(p, q), ev2=commonedge(p, g), ev3=commonedge(q, g), fv1=support1(p), fv2=support1(q), fv3=support1(g);
    36993681  std::vector<std::vector<int> > vecs=supports(h), vs1, evec;
     
    37013683  evec.push_back(ev2);
    37023684  evec.push_back(ev3);
    3703   for(i=0;i<evec.size();i++)
     3685  for(unsigned i=0;i<evec.size();i++)
    37043686  {
    37053687    if(evec[i].size()==2)
     
    37223704int valency(ideal h, poly p)
    37233705{
    3724   int i, val=0;
     3706  int val=0;
    37253707  std::vector<int> ev=support1(pCopy(p));
    37263708  int ver=ev[0];
    37273709//PrintS("the vertex is :\n"); listprint(p);
    37283710  std::vector<std::vector<int> > vecs=supports(idCopy(h));
    3729   for(i=0;i<vecs.size();i++)
     3711  for(unsigned i=0;i<vecs.size();i++)
    37303712  {
    37313713    if(vecs[i].size()==2 && IsinL(ver, vecs[i]))
     
    37723754std::vector<std::vector<int> > vsMinusvs(std::vector<std::vector<int> > vs1, std::vector<std::vector<int> > vs2)
    37733755{
    3774   int i;
    37753756  std::vector<std::vector<int> > vs=vs1;
    3776   for(i=0;i<vs2.size();i++)
     3757  for(unsigned i=0;i<vs2.size();i++)
    37773758  {
    37783759    vs=vsMinusv(vs, vs2[i]);
     
    37853766{
    37863767  std::vector<std::vector<int> >  sset, bv;
    3787   for(int i=0;i<vs.size();i++)
     3768  for(unsigned i=0;i<vs.size();i++)
    37883769  {
    37893770    bv=b_subsets(vs[i]);
     
    38493830ideal c_New(ideal Io, ideal sig)
    38503831{
    3851   poly p, q, g;
    38523832  std::vector<std::vector<int> > vs1=p_constant(Io, sig), vs2=p_change(sig), vs3=p_new(Io, sig), vsig=supports(sig), vs;
    38533833  std::vector<int> ev;
     
    38903870  std::vector<std::vector<int> > ss=supports(idCopy(Sigma)), fvs;
    38913871  std::vector<int> av=support1(a), intvec, vv;
    3892   for(int i=0;i<ss.size();i++)
     3872  for(unsigned i=0;i<ss.size();i++)
    38933873  {
    38943874    intvec=vecIntersection(ss[i], av);
     
    39043884
    39053885
    3906 std::vector<std::vector<int> > phi2(poly a, ideal Xo, ideal Sigma, int vert)
     3886static std::vector<std::vector<int> > phi2(poly a, ideal Xo, ideal Sigma)
    39073887{
    39083888
     
    39323912  {
    39333913    lk1=phi1(a, Sigma);
    3934     lk2=phi2(a, Xo, Sigma, vert);
     3914    lk2=phi2(a, Xo, Sigma);
    39353915    lkn=vsMinusvs(lko, lk1);
    39363916    lkn=vsUnion(lkn, lk2);
     
    39393919  if(ord==3)
    39403920  {
    3941     lkn=phi2(a, Xo, Sigma, vert);
     3921    lkn=phi2(a, Xo, Sigma);
    39423922    return lkn;
    39433923  }
     
    39553935  std::vector<std::vector<int> > xvs=supports(idCopy(Xs)), bs=b_subsets(bv);
    39563936  bs=vsMinusv(bs, bv);
    3957   for(int i=0;i<bs.size();i++)
     3937  for(unsigned i=0;i<bs.size();i++)
    39583938  {
    39593939    if(!vInvsl(bs[i], xvs))
     
    39883968int ifIso(poly p, poly q, poly f, poly g, poly a, poly b)
    39893969{
    3990   int i;
    39913970  std::vector<int> va=support1(a), vb=support1(b), vp=support1(p),  vq=support1(q), vf=support1(f), vg=support1(g);
    39923971  std::vector<int>   v1=phimagel(vp, va, vb), v2=phimagel(vq, va, vb), v3=phimagel(vf, va, vb), v4=phimagel(vg, va, vb);
     
    40043983{
    40053984  ideal h=idInit(1,1);
    4006   int i,j,eq=0;
     3985  int i;
    40073986  for(i=0;i<IDELEMS(I);i++)
    40083987  {
     
    40304009{
    40314010  std::vector<int> vec;
    4032   for(int i=0;i<v1.size();i++)
     4011  for(unsigned i=0;i<v1.size();i++)
    40334012  {
    40344013    vec.push_back(v1[i]-v2[i]);
     
    40404019std::vector<int> gdegree(poly a, poly b)
    40414020{
    4042   int i,j;
     4021  int i;
    40434022  std::vector<int> av,bv;
    40444023  for(i=1;i<=currRing->N;i++)
     
    42734252{
    42744253  intvec *m;
    4275   int i,j;
     4254  int i;
    42764255  std::vector<int> dg=gdegree(a,b);
    42774256  int lg=dg.size();
Note: See TracChangeset for help on using the changeset viewer.