June 2007 - PPL-devel - BUGSENG Mailman

ppl/Watchdog/m4
by Roberto Bagnara 08 Jun '07

08 Jun '07

CVSROOT: /cvs/ppl Module name: ppl Changes by: bagnara(a)cs.unipr.it 2007-06-08 17:12:28 ppl/Watchdog/m4 Update of /cvs/ppl/ppl/Watchdog/m4 In directory spartacus.cs.unipr.it:/tmp/cvs-serv24150/m4 Log Message: Directory /cvs/ppl/ppl/Watchdog/m4 added to the repository

1 0

ppl ppl.spec.in
by Roberto Bagnara 07 Jun '07

07 Jun '07

CVSROOT: /cvs/ppl Module name: ppl Branch: ppl-0_9-branch Changes by: bagnara(a)cs.unipr.it 2007-06-07 16:17:19 Modified files: . : ppl.spec.in Log message: Apply both patches. Patches: http://www.cs.unipr.it/cgi-bin/cvsweb.cgi/ppl/ppl.spec.in.diff?cvsroot=ppl&…

1 0

ppl ppl.spec.in
by Roberto Bagnara 07 Jun '07

07 Jun '07

CVSROOT: /cvs/ppl Module name: ppl Branch: ppl-0_9-branch Changes by: bagnara(a)cs.unipr.it 2007-06-07 15:51:34 Modified files: . : ppl.spec.in Log message: Fixed the command to remove all installed *.la files. Patches: http://www.cs.unipr.it/cgi-bin/cvsweb.cgi/ppl/ppl.spec.in.diff?cvsroot=ppl&…

1 0

ppl ppl.spec.in
by Roberto Bagnara 07 Jun '07

07 Jun '07

CVSROOT: /cvs/ppl Module name: ppl Branch: ppl-0_9-branch Changes by: bagnara(a)cs.unipr.it 2007-06-07 15:43:59 Modified files: . : ppl.spec.in Log message: Remove installed *.la files. Added a `ppl-0.9-configure.patch' to avoid overriding CFLAGS and CXXFLAGS. Patches: http://www.cs.unipr.it/cgi-bin/cvsweb.cgi/ppl/ppl.spec.in.diff?cvsroot=ppl&…

1 0

ppl ppl.spec.in
by Roberto Bagnara 07 Jun '07

07 Jun '07

CVSROOT: /cvs/ppl Module name: ppl Branch: ppl-0_9-branch Changes by: bagnara(a)cs.unipr.it 2007-06-07 12:27:56 Modified files: . : ppl.spec.in Log message: All the static archives are now in `*-static' packages. Packages `ppl-gprolog-devel', `ppl-swiprolog-devel' and `ppl-yap-devel' renamed `ppl-gprolog', `ppl-swiprolog' and `ppl-yap', respectively. As a workaround for a bug in the `pl' package (Bugzilla Bug 243084) `ppl-swiprolog' is now dependent on `readline-devel'. Added `%%dir %%{_datadir}/doc/%%{name}'. The `ppl-user-0.9-html' documentation directory is now properly listed. Patches: http://www.cs.unipr.it/cgi-bin/cvsweb.cgi/ppl/ppl.spec.in.diff?cvsroot=ppl&…

1 0

ppl ppl.spec.in
by Roberto Bagnara 06 Jun '07

06 Jun '07

CVSROOT: /cvs/ppl Module name: ppl Branch: ppl-0_9-branch Changes by: bagnara(a)cs.unipr.it 2007-06-06 22:16:23 Modified files: . : ppl.spec.in Log message: Use `%{_includedir}' and `%{_libdir}' instead of `/usr/include' and `/usr/lib', respectively. Use `%{_datadir}/doc/%{name}' instead of `/usr/share/doc/ppl'. Replaced `%defattr(-,root,root)' with `%defattr(-,root,root,-)'. (Thanks to Mamoru Tasaka.) Patches: http://www.cs.unipr.it/cgi-bin/cvsweb.cgi/ppl/ppl.spec.in.diff?cvsroot=ppl&…

1 0

ppl/interfaces/Java/jni ppl_java_common.hh
by Roberto Bagnara 05 Jun '07

05 Jun '07

CVSROOT: /cvs/ppl Module name: ppl Changes by: bagnara(a)cs.unipr.it 2007-06-05 21:36:52 Modified files: interfaces/Java/jni: ppl_java_common.hh Log message: Shadowing avoided. Patches: http://www.cs.unipr.it/cgi-bin/cvsweb.cgi/ppl/interfaces/Java/jni/ppl_java_…

1 0

