Changeset 19addd1 in git for gfanlib/gfanlib_zcone.cpp


Ignore:
Timestamp:
Oct 25, 2013, 6:22:01 PM (10 years ago)
Author:
Yue Ren <ren@…>
Branches:
(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
Children:
69b2c1754cbeecd7c1b844e71e2101e1b3794062
Parents:
b7b2c516034df879b01c7c8666bbdf0928da91dc
git-author:
Yue Ren <ren@mathematik.uni-kl.de>2013-10-25 18:22:01+02:00
git-committer:
Yue Ren <ren@mathematik.uni-kl.de>2013-10-25 18:25:28+02:00
Message:
fix: warnings
File:
1 edited

Legend:

Unmodified
Added
Removed
  • gfanlib/gfanlib_zcone.cpp

    rb7b2c5 r19addd1  
    4444    dd_colrange d_input,j;
    4545    dd_RepresentationType rep=dd_Inequality;
    46     dd_boolean found=dd_FALSE, newformat=dd_FALSE, successful=dd_FALSE;
    47     char command[dd_linelenmax], comsave[dd_linelenmax];
     46    // dd_boolean found=dd_FALSE, newformat=dd_FALSE, successful=dd_FALSE;
     47    // char command[dd_linelenmax], comsave[dd_linelenmax];
    4848    dd_NumberType NT;
    4949
    5050    (*Error)=dd_NoError;
    5151
    52     rep=dd_Inequality; newformat=dd_TRUE;
     52    rep=dd_Inequality; // newformat=dd_TRUE;
    5353
    5454    m_input=g.getHeight();
     
    7070    }
    7171
    72     successful=dd_TRUE;
     72    // successful=dd_TRUE;
    7373
    7474    return M;
     
    108108  {
    109109    bool ret;
    110     dd_MatrixPtr M=NULL,M2=NULL,M3=NULL;
    111     dd_colrange d;
     110    dd_MatrixPtr M=NULL/*,M2=NULL,M3=NULL*/;
     111    // dd_colrange d;
    112112    dd_ErrorType err=dd_NoError;
    113     dd_rowset redrows,linrows,ignoredrows, basisrows;
    114     dd_colset ignoredcols, basiscols;
    115     dd_DataFileType inputfile;
    116     FILE *reading=NULL;
     113    // dd_rowset redrows,linrows,ignoredrows, basisrows;
     114    // dd_colset ignoredcols, basiscols;
     115    // dd_DataFileType inputfile;
     116    // FILE *reading=NULL;
    117117
    118118    cddinitGmp();
     
    121121    if (err!=dd_NoError) goto _L99;
    122122
    123     d=M->colsize;
     123    // d=M->colsize;
    124124
    125125    static dd_Arow temp;
     
    339339    if(numberOfRows==0)return;//the full space, so description is already irredundant
    340340
    341     dd_rowset r=NULL;
     341    // dd_rowset r=NULL;
    342342    ZMatrix g=inequalities;
    343343    g.append(equations);
    344344
    345     dd_LPSolverType solver=dd_DualSimplex;
     345    // dd_LPSolverType solver=dd_DualSimplex;
    346346    dd_MatrixPtr A=NULL;
    347347    dd_ErrorType err=dd_NoError;
     
    402402    int numberOfRows=numberOfEqualities+numberOfInequalities;
    403403
    404     dd_rowset r=NULL;
     404    // dd_rowset r=NULL;
    405405    ZMatrix g=inequalities;
    406406    g.append(equations);
     
    448448  void dual(ZMatrix const &inequalities, ZMatrix const &equations, ZMatrix &dualInequalities, ZMatrix &dualEquations)
    449449  {
    450     int result;
     450    // int result;
    451451
    452452    dd_MatrixPtr A=NULL;
     
    472472
    473473    return;
    474    _L99:
    475     assert(0);
     474   // _L99:
     475   // assert(0);
    476476  }
    477477  // this procedure is take from cddio.c.
     
    556556    int dim2=inequalities.getHeight();
    557557    if(dim2==0)return std::vector<std::vector<int> >();
    558     int dimension=inequalities.getWidth();
     558    // int dimension=inequalities.getWidth();
    559559
    560560    dd_MatrixPtr A=NULL;
     
    603603
    604604    return ret;
    605    _L99:
    606     assert(0);
    607     return std::vector<std::vector<int> >();
     605   // _L99:
     606   // assert(0);
     607   // return std::vector<std::vector<int> >();
    608608  }
    609609
     
    707707          inequalities=LpSolver::fastNormals(inequalities2);
    708708          goto noFallBack;
    709         fallBack://alternativ (disabled)
    710           lpSolver.removeRedundantRows(inequalities,equations,true);
     709        // fallBack://alternativ (disabled)
     710        //   lpSolver.removeRedundantRows(inequalities,equations,true);
    711711        noFallBack:;
    712712        }
     
    730730}
    731731
    732 std::ostream &operator<<(std::ostream &f, ZCone const &c)
     732void operator<<(std::ostream &f, ZCone const &c)
    733733{
    734734  f<<"Ambient dimension:"<<c.n<<std::endl;
     
    741741
    742742ZCone::ZCone(int ambientDimension):
     743  preassumptions(PCP_impliedEquationsKnown|PCP_facetsKnown),
     744  state(1),
     745  n(ambientDimension),
     746  multiplicity(1),
     747  linearForms(ZMatrix(0,ambientDimension)),
    743748  inequalities(ZMatrix(0,ambientDimension)),
    744749  equations(ZMatrix(0,ambientDimension)),
    745   n(ambientDimension),
    746   state(1),
    747   preassumptions(PCP_impliedEquationsKnown|PCP_facetsKnown),
     750  haveExtremeRaysBeenCached(false)
     751{
     752}
     753
     754
     755ZCone::ZCone(ZMatrix const &inequalities_, ZMatrix const &equations_, int preassumptions_):
     756  preassumptions(preassumptions_),
     757  state(0),
     758  n(inequalities_.getWidth()),
    748759  multiplicity(1),
    749   haveExtremeRaysBeenCached(false),
    750   linearForms(ZMatrix(0,ambientDimension))
    751 {
    752 }
    753 
    754 
    755 ZCone::ZCone(ZMatrix const &inequalities_, ZMatrix const &equations_, int preassumptions_):
     760  linearForms(ZMatrix(0,inequalities_.getWidth())),
    756761  inequalities(inequalities_),
    757762  equations(equations_),
    758   state(0),
    759   preassumptions(preassumptions_),
    760   multiplicity(1),
    761   haveExtremeRaysBeenCached(false),
    762   n(inequalities_.getWidth()),
    763   linearForms(ZMatrix(0,inequalities_.getWidth()))
     763  haveExtremeRaysBeenCached(false)
    764764  {
    765765  assert(preassumptions_<4);//OTHERWISE WE ARE DOING SOMETHING STUPID LIKE SPECIFYING AMBIENT DIMENSION
     
    10421042  std::vector<std::vector<int> > indices=lpSolver.extremeRaysInequalityIndices(inequalities);
    10431043
    1044   for(int i=0;i<indices.size();i++)
     1044  for(unsigned i=0;i<indices.size();i++)
    10451045    {
    10461046      /* At this point we know lineality space, implied equations and
     
    10821082
    10831083          std::vector<int> asVector(inequalities.getHeight());
    1084           for(int j=0;j<indices[i].size();j++){asVector[indices[i][j]]=1;}
     1084          for(unsigned j=0;j<indices[i].size();j++){asVector[indices[i][j]]=1;}
    10851085          ZMatrix equations=this->equations;
    10861086          ZVector theInequality;
    10871087
    1088           for(int j=0;j<asVector.size();j++)
     1088          for(unsigned j=0;j<asVector.size();j++)
    10891089            if(asVector[j])
    10901090              equations.appendRow(inequalities[j]);
     
    12291229ZCone ZCone::faceContaining(ZVector const &v)const
    12301230{
    1231   assert(n==v.size());
     1231  assert(n==(int)v.size());
    12321232  assert(contains(v));
    12331233  ZMatrix newEquations=equations;
Note: See TracChangeset for help on using the changeset viewer.