/* uniqtest.c: Test program to call the cdd library cddlib written by Komei Fukuda, fukuda@ifor.math.ethz.ch Version 0.93c, Nov. 14, 2003 Standard ftp site: ftp.ifor.math.ethz.ch, Directory: pub/fukuda/cdd */ /* This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version. This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. */ #include "setoper.h" #include "cdd.h" #include <stdio.h> #include <stdlib.h> #include <time.h> #include <math.h> #include <string.h> dd_boolean SetInputFile(FILE **f, dd_DataFileType fname) { dd_boolean success=dd_FALSE; success=dd_FALSE; if ( ( *f = fopen(fname, "r") )!= NULL) { printf("input file %s is open\n", fname); success=dd_TRUE; } else{ printf("The input file %s not found\n",fname); } return success; } dd_boolean SetWriteFile(FILE **f, dd_DataFileType fname) { dd_boolean success=dd_FALSE; if ( (*f = fopen(fname, "w")) != NULL){ printf("output file %s is open\n",fname); success=dd_TRUE; } else{ printf("The output file %s cannot be opened\n",fname); } return success; } int main(int argc, char *argv[]) { dd_MatrixPtr M=NULL,M1=NULL,M2=NULL; dd_rowrange i; dd_colrange d; dd_ErrorType err=dd_NoError; dd_rowset redrows,linrows,posset; dd_rowindex newpos=NULL, newpos1=NULL, newpos2=NULL; mytype val; dd_DataFileType inputfile; FILE *reading=NULL; int foi; dd_Arow certificate; dd_set_global_constants(); /* First, this must be called. */ dd_init(val); if (argc>1) strcpy(inputfile,argv[1]); if (argc<=1 || !SetInputFile(&reading,argv[1])){ dd_WriteProgramDescription(stdout); fprintf(stdout,"\ncddlib test program to remove duplicates.\n"); dd_SetInputFile(&reading,inputfile, &err); } if (err==dd_NoError) { M=dd_PolyFile2Matrix(reading, &err); } else { fprintf(stderr,"Input file not found\n"); goto _L99; } if (err!=dd_NoError) goto _L99; d=M->colsize; /* printf("\nInput Matrix.\n"); dd_WriteMatrix(stdout, M); */ M1=dd_MatrixSortedUniqueCopy(M,&newpos1); printf("\nSort and remove duplicates with dd_MatrixSortedUniqueCopy.\n"); printf(" Row index changes -- original:new\n"); for (i=1;i<=M->rowsize; i++){ printf(" %ld:%ld",i,newpos1[i]); } printf("\n"); /* dd_WriteMatrix(stdout, M1); */ dd_InitializeArow(M1->colsize+2, &certificate); fprintf(stdout, "\nCheck whether the system contains an implicit linearity.\n"); foi=dd_FreeOfImplicitLinearity(M1, certificate, &posset, &err); switch (foi) { case 1: fprintf(stdout, " It is free of implicit linearity.\n"); break; case 0: fprintf(stdout, " It is not free of implicit linearity.\n"); break; case -1: fprintf(stdout, " The input system is trivial (i.e. the empty H-polytope or the V-rep of the whole space.\n"); break; } fprintf(stdout, "\nOne can then remove nontrivial redundant rows with dd_RedundantRows.\n Redundant rows:\n"); redrows=dd_RedundantRowsWithPresorting(M1, &err); /* redrows=dd_RedundantRows(M1, &err); */ set_fwrite(stdout, redrows); printf("\n"); M2=dd_MatrixSubmatrix2(M1, redrows, &newpos); dd_FreeMatrix(M1); set_free(redrows); free(newpos); fprintf(stdout, " Implicit linearity (after removal of redundant rows): "); linrows=dd_ImplicitLinearityRows(M2, &err); if (M->representation==dd_Generator) fprintf(stdout," %ld ", set_card(linrows)); else fprintf(stdout," %ld ", set_card(linrows)); set_fwrite(stdout,linrows); set_uni(M2->linset, M2->linset, linrows); /* add the implicit linrows to the given linearity rows */ printf("\nNonredundant representation (except for the linearity part):\n"); dd_WriteMatrix(stdout, M2); dd_WriteLPStats(stdout); set_free(linrows); dd_FreeMatrix(M); dd_FreeMatrix(M2); dd_clear(val); free(newpos1); free(newpos2); _L99:; /* if (err!=dd_NoError) dd_WriteErrorMessages(stderr,err); */ return 0; } /* end of uniqtest.c */