We study geometrical properties of the complete set of solutions of the random 3-satisfiability problem. We show that even for moderate system sizes the number of clusters corresponds surprisingly well with the theoretic asymptotic prediction. We locate the freezing transition in the space of solutions which has been conjectured to be relevant in explaining the onset of computational hardness in random constraint satisfaction problems. PACS numbers: 89.20.Ff,75.10.Nr,89.70.Eg Satisfiability (SAT) is one of the most important problems in theoretical computer science. It was the first problem shown to be NP-complete [1,2], and it is of central relevance in various practical applications, including artificial intelligence, planning, hardware and electronic design, automation, verification and more. It can thus be pictorially thought of as the Ising model of computer science. Ensembles of randomly generated SAT instances emerged in computer science as a way of evaluating algorithmic performance and addressing questions regarding the average case complexity.An instance of random K-SAT problem consists of N Boolean variables and M clauses. Each clause contains a subset of K distinct variables chosen uniformly at random, and each clause forbids one random assignment of the K variables out of the 2 K possible ones. The problem is satisfiable if there exists a variable assignment that simultaneously satisfies all clauses and we call such an assignments a solution to the problem. When the density of constraints α = M/N is increased, the formulas become less likely to be satisfiable. In the thermodynamical limit there is a sharp transition from a phase in which the formulas are almost surely satisfiable to a phase where they are almost surely unsatisfiable. The existence of this transition is partly established rigorously [3]. It is also a well known empirical result that the hardest instances are found near to this threshold [4,5,6].Random K-SAT has attracted interest of statistical physicists because of its equivalence to mean field spin glasses [7]. Indeed, the problem can be rephrased as minimizing a spin glass-like energy function which counts the number of violated clauses. The results and insights coming from this equivalence are remarkable. The satisfiability threshold and other phase transitions in the structure of solutions are described in [8,9,10]. In particular, it was shown that for K ≥ 3 the space of solutions for highly constrained but still satisfiable instances splits into exponentially many clusters and in some cases this clustering has been rigorously confirmed [11,12]. The so-called freezing of variables in clusters is another rich concept studied recently [13,14]. However, a detailed understanding of how the clustering or freezing of solutions affects the average computational hardness is still one of the most interesting open problems in the field.Since the exact statistical physics solution of the random satisfiability problem appeared [9,15] dozens of directly related articles followed. Mathemat...