ppl/interfaces/Java/jni ppl_java_common.cc
by Roberto Bagnara 05 Jun '07

05 Jun '07

CVSROOT: /cvs/ppl Module name: ppl Changes by: bagnara(a)cs.unipr.it 2007-06-05 21:31:41 Modified files: interfaces/Java/jni: ppl_java_common.cc Log message: Useless break's and braces removed. Indentation fixed in several places. Patches: http://www.cs.unipr.it/cgi-bin/cvsweb.cgi/ppl/interfaces/Java/jni/ppl_java_…

1 0

ppl/src Polyhedron.inlines.hh
by Roberto Bagnara 05 Jun '07

05 Jun '07

CVSROOT: /cvs/ppl Module name: ppl Changes by: bagnara(a)cs.unipr.it 2007-06-05 21:31:00 Modified files: src : Polyhedron.inlines.hh Log message: Make sure Polyhedron::space_dimension() is defined before being used. Patches: http://www.cs.unipr.it/cgi-bin/cvsweb.cgi/ppl/src/Polyhedron.inlines.hh.dif…

1 0

Re: [PPL-devel] PPL's handling of floating-point coefficients
by Enea Zaffanella 05 Jun '07

05 Jun '07

Sebastian Nowozin wrote: > > Hello, > > thanks for making the PPL library open-source software. I have been > using it for a few days now and it works quite well. > > For solving an optimization problem I use a polytope to > outer-approximate a convex set incrementally. I want to use PPL to find > the set of new vertices of this polytope when I add a new halfspace to > it (ie. refine the outer-approximation). The polytope lives in R^n, n > being a fixed dimension and I have no prior knowledge of the nature of > the coefficients, that is, they are not integer but just some floating > point numbers. > > For this, my code looks like (given Ax <= b halfspace representation, > what new vertices are created by a new halfspace row'x <= rowB): > > PPL::Constraint_System cs; > for (unsigned int n = 0 ; n < b.size () ; ++n) { > PPL::Linear_Expression e; > > for (unsigned int j = 0 ; j < gmm::mat_ncols (A) ; ++j) { > if (j > 0) > std::cout << " + "; > > std::cout << A(n,j) << "*x_" << j; > e += mpq_class (A(n,j)) * PPL::Variable (j); > } > > cs.insert (e <= mpq_class (b[n])); > std::cout << " <= " << b[n] << std::endl; > } > > /* Add new row'x = rowB constraint > */ > PPL::Linear_Expression e; > for (unsigned int j = 0 ; j < gmm::mat_ncols (A) ; ++j) { > if (j > 0) > std::cout << " + "; > > std::cout << row[j] << "*x_" << j; > e += mpf_class (row[j]) * PPL::Variable (j); > } > std::cout << " == " << rowB << std::endl; > cs.insert (e == mpf_class (rowB)); > > using namespace PPL::IO_Operators; > std::cout << "Constraint system:\n" << cs << std::endl; > > However, to my surprise, PPL seems to round the coefficients to integer > numbers (as printed by the constraint system). I did not find any hint > in the documentation, except the mention of rational polytopes. Is PPL > able to handle arbitrary floating point coefficients or do I always have > to normalize them in some way? If so, what is the easiest way to do so? Hello. The coefficients of constraints and generators are integers (mpz_class in the default configuration). The domain of convex polyhedra implements rational polyhedra, but only allows for integer coefficients as far as the representation and input/output are concerned. Any constraint having rational coefficients can be normalized so as to have integer coefficients (as for generators, the same can be done by exploiting the integer divisor, which by default is 1; but you seem to be most interested to the constraint representation). The translation from a rational-coefficient constraint (i.e., your representation) to an equivalent integer-coefficient constraint (i.e., the PPL representation) is typically done by computing the lcm of all the denominators and normalizing the numerators accordingly. To this end, you can use helper functions computing the lcm (or gcd) of two coefficients: //! Assigns to \p x the least common multiple of \p y and \p z. /*! \relates GMP_Integer */ void lcm_assign(GMP_Integer& x, const GMP_Integer& y, const GMP_Integer& z); //! Assigns to \p x the greatest common divisor of \p y and \p z. /*! \relates GMP_Integer */ void gcd_assign(GMP_Integer& x, const GMP_Integer& y, const GMP_Integer& z); [NOTE: GMP_Integer is just an alias for mpz_class.] I hope this is enough to bring you to the right track. I am also forwarding this reply to the PPL-devel mailing list, since other users/developers may want to provide more details. Cheers, Enea Zaffanella. > > Thanks, > Sebastian Nowozin > >

1 0