SUMMARYA self-adaptive method for the solution of elliptic and parabolic problems defined over 3D general multiconnected regions by generalized finite-difference formulae, with boundary approximation by means of polyhedra, is proposed. Different approaches to space refinements are studied and implemented. Numerical examples are given.