: In computer-aided design, geometric modeling by constraints enables users to describe shapes by relationships called constraints between geometric elements. The problem is to derive automatically these geometric elements. Many resolution methods have been proposed for solving systems of geometric constraints. Geometric methods can be very efficient but are only applicable to particular kinds of problems (typically rule and compass constructive problems). Some schemes can't be solved by any geometric method. A numerical method is then necessary. We use the bisection method to solve these cases. The bisection method enables users to reliably find all solutions to a system of non-linear equations within a region defined by bounds on each individual coordinate of the geometric objects. But bisection is very time consuming. To overcome this problem, we present first a method that speed-up the resolution process. This improvement is still insufficient; we propose then to distribute th...