This is a follow-up (but self-contained) question to my previous one. There I asked about state-of-the-art methods to solve multivariate polynomials systems over non-algebraically closed fields in general.

I learned that the theory is more involved that I thought (I'm not working in algorithmic algebraic geometry, so I'm only familiar with the very basics, like Buchberger's algorithm, or the definition of the dimension of an ideal). Therefore it is necessary to ask a more specific questions that the previous, general one, which is more tailored to my needs.

My setup is the following:

*Regarding complexity:* I'm interested in solving a large number of polynomial systems (on commodity hardware), on the order of $10^4$. But each of the systems is of relatively small size - my baseline consists of least 6 different variable and 4 equations. If I could tackle this, I'd already be happy. Going further, I don't expect the systems to grow beyond about 20 different variables and 20 equations.

So perhaps I don't actually need the fastest possible algorithm and can make do with simpler, older ones - but I will let you be the judge of that.

*Regarding the polynomials:* There are no restrictions their coefficients, so, depending on the field I'm working in they can take any number.

*Regarding the field:* Regarding the field I'm working in, my baseline is $\mathbb{R}$, but I'd also be interested in $\mathbb{Q}$ and $\mathbb{Z}$. If there are methods that are much easier for one field than another, than I will the choice of the field to study be influence by the time I need to invest to learn that method, i.e. the easiest one wins.

*Regardin the dimension of the ideal spanned by the polynomials:* The ideal has dimension $2$ or $3$ over the complex numbers, in most cases I tested so far with the help of CAS.

*What I'm looking for*: I'm interested is learning about methods (I'm happy with specific references) that tell me

1) *whether* the system has a solution at all or not. Working over, e.g., $\mathbb{C}$, this would be easy (e.g. compute a Gröbner basis: If it contain the $1$, if and only if the solution variety is empty). But this doesn't work unfortunately for non-algebraically closed fields. Given the answers from my previous question, I'm inclined to think that answer this question shouldn't be too hard (perhaps even trivial for the expert computational geometer, which I'm not unfortunately).

2) if it has an infinite number of solutions (if the variety is zero-dimensional, things are easy of course), I would like to pick out one single variable, say $n_0\in \{1,\ldots,n\}$, project the solution variety $V(f_1,\ldots,f_s)\subseteq \mathbb{ R}^n$ (supposing we work over the field $\mathbb{R}$) along this variable onto $\mathbb{R}$ to investigate whether there exists an interval $[-\alpha,\alpha]$ around $0$ which is contained in this projected set (I don't need to understand the projected set fully). That is what I menat by "rough shape" in the title.