Ask Your Question

Reformulation of a semialgebraic set, without quantifiers

asked 2015-10-27 04:26:42 -0500

pang gravatar image

According to wikipedia:

The Tarski–Seidenberg theorem states that a set in (n + 1)-dimensional space defined by polynomial equations and inequalities can be projected down onto n-dimensional space, and the resulting set is still definable in terms of polynomial identities and inequalities.

That's great, and my question is if it can be done in practice, in a particular case. I have read that the original Tarski–Seidenberg algorithm is of little practical use, but maybe my problem is tractable with another algorithm:

We have a finite set of points in C^2 (z_i,w_i), and I want to project onto the first C coordinate the intersection of the cones {(z,w): |z-z_i|>|w-w_i|}. It would be awesome, for example for plotting the set, if I could actually get the polynomial inequalities that define the projection. Getting only the higher dimensional strata is fine (an open set in C).

Any idea? Thanks in advance.

edit retag flag offensive close merge delete


For the record: I couldn't get qepcad to make a projection of codimension 2 using even 4 points, but I did get some significant results by restricting to a three dimensional space, so that qepcad only projects along lines.

pang gravatar imagepang ( 2015-10-28 06:27:21 -0500 )edit

2 answers

Sort by » oldest newest most voted

answered 2015-10-27 07:20:23 -0500

tmonteil gravatar image

updated 2015-10-27 07:20:48 -0500

You should install and use qepcad. It has been updated recenty, so it should work, just install the latest version of Sage and run sage -i qepcad, see this page for more information.

edit flag offensive delete link more


I downloaded the program and I'm calling it from the command line. I can't get it to work on any non-trivial example (it does work with all the examples from its documentation). The output is:

Failure occurred in: GCSI (final check) Reason for the failure: Too few cells reclaimed.

pang gravatar imagepang ( 2015-10-27 10:19:49 -0500 )edit

It just finished installing. Thanks for taking care of this packege, Thierry!

pang gravatar imagepang ( 2015-10-27 10:31:11 -0500 )edit

Great. Please do not hesitate to tell if there is some issue with this package.

tmonteil gravatar imagetmonteil ( 2015-10-27 10:38:56 -0500 )edit

I get a "RuntimeError: unable to start QEPCAD". I upgraded to sage 6.9 with

./sage -upgrade

and I installed qepcad with

./sage -i qepcad

The first upgrade reported errors when building the html docs, but nothing else, and everything seems to work fine. The installation of qepcad produced no errors (I think), and there are qepcad and saclib folders in sage/pkgs.

pang gravatar imagepang ( 2015-10-28 04:26:30 -0500 )edit

I have bug traces etc: should I open a ticket?

pang gravatar imagepang ( 2015-10-28 04:36:28 -0500 )edit

answered 2015-10-27 04:30:09 -0500

pang gravatar image

updated 2015-10-27 10:31:54 -0500

I found the optional package qepcad:

mentioned in this site:

edit flag offensive delete link more

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

1 follower


Asked: 2015-10-27 04:26:42 -0500

Seen: 214 times

Last updated: Oct 27 '15